MATLAB 上机作业1对以下问题,编写M 文件：(1) 用起泡法对10个数由小到大排序。

即将相邻两个数比较,将小的调到前头。

function f=qipaofa(x)for j=9:-1:1for i=1:jif(x(i)>x(i+1))t=x(i);x(i)=x(i+1);x(i+1)=t;endendendf=xx=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:4for j=1:5if x(i,j)>ca=i;b=j;c=x(i,j);endendendf=[c,a,b]x=rand(4,5)zuidazhi(x)(3) 编程求∑=201!n n 。

function f=qiuhe(x)b=0;for i=1:xa=prod(1:i);b=b+a;endf=bqiuhe(20)(4)一球从100米高度自由落下,每次落地后反跳回原高度的一半，再落下。

求它在第10次落地时，共经过多少米?第10次反弹有多高?function f=gao(x)b=x;for i=2:10x=x/2;a=x*2;b=b+a;endf=[b x/2]gao(100)(5)有一函数y xy x y x f 2sin ),(2++=，写一程序，输入自变量的值,输出函数值。

Function f=fun(x)f=x(1)^2+sin(x(1)*x(2))+2*x(2)MATLAB 上机作业21. 求和 4024441++++=ΛY 。

syms ks=4^k;S=symsum(s,k,0,40)2. 求函数71862)(23+--=x x x x f 的极值，并作图。

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. 设25cos 7sin x x x e y x +-=，求??,22==dx y d dx dy y=exp(x)*sin(x)-7*cos(x)+5*x^2;dy=diff(y)d2y=diff(y,2)4. 求积分 dx x I ⎰-=2/012sin 21π。

syms xy=sqrt(1-2*sin(2*x));F=int(y,0,pi/2)5. 求积分 ⎰⎰=D dxdy y x I 222，其中D 为直线x y x y x y -===12,2/,2围成的区域。

syms x yf=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)7. 求解方程组的通解 ⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=-++-+-=+++--=+++-=-+++=-++-+=+-+++246635554254843219526335544242364265432165432642165321654321654321x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x xA=[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 ratif R_A==R_B&R_A==nX=A\belseif R_A==R_B&R_A<nX=A\bC=null(A,'r')else X='equation has no solves'endsyms k1 k2X=k1*C+Xpretty(X)8. 求符号微分方程 x x y dxdy cos tan =+的通解和当y (0)=2的特解。

eq='Dy+y*tan(x)=cos(x)';y1=dsolve(eq,'x')y2=dsolve(eq,'y(0)=2','x')MATLAB 上机作业31、某校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 551)计算均值、标准差、极差、偏度、峰度，画出直方图；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 1251）分别用两个月的数据验证这种说法的可靠性；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 on2、以不同的视角观察球面 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 onezplot('x^2+y^2=x')（注：专业文档是经验性极强的领域，无法思考和涵盖全面，素材和资料部分来自网络，供参考。

可复制、编制，期待你的好评与关注）。