(完整word版)MATLAB上机作业

1、MATLAB 上机作业 1对以下问题编写 M 文件：(1) 用起泡法对 10 个数由小到大排序。即将相邻两个数比较fun cti on 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=rou nd(10*ra nd(1,10)qipaofa(x);(2) 有一个 4X5 矩阵，编程求出其最大值及其所处的位置。fun cti on 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(

2、i,j);endendendf=c,a,bx=ra nd(4,5)zuidazhi(x)20编程求n!。n =1fun cti on f=qiuhe(x)b=0;for i=1:xa=prod(1:i);b=b+a;endf=bqiuhe(20)(4)一球从 100 米高度自由落下，每次落地后反跳回原高度的一 在第 10 次落地时，共经过多少米？第 10 次反弹有多高？fun cti on f=gao(x)，将小的调到前再落下。求它b=x;for i=2:10 x=x/2; a=x*2;b=b+a; end f=b x/2gao(100)有一函数f(x, y)=x2 sin xy 2y，写一程

3、序，输入自变量的值，输出函数值。Function f=fun(x)f=x(1)A2+si n(x(1)*x(2)+2*x(2)MATLAB 上机作业 21.求和 Y =1 4 42 440。syms k s=4Ak;S=symsum(s,k,0,40)2.求函数f (x) =2x3-6x2-18x 7的极值，并作图。y=2*xA3-6*xA2-18*x+7; y_=-2*xA3+6*xA2+18*x-7;x_mi n, y_mi n=fmi nbn d(y,-7,7) x_max,y_max=fmi nbn d(y_,-7,7) ezplot(y) y=exp(x)*sin(x)-7*cos(

4、x)+5*xA2;dy=diff(y)d2y=diff(y,2)4.求积分 h = o d - 2 sin2xdx。syms xy=sqrt(1-2*si n(2*x);F=in t(y,0,pi/2)x5. 求积分丨丨22dxdy，其中 D 为直线y = 2x, y =x/2, y =12-x围成的区域。ysyms x yf=xA2/yA2;in t(i nt(f,y,x/2,x*2),x,0,4)+in t(i nt(f,y,x/2,12-x),x,4,8)6. 分别随机产生一个 6X6 的整数矩阵(元素可在-20, 20之间)，求该随机阵 的行最简形，秩，行列式，逆矩阵，列向量组的一个极

5、大线性无关组，特征值和 特征向量。a=rou nd(-20+40*ra nd)rref(a)ran k(a)det(a)in v(a)3.2sin x - 7 cosx 5x，V,D=eig(a)7.求解方程组的通解捲+2x2+4x3+6x4_3x5+2x6=42 x|+ 4 X2 4 X3+5X4+X5 5 x 3 +6x2+ 2x3+5x5-9x6=-12% +3x2+4x4+x6=84x25x3+ 2x4+ x5+ 4x6= -55% +5x2_3x3+6X4+6X5_4x6= 2A=1 2 4 6 -3 2;2 4 -4 5 1b=4 3 -1 8 -5 2;B=A b;-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 ;C=rref(B)A=1 2 4 6 -3 2;2 4 -4 5 1b=4 3 -1 8 -5 2;B=A b;n=6;R_A=ra nk(A)R_B=ra nk(B)-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 ;format ratif R_