Matlab上機(jī)實(shí)驗(yàn)答案

1、Matlab上機(jī)實(shí)驗(yàn)答案實(shí)驗(yàn)一 MATLAB運(yùn)算基礎(chǔ)1. 先求下列表達(dá)式的值，然后顯示MATLAB工作空間的使用情況并保存全部變量。(1) >> z1=2*sin(85*pi/180)/(1+exp(2)z1 = 0.2375(2) ，其中>> x=2 1+2i;-0.45 5;>> z2=1/2*log(x+sqrt(1+x2)z2 = 0.7114 - 0.0253i 0.8968 + 0.3658i 0.2139 + 0.9343i 1.1541 - 0.0044i(3) >> a=-3.0:0.1:3.0;>> z3=(exp

2、(0.3.*a)-exp(-0.3.*a)./2.*sin(a+0.3)+log(0.3+a)./2)(>> z33=(exp(0.3*a)-exp(-0.3*a)/2.*sin(a+0.3)+log(0.3+a)/2)可以驗(yàn)證z3=z33，是否都為1）z3 = Columns 1 through 5 0.7388 + 3.1416i 0.7696 + 3.1416i 0.7871 + 3.1416i 0.7913 + 3.1416i 0.7822 + 3.1416i Columns 6 through 10 0.7602 + 3.1416i 0.7254 + 3.1416i 0.

3、6784 + 3.1416i 0.6196 + 3.1416i 0.5496 + 3.1416i Columns 11 through 15 0.4688 + 3.1416i 0.3780 + 3.1416i 0.2775 + 3.1416i 0.1680 + 3.1416i 0.0497 + 3.1416i Columns 16 through 20 -0.0771 + 3.1416i -0.2124 + 3.1416i -0.3566 + 3.1416i -0.5104 + 3.1416i -0.6752 + 3.1416i Columns 21 through 25 -0.8536 +

4、3.1416i -1.0497 + 3.1416i -1.2701 + 3.1416i -1.5271 + 3.1416i -1.8436 + 3.1416i Columns 26 through 30 -2.2727 + 3.1416i -2.9837 + 3.1416i -37.0245 + 0.0000i -3.0017 + 0.0000i -2.3085 + 0.0000i Columns 31 through 35 -1.8971 + 0.0000i -1.5978 + 0.0000i -1.3575 + 0.0000i -1.1531 + 0.0000i -0.9723 + 0.0

5、000i Columns 36 through 40 -0.8083 + 0.0000i -0.6567 + 0.0000i -0.5151 + 0.0000i -0.3819 + 0.0000i -0.2561 + 0.0000i Columns 41 through 45 -0.1374 + 0.0000i -0.0255 + 0.0000i 0.0792 + 0.0000i 0.1766 + 0.0000i 0.2663 + 0.0000i Columns 46 through 50 0.3478 + 0.0000i 0.4206 + 0.0000i 0.4841 + 0.0000i 0

6、.5379 + 0.0000i 0.5815 + 0.0000i Columns 51 through 55 0.6145 + 0.0000i 0.6366 + 0.0000i 0.6474 + 0.0000i 0.6470 + 0.0000i 0.6351 + 0.0000i Columns 56 through 60 0.6119 + 0.0000i 0.5777 + 0.0000i 0.5327 + 0.0000i 0.4774 + 0.0000i 0.4126 + 0.0000i Column 61 0.3388 + 0.0000i(4) ，其中t=0:0.5:2.5>>

7、t=0:0.5:2.5;>> z4=(t>=0&t<1).*(t.2)+(t>=1&t<2).*(t.2-1)+(t>=2&t<3).*(t.2-2.*t+1)z4 = 0 0.2500 0 1.2500 1.0000 2.25002. 已知：求下列表達(dá)式的值：(1) A+6*B和A-B+I（其中I為單位矩陣）(2) A*B和A.*B(3) A3和A.3(4) A/B及BA(5) A,B和A(1,3,:);B2>> A=12 34 -4;34 7 87;3 65 7;>> B=1 3 -1; 2 0

8、 3;3 -2 7;>> A+6*Bans = 18 52 -10 46 7 105 21 53 49>> I=eye(3);>> A-B+Ians = 12 31 -3 32 8 84 0 67 1(2)>> A*Bans = 68 44 62 309 -72 596 154 -5 241>> A.*Bans = 12 102 4 68 0 261 9 -130 49(3)>> A3ans = 37226 233824 48604 247370 149188 600766 78688 454142 118820>&

9、gt; A.3ans = 1728 39304 -64 39304 343 658503 27 274625 343(4)>> A/Bans = 16.4000 -13.6000 7.6000 35.8000 -76.2000 50.2000 67.0000 -134.0000 68.0000>> BAans = 109.4000 -131.2000 322.8000 -53.0000 85.0000 -171.0000 -61.6000 89.8000 -186.2000(5)>> A,Bans = 12 34 -4 1 3 -1 34 7 87 2 0

10、3 3 65 7 3 -2 7>> A(1,3,:);B2ans = 12 34 -4 3 65 7 4 5 1 11 0 19 20 -5 403. 設(shè)有矩陣A和B(1) 求它們的乘積C。(2) 將矩陣C的右下角3×2子矩陣賦給D。(3) 查看MATLAB工作空間的使用情況。>> A=(reshape(1:25,5,5)'>> B=3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11;>> C=A*BC = 93 150 77 258 335 237 423 520 397 588 705 557 753

11、890 717>> D=C(3:5,2:3)D = 520 397 705 557 890 717>> whos Name Size Bytes Class Attributes A 5x5 200 double B 5x3 120 double C 5x3 120 double D 3x2 48 double 4. 完成下列操作：(1) 求100,999之間能被21整除的數(shù)的個(gè)數(shù)。(2) 建立一個(gè)字符串向量，刪除其中的大寫字母。>> n=100:999;>> l=find(rem(n,21)=0); >> length(l)ans

12、=43>> ch='aegbBOIEG0je23RGnc'>> wz=find(ch>='A'&ch<='Z');>> ch(wz)=ch =aegb0je23nc實(shí)驗(yàn)二 MATLAB矩陣分析與處理1. 設(shè)有分塊矩陣，其中E、R、O、S分別為單位矩陣、隨機(jī)矩陣、零矩陣和對(duì)角陣，試通過(guò)數(shù)值計(jì)算驗(yàn)證。>> E=eye(3);>> R=rand(3,2);>> O=zeros(2,3);>> S=diag(1:2);>> A=E,R;O,

13、SA = 1.0000 0 0 0.8147 0.9134 0 1.0000 0 0.9058 0.6324 0 0 1.0000 0.1270 0.0975 0 0 0 1.0000 0 0 0 0 0 2.0000>> A2ans = 1.0000 0 0 1.6294 2.7401 0 1.0000 0 1.8116 1.8971 0 0 1.0000 0.2540 0.2926 0 0 0 1.0000 0 0 0 0 0 4.0000>> E,R+R*S;O,S2ans = 1.0000 0 0 1.6294 2.7401 0 1.0000 0 1.8116

14、1.8971 0 0 1.0000 0.2540 0.2926 0 0 0 1.0000 0 0 0 0 0 4.0000>> A2=E,R+R*S;O,S2ans = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1通過(guò)驗(yàn)證，矩陣成立。2. 產(chǎn)生5階希爾伯特矩陣H和5階帕斯卡矩陣P，且求其行列式的值Hh和Hp以及它們的條件數(shù)Th和Tp，判斷哪個(gè)矩陣性能更好。為什么？>> H=hilb(5)H = 1.0000 0.5000 0.3333 0.2500 0.2000 0.5000 0.3333 0.2500 0.2000

15、 0.1667 0.3333 0.2500 0.2000 0.1667 0.1429 0.2500 0.2000 0.1667 0.1429 0.1250 0.2000 0.1667 0.1429 0.1250 0.1111>> P=pascal(5)P = 1 1 1 1 1 1 2 3 4 5 1 3 6 10 15 1 4 10 20 35 1 5 15 35 70>> Hh=det(H)Hh = 3.7493e-12>> Hp=det(P)Hp = 1>> Th=cond(H)Th = 4.7661e+05>> Tp=cond

16、(P)Tp = 8.5175e+03答：5階帕斯卡矩陣P的性能好。矩陣的性能是由條件數(shù)決定的，條件數(shù)越接近于1其性能就越好。由上機(jī)操作求得Th=4.7661e+005，Tp=8.5175e+003。Tp的值更接近于1則其性能要好。所以5階帕斯卡矩陣P的性能好。3. 建立一個(gè)5×5矩陣，求它的行列式值、跡、秩和范數(shù)。>> A=rand(5)A = 0.2785 0.9706 0.4218 0.0357 0.7431 0.5469 0.9572 0.9157 0.8491 0.3922 0.9575 0.4854 0.7922 0.9340 0.6555 0.9649 0.8

17、003 0.9595 0.6787 0.1712 0.1576 0.1419 0.6557 0.7577 0.7060>> det(A)ans = -0.1322>> trace(A)ans = 3.4127>> rank(A)ans = 5>> norm(A)ans = 3.27594. 已知求A的特征值及特征向量，并分析其數(shù)學(xué)意義。>> A=-29,6,18;20,5,12;-8,8,5A = -29 6 18 20 5 12 -8 8 5>> V,D=eig(A)V = 0.7130 0.2803 0.2733 -0

18、.6084 -0.7867 0.8725 0.3487 0.5501 0.4050D = -25.3169 0 0 0 -10.5182 0 0 0 16.8351在數(shù)學(xué)上，特別是線性代數(shù)中，對(duì)于一個(gè)給定的線性變換，它的特征向量（本征向量或稱正規(guī)正交向量）是這樣一個(gè)非零的向量v：當(dāng)v 經(jīng)過(guò)這個(gè)線性變換的作用之后，得到的新向量（長(zhǎng)度也許改變）仍然與原來(lái)的v 保持在同一條線上。一個(gè)特征向量的長(zhǎng)度在該線性變換下縮放的比例稱為其特征值（本征值）。如果特征值為正，則表示v 在經(jīng)過(guò)線性變換的作用后方向也不變；如果特征值為負(fù)，說(shuō)明方向會(huì)反轉(zhuǎn)；如果特征值為0，則是表示縮回零點(diǎn)。但無(wú)論怎樣，仍在同一條直線上。5

19、. 下面是一個(gè)線性方程組：(1) 求方程的解。(2) 將方程右邊向量元素b3改為0.53再求解，并比較b3的變化和解的相對(duì)變化。(3) 計(jì)算系數(shù)矩陣A的條件數(shù)并分析結(jié)論。>> format rat %用分?jǐn)?shù)格式顯示>> A=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6A = 1/2 1/3 1/4 1/3 1/4 1/5 1/4 1/5 1/6 >> format %恢復(fù)默認(rèn)格式>> b=0.95;0.67;0.52b = 0.9500 0.6700 0.5200>> X=AbX = 1.2000 0.600

20、0 0.6000>> b2=0.95;0.67;0.53b2 = 0.9500 0.6700 0.5300>> X2=Ab2X2 = 3.0000 -6.6000 6.6000>> D=cond(A)D = 1.3533e+03矩陣的條件數(shù)決定矩陣的性能，條件數(shù)越接近于1其性能越好，系數(shù)矩陣A的條件數(shù)為1.3533e+003，和1相差很大，則其性能不好。因此b矩陣個(gè)別元素的微小變動(dòng)，對(duì)方程的解（X值）影響很大。6. 建立A矩陣，試比較sqrtm(A)和sqrt(A)，分析它們的區(qū)別。實(shí)驗(yàn)三 選擇結(jié)構(gòu)程序設(shè)計(jì)1. 求分段函數(shù)的值。用if語(yǔ)句實(shí)現(xiàn)，分別輸出x=-

21、5.0,-3.0,1.0,2.0,2.5,3.0,5.0時(shí)的y值。程序設(shè)計(jì)：clear;x=-5.0 -3.0 1.0 2.0 2.5 3.0 5.0;if x<0&x=-3 y=x.2+x-6;elseif x>=0&x<5&x=2&x=3 y=x.2-5*x+6;else y=x.2-x-1;endy運(yùn)行結(jié)果：y = 29.0000 11.0000 -1.0000 1.0000 2.7500 5.0000 19.00002. 輸入一個(gè)百分制成績(jī)，要求輸出成績(jī)等級(jí)A、B、C、D、E。其中90分100分為A，80分89分為B，79分79分為C，

22、60分69分為D，60分以下為E。要求：(1) 分別用if語(yǔ)句和switch語(yǔ)句實(shí)現(xiàn)。程序設(shè)計(jì)：（一） 用if語(yǔ)句clear;score=input('請(qǐng)輸入百分制成績(jī)：');a=fix(score/10);if a=9 | a=10 grade='A'elseif a=8 grade='B'elseif a=7 grade='C'elseif a=6 grade='D'else grade='E'endgrade（二） 用if語(yǔ)句clear;score=input('請(qǐng)輸入百分制成績(jī)：&#

23、39;);if score>=90&score<=100 grade='A'elseif score>=80&score<90 grade='B'elseif score>=70&score<80 grade='C'elseif score>=60&score<70 grade='D'else grade='E'endgrade（三） 用switch語(yǔ)句clear;score=input('請(qǐng)輸入百分制成績(jī)：');swi

24、tch fix(score/10) case 9,10 grade='A' case 8 grade='B' case 7 grade='C' case 6 grade='D' otherwise grade='E'endgrade (2) 輸入百分制成績(jī)后要判斷該成績(jī)的合理性，對(duì)不合理的成績(jī)應(yīng)輸出出錯(cuò)信息。clear;score=input('請(qǐng)輸入百分制成績(jī)：');if score>=0&score<=100 %判斷成績(jī)合理性 switch fix(score/10) cas

25、e 9,10 grade='A' case 8 grade='B' case 7 grade='C' case 6 grade='D' otherwise grade='E' end gradeelse disp('error') %成績(jī)不合理時(shí)輸出出錯(cuò)信息end3. 硅谷公司員工的工資計(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ā)工資。程序設(shè)計(jì)：clea

26、r;num=input('請(qǐng)輸入員工工號(hào)：');time=input('請(qǐng)輸入員工工時(shí)數(shù)：');if time<60 wages=time*84-700;elseif time>120 wages=84*120+(time-120)*84*(1+0.15);else wages=time*84;enddisp=('員工工號(hào)：' num2str(num) ' 應(yīng)發(fā)工資：' num2str(wages)運(yùn)行結(jié)果：請(qǐng)輸入員工工號(hào)：4請(qǐng)輸入員工工時(shí)數(shù)：35disp =員工工號(hào)：4 應(yīng)發(fā)工資：22404. 設(shè)計(jì)程序，完成兩位數(shù)的

27、加、減、乘、除四則運(yùn)算，即產(chǎn)生兩個(gè)兩位隨機(jī)整數(shù)，再輸入一個(gè)運(yùn)算符號(hào)，做相應(yīng)的運(yùn)算，并顯示相應(yīng)的結(jié)果。解：程序設(shè)計(jì)：clear;a=floor(rand(1)*90+10)b=floor(rand(1)*90+10)c=input('請(qǐng)輸入運(yùn)算符號(hào)：','s');if c='+' s=a+belseif c='-' s=a-belseif c='*' s=a*belseif c='' s=abelse disp('error')end運(yùn)行結(jié)果：a = 30b = 77請(qǐng)輸入運(yùn)算符號(hào)：-

28、s = -475. 建立5×6矩陣，要求輸出矩陣第n行元素。當(dāng)n值超過(guò)矩陣的行數(shù)時(shí)，自動(dòng)轉(zhuǎn)為輸出矩陣最后一行元素，并給出出錯(cuò)信息。程序設(shè)計(jì)：clear;A=rand(5,6)n=input('n=?')if n>=1&n<=5 B=A(n,1:6);else disp('erro!') B=A(5,1:6);endB運(yùn)行結(jié)果：A = 0.7792 0.0119 0.5285 0.6892 0.9133 0.0782 0.9340 0.3371 0.1656 0.7482 0.1524 0.4427 0.1299 0.1622 0.6

29、020 0.4505 0.8258 0.1067 0.5688 0.7943 0.2630 0.0838 0.5383 0.9619 0.4694 0.3112 0.6541 0.2290 0.9961 0.0046n=?3n = 3B = 0.1299 0.1622 0.6020 0.4505 0.8258 0.1067實(shí)驗(yàn)四 循環(huán)結(jié)構(gòu)程序設(shè)計(jì)1. 根據(jù)，求的近似值。當(dāng)n分別取100、1000、10000時(shí)，結(jié)果是多少？要求：分別用循環(huán)結(jié)構(gòu)和向量運(yùn)算（使用sum函數(shù)）來(lái)實(shí)現(xiàn)。程序設(shè)計(jì)（一）：clear;for n=100,1000,10000 sum=0; for x=1:n sum=sum

30、+1/(x.2); end result=sqrt(6*sum)end運(yùn)行結(jié)果：result = 3.1321result = 3.1406result = 3.1415程序設(shè)計(jì)（二）：clear;for n=100,1000,10000 x=1:n; result=sqrt(6*sum(1./x./x)end運(yùn)行結(jié)果：result = 3.1321result = 3.1406result =3.14152. 根據(jù)，求：(1) y<3時(shí)的最大n值。(2) 與(1)的n值對(duì)應(yīng)的y值。程序設(shè)計(jì)：y=0;n=1;while y<3 y=y+1/(2*n-1); n=n+1;endn=n

31、-1;y=y-1/(2*n-1)n=n-1運(yùn)行結(jié)果：y = 2.9944n = 56【%驗(yàn)證n=56時(shí)的y值clear all;clc;n=56;i=1:n;f=1./(2*i-1);y=sum(f)】3. 考慮以下迭代公式：其中a、b為正的學(xué)數(shù)。(1) 編寫程序求迭代的結(jié)果，迭代的終止條件為|xn+1-xn|10-5，迭代初值x0=1.0，迭代次數(shù)不超過(guò)500次。(2) 如果迭代過(guò)程收斂于r，那么r的準(zhǔn)確值是，當(dāng)(a,b)的值取(1,1)、(8,3)、(10,0.1)時(shí)，分別對(duì)迭代結(jié)果和準(zhǔn)確值進(jìn)行比較。程序設(shè)計(jì)：clear;a=input('a=?');b=input(

32、9;b=?');x0=1;x1=a/(b+x0);y=abs(x1-x0);n=1;while y>10(-5)&n<=499 x0=x1; x1=a/(b+x0); y=abs(x1-x0); n=n+1;endx1 %比較收斂準(zhǔn)確值和迭代結(jié)果r1=(-b+sqrt(b2+4*a)/2r2=(-b-sqrt(b2+4*a)/2運(yùn)行結(jié)果：a=?1b=?1x1 = 0.6180r1 = 0.6180r2 = -1.6180a=?8b=?3x1 = 1.7016r1 = 1.7016r2 = -4.7016a=?10b=?0.1x1 = 3.1127r1 = 3.112

33、7r2 = -3.21274. 已知求f1f100中：(1) 最大值、最小值、各數(shù)之和。(2) 正數(shù)、零、負(fù)數(shù)的個(gè)數(shù)。程序設(shè)計(jì)：clear;f1=1;f2=0;f3=1;f=f1,f2,f3;for n=4:100 f4=f3-2*f2+f1; f=f f4; f1=f2; f2=f3; f3=f4;end maxf=max(f)minf=min(f)sumf=sum(f) positivenum=length(find(f>0)negtivenum=length(find(f<0)zeronum=length(find(f=0)運(yùn)行結(jié)果：maxf = 4.3776e+11minf

34、 = -8.9941e+11sumf = -7.4275e+11positivenum = 49negtivenum = 49zeronum = 25. 若兩個(gè)連續(xù)自然數(shù)的乘積減1是素?cái)?shù)，則稱這兩個(gè)邊疆自然數(shù)是親密數(shù)對(duì)，該素?cái)?shù)是親密素?cái)?shù)。例如，2×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ù)之和。程序設(shè)計(jì)（一）：clearpair=;for m=2:49 f=(m+1)*m-1; signf=0; %標(biāo)記為0表示f為素?cái)?shù) for k=2:sqrt(f) if (rem(f,k)=0) s

35、ignf=1; %標(biāo)記為1表示f不是素?cái)?shù) end end if signf=0 pair=pair;m m+1; %親密數(shù)對(duì)矩陣 endendnum=size(pair);susum=sum(pair(:,1).*pair(:,2)-1);disp('親密數(shù)對(duì)的對(duì)數(shù)',num2str(num(1)disp('親密素?cái)?shù)之和',num2str(susum)運(yùn)行結(jié)果：親密數(shù)對(duì)的對(duì)數(shù)28親密素?cái)?shù)之和21066程序設(shè)計(jì)：clearsusum=0;num=0;for m=2:49 f=(m+1)*m-1; n=fix(sqrt(f); for k=2:n if (rem(f

36、,k)=0) break; %f不是素?cái)?shù)時(shí)跳出循環(huán) end end if k=n num=num+1; %親密數(shù)對(duì)個(gè)數(shù) susum=susum+f; %親密素?cái)?shù)和累加 endend disp('親密數(shù)對(duì)的對(duì)數(shù)',num2str(num)disp('親密素?cái)?shù)之和',num2str(susum)運(yùn)行結(jié)果：親密數(shù)對(duì)的對(duì)數(shù)28親密素?cái)?shù)之和21066實(shí)驗(yàn)五 函數(shù)文件1. 定義一個(gè)函數(shù)文件，求給定復(fù)數(shù)的指數(shù)、對(duì)數(shù)、正弦和余弦，并在命令文件中調(diào)用該函數(shù)文件。建立函數(shù)文件fushu.m：function e,ln,s,c=fushu(x)e=exp(x);ln=log(x);s

37、=sin(x);c=cos(x);建立命令文件fushu2.mA=input('請(qǐng)輸入任意一個(gè)復(fù)數(shù)A');a,b,c,d=fushu(A); %調(diào)用fushu.m的函數(shù)文件abcdfushu2.m運(yùn)行結(jié)果：>> fushu2請(qǐng)輸入任意一個(gè)復(fù)數(shù)A：4ia = -0.6536 - 0.7568ib = 1.3863 + 1.5708ic = 0.0000 +27.2899id = 27.30822. 一物理系統(tǒng)可用下列方程組來(lái)表示：從鍵盤輸入m1、m2和的值，求a1、a2、N1和N2的值。其中g(shù)取9.8，輸入時(shí)以角度為單位。要求：定義一個(gè)求解線性方程組AX=B的函數(shù)文件

38、，然后在命令文件中調(diào)用該函數(shù)文件。建立函數(shù)文件xxfcz.mfunction a1,a2,N1,N2=xxfcz(m1,m2,t)t=t*pi/180;A=m1*cos(t) -m1 -sin(t) 0; m1*sin(t) 0 cos(t) 0; 0 m2 -sin(t) 0; 0 0 -cos(t) 1;g=9.8;B=0;m1*g;0;m2*g;X=inv(A)*B;a1=X(1);a2=X(2);N1=X(3);N2=X(4);建立命令文件xxfcz2.mclearm1=input('請(qǐng)輸入m1：');m2=input('請(qǐng)輸入m2：');t=input

39、('請(qǐng)輸入角度t：');a1,a2,N1,N2=xxfcz(m1,m2,t);a1a2N1N2運(yùn)行結(jié)果：>> xxfcz2請(qǐng)輸入m1：5請(qǐng)輸入m2：3請(qǐng)輸入角度t：90a1 = 9.8000a2 = 3.7505e-16N1 = 1.1251e-15N2 = 29.40003. 一個(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ù)文件。建立判斷素?cái)?shù)的函數(shù)文件sushu.mfunction x=sushu(m)%判斷數(shù)m是否為素?cái)?shù)for M=2:(m/2) if rem(m,M)=0 b=0;

40、break else b=1; endendif b=1 x=m; %判斷數(shù)m為素?cái)?shù)，并返回else x=; %數(shù)m不為素?cái)?shù)endend建立命令文件sushu2.mclear;js=; %初始絕對(duì)素?cái)?shù)for m=10:99 x=sushu(m); if x x=fix(x/10)+rem(x,10)*10; y=sushu(x); js=js,y; %絕對(duì)素?cái)?shù) else continue endendjs執(zhí)行命令文件sushu2.m>> sushu2js = 11 31 71 13 73 17 37 97 79法二：建立函數(shù)文件sushu3.mfunction p=sushu3(m

41、)b=m(end); p=m;for i=2:sqrt(b) n=find(rem(p,i)=0&p=i);%將m中能被i整除，且不等于i的元素的下標(biāo)找出，即n p(n)=; %將下標(biāo)為n的元素剔除，其余的即為素?cái)?shù)end建立命令文件sushu4.mclear;m=10:99;p=sushu3(m); %找出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=sushu3(p) %再對(duì)調(diào)換順序后的素?cái)?shù)矩陣找出其中的素?cái)?shù)執(zhí)行命令文件sushu4.m>> sushu4p = 11 31 71 13 73 1

42、7 37 97 794. 設(shè)，編寫一個(gè)MATLAB函數(shù)文件fx.m，使得調(diào)用f(x)時(shí)，x可用矩陣代入，得出的f(x)為同階矩陣。建立函數(shù)文件fx.mfunction f=fx(x)f=1./(x-2).2+0.1)+1./(x-3).4+0.01);end5. 已知(1) 當(dāng)f(n)=n+10ln(n2+5)時(shí)，求y的值。(2) 當(dāng)f(n)=1×2+2×3+3×4+.+n×(n+1)時(shí)，求y的值。建立函數(shù)文件fn1.mfunction fn=fn1(n)fn=n+10*log(n.2+5);end建立函數(shù)文件fn2.mfunction fn=fn2(n

43、)fn=0;for x=1:n fn=fn+x*(x+1); end求（1）（2）中y的值：y1=fn1(40)/(fn1(30)+fn1(20);y2=fn2(40)/(fn2(30)+fn2(20);運(yùn)行結(jié)果：>> yfny1 = 0.6390y2 =1.7662實(shí)驗(yàn)六 高層繪圖操作1. 設(shè)，在x=02區(qū)間取101點(diǎn)，繪制函數(shù)的曲線。>> x=linspace(0,2*pi,101);>> y=(0.5+3.*sin(x)./(1+x.2).*cos(x);>> plot(x,y);2. 已知y1=x2，y2=cos(2x)，y3=y1

44、15;y2，完成下列操作：(1) 在同一坐標(biāo)系下用不同的顏色和線型繪制三條曲線。>> x=linspace(0,2*pi,101);>> y1=x.2;>> y2=cos(2*x);>> y3=y1.*y2;>> plot(x,y1,'k:',x,y2,'b-',x,y3,'m-.');(2) 以子圖形式繪制三條曲線。x=linspace(0,2*pi,101);y1=x.2;y2=cos(2*x);y3=y1.*y2;subplot(3,1,1); plot(x,y1); title(

45、'y1'); subplot(3,1,2); plot(x,y2); title('y2'); subplot(3,1,3); plot(x,y3); title('y3');(3) 分別用條形圖、階梯圖、桿圖和填充圖繪制三條曲線。subplot(3,4,1);bar(x,y1,'y');title('bar(x,y1,"y")'); subplot(3,4,2);stairs(x,y1,'m');title('stairs(x,y1,"m")'); subplot(3,4,3);stem(x,y