Matlab上機實驗答案重點講義

1、Matlab上機實驗答案實驗一 MATLAB運算基礎1. 先求下列表達式的值，然后顯示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)可以驗證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. 已知：求下列表達式的值：(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. 設有矩陣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ù)的個數(shù)。(2) 建立一個字符串向量，刪除其中的大寫字母。>> 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實驗二 MATLAB矩陣分析與處理1. 設有分塊矩陣，其中E、R、O、S分別為單位矩陣、隨機矩陣、零矩陣和對角陣，試通過數(shù)值計算驗證。>> 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通過驗證，矩陣成立。2. 產(chǎn)生5階希爾伯特矩陣H和5階帕斯卡矩陣P，且求其行列式的值Hh和Hp以及它們的條件數(shù)Th和Tp，判斷哪個矩陣性能更好。為什么？>> 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其性能就越好。由上機操作求得Th=4.7661e+005，Tp=8.5175e+003。Tp的值更接近于1則其性能要好。所以5階帕斯卡矩陣P的性能好。3. 建立一個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ù)學意義。>> 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ù)學上，特別是線性代數(shù)中，對于一個給定的線性變換，它的特征向量（本征向量或稱正規(guī)正交向量）是這樣一個非零的向量v：當v 經(jīng)過這個線性變換的作用之后，得到的新向量（長度也許改變）仍然與原來的v 保持在同一條線上。一個特征向量的長度在該線性變換下縮放的比例稱為其特征值（本征值）。如果特征值為正，則表示v 在經(jīng)過線性變換的作用后方向也不變；如果特征值為負，說明方向會反轉(zhuǎn)；如果特征值為0，則是表示縮回零點。但無論怎樣，仍在同一條直線上。5

19、. 下面是一個線性方程組：(1) 求方程的解。(2) 將方程右邊向量元素b3改為0.53再求解，并比較b3的變化和解的相對變化。(3) 計算系數(shù)矩陣A的條件數(shù)并分析結(jié)論。>> format rat %用分數(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 %恢復默認格式>> 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矩陣個別元素的微小變動，對方程的解（X值）影響很大。6. 建立A矩陣，試比較sqrtm(A)和sqrt(A)，分析它們的區(qū)別。實驗三 選擇結(jié)構(gòu)程序設計1. 求分段函數(shù)的值。用if語句實現(xiàn)，分別輸出x=-

21、5.0,-3.0,1.0,2.0,2.5,3.0,5.0時的y值。程序設計：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運行結(jié)果：y = 29.0000 11.0000 -1.0000 1.0000 2.7500 5.0000 19.00002. 輸入一個百分制成績，要求輸出成績等級A、B、C、D、E。其中90分100分為A，80分89分為B，79分79分為C，

22、60分69分為D，60分以下為E。要求：(1) 分別用if語句和switch語句實現(xiàn)。程序設計：（一） 用if語句clear;score=input('請輸入百分制成績：');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語句clear;score=input('請輸入百分制成績：&#

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語句clear;score=input('請輸入百分制成績：');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) 輸入百分制成績后要判斷該成績的合理性，對不合理的成績應輸出出錯信息。clear;score=input('請輸入百分制成績：');if score>=0&score<=100 %判斷成績合理性 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') %成績不合理時輸出出錯信息end3. 硅谷公司員工的工資計算方法如下：(1) 工作時數(shù)超過120小時者，超過部分加發(fā)15%。(2) 工作時數(shù)低于60小時者，扣發(fā)700元。(3) 其余按每小時84元計發(fā)。試編程按輸入的工號和該號員工的工時數(shù)，計算應發(fā)工資。程序設計：clea

26、r;num=input('請輸入員工工號：');time=input('請輸入員工工時數(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=('員工工號：' num2str(num) ' 應發(fā)工資：' num2str(wages)運行結(jié)果：請輸入員工工號：4請輸入員工工時數(shù)：35disp =員工工號：4 應發(fā)工資：22404. 設計程序，完成兩位數(shù)的

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

28、s = -475. 建立5×6矩陣，要求輸出矩陣第n行元素。當n值超過矩陣的行數(shù)時，自動轉(zhuǎn)為輸出矩陣最后一行元素，并給出出錯信息。程序設計：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運行結(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實驗四 循環(huán)結(jié)構(gòu)程序設計1. 根據(jù)，求的近似值。當n分別取100、1000、10000時，結(jié)果是多少？要求：分別用循環(huán)結(jié)構(gòu)和向量運算（使用sum函數(shù)）來實現(xiàn)。程序設計（一）：clear;for n=100,1000,10000 sum=0; for x=1:n sum=sum

30、+1/(x.2); end result=sqrt(6*sum)end運行結(jié)果：result = 3.1321result = 3.1406result = 3.1415程序設計（二）：clear;for n=100,1000,10000 x=1:n; result=sqrt(6*sum(1./x./x)end運行結(jié)果：result = 3.1321result = 3.1406result =3.14152. 根據(jù)，求：(1) y<3時的最大n值。(2) 與(1)的n值對應的y值。程序設計：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運行結(jié)果：y = 2.9944n = 56【%驗證n=56時的y值clear all;clc;n=56;i=1:n;f=1./(2*i-1);y=sum(f)】3. 考慮以下迭代公式：其中a、b為正的學數(shù)。(1) 編寫程序求迭代的結(jié)果，迭代的終止條件為|xn+1-xn|10-5，迭代初值x0=1.0，迭代次數(shù)不超過500次。(2) 如果迭代過程收斂于r，那么r的準確值是，當(a,b)的值取(1,1)、(8,3)、(10,0.1)時，分別對迭代結(jié)果和準確值進行比較。程序設計：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 %比較收斂準確值和迭代結(jié)果r1=(-b+sqrt(b2+4*a)/2r2=(-b-sqrt(b2+4*a)/2運行結(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ù)、零、負數(shù)的個數(shù)。程序設計：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)運行結(jié)果：maxf = 4.3776e+11minf

34、 = -8.9941e+11sumf = -7.4275e+11positivenum = 49negtivenum = 49zeronum = 25. 若兩個連續(xù)自然數(shù)的乘積減1是素數(shù)，則稱這兩個邊疆自然數(shù)是親密數(shù)對，該素數(shù)是親密素數(shù)。例如，2×3-1=5，由于5是素數(shù)，所以2和3是親密數(shù)，5是親密素數(shù)。求2,50區(qū)間內(nèi)：(1) 親密數(shù)對的對數(shù)。(2) 與上述親密數(shù)對對應的所有親密素數(shù)之和。程序設計（一）：clearpair=;for m=2:49 f=(m+1)*m-1; signf=0; %標記為0表示f為素數(shù) for k=2:sqrt(f) if (rem(f,k)=0) s

35、ignf=1; %標記為1表示f不是素數(shù) end end if signf=0 pair=pair;m m+1; %親密數(shù)對矩陣 endendnum=size(pair);susum=sum(pair(:,1).*pair(:,2)-1);disp('親密數(shù)對的對數(shù)',num2str(num(1)disp('親密素數(shù)之和',num2str(susum)運行結(jié)果：親密數(shù)對的對數(shù)28親密素數(shù)之和21066程序設計：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不是素數(shù)時跳出循環(huán) end end if k=n num=num+1; %親密數(shù)對個數(shù) susum=susum+f; %親密素數(shù)和累加 endend disp('親密數(shù)對的對數(shù)',num2str(num)disp('親密素數(shù)之和',num2str(susum)運行結(jié)果：親密數(shù)對的對數(shù)28親密素數(shù)之和21066實驗五 函數(shù)文件1. 定義一個函數(shù)文件，求給定復數(shù)的指數(shù)、對數(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('請輸入任意一個復數(shù)A');a,b,c,d=fushu(A); %調(diào)用fushu.m的函數(shù)文件abcdfushu2.m運行結(jié)果：>> fushu2請輸入任意一個復數(shù)A：4ia = -0.6536 - 0.7568ib = 1.3863 + 1.5708ic = 0.0000 +27.2899id = 27.30822. 一物理系統(tǒng)可用下列方程組來表示：從鍵盤輸入m1、m2和的值，求a1、a2、N1和N2的值。其中g(shù)取9.8，輸入時以角度為單位。要求：定義一個求解線性方程組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('請輸入m1：');m2=input('請輸入m2：');t=input

39、('請輸入角度t：');a1,a2,N1,N2=xxfcz(m1,m2,t);a1a2N1N2運行結(jié)果：>> xxfcz2請輸入m1：5請輸入m2：3請輸入角度t：90a1 = 9.8000a2 = 3.7505e-16N1 = 1.1251e-15N2 = 29.40003. 一個自然數(shù)是素數(shù)，且它的數(shù)字位置經(jīng)過任意對換后仍為素數(shù)。例如13是絕對素數(shù)。試求所有兩位絕對素數(shù)。要求：定義一個判斷素數(shù)的函數(shù)文件。建立判斷素數(shù)的函數(shù)文件sushu.mfunction x=sushu(m)%判斷數(shù)m是否為素數(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為素數(shù)，并返回else x=; %數(shù)m不為素數(shù)endend建立命令文件sushu2.mclear;js=; %初始絕對素數(shù)for m=10:99 x=sushu(m); if x x=fix(x/10)+rem(x,10)*10; y=sushu(x); js=js,y; %絕對素數(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的元素的下標找出，即n p(n)=; %將下標為n的元素剔除，其余的即為素數(shù)end建立命令文件sushu4.mclear;m=10:99;p=sushu3(m); %找出10到99內(nèi)所有素數(shù)p=10*rem(p,10)+(p-rem(p,10)/10; %將p素數(shù)矩陣每個元素個位十位調(diào)換順序p=sushu3(p) %再對調(diào)換順序后的素數(shù)矩陣找出其中的素數(shù)執(zhí)行命令文件sushu4.m>> sushu4p = 11 31 71 13 73 1

42、7 37 97 794. 設，編寫一個MATLAB函數(shù)文件fx.m，使得調(diào)用f(x)時，x可用矩陣代入，得出的f(x)為同階矩陣。建立函數(shù)文件fx.mfunction f=fx(x)f=1./(x-2).2+0.1)+1./(x-3).4+0.01);end5. 已知(1) 當f(n)=n+10ln(n2+5)時，求y的值。(2) 當f(n)=1×2+2×3+3×4+.+n×(n+1)時，求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);運行結(jié)果：>> yfny1 = 0.6390y2 =1.7662實驗六 高層繪圖操作1. 設，在x=02區(qū)間取101點，繪制函數(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) 在同一坐標系下用不同的顏色和線型繪制三條曲線。>> 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