MATLAB實(shí)驗(yàn)題答案

1、1、求以下變量的值，并在MATLAB中驗(yàn)證。( 1 ) a = 1 : 2 : 5a =135( 2 ) b = a' , a' , a' ； a b =111333555135( 3 ) c = a + b ( 2 , : )c =4682、下列運(yùn)算是否合法，為什么？如合法，結(jié)果是多少？>> result2=a*bError using*Inner matrix dimensions must agree.>> result3=a+b result3 =3625811>> result4=b*d result4 =312222404

2、913>> result5=b;c'*d result5 =312222404913-5-87>> result6=a.*b result6 =28-341530>> result7=a./b result7 =0.50000.5000-3.00004.00001.66671.2000>> result8=a.cAttempt to reference field of non-structure array.>> result9=a.b result9 =2.00002.0000-0.33330.25000.60000.83

3、33>> result10=a.2 result10 =149162536>> result11=2.a result11 =2481632643、用 MATLAB求解下面的的方程組。7212x1491532x27(1)2115x31213213x40>> A=7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13>> B=4 7 -1 0>> B=B'>> x=inv(A)*B>> A1=1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6>> B2

4、=1;8;3;5>> x2=inv(A1)*B2721291532A211524、已知13213(1) 求矩陣 A 的秩 (rank)(2) 求矩陣 A 的行列式 (determinant)(3) 求矩陣 A 的逆 (inverse)(4) 求矩陣 A 的特征值及特征向量(eigenvalue and eigenvector)>> A3=7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13>> r=rank(A3)>> b=inv(A3)>> a=det(A3)>> V,D=eig(A3)5、102

5、 102 929210y2nn10, 求 y=？（運(yùn)行 format long g命令后，查看 y 的值）m1=0;for m=-10:10m1=m1+2m;endm1m1 =6、求分段函數(shù)的值。用 if語句實(shí)現(xiàn)，算出下列表中x 對應(yīng)的 y值。x=input('enter x=');if x<0y=x2+x-6;elseif x>=0&&x<5y=x2-5*x+6;elsey=x2-x-1;endy7、分別用 if 和 switch 語句實(shí)現(xiàn)，將百分制成績轉(zhuǎn)換為成績等級 A、B、C、D、E。其中 90100 分為 A，8089 分為 B，707

6、9分為 C，6069 分為 D，60 分以下為 E。對超出百分制范圍的成績， 給出錯誤提示信息。if結(jié)構(gòu)程序：x=input('please enter score=');if x>=90&&x<=100disp('A')elseif x<90&&x>=80disp('B')elseif x<80&&x>=70disp('C')elseif x<70&&x>=60disp('D')elseif x<

7、60&&x>=0disp('E')elsedisp('error')endswitch結(jié)構(gòu)程序：x=input('please enter sco re=');switch fix(x/10)case10,9if x>100disp('error')elsedisp('A')endcase8disp('B')case7disp('C')case6disp('D')case0,1,2,3,4,5disp('E')otherwi

8、sedisp('error')end8、思考題設(shè)計(jì)程序，完成成兩位數(shù)的加、減、乘、除四則運(yùn)算，即產(chǎn)生兩個兩位隨機(jī)整數(shù)，再輸入一個運(yùn)算符號，做相應(yīng)的運(yùn)算，并顯示相應(yīng)的結(jié)果。x=input(' 請輸入運(yùn)算符')a=num2str(floor(rand(1)*90+10);ab=num2str(floor(rand(1)*90+10);bif x='+'y=a+b;elseif x='-'y=a-b;elseif x='*'y=a*b;elseif x='/'y=a/b;elsedisp('err

9、or')endy9、啟動 MATLAB后，點(diǎn)擊 File|New|M-File ， 啟動 MATLAB的程序編輯及調(diào)試器（ Editor/Debugger ），編輯以下程序，點(diǎn)擊 File|Save 保存程序， 注意文件名最好用英文字符。 點(diǎn)擊 Debug|Run 運(yùn)行程序，在命令窗口查看運(yùn)行結(jié)果， 程序如有錯誤則改正。注：數(shù)論中一個有趣的題目： 任意一個正整數(shù)，若為偶數(shù)， 則用 2 除之，若為奇數(shù)，則與 3 相乘再加上 1。重復(fù)此過程，最終得到的結(jié)果為 1。n=input(' 請輸入 n 值:');a=n;while n>1if rem(n,2)=0n=n/2;

10、elsen=3*n+1;enda=a,n;enda10、根據(jù)x211116122232n2，當(dāng) n 分別取 100、 1000、 10000 時，求 x 的值分別是多少？a=input(' 請輸入數(shù)值 ')for m=1:100n=n+1/(m*m);endn=6*n;x=sqrt(n);x11、編程求滿足km1 2k10000的最小 m值。sum=0;m=2;a=1;while afor i=1:msum=sum+2i;if sum>10000a=0;endendm=m+1;endm12、思考題已知y和t的函數(shù)關(guān)系：y ( t )1tt 2 / 2!t 3 / 3!.t

11、 20 / 20 !求下面表格中與t 對應(yīng)的 y 值n=0;t=input(' 請輸入 t 值： ')sum=0;for i=1:20b=factorial(i);sum=sum+ti/b;endsum=sum+1;sum13、編寫一個函數(shù)，計(jì)算下面函數(shù)的值，給出標(biāo)量 x 的值，調(diào)用該函數(shù)后，返回 y 的值。function y=myfun1(x)選擇一些數(shù)據(jù)測試你編寫的函數(shù)。functiony=myfun1(x)if x<=0y=sin(x);elseif x>0&x<=3y=x;elseif x>3y=-x+6;end14、編寫一個函數(shù)求向量

12、 x 中元素的平均值、最大值、最小值、均方根值。functionm_x,max_x,min_x,rms_x=myfun2(x)方均根值 (Root Mean Square)的計(jì)算公式為：用下面數(shù)據(jù)測試你寫的函數(shù)：（ 1） x=sin(0:0.01:6*pi)t0.20.40.60.81.0y（ 2）x=rand(1,200) ，得到的 x 為 200 個（0，1）之間均勻分布的隨機(jī)數(shù)。functionm_x,max_x,min_x,rms_x=myfun2(x)%求平均值sum_x=sum(x); %向量元素求和m,n=size(x);%最好用 n=length(x);m_x=sum_x/n;

13、%求最大值采用逐個比較方式if x(1)>x(2)max_x=x(1);elsemax_x=x(2);endfor k=3:nif max_x<x(k)max_x=x(k);elsemax_x=max_x;%可省略endend%求最小值if x(1)<x(2)min_x=x(1);elsemin_x=x(2);endfor k=3:nif min_x>x(k)min_x=x(k);elsemin_x=min_x;%可省略endend%求均方根值sum_x2=0;for k=1:nsum_x2=sum_x2+x(k).2;rms_x=sqrt(sum_x2/n);endm

14、_x;max_x;min_x;rms_x;%按照函數(shù)值行參順序輸出結(jié)果15、編寫一個函數(shù)，給出一個向量x x1, x2,xn ，生成如下范德蒙矩陣。function v=myvander(x)例如： >>v=myvander(2 3 4 5)得 v=function v=myvander(x)v1=vander(x);%生成范德蒙矩陣v2=v1'v=flipud(v2);%實(shí)現(xiàn)矩陣上下翻轉(zhuǎn)16、思考題編寫程序， 用如下迭代公式求a ，a 的值分別為： 3，17，113。迭代的終止條件xx 10 5為n 1n，迭代初值x0 1.0 , 迭代次數(shù)不超過 100 次。分別對迭代結(jié)

15、果和準(zhǔn)確值進(jìn)行比較， 并xn 1xna22xn統(tǒng)計(jì)迭代次數(shù)。functionx,n=sqrt_a(a)x=1.0;for k=1:100m=x;x=x/2+a/(2*x);if abs(x-m)<=10(-5)breakendendx;n=k;s=(x-sqrt(a);if s<=10(-5)disp(' 正確 ');elsedisp(' 錯誤 ');enduc=cos(2*pi/3)+sin(2*pi/3)*i;17、在同一坐標(biāo)系下繪制下面三個函數(shù)在compass(ua,ub,uc,ua-ub,ub-uc,uc-ua)t0 ，4 的圖象。>&

16、gt; ua=1;>> t=0:0.1:4*pi;>> ob=cos(-2*pi/ 3)>> y1=t;>> ob=cos(-2*pi/ 3)+sin(-2*pi/ 3)*i;>> y2=t(0.5);>> ub=cos(-2*pi/ 3)+sin(-2*pi/ 3)*i;>> y2=sqrt(t);>> uc=cos(2*pi/ 3)+sin(2*pi/ 3)*i;>> y3=4*pi*eps(-0.1*t).*sin(t);>> compass(ua,ub,uc,ua-u

17、b,ub-uc,uc-ua)>> plot(t,y1,t,y2,t,y3)20、三維空間曲線繪制18、編寫程序，選擇合適的步距，繪制下面z=0:0.1:4*pi;函數(shù)在區(qū)間 -6， 6 中的圖象。function y=fun1(x)x=cos(z);y=;y=sin(z);for x0=xplot3(x,y,z)if x0<=0>> z=0:0.1:4*pi;y=y,sin(x0);>> x=cos(z);elseif x0>0&&x0<=3>> y=sin(z);y=y,x0;>> plot3(x,

18、y,z)elseif x0>321、用 mesh或 surf 函數(shù)，繪制下面方程所表示的三維空間曲面， x 和 y 的取值范y=y,6-x0;圍設(shè)為 -3 ，3 。end>> x=-3:0.01:3;end>> y=-3:0.01:3;>> x=-6:0.1:6;>> x,y=meshgrid(x,y);>> y=fun1(x);>> z=-x.*x/10+y.*y/10;>> plot(x,y)>> mesh(x,y,z)19、用 compass 函數(shù)畫下面相量圖>> z=x.*

19、x/10+y.*y/10;ua=1;ub=>> mesh(x,y,z)cos(-2*pi/3)+sin(-2*pi/3)*i;22、思考題用對分法求解方程2e xsin x 在 0 ，1 內(nèi)的解， 并驗(yàn)證， 在程序中統(tǒng)計(jì)出對分次數(shù)。提示：先將原方程轉(zhuǎn)化成f (x) 2e xsin x 0的 形式。對分法的基本思想是：一個一元方程 f(x)=0 ，若 f(x1)*f(x2)<0 ，則在 x1,x2 區(qū)間內(nèi)有實(shí)數(shù)解。 取該區(qū) 間的中點(diǎn) xm=(x1+x2)/2 ，判定 f(x1) 和 f(x2) 二者 中 哪 一 個 與 f(xm) 異 號 ， 若 f(x1)*f(xm)<

20、0 ，則解存在的區(qū)間縮小為 x1,xm ，否則解存在的區(qū)間縮小為xm,x2 。重復(fù)這樣的步驟，直到區(qū)間的長度小于一個可以接受的小數(shù) （比如 1e-10 ） , 則認(rèn)為中點(diǎn)即是原方程的解。>> x0=0;x1=1;y=1;i=0;while y=0i=i+1;x3=(x0+x1)/2;y=2*eps(-x3)-sin(x3);y0=2*eps(-x0)-sin(x0);y1=2*eps(-x1)-sin(x1);if y*y0<0x0=x3;elseif y*y1<0x1=x3;elseif y=0x3=x3;endif i=1000000y=0;endEnd23、求代數(shù)

21、方程3x54x 47 x32x 29x120 的 5個根，并將其用星號（* ）標(biāo)記在復(fù)平面圖上。（用 roots和 plot函數(shù)）>> a=3,4,7,2,9,12;>> y=roots(a)>> plot(y,'*')求代數(shù)方程 x510 的 5 個根，并將其用星號（ * ）標(biāo)記在復(fù)平面圖上。（用 roots和 plot 函數(shù)）程序：>> a1=1,0,0,0,0,-1;>> y1=roots(a1)>> plot(y1,'*')24、求下面函數(shù)在0.5,4區(qū)間內(nèi)的過零點(diǎn)。（用 fzero函）>> x=0.5:0.01:4;>> y2=x.3-2*x.2.*sin(x)+5*x.*cos(x)+1./x;>> plot(x,y2)函數(shù)：function y2=fz(x)y2=x.3-2*x.2.*sin(x)+5*x.*cos(x)+1./x;分析：根據(jù)圖像可知零點(diǎn)在x=1 和 x=3 附近所以>> fzero('