(完整word版)MATLAB上機(jī)作業(yè)

1、MATLAB 上機(jī)作業(yè) 1對(duì)以下問(wèn)題編寫(xiě) M 文件：(1) 用起泡法對(duì) 10 個(gè)數(shù)由小到大排序。即將相鄰兩個(gè)數(shù)比較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) 有一個(gè) 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 次落地時(shí)，共經(jīng)過(guò)多少米？第 10 次反彈有多高？fun cti on f=gao(x)，將小的調(diào)到前再落下。求它b=x;for i=2:10 x=x/2; a=x*2;b=b+a; end f=b x/2gao(100)有一函數(shù)f(x, y)=x2 sin xy 2y，寫(xiě)一程

3、序，輸入自變量的值，輸出函數(shù)值。Function f=fun(x)f=x(1)A2+si n(x(1)*x(2)+2*x(2)MATLAB 上機(jī)作業(yè) 21.求和 Y =1 4 42 440。syms k s=4Ak;S=symsum(s,k,0,40)2.求函數(shù)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圍成的區(qū)域。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. 分別隨機(jī)產(chǎn)生一個(gè) 6X6 的整數(shù)矩陣(元素可在-20, 20之間)，求該隨機(jī)陣 的行最簡(jiǎn)形，秩，行列式，逆矩陣，列向量組的一個(gè)極

5、大線性無(wú)關(guān)組，特征值和 特征向量。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_