matlab期中程序

1、測(cè)試(第2,3,4題選做，其余題必做)1.解答題(1) 寫(xiě)出執(zhí)行以下代碼后，MATLAB命令窗口上顯示的x矩陣的值 x=0,2,0,2,0,3,0,4;for k=1:8if x(k)=0x(k)=k;elsex(k)=2*k+1;endend dis p(x);結(jié)果:39513717(2)寫(xiě)出執(zhí)行以下代碼后C, D , E的值A(chǔ)=1,2,3;4,5,6;7,8,9;C=A;10,17,12D=C(1:3,2 3)E=C(2,1 2)結(jié)果: C1012 D E2.MATLAB7.0的符號(hào)計(jì)算功能計(jì)算下列題目x(1 )求 exsin x 在3,18上的定積分;(2) 求解非線性方程組sin(x)

2、 2cos(y) 02cos(x) 3sin(y) 0 (3 )計(jì)算級(jí)數(shù)501S；n 1 n(4 )設(shè) f Xyey 2t,求復(fù)合函數(shù)f (g(z)；(5)計(jì)算zxy32e xy 12x的二.x3.利用MATLAB7.0的符號(hào)計(jì)算功能計(jì)算下列題目(1 )求定積分 ：e dx ；求解非線性方程組x 0.7sin(x) 0.2cos(y)y 0.7cos(x) 0.2sin(y)計(jì)算級(jí)數(shù)S nS ；求極限 lim tanx 3sinxX 0sin X223Z(5)計(jì)算Z yx 3y x 3y的導(dǎo)數(shù)一。y4.利用matlab7.0 計(jì)算下列題目2小時(shí)的室外溫度(C )(1 )某觀測(cè)站測(cè)得某日6:00

3、時(shí)至18:00時(shí)之間每隔t=15,19,24,28,34,32,29用3次樣條插值求該日室外6:30至17:30之間每隔2小時(shí)各點(diǎn)的近似溫度(C )；(2)采用8階多項(xiàng)式對(duì)函數(shù) y x3 5sin(x),x0,10進(jìn)行擬合;(3)求f(x) x3 2x 5在0,10內(nèi)的最小值點(diǎn)。5.利用MATLAB7.0 計(jì)算下列題目:建立一個(gè)6階魔方矩陣A,(1 )求A的逆;(2)求A的行列式;(3 )求A的秩；(4)求A的跡；(5)求A的特征值。解:6階魔方矩陣A: A=magic(6)351626192433272123253192222720828331710153053412141643629131

4、8111.60840.0000-1.6084-1.60840.00001.60843.21690.0000-3.2169-3.21690.00003.21693.21690.0000-3.2169-3.21690.00003.21691.60840.00001.60841.60840.0000-1.6084(1 ) A 的逆:inv(A)結(jié)果: inv(A)ans =1.0e+014 *-3.2169-0.00003.21693.2169-0.0000-3.2169-3.2169-0.00003.21693.2169-0.0000-3.2169(2 ) A 的行列式：d=det(A) d=de

5、-t (A結(jié)果: d=det(A)0(3) A 的秩：rank(A)結(jié)果: ran k(A)ans =(4) A 的跡：trace(A)結(jié)果: trace(A)ans =111（5） A的特征值：eig（A）結(jié)果: eig(A)ans =111.000027.0000-27.00009.7980-0.0000-9.79806.利用matlab7.0 計(jì)算下列題目（1 ）某實(shí)驗(yàn)對(duì)一根長(zhǎng)10米的鋼軌進(jìn)行熱源的溫度傳播測(cè)試。用X表示測(cè)量點(diǎn)025:10（米），T=95,14,0,0,0;88,48,32,12,6;67,64,54,48,42( 秒、鋼軌每隔1米處的溫度TI;用h表示測(cè)量時(shí)間0:30:

6、60（秒），用T表示測(cè)試所得各點(diǎn)的溫度C ）。試用線性插值求出在一分鐘內(nèi)每隔20（2 ）隨機(jī)產(chǎn)生一向量 X,求向量x的最大值;解:(1)命令如下:x=0:2.5:10;h=0:30:60;T=95,14,0,0,0;88,48,32,12,4;47,44,54,48,41;xi=0:10;hi=0:20:60;TI=i nterp2(x,h,T,xi,hi)結(jié)果:TI =Columns 1 through 795.000062.600030.200011.20005.60000090.333368.866747.400033.600027.466721.333316.000074.333363

7、.266752.200045.200042.266739.333333.200047.000045.800044.600046.000050.000054.000051.6000Colu mns 8 through 11000010.66676.93334.80002.666727.066722.466719.400016.333349.200046.600043.800041.0000(2 ) x=rand(); m=max(x)0.4001運(yùn)用函數(shù)rand生成了一組10個(gè)在01之間的隨機(jī)數(shù)X中最大值：max(x)結(jié)果: 最大值為0.9830(3) 命令如下: syms x fun=in l

8、i ne(x -0x + 2 ,x) x,f=fsolve(fu n,0.55)結(jié)果:0.3758-9.2419e-0117.利用matlab7.0 繪制以下圖形解:(1)( 1)在同一直角坐標(biāo)系內(nèi)畫(huà)出y xsinx和y xsin1在區(qū)間1,3 的圖形，添加xx軸和y軸的簡(jiǎn)單注解，對(duì)圖形進(jìn)行文本標(biāo)注，并嘗試不同的線形、顏色;命令如下x=li nsp ace(1,3* pi,1000);y1=x.*s in( x);y2=x.*si n( 1./x);plot(x,y1,-r,x,y2,-g)xlabel(x);ylabel(y=xs in x,y=xs in (1/x);(2)分別在同一圖形窗

9、口的不同子圖繪制y1 si nt, y2 cost, y3 si n 2t, y4 si nt在區(qū)間0,3 上的圖像。命令如下:t=0: pi/10:3* pi;sub plot(2,2,1)plot(t,s in (t)sub plot(2,2,2)p lot(t,cos(t) sub plot(2,2,3) p lot(t,si n( 2*t) sub plot(2,2,4)plot(t,abs(si n( t)8.利用matlab7.0 繪制以下圖形:解:(1)在同一窗口分別用紅色實(shí)線和綠色虛線繪制y1 sinx和y2 cosx在區(qū)間0,2上的曲線，并用星號(hào)*標(biāo)出兩條曲線的交點(diǎn)以及建立圖例命令如下:x=li nsp ace(0,2* pi,1000);y1=s in( x);y2=cos(x);k=fi nd(abs(y1-y2)=90dis p(A);elseif p=80dis p(B);elseif p=70dis p(C);else dis p(D);end結(jié)果:請(qǐng)輸入分埶BEA請(qǐng)輸入分?jǐn)?shù)怕rB請(qǐng)輸入分?jǐn)?shù)M目C請(qǐng)輸人分埶曲HD6=1+2+3 ，所以6是完數(shù)。10若一個(gè)數(shù)等于它的各個(gè)真因子之和，則稱該數(shù)為完數(shù)，如求1,1000之間的全部完數(shù)。解:function fun cti on5()result=;for n=1:1