MATLAB上机作业

MATLAB上机作业1

对以下问题,编写M文件：

(1) 用起泡法对10个数由小到大排序。即将相邻两个数比较,将小的调到前头。 function f=qipaofa(x) for j=9:-1:1 for i=1:j

if(x(i)>x(i+1))

t=x(i);x(i)=x(i+1);x(i+1)=t; end end end f=x

x=round(10*rand(1,10)) qipaofa(x);

(2) 有一个4×5矩阵，编程求出其最大值及其所处的位置。 function f=zuidazhi(x) a=1;b=1;c=x(1,1); for i=1:4 for j=1:5

if x(i,j)>c

a=i;b=j;c=x(i,j); end end end

f=[c,a,b]

x=rand(4,5) zuidazhi(x)

20(3) 编程求?n!。

n?1function f=qiuhe(x) b=0; for i=1:x a=prod(1:i); b=b+a; end f=b

qiuhe(20)

(4)一球从100米高度自由落下,每次落地后反跳回原高度的一半，再落下。在第10次落地时，共经过多少米?第10次反弹有多高? function f=gao(x)

求它 b=x; for i=2:10 x=x/2; a=x*2; b=b+a; end f=[b x/2]

gao(100)

(5)有一函数

f(x,y)?x2?sinxy?2y，写一程序，输入自变量的值,输出函数值。

Function f=fun(x)

f=x(1)^2+sin(x(1)*x(2))+2*x(2)

MATLAB上机作业2

1. 求和

Y?1?4?4???4240。

syms k s=4^k;

S=symsum(s,k,0,40)

2. 求函数

f(x)?2x?6x32?18x?7的极值，并作图。

y='2*x^3-6*x^2-18*x+7'; y_='-2*x^3+6*x^2+18*x-7'; [x_min,y_min]=fminbnd(y,-7,7) [x_max,y_max]=fminbnd(y_,-7,7) ezplot(y)

3. 设y?esinx?7cosx?5xx2，求

dydx??,dydx22??

y=exp(x)*sin(x)-7*cos(x)+5*x^2; dy=diff(y) d2y=diff(y,2)

4. 求积分

I1???/201?2sin2xdx。

syms x

y=sqrt(1-2*sin(2*x)); F=int(y,0,pi/2)

5. 求积分

I2???Dxy22dxdy，其中D为直线y?2x,y?x/2,y?12?x围成的区域。

syms x y

f=x^2/y^2;

int(int(f,y,x/2,x*2),x,0,4)+int(int(f,y,x/2,12-x),x,4,8)

6. 分别随机产生一个6×6的整数矩阵（元素可在[?20,20]之间），求该随机阵

的行最简形，秩，行列式，逆矩阵，列向量组的一个极大线性无关组，特征值和特征向量。

a=round(-20+40*rand(6)) rref(a) rank(a) det(a) inv(a)

[V,D]=eig(a)

?x1?2x2?4x3?6x4?3x5?2x6?4?2x?4x2?4x3?5x4?x5?5x6?3?1??3x1?6x2?2x3?5x5?9x6??1 ?2x?3x?4x?x?8246?1??4x?5x?2x?x?4x??523456???5x1?5x2?3x3?6x4?6x5?4x6?27. 求解方程组的通解

A=[1 2 4 6 -3 2;2 4 -4 5 1 -5;3 6 2 0 5 -9;2 3 0 4 0 1;0 -4 -5 2 1 4 ;5 5 -3 6 6 -4 ];

b=[4 3 -1 8 -5 2]'; B=[A b]; C=rref(B)

A=[1 2 4 6 -3 2;2 4 -4 5 1 -5;3 6 2 0 5 -9;2 3 0 4 0 1;0 -4 -5 2 1 4 ;5 5 -3 6 6 -4 ]; b=[4 3 -1 8 -5 2]'; B=[A b]; n=6;

R_A=rank(A) R_B=rank(B) format rat

if R_A==R_B&R_A==n X=A\\b

elseif R_A==R_B&R_A

C=null(A,'r')

else X='equation has no solves' end

syms k1 k2 X=k1*C+X pretty(X)

8. 求符号微分方程

dydx?ytanx?cosx的通解和当y(0)=2的特解。

eq='Dy+y*tan(x)=cos(x)'; y1=dsolve(eq,'x')

y2=dsolve(eq,'y(0)=2','x')

MATLAB上机作业3

1、某校60名学生的一次考试成绩如下:

93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55

1)计算均值、标准差、极差、偏度、峰度，画出直方图； 2)检验分布的正态性；

3)若检验符合正态分布，估计正态分布的参数并检验参数。

x=[93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55]; mean(x)

std(x) range(x) skewness(x) kurtosis(x) hist(x)

h=normplot(x)

[muhat,sigmahat,muci,sigmaci]=normfit(x) [H,sig,ci]=ttest(x,80.1)

2、据说某地汽油的价格是每加仑115美分，为了验证这种说法，一位学者开车随机选择了一些加油站，得到某年一月和二月的数据如下：

一月：119 117 115 116 112 121 115 122 116 118 109 112 119 112 117 113 114 109 109 118

二月：118 119 115 122 118 121 120 122 128 116 120 123 121 119 117 119 128 126 118 125

1）分别用两个月的数据验证这种说法的可靠性； 2）分别给出1月和2月汽油价格的置信区间； 3）给出1月和2月汽油价格差的置信区间。

x=[119 117 115 116 112 121 115 122 116 118 109 112 119 112 117 113 114 109 109 118]; y=[118 119 115 122 118 121 120 122 128 116 120 123 121 119 117 119 128 126 118 125]; [H,sig,ci]=ttest(x,115) [H,sig,ci]=ttest(y,115) [h,sig,ci]=ttest2(x,y)

1、在同一平面中的两个窗口分别画出心形线和马鞍面。 要求：1、在图形上加格栅、图例和标注

2、定制坐标

3、以不同角度观察马鞍面 a=3;b=1;

theta=linspace(-2*pi,2*pi); rho=a*(1-cos(theta));

subplot(1,2,1);polar(theta,rho)

gtext('rho');gtext('theta');title('心形线'); axis([-10 10 -10 10]);grid on;

x=-a:0.1:a;y=-b:0.1:b; [x,y]=meshgrid(x,y);

z=x.^2/(a^2*2)-y.^2/(b^2*2); subplot(1,2,2);surf(z)

xlabel('x轴');ylabel('y轴');zlabel('z轴');title('马鞍面'); axis([0 100 0 30 -1 1]);grid on

2、以不同的视角观察球面 x2+y2+z2=r2 和圆柱面 x2+y2=rx所围区域。

[x,y,z]=sphere(200); i=find(x.^2+y.^2-x>=0); z(i)=NaN; mesh(x,y,z) hold on

ezplot('x^2+y^2=x')

本文来源：https://www.bwwdw.com/article/sg7x.html