MATLAB試驗題答案

1、>> result5=b;c'*d( 1 ) a = 1: 2 :5result5 =a =1 3531 2222( 2 ) b = a' , a' , a'； a 40 4913b =-5 -871 11>> result6=a.*b333result6 =5551352 8-3( 3 ) c = a+ b ( 2 , : )4 1530c =>> result7=a./b468result7 =1、求以下變量的值，并在MATLAB 驗證。2、下列運算是否合法，為什么如合法，結(jié)果是多少>> result2=a*b

2、>> result8=result10 =40 49 131 4 9Error usingAttemptto referencefieldofInner matrix dimensions must agree.non-structure array.>> result3=a+b>> result9=a.bresult3 =result9 =5 8 11>> result4=b*dresult4 =>> result10=a.A231 22 221516 25 36>> result11=2.Aaresult11 =114

3、、已知13(1)求矩陣A 的秩 (rank)16 32 64(2)求矩陣A 的行列式 (determinant)3、用MATLAB求解下面的的方程組。x1(3)(4)求矩陣A 的逆 (inverse)求 矩陣 A 的特征值 及特征向 量 (eigenvalue and eigenvector)9(1) 215x211x3>> A3=7 2 1 -2;9 15 3 -2;-2 -2 11 5;13 2 1313x4>> r=rank(A3)>> A=7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13>> b=inv(A3)

4、>> a=det(A3)>> B=4 7 -1 0>> V,D=eig(A3)>> B=B'>> x=inv(A)*Bz1x(2) 2xy2y z w 83w 33xy3y 5z 6w 5102n2 10 2 9n 10, 求 y=( 運行 format long g 命令后， 查看 y 的值)29 210>> A1=1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6m1=0;for m=-10:10>> B2=1;8;3;5m1=m1+2m;>> x2=inv(A1)*

5、B2endm1m1 =6、求分段函數(shù)的值。elseif x<80&&x>=70x2x2x2disp('C')elseif x<70&&x>=605xx5disp('D')elseif x<60&&x>=0if 語句實現(xiàn),算出下列表中 x 對應(yīng)的 y 值。disp('E')elsex=input('enter x=');disp('error')if x<0endy=xA2+x-6;switch 結(jié)構(gòu)程序：elseif x>

6、;=0&&x<5x=input('please enter score=');y=x2-5*x+6;switch fix(x/10)elsecase10,9y=x2-x-1;if x>100enddisp('error')else7、分別用 if 和 switch 語句實現(xiàn)，將百分 制成績轉(zhuǎn)換為成績等級 A、B、C、D、E。其中 90100 分為 A,8089 分為 B, 7079分為 C，6069 分為 D，60 分以下為 E。對超出百分制范圍的成績, 給出錯誤提 示信息。disp('A')endcase8disp(

7、'B')if 結(jié)構(gòu)程序：case7x=input('please enter score=');disp('C')if x>=90&&x<=100disp('A')disp('D')elseif x<90&&x>=80case0,1,2,3,4,5disp('B')disp('E')case6otherwisedis p('eiTO門end（Editor/Debugger ），編輯以下程序， 點擊FilelSave 保存

8、程序，注意文件名 最好用英文字符。點擊Debug|Run運行 程序，在命令窗口查看運行結(jié)果，程序如有錯誤則改正。n=0;8、思考題設(shè)計程序，完成成兩位數(shù)的加、減、乘、除 四則運算，即產(chǎn)生兩個兩位隨機整數(shù)， 再輸入一個運算符號，做相應(yīng)的運算， 并顯示相應(yīng)的結(jié)果。注：數(shù)論中一個有趣的題目： 任意一個正整 數(shù)，若為偶數(shù)，則用2除之，若為奇數(shù), 則與3相乘再加上1。重復(fù)此過程，最 終得到的結(jié)果為1。n=inpu t(' 請輸入a=n;x=input（'請輸入運算符'）while n>1a=nu m2str(floor(ra nd(1)*90+10);if rem( n,2

9、)=0n=n/2;b=nu m2str(floor(ra nd(1)*90+10);elsen=3*n+1;if x='+'endy=a+b;a=a, n;elseif x='-'endy=a-b;elseif x='*'y=a*b;10elseif x='/'x2y=a/b;丄22else6，當(dāng)求x的值分別是多少分別取100、1n2丄321000、 10000 時,dis p('error')enda=input（'請輸入數(shù)值'）endfor m=1:100 n=n+1/(m*m);9、啟動 MA

10、TLAB，點擊 File|New|M-File ，啟動MATLAB的程序編輯及調(diào)試器n=6* n;b=factonal(i);x=sqrt( n);sum=sum+ti/b;endsum=sum+1;11、編程求滿足m k2k10000sum最小m值。sum=0;13、編寫一個函數(shù)，計算下面函數(shù)的值，給 出標量x的值，調(diào)用該函數(shù)后， 的值。返回 ym=2;fun cti on y=myfun 1(x)a=1;while asin x,for i=1:my(x)x,sum=sum+2i;x 6,if sum>10000選擇一些數(shù)據(jù)測試你編寫的函數(shù)。a=0;fun cti on y=myfu

11、 n1(x)endif x<=0endy=s in(X);m=m+1;elseif x>0& x<=3endy=x;elseif x>3y=-x+6;12、思考題end已知 y 和 t數(shù)關(guān)系：y(t) 1t2/2! t3/3!求下面表格中與對應(yīng)的y值t=input('請輸入t值:')sum=0;tyt20 / 20!14、編寫一個函數(shù)求向量x中元素的平均for i=1:20值、最大值、最小值、均方根值。max_x=max_x; % 可省略endendfun ctio nm_x,max_x,min_x,rms_x=myfu n2(x)if max_

12、x<x(k)min_x;%求最小值if x(1)<x(2)方均根值(Root Mean Square)的計算公式為：min_x=x(1);rms1 NNiixi2elsemin_x=x(2);用下面數(shù)據(jù)測試你寫的函數(shù):end(1)x=sin(0:6*pi)for k=3:n(2)x=rand(1,200)，得到的 x 為 200 個(0,1)之間均勻分布的隨機數(shù)。if min_x>x(k)min_x=x(k);fun ctio nm_x,max_x,min_x,rms_x=myfu n2(x)%求平均值elsemin_x=min_x; %可省略sum_x=sum(x); %向

13、量元素求禾口endendm,n=size(x); % 最好用 n=length(x);%求均方根值m_x=sum_x/n;%求最大值 采用逐個比較方式sum_x2=0;for k=1: nif x(1)>x(2)max_x=x(1);elsesum_x2=sum_x2+x(k).A2;rms_x=sqrt(sum_x2/n);max_x=x(2);endendm_x;for k=3:nmax_x;max_x=x(k);rms_x; %按照函數(shù)值行參順序輸出結(jié)果else15、編寫一個函數(shù)，給出一個向量XX1*,Xn，生成如下范德蒙矩陣。Xn統(tǒng)計迭代次數(shù)。Xn22Xnfun ctio n v

14、=myva nder(x)注：X2 a 2x2x2X aX 2 2xxXn 1-X12X1X2X；Xn2Xnfun cti on x,n=sqrt_a(a)n 1X1n 1Xnx=;例如:>>v=m yvan der(2 3 4 5)for k=1:100m=x;end得v=x=x/2+a/(2*x);2345if abs(x-m)<=10A(-5)491625break827 64125endfunction v=myva nder(x)endv1=va nder(x); %生成范德蒙矩陣x;v2=v1'n=k;v=flip ud(v2); %實現(xiàn)矩陣上下翻轉(zhuǎn)s=(

15、x-sqrt(a);if s<=10A(-5)16、思考題disp('正確');elsedisp('錯誤');Ja編寫程序，用如下迭代公式求wa，a的值 分別為：3,17, 113。迭代的終止條件Xn 1Xn10 5，迭代初值Xo1.0 ,迭代次數(shù)不超過100次。分別對迭代結(jié)果和準確值進行比較，并17、在同一坐標系下繪制下面三個函數(shù)在 t 0，4的圖象。y 1 t>> y=fun 1(x);>> plot(x,y)y 2 VTy 34 e 0.1t sin( t)y=y,6-x0;>> y=s in(z);>>

16、; t=0:4* pi;19、用comp ass函數(shù)畫下面相量圖>> y1=t;uaub>> y2=t;cos(-2* pi/3)+si n(-2* pi/3)*iuc=cos(2* pi/3)+si n(2* pi/3)*i;>> y2=sqrt(t);comp ass(ua,ub,uc,ua-ub,ub-uc,uc-ua)>> y3=4* pi*e ps*t).*si n(t);>> plot(t,y1,t,y2,t,y3)>> ua=1;>> ob=cos(-2* pi/3)>> ob=cos

17、(-2* pi/3)+si n(-2* pi/3)*i;18、編寫程序，選擇合適的步距，繪制下面函數(shù)在區(qū)間-6，6中的圖象。>> ub=cos(-2* pi/3)+si n(-2* pi/3)*i;y(x)sin x,x,>> uc=cos(2* p“3)+si n(2* pi/3)*i;>>x 6,comp ass(ua,ub,uc,ua-ub,ub-uc,uc-ua)function y=fun 1(x)y=;20、三維空間曲線繪制for x0=xz=0:4* pi;if x0<=0x=cos(z);y=y,s in (x0);y=s in(z);

18、elseif x0>0&&x0<=3plot3(x,y,z)y=y,x0;>> z=0:4* pi;elseif x0>3>> x=cos(z);end>> p Iot3(x,y,z)end>> x=-6:6;21、用mesh或surf函數(shù)，繪制下面方程所表示的三維空間曲面， x 和 y 的取值范 圍設(shè)為 -3 ，3 。while y=0i=i+1;>> x=-3:3;x3=(x0+x1)/2;>> y=-3:3;y=2*eps(-x3)-sin(x3);>> x,y=mesh

19、grid(x,y);y0=2*eps(-x0)-sin(x0);>> z=-x.*x/10+y.*y/10;y1=2*eps(-x1)-sin(x1);>> mesh(x,y,z)if y*y0<0>> z=x.*x/10+y.*y/10;x0=x3;>> mesh(x,y,z)elseif y*y1<0x1=x3;22、思考題elseif y=0用對分法求解方程 2e1 內(nèi)的解， 分次數(shù)。并驗證，sinx在0,在程序中統(tǒng)計出對x3=x3;endif i=1000000提示：先將原方程轉(zhuǎn)化成f (x)2e xsinx 0 的形y=0;

20、式。end對分法的基本思想是：一個一元方程 f(x)=O ,若 f(x1)*f(x2)<0,則在x1,x2區(qū) 間內(nèi) 有實 數(shù)解。 取該區(qū) 間的 中 點 xm=(x1+x2)/2 ，判定 f(x1) 和 f(x2) 二 者中哪一 個與 f(xm) 異 號， 若 f(x1)*f(xm)<0 ，則解存在的區(qū)間縮小 為x1,xm，否則解存在的區(qū)間縮小為 xm,x2 。重復(fù)這樣的步驟，直到區(qū)間 的長度小于一個可以接受的小數(shù) (比如 1e-10 ), 則認為中點即是原方程的解。End233x5求代數(shù)4x4 7x3 2x2 9x方程12 0 的 5并將其用星號( * )標記在復(fù)平面圖個根，上。(

21、用 roots 和 plot 函數(shù))>> x0=0;>> a=3,4,7,2,9,12;y=1;>> plot(y,'*')i=0;x1=1;>> y=roots(a)求代數(shù)方程X510的5個根，并將其用25、星號（*）標記在復(fù)平面圖上。（用roots 和plot函數(shù)）程序:>> a1=1,0,0,0,0,-1;>> y1=roots(a1)>> plot(y1,'*')已知R=50歐姆，U=4V二極管 D正向電流 與電壓的關(guān)系為：24、求下面函數(shù)在,4區(qū)間內(nèi)的過零點。fzero 函）（用Udq 1IIkt1 d1 se>>y2=x.A3-2*x.A2.*si n(x)+5*x.*cos(x )+1./x;求此電路中的電流Id和二極管正向電壓 Ud （要求用fsolve函數(shù)求解）程序: