MATLAB程序設(shè)計(jì)與應(yīng)用(第二版)課后實(shí)驗(yàn)答案

1、Matlab課后實(shí)驗(yàn)題答案實(shí)驗(yàn)一 MATLAB運(yùn)算基礎(chǔ)1 .先求下列表達(dá)式的值，然后顯示MATLAB工作空間的使用情況并保存全部變量。2sin85(2)z2 = ln( x + M +x2),其中 x =2lt-0.451 2i5_0.3asin(a 0.3) ln0.3 aa = -3.0, -2.9, HI, 2.9, 3.00 Mt ； 1(4)4=t2 11 Mt <2 ,其中 t=0:0.5:2.52 _. 一 一t -2t+1 2 Mt <3解：M文件：z1=2*sin(85*pi/180)/(1+exp(2)x=2 1+2*i;-.45 5;z2=1/2*log(x+

2、sqrt(1+xA2)a=-3.0:0.1:3.0;z3=(exp(0.3.*a)-exp(-0.3.*a)./2.*sin(a+0.3)+log(0.3+a)./2)t=0:0.5:2.5;z4=(t>=0&t<1).*(t.A2)+(t>=1&t<2).*(t.A2-1)+(t>=2&t<3) .*(t.A2-2*t+1)2.已知:1234-4一13A =34787,B= 2033657 一1?-27求下列表達(dá)式的值：(1) A+6*B和A-B+I (其中I為單位矩陣)(2) A*B 和 A.*B(3) AA3 和 A.A3(4)

3、 A/B 及 BA(5) A,B和A(1,3,:);BA2解：M文件：A=12 34 -4;34 7 87;3 65 7；B=1 3 -1;2 0 3;3 -2 7;A+6.*BA-B+eye(3)A*BA.*BAA3A.A3A/BBAA，BA(1,3,:);BA2(6) 有矩陣A和B45 13091017_61415,B =0231920972425 14131619-4011123678A = 11 12 1316 17 18 -21 22 23(1)求它們的乘積 Co(2)將矩陣C的右下角3X2子矩陣賦給 Do(3)查看MATLAB工作空間的使用情況。解：.運(yùn)算結(jié)果：E=(reshape

4、(1:1:25,5,5)'F=3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11;C= E*FH=C(3:5,2:3)C =9315077258335237423520397588705557753890717H =5203977055578907174.完成下列操作：(1) 求100,999之間能被21整除的數(shù)的個(gè)數(shù)。(2)建立一個(gè)字符串向量，刪除其中的大寫字母。解：(1)結(jié)果：m=100:999;n=find(mod(m,21)=0);length(n)ans =43(2) .建立一個(gè)字符串向量例如：ch='ABC123d4e56Fg9'則要求

5、結(jié)果是:ch='ABC123d4e56Fg9'k=find(ch>='A'&ch<='Z');ch(k戶口ch =123d4e56g9實(shí)驗(yàn)二MATLAB矩陣分析與處理1 .設(shè)有分塊矩陣A = ,E3>3 R3* I其中E、R、O、S分別為單位矩陣、隨機(jī)矩陣、零矩 02 3 S22陣和對(duì)角陣，試通過(guò)數(shù)值計(jì)算驗(yàn)證A2E R RSI2。PS2 1解：M文件如下;由ans,所以A2E R+RS1|OS -2 .產(chǎn)生5階希爾伯特矩陣 H和5階帕斯卡矩陣 P,且求其行列式的值 Hh和Hp以及它們 的條件數(shù)Th和Tp,判斷哪個(gè)矩陣性能

6、更好。為什么？解：M文件如下：因?yàn)樗鼈兊臈l件數(shù) Th>>Tp,所以pascal矩陣性能更好。3 .建立一個(gè)5X5矩陣，求它的行列式值、跡、秩和范數(shù)。解： M文件如下：4 .已知-29 6 18A= 20 5 12-8 85_求A的特征值及特征向量，并分析其數(shù)學(xué)意義。 解：M文件如圖:Edit or - Dnt i±led*百叵區(qū)IFile Ed.i t Text Go CeU TodIs. Debu; Desktop Window Help唱唱 |i o_ + _ 1. J « 滋滋 0<clcA=-29 6 18 ; 20 5 12.-6 8 51 .V

7、,D=eiff(A)因全部特征值構(gòu)成對(duì)用陣D，¥是特征向量做為列向量構(gòu)成數(shù)學(xué)意義：V的3個(gè)列向量是A的特征向量，D的主對(duì)角線上3個(gè)是A的特征值，特另IJ的,V的3個(gè)列向量分別是 D的3個(gè)特征值的特征向量。5.下面是一個(gè)線性方程組：111r 3 4 -方1-0.9511111I - x2 = 0.673 4 5111p3P52.4 5 6(1)求方程的解。(2)將方程右邊向量元素 b3改為0.53再求解，并比較 b3的變化和解的相對(duì)變化。(3)計(jì)算系數(shù)矩陣A的條件數(shù)并分析結(jié)論。解：M文件如下：IBB Editor - Tintitlfd2»輸出結(jié)果：由結(jié)果，X和X2的值一樣，

8、這表示b的微小變化對(duì)方程解也影響較小，而 A的條件數(shù)算得較小，所以數(shù)值穩(wěn)定性較好，A是較好的矩陣。6.建立A矩陣，1式比較 sqrtm(A)和sqrt(A),分析它們的區(qū)別。解：M文件如下：分析結(jié)果知：sqrtm(A)是類似A的數(shù)值平方根(這可由 b1*b1=A的結(jié)果看出)，而sqrt(A) 則是對(duì)A中的每個(gè)元素開根號(hào)，兩則區(qū)別就在于此。實(shí)驗(yàn)三選擇結(jié)構(gòu)程序設(shè)計(jì)1.求分段函數(shù)的值。x2+x_6x<0且 x #-32y=1x -5x+6 0 < x < 5H x ¥ 2及x # 3x2 -x -1其他I用if語(yǔ)句實(shí)現(xiàn)，分別輸出x=-5.0,-3.0,1.0,2.0,2.

9、5,3.0,5.0 時(shí)的y值。解：M文件如下：2.輸入一個(gè)百分制成績(jī)，要求輸出成績(jī)等級(jí)A、B、C、D、E。其中90分100分為A,80分89分為B, 79分79分為C, 60分69分為D, 60分以下為 E。要求：(1)分別用if語(yǔ)句和switch語(yǔ)句實(shí)現(xiàn)。(2)輸入百分制成績(jī)后要判斷該成績(jī)的合理性，對(duì)不合理的成績(jī)應(yīng)輸出出錯(cuò)信息。解：M文件如下 Editor 口：、Backup、我的交檔AILKB'.c.re. File Edi t Text Go Cell Tools Debug Desktop Window Helpw X:口妗旦|¥ 喑句牙“胸也|同目相三融回三 :唱唱

10、-To + 口一 *滋癖41 -clc口2 -sc or e=inputscore-'');3 -if score<=100fecore>=04 -switch fix (score/10)5 -case瓦10G -grades'T -case88 -grades，E”;9 -case710 -grad&=' C"11 -caseS12 -grade='D',13 -otherwise14 -grades'E' ；15 -end16 -grade17 - else18 -di=pC錯(cuò)誤士輸入的成績(jī)不是

11、百分制成績(jī))；Teri ptLn 19 Col 4score=88 grade =B score=123錯(cuò)誤：輸入的成績(jī)不是百分制成績(jī)3 .硅谷公司員工的工資計(jì)算方法如下：(1)工作時(shí)數(shù)超過(guò)120小時(shí)者，超過(guò)部分加發(fā)15%。(2)工作時(shí)數(shù)低于60小時(shí)者，扣發(fā)700元。(3)其余按每小時(shí)84元計(jì)發(fā)。試編程按輸入的工號(hào)和該號(hào)員工的工時(shí)數(shù)，計(jì)算應(yīng)發(fā)工資。解：M文件下4 .設(shè)計(jì)程序，完成兩位數(shù)的加、減、乘、除四則運(yùn)算，即產(chǎn)生兩個(gè)兩位隨機(jī)整數(shù)，再輸入 一個(gè)運(yùn)算符號(hào)，做相應(yīng)的運(yùn)算，并顯示相應(yīng)的結(jié)果。解：M文件如下；5 .建立5X6矩陣，要求輸出矩陣第 n行元素。當(dāng)n值超過(guò)矩陣的行數(shù)時(shí)，自動(dòng)轉(zhuǎn)為輸出 矩陣最

12、后一行元素，并給出出錯(cuò)信息。解：M文件如下:實(shí)驗(yàn)四循環(huán)結(jié)構(gòu)程序設(shè)計(jì)211111.根據(jù) = + +| +,求兀的近似值。當(dāng) n分別取100、1000、10000 612 22 32n2時(shí)，結(jié)果是多少？要求：分別用循環(huán)結(jié)構(gòu)和向量運(yùn)算（使用 sum函數(shù)）來(lái)實(shí)現(xiàn)。解：M文件如下:B 編輯器-Untitled9*SB®11. i2.根據(jù) y=1+| +,求:3 5 2n -1y<3時(shí)的最大n值。(2)與(1)的n值對(duì)應(yīng)的y值。解：M一文件如下：3.考慮以下迭代公式:abXn其中a、b為正的學(xué)數(shù)。迭代初值X0=1.0,迭代(1)編寫程序求迭代的結(jié)果，迭代的終止條件為|Xn+i-Xn|Wl0

13、-5,次數(shù)不超過(guò)500次。(2)如果迭代過(guò)程收斂于r,那么r的準(zhǔn)確值是一b 一 . b2 4a2,當(dāng)(a,b)的值取(1,1)、(8,3)、(10,0.1)時(shí)，分別對(duì)迭代結(jié)果和準(zhǔn)確值進(jìn)行比較。解：M文件如下：4.已知f1 =1n = 1j f2 =0n = 21f3 =1n = 3fn =fn 廠 2fnn 3求 fifi00 中：(1)最大值、最小值、各數(shù)之和。(2)正數(shù)、零、負(fù)數(shù)的個(gè)數(shù)。解：M一文件-|D| X|u r x1011121314151617IS1920clcf(l)=l;f(2)=D;f(3)= l；for n=4:100f(n)=f(n-2)+f(n-3);enddiwp

14、( max (f) = ? j num2str (naxCf) disp (r min(f)=3 num2str (min (f) disp ( sujuff) = , num2str (sun(f) disp ( c l-b j nun23tr (length(f ind (f>0) disp (' c2='j nun2str (1 erigth(f ind(£=0) disp ( c3=nun2str (length(find (f<0)蜒1到1 0口最大值圖fl到1QQ最小值 枇1到10。總和Kcl為正數(shù)的個(gè)數(shù)*cZ為。的個(gè)數(shù)骯3為f負(fù)數(shù)的個(gè)數(shù)Ldi

15、ptLn 9 Col 34OVR Editor - VntHled*File： 旦di t 工電:st Go! Cell1電池口/ 口ktgp Window Ulp酶 同三I冷燈三1.11.0以下是運(yùn)算結(jié)果:max=437763282635 min(f)=-899412113528 sum(f)=-742745601951 c1=49c2=2 c3=495.若兩個(gè)連續(xù)自然數(shù)的乘積減1是素?cái)?shù)，則稱這兩個(gè)邊疆自然數(shù)是親密數(shù)對(duì)，該素?cái)?shù)是親密素?cái)?shù)。例如，2X 3-1=5 ,由于5是素?cái)?shù)，所以2和3是親密數(shù)，5是親密素?cái)?shù)。求2,50 區(qū)間內(nèi)：(1)親密數(shù)對(duì)的對(duì)數(shù)。(2)與上述親密數(shù)對(duì)對(duì)應(yīng)的所有親密素?cái)?shù)之

16、和。解：M文件：運(yùn)算結(jié)果為:29s =23615實(shí)驗(yàn)五函數(shù)文件、實(shí)驗(yàn)?zāi)康? .理解函數(shù)文件的概念。2 .掌握定義和調(diào)用 MATLAB函數(shù)的方法。、實(shí)驗(yàn)內(nèi)容1. 定義一個(gè)函數(shù)文件，求給定復(fù)數(shù)的指數(shù)、對(duì)數(shù)、正弦和余弦，并在命令文件中調(diào)用 該函數(shù)文件。解：M文件如下：函數(shù)fushu.M文件：function e,l,s,c = fushu(z)%fushu 復(fù)數(shù)的指數(shù)，對(duì)數(shù)，正弦，余弦的計(jì)算%e復(fù)數(shù)的指數(shù)函數(shù)值%1 復(fù)數(shù)的對(duì)數(shù)函數(shù)值%s復(fù)數(shù)的正弦函數(shù)值%c復(fù)數(shù)的余弦函數(shù)值e=exp(z);l=log(z)；s=sin(z);c=cos(z);命令文件M :z=input('請(qǐng)輸入一個(gè)復(fù)數(shù)z=

17、');a,b,c,d=fushu(z)運(yùn)算結(jié)果如下：z=input('請(qǐng)輸入一個(gè)復(fù)數(shù)z=');a,b,c,d=fushu(z)請(qǐng)輸入一個(gè)復(fù)數(shù)z=1+ia =1.4687 + 2.2874ib =0.3466 + 0.7854c1.2985 + 0.6350id = 0.8337 - 0.9889i2. 一物理系統(tǒng)可用下列方程組來(lái)表示:從鍵盤輸入m一mi cosQ-m1-sin 001-ai 1f-0m1 sin 90cos 60a2mig0m2-sin 00Ni0-00-cos 61%m2 g1、m2和。的值,求 a1、a2、Ni和N2的值。其中g(shù)取9.8,輸入0時(shí)以角

18、度為單位。要求：定義一個(gè)求解線性方程組AX=B的函數(shù)文件，然后在命令文件中調(diào)用該函數(shù)文件。解：M文件函數(shù)fc.M文件：function X= fc(A,B)%fcfc是求解線性方程的函數(shù)%AA是未知矩陣的系數(shù)矩陣X=AB ;命令M文件：clc;m1=input('輸入 m1=');m2=input('輸入 m2=');theta=input('輸入 theta=');x=theta*pi/180;g=9.8;A=m1*cos(x) -ml -sin(x) 0m1*sin(x) 0 cos(x) 00 m2 -sin(x) 00 0 -cos(x)

19、 1;B=0;m1*g;0;m2*g;X=fc(A,B)運(yùn)算結(jié)果：輸入m1=1輸入m2=1輸入 theta=30X =7.84003.39486.789615.68003. 一個(gè)自然數(shù)是素?cái)?shù)，且它的數(shù)字位置經(jīng)過(guò)任意對(duì)換后仍為素?cái)?shù)。例如13是絕對(duì)素?cái)?shù)。試求所有兩位絕對(duì)素?cái)?shù)。要求：定義一個(gè)判斷素?cái)?shù)的函數(shù)文件。解：M文件：函數(shù)prime.m文件function p = prime(p)%輸入p的范圍，找出其中的素?cái)?shù)m=p(length(p);for i=2:sqrt(m)n=find(rem(p,i)=0&p=i);p(n)=；%將p中能被i整除，而卻不等于i的元素，即下標(biāo)為n的元素剔除，其

20、余的即為素?cái)?shù)endp；命令文件：clc;p=10:99;p=prime(p); %找出10到99內(nèi)的所有素?cái)?shù)p=10*rem(p,10)+(p-rem(p,10)/10;%將p素?cái)?shù)矩陣每個(gè)元素個(gè)位十位調(diào)換順序p=prime(p)%再對(duì)對(duì)換后的素?cái)?shù)矩陣找出所有的素?cái)?shù)運(yùn)算結(jié)果:p =113171137317379779一114.設(shè)f(x)=2+4,編寫一個(gè) MATLAB函數(shù)文件 僅巾，使得(x -2)0.1 (x -3)0.01調(diào)用f(x)時(shí)，x可用矩陣代入，得出的 f(x)為同階矩陣。解：函數(shù)僅.m文件：function f=僅(x)%fxfx求算x矩陣下的f(x)的函數(shù)值A(chǔ)=0.1+(x-2)

21、.A2;B=0.01+(x-3).A4;f=1./A+1./B;命令文件：clc;x=input('輸入矩陣 x=');f=fx(x)運(yùn)算結(jié)果:>> x=input('輸入矩陣 x='); f=fx(x)輸入矩陣x=7 2;12 5f =0.043710.99010.01010.1724當(dāng)f(n尸n+101n(n 2+5)時(shí),求y的值。(2)當(dāng) f(n)=1 X2+2 X3+3X4+.+n X (n+1)時(shí)，求 y 的值。解： 函數(shù)f.m文件： function f=f(x) f=x+10*1og(xA2+5);命令文件：clc;n1=input(&

22、#39;n1=');n2=input('n2=');n3=input('n3=');y1=f(n1);y2=f(n2);y3=f(n3);y=y1/(y2+y3)運(yùn)算結(jié)果如下:n1=40n2=30n3=20y =0.6390(2).函數(shù)g.m文件functions= g(n)for i=1:ng(i)=i*(i+1);ends=sum(g);命令文件：clc;n1=input('n1=' );n2=input('n2=');n3=input( y1=g(n1); y2=g(n2); y3=g(n3);'n3=

23、9;);y=y1/(y2+y3)運(yùn)算結(jié)果如下:n1=40n2=30n3=20 y =1.7662實(shí)驗(yàn)六高層繪圖操作、實(shí)驗(yàn)?zāi)康? .掌握繪制二維圖形的常用函數(shù)。2 .掌握繪制三維圖形的常用函數(shù)。3 .掌握繪制圖形的輔助操作。、實(shí)驗(yàn)內(nèi)容3sin xI .設(shè)y = 0.5 +3sin-x icosx,在x=02兀區(qū)間取101點(diǎn)，繪制函數(shù)的曲線。II 1 x2解：M文件如下:clc;x=linspace(0,2*pi,101);y=(0.5+3*sin(x)./(1+x，2); plot(x,y)運(yùn)行結(jié)果有:2.已知 y1 = x2, y2=cos(2 x), y3=y1Xy2,完成下列操作:(1)在

24、同一坐標(biāo)系下用不同的顏色和線型繪制三條曲線。(2)以子圖形式繪制三條曲線。(3)分別用條形圖、階梯圖、桿圖和填充圖繪制三條曲線。解：(1) M文件：clc;x=-pi:pi/100:pi;y1=x.A2;y2=cos(2*x);y3=y1.*y2;plot(x,y1,'b-',x,y2,'r:',x,y3,'k-')運(yùn)行結(jié)果:(2) M文件：clc;x=-pi:pi/100:pi;y1=x.A2;y2=cos(2*x);y3=y1.*y2;subplot(1,3,1);plot(x,y1, 'b-');title( 'y1

25、=xA2');subplot(1,3,2);plot(x,y2, 'r:');title( 'y2=cos(2x)');subplot(1,3,3);plot(x,y3, 'k-');title( 'y3=y1*y2');.運(yùn)行結(jié)果:(3) M文件:clc;x=-pi:pi/100:pi;y1=x.A2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,'b-',x,y2, 'rsubplot(2,2,2);bar(x,y1,'b');ti

26、tle( 'y1=xA2');subplot(2,2,3);bar(x,y2, 'r');title( 'y2=cos(2x)');subplot(2,2,4);bar(x,y3, 'k');title( 'y3=y1*y2');：',x,y3,'k-')；由上面的M文件，只要依次將“ bar”改為“ stairs "、“stem ”、“fill”，再適當(dāng)更改區(qū)間取的 點(diǎn)數(shù)，運(yùn)行程序即可，即有下面的結(jié)果：. J J i &工工門©收/ 一崎 口國(guó)Q 3.已知x&l

27、t; 0x 0x餐|2y =e _工 ln(x 、1 x2)2在-5W xw 5區(qū)間繪制函數(shù)曲線。 解：M文件：clc;x=-5:0.01:5; y=(x+sqrt(pi)/(exp(2).*(x<=0)+0.5*log(x+sqrt(1+x.A2).*(x>0); plot(x,y)運(yùn)行結(jié)果:由圖可看出，函數(shù)在零點(diǎn)不連續(xù)。4 .繪制極坐標(biāo)曲線p =asin(b+n 0 ),并分析參數(shù)a、b、n對(duì)曲線形狀的影響。解：M文件如下：clc;theta=0:pi/100:2*pi;a=input('輸入 a=');b=input('輸入 b=');n=in

28、put('輸入 n=');rho=a*sin(b+n*theta);polar(theta,rho,'m')采用控制變量法的辦法，固定兩個(gè)參數(shù)，變動(dòng)第三個(gè)參數(shù)觀察輸出圖象的變化。分析結(jié)果：由這8個(gè)圖知道，當(dāng)a,n固定時(shí)，圖形的形狀也就固定了，b只影響圖形的旋轉(zhuǎn)的角度；當(dāng)a,b固定時(shí)，n只影響圖形的扇形數(shù)，特別地，當(dāng) n是奇數(shù)時(shí)，扇葉數(shù)就是n,當(dāng)是偶數(shù)時(shí)，扇葉數(shù)則是2n個(gè)；當(dāng)b,n固定時(shí)，a影響的是圖形大小，特別地，當(dāng) a是整數(shù)時(shí)，圖形半徑大小就是a。5 .繪制函數(shù)的曲線圖和等高線。22x y_ 一 4-z = cos x cos ye其中x的21個(gè)值均勻分布-5

29、, 5范圍，y的31個(gè)值均勻分布在0, 10,要求使用 subplot(2,1,1)和subplot(2,1,2)將產(chǎn)生的曲面圖和等高線圖畫在同一個(gè)窗口上。解：M文件：clc;x=linspace(-5,5,21);y=linspace(0,10,31);x,y=meshgrid(x,y);z=cos(x).*cos(y).*exp(-sqrt(x.A2+y.A2)/4);subplot(2,1,1);surf(x,y,z);title('曲面圖');subplot(2,1,2);surfc(x,y,z);title('等高線圖');運(yùn)行結(jié)果:6 .繪制曲面圖形

30、，并進(jìn)行插值著色處理。x = cosscosty = cosssint八二八.3二0 - s - -, 0 -1 -22z =sin s解：M文件：clc;s=0:pi/100:pi/2;t=0:pi/100:3*pi/2;s,t=meshgrid(s,t);x=cos(s).*cos(t);y=cos(s).*sin；z=sin(s);subplot(2,2,1);mesh(x,y,z);title('未著色的圖形');subplot(2,2,2);surf(x,y,z);title('shading faceted (缺省)');subplot(2,2,3)

31、;surf(x,y,z);shading flat; title('shading flat');subplot(2,2,4);surf(x,y,z);shading interp;%插值著色title('shading interp');運(yùn)行結(jié)果有:-Ipl.xlFile Elit Vi ew Insert T»ols Desktop Window Help3二四痘圖因向國(guó)未著色處理的圖形I.超常,H!-G.1設(shè)！'JXK' J!'/'.,«.rl,11S9 感，潴shading intenp實(shí)驗(yàn)八數(shù)據(jù)處理與

32、多項(xiàng)式計(jì)算二、實(shí)驗(yàn)內(nèi)容1 .利用MATLAB提供的rand函數(shù)生成30000個(gè)符合均勻分布的隨機(jī)數(shù)，然后檢驗(yàn)隨 機(jī)數(shù)的性質(zhì)：(1)均值和標(biāo)準(zhǔn)方差。(2)最大元素和最小元素。(3)大于0.5的隨機(jī)數(shù)個(gè)數(shù)占總數(shù)的百分比。解：M文件：clc;x=rand(1,30000);mu=mean(x) %求這30000個(gè)均勻分布隨機(jī)數(shù)的平均值sig=std(x) %求其標(biāo)準(zhǔn)差o- iy=length(find(x>0.5); % 找出大于 0.5 數(shù)的個(gè)數(shù)p=y/30000%大于0.5的所占百分比運(yùn)行結(jié)果:mu =0.499488553231043 sig =0.288599933559786P =0

33、.4994000000000002 .將100個(gè)學(xué)生5門功課的成績(jī)存入矩陣P中，進(jìn)行如下處理：(1)分別求每門課的最高分、最低分及相應(yīng)學(xué)生序號(hào)。(2)分別求每門課的平均分和標(biāo)準(zhǔn)方差。3 3) 5門課總分的最高分、最低分及相應(yīng)學(xué)生序號(hào)。(4)將5門課總分按從大到小順序存入zcj中，相應(yīng)學(xué)生序號(hào)存入xsxh。提示：上機(jī)調(diào)試時(shí)，為避免輸入學(xué)生成績(jī)的麻煩，可用取值范圍在45,95之間的隨機(jī)矩陣來(lái)表示學(xué)生成績(jī)。解：M文件：clc;t=45+50*rand(100,5);P=fix(t); %生成100個(gè)學(xué)生5門功課成績(jī)x,l=max(P)%x為每門課最高分行向量，1為相應(yīng)學(xué)生序號(hào)y,k尸min(P)%y

34、為每門課最低分行向列,k為相應(yīng)學(xué)生序號(hào)mu=mean(P)%每門課的平均值行向量sig=std(P)%每門課的標(biāo)準(zhǔn)差行向量s=sum(P ,2) %5門課總分的列向量X,m=max(s)%5門課總分的最高分 X與相應(yīng)學(xué)生序號(hào)mY,n=min(s)%5門課總分的最低分 Y與相應(yīng)學(xué)生序號(hào)nzcj,xsxh=sort(s)%zcj為5門課總分從大到小排序，相應(yīng)學(xué)生序號(hào)xsxh運(yùn)行結(jié)果:3 .某氣象觀測(cè)得某日6:0018:00之間每隔2h的室內(nèi)外溫度（0C）如實(shí)驗(yàn)表1所示。實(shí)驗(yàn)表1室內(nèi)外溫度觀測(cè)結(jié)果（°C）h681012141618室內(nèi)溫度t118.020.022.025.030.028.0

35、24.0室外溫度t215.019.024.028.034.032.030.0試用三次樣條插值分另J求出該日室內(nèi)外6:3018:30之間每隔2h各點(diǎn)的近似溫度(0C)。解：M文件：clc;h=6:2:18;t1=18.0 20.0 22.0 25.0 30.0 28.0 24.0;t2=15.0 19.0 24.0 28.0 34.0 32.0 30.0;T1=interp1(h,t1,'spline')% 室內(nèi)的3次樣條插值溫度T2=interp1(h,t2,'spline')% 室外的3次樣條插值溫度運(yùn)行結(jié)果：T1 =Columns 1 through 340

36、.000000000000703Columns 4 through 654.000000000002885Column 752.000000000002444T2 =Columns 1 through 334.000000000000284Columns 4 through 660.00000000000451244.00000000000113064.00000000000588342.00000000000090272.00000000000940848.00000000000170560.00000000000451252.00000000000244468.000000000007503Column 764.0000000