MATLAB程序設(shè)計與應(yīng)用(劉衛(wèi)國編)課后實驗答案

1、MATLAB程序設(shè)計與應(yīng)用(劉衛(wèi)國編)課后實驗答案真驗一 MATLAB 運算基本1. 先供以下抒發(fā)式的值，而后隱示MATLAB 事情空間的利用情形并保留齊部變量。 (1) 0122sin851z e =+ (2) 221ln(1)2z x x =+，個中2120.455i x +=- (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,2.9,3.022a a e e a z a a -+=+=- (4) 2242011122123t t z t t t t t ，個中t =0:0.5:2.5 解： M 文件:z1=2*sin(85*pi/180)/(1+exp(2)x=2

2、 1+2*i;-0.45 5;z2=1/2*log(x+sqrt(1+x2)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運算了局：z1=2*sin(85*pi/180)/(1+exp(2)x=2 1+2*i;-.45 5;z2=1/2*log(x+sqrt(1+x2)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&

3、t0.2375z2 =0.7114 - 0.0253i 0.8968 + 0.3658i0.2139 + 0.9343i 1.1541 - 0.0044iz3 =Columns 1 through 40.7388 + 3.1416i 0.7696 + 3.1416i 0.7871 + 3.1416i 0.7913 + 3.1416iColumns 5 through 80.7822 + 3.1416i 0.7602 + 3.1416i 0.7254 + 3.1416i 0.6784 + 3.1416iColumns 9 through 120.6196 + 3.1416i 0.5496 + 3

4、.1416i 0.4688 + 3.1416i 0.3780 + 3.1416iColumns 13 through 160.2775 + 3.1416i 0.1680 + 3.1416i 0.0497 + 3.1416i -0.0771 + 3.1416iColumns 17 through 20-0.2124 + 3.1416i -0.3566 + 3.1416i -0.5104 + 3.1416i -0.6752 + 3.1416iColumns 21 through 24-0.8536 + 3.1416i -1.0497 + 3.1416i -1.2701 + 3.1416i -1.5

5、271 + 3.1416iColumns 25 through 28-1.8436 + 3.1416i -2.2727 + 3.1416i -2.9837 + 3.1416i -37.0245Columns 29 through 32-3.0017 -2.3085 -1.8971 -1.5978Columns 33 through 36-1.3575 -1.1531 -0.9723 -0.8083Columns 37 through 40-0.6567 -0.5151 -0.3819 -0.2561Columns 41 through 44-0.1374 -0.0255 0.0792 0.17

6、66Columns 45 through 480.2663 0.3478 0.4206 0.4841Columns 49 through 520.5379 0.5815 0.6145 0.6366Columns 53 through 560.6474 0.6470 0.63510.6119Columns 57 through 600.5777 0.5327 0.4774 0.4126Column 610.3388z4 =0 0.2500 0 1.2500 1.0000 2.25002. 已經(jīng)知：1234413134787,2033657327A B -=-供以下抒發(fā)式的值：(1) A+6*B

7、以及A-B+I （個中I 為單元矩陣）(2) A*B 以及A.*B(3) A3以及A.3(4) A/B 及BA(5) A,B以及A(1,3,:);B2解： 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.*BA3A.3A/BBAA,BA(1,3,:);B2運算了局：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.*BA3A.3A/BBAA,BA(1,3,:);B2ans =18 52 -1046 7

8、10521 53 49ans =12 31 -332 8 840 67 1ans =68 44 62309 -72 596154 -5 241ans =12 102 468 0 2619 -130 49ans =37226 233824 48604247370 149188 60076678688 454142 118820 ans =1728 39304 -6439304 343 65850327 274625 343 ans =16.4000 -13.6000 7.600035.8000 -76.2000 50.200067.0000 -134.0000 68.0000ans =109.4

9、000 -131.2000 322.8000-53.0000 85.0000 -171.0000-61.6000 89.8000 -186.2000 ans = 12 34 -4 1 3 -1 34 7 87 2 0 3 3 65 7 3 -2 7 ans = 12 34 -43 65 74 5 111 0 1920 -5 40 3. 設(shè)有矩陣A 以及B1234530166789101769,111213141502341617181920970212223242541311A B -=-(1) 供它們的乘積C。(2) 將矩陣C的左下角32子矩陣賦給D。(3) 檢察MATLAB事情空間的利用情

10、形。解：. 運算了局：E=(reshape(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 =93 150 77258 335 237423 520 397588 705 557753 890 717H =520 397705 557890 7174. 實現(xiàn)以下操縱：(1) 供100,999之間能被21整除了的數(shù)的個數(shù)。(2) 創(chuàng)建一個字符串背量，刪除了個中的年夜寫字母。解：(1) 了局：m=100:999;n=find(mod(m,21)=0);length(n)ans =43(2). 創(chuàng)建一個字

11、符串背量比方：ch=ABC123d4e56Fg9;則請求了局是：ch=ABC123d4e56Fg9;k=find(ch=A&chch(k)= 真驗2 MATLAB 矩陣剖析取處置1. 設(shè)有分塊矩陣33322322E R A O S =，個中E 、R 、O 、S 分手為單元矩陣、隨機矩陣、整矩陣以及對于角陣，試經(jīng)由過程數(shù)值盤算考證22E R RS A O S +=。 解: M 文件以下； 輸入了局：S = 1 00 2 ch = 123d4e56g9 A = 1.0000 0 0 0.5383 0.44270 1.0000 0 0.9961 0.10670 0 1.0000 0.0782 0.9

12、6190 0 0 1.0000 00 0 0 0 2.0000 a = 1.0000 0 0 1.0767 1.32800 1.0000 0 1.9923 0.32000 0 1.0000 0.1564 2.88570 0 0 1.0000 00 0 0 0 4.0000 ans = 0 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0由ans,以是22E R RS A O S += 2. 發(fā)生5階希我伯特矩陣H 以及5階帕斯卡矩陣P ，且供其止列式的值Hh 以及Hp 和它們的前提數(shù)Th 以及Tp ，判別哪一個矩陣功能更好。為何？解：M 文件以下：輸入了局

13、：H =1.0000 0.5000 0.3333 0.2500 0.20000.5000 0.3333 0.2500 0.2000 0.16670.3333 0.2500 0.2000 0.1667 0.14290.2500 0.2000 0.1667 0.1429 0.12500.2000 0.1667 0.1429 0.1250 0.1111 P =1 1 1 1 11 2 3 4 51 3 6 10 151 4 10 20 351 5 15 35 70Hh =3.7493e-012Hp =1Th =4.7661e+005Tp =8.5175e+003果為它們的前提數(shù)ThTp,以是pasc

14、al矩陣功能更好。3. 創(chuàng)建一個55矩陣，供它的止列式值、跡、秩以及范數(shù)。解：M文件以下：輸入了局為：A =17 24 1 8 1523 5 7 14 164 6 13 20 2210 12 19 21 311 18 25 2 9d =5070000t = 65 c1 = 6.8500 c2 = 5.4618 cinf = 6.8500 4. 已經(jīng)知2961820512885A -=-供A 的特性值及特性背量，并剖析其數(shù)教意思。 解： M 文件如圖： 輸入了局為： V = 0.7130 0.2803 0.2733-0.6084 -0.7867 0.8725 0.3487 0.5501 0.40

15、50 D = -25.3169 0 00 -10.5182 00 0 16.8351數(shù)教意思：V 的3個列背量是A 的特性背量，D 的主對于角線上3個是A 的特性值，出格的，V 的3個列背量分手是D 的3個特性值的特性背量。5. 上面是一個線性圓程組：1231112340.951110.673450.52111456x x x =(1) 供圓程的解。(2) 將圓程左邊背量元素b 3改成0.53再供解，并對比b 3的變動以及解的相對于變動。(3) 盤算系數(shù)矩陣A 的前提數(shù)并剖析論斷。解： M 文件以下：輸入了局：X =1.20000.60000.6000X2 =1.20000.60000.600

16、0C =1.3533e+003由了局，X以及X2的值同樣，那暗示b的巨大變動對于圓程解也影響較小，而A的前提數(shù)算患上較小，以是數(shù)值不亂性較好，A是較好的矩陣。6. 創(chuàng)建A矩陣，試對比sqrtm(A)以及sqrt(A)，剖析它們的區(qū)分。解：M文件以下：運轉(zhuǎn)了局有：A =16 6 1820 5 129 8 5b1 =3.8891 -0.1102 3.21033.2917 2.1436 0.36980.3855 2.0760 1.7305b2 =4.0000 2.4495 4.24264.4721 2.2361 3.46413.0000 2.8284 2.2361b =16.0000 6.0000

17、18.000020.0000 5.0000 12.00009.0000 8.0000 5.0000剖析了局知：sqrtm(A)是相似A的數(shù)值仄圓根（那可由b1*b1=A的了局瞧出），而sqrt(A)則是對于A中的每一個元素開根號，兩則區(qū)分便正在于此。真驗3 取舍布局步伐計劃一、真驗?zāi)繕?biāo)1. 把握創(chuàng)建以及實行M 文件的圓法。2. 把握使用if 語句真現(xiàn)取舍布局的圓法。3. 把握使用switch 語句真現(xiàn)多分收取舍布局的圓法。4. 把握try 語句的利用。2、真驗內(nèi)容1. 供分段函數(shù)的值。2226035605231x x x x y x x x x x x x +-且且及其余用if 語句真現(xiàn)，分手

18、輸入x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0時的y 值。 解：M 文件以下： 運算了局有： f(-5)y =14 f(-3)y =11 f(1)y =2 f(2)y =1 f(2.5)y =-0.2500 f(3)y =5 f(5)y =192. 輸出一個百分造成就，請求輸入成就品級A、B、C、D、E。個中90分100分為A，80分89分為B，79分79分為C，60分69分為D，60分下列為E。請求：(1) 分手用if語句以及switch語句真現(xiàn)。(2) 輸出百分造成就后要判別該成就的開感性，對于沒有開理的成就應(yīng)輸入墮落疑息。解：M文件以下試算了局：score=88gra

19、de =Bscore=123同伴：輸出的成就沒有是百分造成就3. 硅谷公司員工的人為盤算圓法以下：(1) 事情時數(shù)凌駕120小時者，凌駕全體減收15%。(2) 事情時數(shù)低于60小時者，扣收700元。(3) 其他按每一小時84元計收。試編程按輸出的工號以及該號員工的工時數(shù)，盤算應(yīng)收人為。解：M文件下4. 計劃步伐，實現(xiàn)兩位數(shù)的減、加、乘、除了4則運算，即發(fā)生兩個兩位隨機整數(shù)，再輸出一個運算標(biāo)記，做響應(yīng)的運算，并隱示響應(yīng)的了局。解：M文件以下；運算了局例：a =38b =33輸出一個運算符: c =falsea =92b =40輸出一個運算符:+c =1325. 創(chuàng)建56矩陣，請求輸入矩陣第n止元

20、素。當(dāng)n值凌駕矩陣的止數(shù)時，主動轉(zhuǎn)為輸入矩陣最初一止元素，并給出墮落疑息。解：M文件以下：運算了局以下：輸出一個5止6列矩陣A=1 2 3 4 5 5;2 3 4 5 7 6;2 2 2 2 2 3;11 2 3 9 7 3;2 3 4 5 6 7 輸出一正整數(shù)n=411 2 3 9 7 3輸出一個5止6列矩陣A=1 2 3 4 5 5;2 3 4 5 7 6;2 2 2 2 2 3;11 2 3 9 7 3;2 3 4 5 6 7 輸出一正整數(shù)n=62 3 4 5 6 7ans =Error using = dispToo many input arguments.真驗4輪回布局步伐計劃一、

21、真驗?zāi)繕?biāo)1. 把握使用for語句真現(xiàn)輪回布局的圓法。2. 把握使用while語句真現(xiàn)輪回布局的圓法。3. 生悉使用背量運算去取代輪回操縱的圓法。2、真驗內(nèi)容1. 依據(jù)2222211116123n=+，供的遠(yuǎn)似值。當(dāng)n分手與100、1000、10000時，了局是幾？請求：分手用輪回布局以及背量運算（利用sum函數(shù)）去真現(xiàn)。解：M文件以下：運轉(zhuǎn)了局以下：K %輪回布局盤算pi值y=0;n=input(n=);for i=1:ny=y+1/i/i;endpi=sqrt(6*y)n=100pi =3.1321n=1000pi =3.1406n=10000pi =3.1415%背量圓法盤算Pi值n=in

22、put(n=);i=1./(1:n).2;s=sum(i);pi=sqrt(6*s)n=100pi =3.1321n=1000pi =3.1406n=10000pi =3.14152. 依據(jù)11113521yn=+-，供：(1) y(2) 取(1)的n值對于應(yīng)的y值。解：M文件以下：運轉(zhuǎn)了局以下：K y=0;n=0;while yn=n+1;y=y+1/(2*n-1);endynif y3n=n-1;endny =3.0033 n = 57 n = 56 3. 思索下列迭代公式：1n na xb x +=+ 個中a 、b 為正的教數(shù)。(1) 編寫步伐供迭代的了局，迭代的末行前提為|x n+1-

23、x n |10-5，迭代初值x 0=1.0，迭代次數(shù)沒有凌駕500次。(2) 假如迭代歷程支斂于r ，那末r 的正確值是242b b a -+，當(dāng)(a,b)的值與(1,1)、(8,3)、(10,0.1)時，分手對于迭代了局以及正確值舉行對比。解：M 文件以下：運算了局以下；請輸出負(fù)數(shù)a=1請輸出負(fù)數(shù)b=1x =0.6180r =0.6180 -4.7016 r =0.6180 -1.6180 s =-0.0000 -2.2361請輸出負(fù)數(shù)a=8請輸出負(fù)數(shù)b=3x =1.7016r =1.7016 -1.6180 r =1.7016 -4.7016 s =0.0 -6.4031請輸出負(fù)數(shù)a=10

24、請輸出負(fù)數(shù)b=0.1x =3.1127r =3.1127 -4.7016 r =3.1127 -3.2127 s =-0.0000 -6.32544. 已經(jīng)知12312311021323n n n n f n f n f n f f f f n -=-+供f 1f 100中：(1) 最年夜值、最小值、各數(shù)之以及。(2) 負(fù)數(shù)、整、背數(shù)的個數(shù)。 解：M 文件下列是運算了局：max(f)=437763282635 min(f)=-899412113528 sum(f)=-742745601951 c1=49c2=2c3=49 5. 若兩個一連做作數(shù)的乘積加1是素數(shù)，則稱那兩個內(nèi)地做作數(shù)是親切數(shù)對于

25、，該素數(shù)是親切素數(shù)。比方，23-1=5，因為5是素數(shù)，以是2以及3是親切數(shù)，5是親切素數(shù)。供2,50區(qū)間內(nèi)：(1) 親切數(shù)對于的對于數(shù)。(2) 取上述親切數(shù)對于對于應(yīng)的一切親切素數(shù)之以及。解：M文件：運算了局為：j =29s =23615真驗5函數(shù)文件一、真驗?zāi)繕?biāo)1. 了解函數(shù)文件的觀點。2. 把握界說以及挪用MATLAB函數(shù)的圓法。2、真驗內(nèi)容1. 界說一個函數(shù)文件，供給定單數(shù)的指數(shù)、對于數(shù)、正弦以及余弦，并正在下令文件中挪用該函數(shù)文件。解：M文件以下：函數(shù)fushu.M文件：function e,l,s,c = fushu(z)%fushu 單數(shù)的指數(shù)，對于數(shù)，正弦，余弦的盤算%e 單數(shù)的

26、指數(shù)函數(shù)值%l 單數(shù)的對于數(shù)函數(shù)值%s 單數(shù)的正弦函數(shù)值%c 單數(shù)的余弦函數(shù)值e=exp(z);l=log(z);s=sin(z);c=cos(z);下令文件M：z=input(請輸出一個單數(shù)z=);a,b,c,d=fushu(z)運算了局以下：z=input(請輸出一個單數(shù)z=);a,b,c,d=fushu(z)請輸出一個單數(shù)z=1+ia =1.4687 +2.2874i b = 0.3466 + 0.7854i c = 1.2985 + 0.6350i d = 0.8337 - 0.9889i 2. 一物理體系可用以下圓程組去暗示：11121112220cos sin 0sin 0cos

27、000sin 000cos 1a m m a m g m N m N m g -=-從鍵盤輸出m 1、m 2以及的值，供a 1、a 2、N 1以及N 2的值。個中g(shù) 與9.8，輸出時以角度為單元。請求：界說一個供解線性圓程組AX=B 的函數(shù)文件，而后正在下令文件中挪用該函數(shù)文件。解： M 文件 函數(shù)fc.M 文件:function X= fc(A,B)%fc fc 是供解線性圓程的函數(shù)%A A 是已知矩陣的系數(shù)矩陣X=AB ； 下令M 文件：clc;m1=input(輸出m1=);m2=input(輸出m2=);theta=input(輸出theta=);x=theta*pi/180;g=9.

28、8;A=m1*cos(x) -m1 -sin(x) 0m1*sin(x) 0 cos(x) 00 m2 -sin(x) 00 0 -cos(x) 1;B=0;m1*g;0;m2*g;X=fc(A,B)運算了局：輸出m1=1輸出m2=1輸出theta=30X =7.84003.39486.789615.68003. 一個做作數(shù)是素數(shù)，且它的數(shù)字地位經(jīng)由恣意對于換后仍為素數(shù)。比方13是盡對于素數(shù)。試供一切兩位盡對于素數(shù)。請求：界說一個判別素數(shù)的函數(shù)文件。解：M文件：函數(shù)prime.m文件function p = prime(p)% 輸出p的局限，尋出個中的素數(shù)m=p(length(p);for i

29、=2:sqrt(m)n=find(rem(p,i)=0&p=i);p(n)=;%將p中能被i整除了，而卻沒有即是i的元素，即下標(biāo)為n的元素剔除了，其他的即為素數(shù)endp;下令文件：clc;p=10:99;p=prime(p); %尋出10到99內(nèi)的一切素數(shù)p=10*rem(p,10)+(p-rem(p,10)/10;%將p 素數(shù)矩陣每一個元素個位10位互換逆序p=prime(p)%再對于對于換后的素數(shù)矩陣尋出一切的素數(shù)運算了局：p = 11 31 71 13 73 17 37 97 79 4. 設(shè)2411()(2)0.1(3)0.01f x x x =+-+-+，編寫一個MATLAB 函數(shù)文件

30、fx.m ，使患上挪用f(x)時，x 可用矩陣代進，患上出的f(x)為同階矩陣。解：函數(shù)fx.m 文件：function f= fx(x)%fx fx 供算x 矩陣下的f(x)的函數(shù)值A(chǔ)=0.1+(x-2).2;B=0.01+(x-3).4;f=1./A+1./B; 下令文件：clc;x=input(輸出矩陣x=); f=fx(x)運算了局： x=input(輸出矩陣x=);f=fx(x)輸出矩陣x=7 2;12 5 f = 0.0437 10.99010.0101 0.1724 5. 已經(jīng)知(40)(30)(20)f y f f =+ (1) 當(dāng)f(n)=n+10ln(n 2+5)時，供y

31、的值。(2) 當(dāng)f(n)=12+23+34+.+n (n+1)時，供y 的值。解：(1)函數(shù)f.m文件: function f=f(x)f=x+10*log(x2+5);下令文件：clc;n1=input(n1=);n2=input(n2=);n3=input(n3=);y1=f(n1);y2=f(n2);y3=f(n3);y=y1/(y2+y3)運算了局以下：n1=40n2=30n3=20y =0.6390(2).函數(shù)g.m文件function s= g(n) for i=1:ng(i)=i*(i+1);ends=sum(g);下令文件：clc;n1=input(n1=);n2=input(

32、n2=);n3=input(n3=);y1=g(n1);y2=g(n2);y3=g(n3);y=y1/(y2+y3) 運算了局以下：n1=40n2=30n3=20 y = 1.7662 真驗6 下層畫圖操縱一、真驗?zāi)繕?biāo)1. 把握畫造2維圖形的經(jīng)常使用函數(shù)。2. 把握畫造3維圖形的經(jīng)常使用函數(shù)。3. 把握畫造圖形的幫助操縱。2、真驗內(nèi)容1. 設(shè)23sin 0.5cos 1x y x x =+，正在x=02區(qū)間與101面，畫造函數(shù)的直線。 解：M 文件以下：clc;x=linspace(0,2*pi,101);y=(0.5+3*sin(x)./(1+x.2); plot(x,y)運轉(zhuǎn)了局有：2.

33、已經(jīng)知y1=x2，y2=cos(2x)，y3=y1y2，實現(xiàn)以下操縱：(1) 正在統(tǒng)一坐標(biāo)系下用沒有同的色彩以及線型畫造3條直線。(2) 以子圖情勢畫造3條直線。(3) 分手用條形圖、門路圖、桿圖以及挖充圖畫造3條直線。解：（1）M文件：clc;x=-pi:pi/100:pi;y1=x.2;y2=cos(2*x);y3=y1.*y2;plot(x,y1,b-,x,y2,r:,x,y3,k-)運轉(zhuǎn)了局：（2）M文件：clc;x=-pi:pi/100:pi;y1=x.2;y2=cos(2*x);y3=y1.*y2;subplot(1,3,1);plot(x,y1,b-);title(y1=x2);

34、 subplot(1,3,2);plot(x,y2,r:);title(y2=cos(2x); subplot(1,3,3);plot(x,y3,k-); title(y3=y1*y2);.運轉(zhuǎn)了局：（3）M文件：clc;x=-pi:pi/100:pi;y1=x.2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,b-,x,y2,r:,x,y3,k-);subplot(2,2,2);bar(x,y1,b);title(y1=x2);subplot(2,2,3);bar(x,y2,r);title(y2=cos(2x);subplot(2,2,4);

35、bar(x,y3,k);title(y3=y1*y2);由下面的M文件，只有挨次將“bar”改成“stairs”、“stem”、“fill”,再得當(dāng)變動區(qū)間與的面數(shù)，運轉(zhuǎn)步伐便可，即有上面的了局： 3. 已經(jīng)知2201ln(1)02x x e y x x x +=+ 正在-5x 5區(qū)間畫造函數(shù)直線。 解：M 文件：clc;x=-5:0.01:5;y=(x+sqrt(pi)/(exp(2).*(x運轉(zhuǎn)了局：由圖可瞧出，函數(shù)正在整面沒有一連。4. 畫造極坐標(biāo)直線=asin(b+n)，并剖析參數(shù)a、b、n對于直線外形的影響。解：M文件以下：clc;theta=0:pi/100:2*pi;a=inpu

36、t(輸出a=);b=input(輸出b=);n=input(輸出n=);rho=a*sin(b+n*theta);polar(theta,rho,m)接納把持變量法的舉措，流動兩個參數(shù)，變化第3個參數(shù)不雅察輸入圖像的變動。 剖析了局：由那8個圖明白，當(dāng)a,n 流動時，圖形的外形也便流動了，b 只影響圖形的扭轉(zhuǎn)的角度;當(dāng)a,b 流動時，n 只影響圖形的扇形數(shù)，出格天，當(dāng)n 是偶數(shù)時，扇葉數(shù)便是n,當(dāng)是奇數(shù)時，扇葉數(shù)則是2n 個;當(dāng)b,n 流動時，a 影響的是圖形年夜小，出格天，當(dāng)a 是整數(shù)時，圖形半徑年夜小便是a 。5. 畫造函數(shù)的直線圖以及等下線。224cos cos x y z x ye +

37、-=個中x 的21個值勻稱散布-5，5局限，y 的31個值勻稱散布正在0，10，請求利用subplot(2,1,1)以及subplot(2,1,2)將發(fā)生的直里圖以及等下線丹青正在統(tǒng)一個窗心上。 解：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.2+y.2)/4);subplot(2,1,1);surf(x,y,z);title(直里圖);subplot(2,1,2);surfc(x,y,z);title(等下線圖); 運轉(zhuǎn)了局： 6. 畫造直里

38、圖形，并舉行插值著色處置。cos cos 3cos sin 0,022sin x s t y s t s t z 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(t);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);surf(x,y,z);shading fl

39、at;title(shading flat);subplot(2,2,4);surf(x,y,z);shading interp; title(shading interp);運轉(zhuǎn)了局有：真驗7低層畫圖操縱2、真驗內(nèi)容1. 創(chuàng)建一個圖形窗心，使之后臺色彩為白色，并正在窗心上保存本有的菜單項，并且正在按下鼠標(biāo)器的左鍵以后隱示出Left Button Pressed字樣。解：M文件以下：clc;hf=figure(color,1 0 0,.WindowButtonDownFcn,disp(Left Button Pressed.);運轉(zhuǎn)了局：左擊鼠標(biāo)后：2. 先使用默許屬性畫造直線y=x2e2x，

40、而后經(jīng)由過程圖形句柄操縱去扭轉(zhuǎn)直線的色彩、線型以及線寬，并使用文件對于象給直線加減筆墨標(biāo)注。解：M文件：clc;x=-2:0.01:2;y=x.2.*exp(2*x);h=plot(x,y);set(h,color,0.4,0.2,0.5,linestyle,-,.linewidth,2);text(1.5,1.52*exp(2*1.5),leftarrow x2exp(2x),fontsize,9);運轉(zhuǎn)了局:3. 使用直里對于象畫造直里v(x,t)=10e-0.01xsin(2000t-0.2x+)。解：M文件：clc;x=0:0.1:2*pi;x,t=meshgrid(x);v=10*e

41、xp(-0.01*x).*sin(2000*pi*t-0.2*x+pi);axes(view,-37,30);hs=surface(x,t,v,facecolor,.0.2,0.3,0.3,edgecolor,flat);grid on;xlabel(x-axis); ylabel(y-axis);zlabel(z-axis);title(mesh-surf);pause %按恣意鍵持續(xù)set(hs,FaceColor,flat); text(0,0,0,直里);運轉(zhuǎn)了局：按恣意鍵持續(xù)：4. 以恣意地位子圖情勢畫造出正弦、余弦、正切以及余切函數(shù)直線。5. 死成一個圓柱體，并舉行光照以及材量處置

42、。解：M文件：x,y,z=cylinder(3,500); %cylinder是死成柱體的函數(shù)surf(x,y,z);title(圓柱體的光照以及質(zhì)料處置);Xlabel(X-axis);Ylabel(Y-axis);Zlabel(Z-axis);axis(-5,5,-5,5,0,1)grid off;light(Color,r,Position,-4,0,0,style,infinite);shading interp;material shiny;view(0,10);lighting phong;axis off;運轉(zhuǎn)了局：真驗8數(shù)據(jù)處置取多項式盤算一、真驗?zāi)繕?biāo)1. 把握數(shù)據(jù)統(tǒng)計以及剖析

43、的圓法。2. 把握數(shù)值插值取直線擬開的圓法及其使用。3. 把握多項式的經(jīng)常使用運算。2、真驗內(nèi)容1. 使用MATLAB供應(yīng)的rand函數(shù)死成30000個切合勻稱散布的隨機數(shù)，而后查驗隨機數(shù)的性子：(1) 均值以及尺度圓好。(2) 最年夜元素以及最小元素。(3) 年夜于0.5的隨機數(shù)個數(shù)占總數(shù)的百分比。解：M文件:clc;x=rand(1,30000);mu=mean(x) %供那30000個勻稱散布隨機數(shù)的仄均值sig=std(x) %供其尺度好1y=length(find(x0.5); %尋出年夜于0.5數(shù)的個數(shù)p=y/30000 %年夜于0.5的所占百分比運轉(zhuǎn)了局：mu =0.499488

44、553231043sig =0.288599933559786p =0.4994000000000002. 將100個教死5門作業(yè)的成就存進矩陣P中，舉行以下處置：(1) 分手供每一門課的最下分、最低分及響應(yīng)教死序號。(2) 分手供每一門課的仄均分以及尺度圓好。(3) 5門課總分的最下分、最低分及響應(yīng)教死序號。(4) 將5門課總分按從年夜到小逆序存進zcj中，響應(yīng)教死序號存進xsxh。提醒：上機調(diào)試時，為躲免輸出教死成就的貧苦，可用與值局限正在45,95之間的隨機矩陣去暗示教死成就。解：M文件：clc;t=45+50*rand(100,5);P=fix(t); %死成100個教死5門作業(yè)成就x

45、,l=max(P)%x為每一門課最下分止背量,l為響應(yīng)教死序號y,k=min(P)%y為每一門課最低分止背列,k為響應(yīng)教死序號mu=mean(P) %每一門課的仄均值止背量sig=std(P) %每一門課的尺度好止背量s=sum(P,2) %5門課總分的列背量X,m=max(s)%5門課總分的最下分X取響應(yīng)教死序號mY,n=min(s)%5門課總分的最低分Y取響應(yīng)教死序號nzcj,xsxh=sort(s)%zcj為5門課總分從年夜到小排序，響應(yīng)教死序號xsxh運轉(zhuǎn)了局：3. 某天氣不雅測患上某日6:0018:00之間每一隔2h的室表里溫度（0C）照實驗表1所示。真驗表1 室表里溫度不雅測了局（

46、0C）光陰h 6 8 10 12 14 16 18 室內(nèi)溫度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 試用3次樣條插值分手供出該日室表里6:3018:30之間每一隔2h各面的遠(yuǎ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,

47、spline)%室中的3次樣條插值溫度運轉(zhuǎn)了局：T1 =Columns 1 through 340.000000000000703 44.000000000001130 48.000000000001705Columns 4 through 654.000000000002885 64.000000000005883 60.000000000004512Column 752.000000000002444T2 =Columns 1 through 334.000000000000284 42.000000000000902 52.000000000002444Columns 4 through

48、 660.000000000004512 72.000000000009408 68.000000000007503Column 764.0000000000058834. 已經(jīng)知lgx正在1,101區(qū)間10個整數(shù)采樣面的函數(shù)值照實驗表2所示。真驗表2 lgx正在10個采樣面的函數(shù)值x 1 11 21 31 41 51 61 71 81 91 101lgx 0 1.0414 1.3222 1.4914 1.6128 1.7076 1.7853 1.8513 1.9085 1.9510 2.0043試供lgx的5次擬開多項式p(x)，并畫造出lgx以及p(x)正在1,101區(qū)間的函數(shù)直線。解：M文件：x=1:10:101;y=lg10(x);P=polyfit(x,y,5)y1=polyval(P,x);plot(x,y,:o,x,y1,-*)運轉(zhuǎn)了局：Warning: Polynomial is badly conditioned. Add points with distinct Xvalues, reduce the degree of the polynomial, or try centering