MATLAB實驗題答案

1、1、求以下變量的值，并在 MATLAB中驗證。 result8=a.c Attempt to reference field of non -structure array. result9=a.b result9 = ( 2 ) b = a , a , a ；a b = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 result10=a.A2 result10 = 16 25 36 result11=2.a result11 = 2、下列運算是否合法，為什么？如合法， 結(jié)果是多少？ result2=a*b Error using * Inner mat

2、rix dimensions must agree. 16 32 64 result3=a+b 3、 用MATLAB求解下面的的方程組。 x1 result3 = 15 11 x2 result4=b*d result4 = 31 22 22 40 49 13 result5=b;c*d result5 = (1) 11 x3 13 x4 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 yz1 x 2y z w 8 31 22 22 (2) 2x y 3w 3 40 49 13 J -5 -8 7 3x 3y

3、 5z 6w 5 result6=a.*b result6 = A1=1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6 B2=1;8;3;5 x2=inv(A1)*B2 -3 15 30 result7=a./b result7 = 15 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000 4、已知 11 13 (1)求矩陣A的秩(rank) 7、分別用 if 和 switch 語句實現(xiàn)，將百分 (2) 求矩陣 A 的行列式 (determinant) 制成績轉(zhuǎn)換為成績等級 A、B、C、D、E。 求矩陣A的逆(inverse) 其中 90

4、100 分為 A,8089 分為 B, 7079 (4) 求 矩 陣 A 的 特 征 值 及 特 征 向 量 分為 C，6069 分為 D，60 分以下為 E。 (eigenvalue and eigenvector) 對超出百分制范圍的成績, 給出錯誤提 A3=7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13 示信息。 r=rank(A3) if 結(jié)構(gòu)程序： b=inv(A3) x=input(please enter score=); a=det(A3) if x=90 disp(E) for m=-10:10 else m1=m1+2m; disp(error

5、) end end m1 switch 結(jié)構(gòu)程序： m1 = x=input(please enter score=); 2047.9990234375 6、 求分段函數(shù)的值。 x2 x2 5x x2 if 語句實現(xiàn)，算出下列表中 值。 x 對應(yīng)的 y x=input(enter x=); if x=0 a b=nu m2str(floor(ra nd(1)*90+10); b if x=+ 10 x2 分別取 x的值分別是多少? a=i npu t(請輸入數(shù)值 n=0; for m=1:100 100、 1000、 10000時，求 ) y=a+b; elseif x=- y=a-b; el

6、seif x=* y=a*b; elseif x=/ y=a/b; else n=n+1/(m*m); end n=6* n; x=sqrt (n); x dis p( error) end 11、編程求滿足 最小m值。 9、啟動 MATLAB后，點擊 File|New|M-File ， 啟動MATLAB的程序編輯及調(diào)試器 (Editor/Debugger )，編輯以下程序， 點擊File|Save保存程序，注意文件名 最好用英文字符。點擊Debug|Run運行 程序，在命令窗口查看運行結(jié)果，程序 如有錯誤則改正。 注：數(shù)論中一個有趣的題目：任意一個正整 數(shù)，若為偶數(shù)，則用2除之，若為奇數(shù)，

7、則與3相乘再加上1。重復(fù)此過程，最 終得到的結(jié)果為1。 n=input(請輸入 n 值:); a=n; while n1 if rem( n,2)=0 n=n/2; else n=3*n+1; end a=a ,n ; end m k k12k sum=0; m=2; a=1; while a for i=1:m sum=sum+2i; if sum10000 a=0; end end m=m+1; end 12、思考題 已知 y 和 t y(t) 1 t2 10000的 數(shù)關(guān)系： /2! t3/3! t20 / 20! t 0.2 0.4 0.6 0.8 1.0 y 求下面表格中與 t對應(yīng)的

8、y值 t=input(請輸入 t 值：) sum=0; for i=1:20 b=factorial(i); end %求最大值： sum=sum+1; if x(1)x(2) sum max_x=x(1); else 13、編寫一個函數(shù)，計算下面函數(shù)的值，給 max_x=x(2); 出標(biāo)量x的值，調(diào)用該函數(shù)后，返回 y end 的值。 for k=3: n fun ctio n y=myfu n1(x) if max_xx(k) max_x=x(k); sum=sum+tAi/b; else m_x=sum_x/ n; 采用逐個比較方式 Sin x, y(x) x, x 6, 選擇一些數(shù)據(jù)測

9、試你編寫的函數(shù)。 fun ctio ny=myfu n1(x) if x0 end max_x=max_x;%可省略 end end %求最小值 if x(1)x(k) min_x=x(k); else min_x=min_x;%可省略 rms 1 N N i ixi 14、編寫一個函數(shù)求向量 x中元素的平均值、 最大值、最小值、均方根值。 fun ctio n m_x,max_x,min_x,rms_x=myfu n2(x) 方均根值(Root Mean Square)的計算公式為: end end %求均方根值 sum_x2=0; for k=1: n sum_x2=sum_x2+x(k)

10、.A2; rms_x=sqrt(sum_x2/n); end m_x; 用下面數(shù)據(jù)測試你寫的函數(shù)： (1) x=sin(0:0.01:6*pi) (2) x=rand(1,200)，得到的 x 為 200 個(0， 1)之間均勻分布的隨機數(shù)。 fun ctio n m_x,max_x,min_x,rms_x=myfu n2(x)% 求平均值 sum_x=sum(x); %向量元素求和 m,n=size(x);%最好用 n=length(x); max_x; min_x; rms_x;%按照函數(shù)值行參順序輸出結(jié)果 15、編寫一個函數(shù)，給出一個向量 x XX2，Xn，生成如下范德 蒙矩陣。 fun

11、 cti on v=myva nder(x) end X1 2 X1 X2 2 X2 Xn 2 Xn n 1 X1 n 1 X2 n 1 Xn 例如: 得v= v=m yvan der(2 3 4 5) x; n=k; s=(x-sqrt(a); if s=10A(-5) dis p(正確); else disp(錯誤); end 16 25 17、在同一坐標(biāo)系下繪制下面三個函數(shù)在 t 0，4 的圖象。 27 64 125 function v=myva nder(x) v1=va nder(x);%生成范德蒙矩陣 v2=v1; v=flip ud(v2);%實現(xiàn)矩陣上下翻轉(zhuǎn) t 4 e 0.

12、1t sin( t) 16、思考題 編寫程序，用如下迭代公式求嘉，a的值 分別為：3，17, 113。迭代的終止條件 10 5 ，迭代初值 Xn 1Xn Xo 1.0 ,迭代次數(shù)不超過100 分別對迭代結(jié)果和準(zhǔn)確值進行比較， 次。 統(tǒng)計迭代次數(shù)。 Xn Xn 2 2Xn 注：x2 a 2x2 X2 fun ctio nx, n=sqrt_a(a) x=1.0; for k=1:100 m=X； x=x/2+a/(2*x); if abs(x-m) t=0:0.1:4* pi; y1=t; y2=tA(0.5); y2=sqrt(t); y3=4*pi*eps( -0.1*t).*sin(t);

13、 plot(t,y1,t,y2,t,y3) 18、編寫程序，選擇合適的步距，繪制下面 函數(shù)在區(qū)間-6，6中的圖象。 Sin x, y(x) x, X 6, fun cti on y=fun 1(x) a=; fr x0=X n 1 2x if x00 end end x=-6:0.1:6; y=fun 1(x); 一個一元方程 f(x)=0 ， ，則在 x1,x2 區(qū)間 取該區(qū)間的中 23 求代數(shù)方程 plot(x,y) 19、用 compass 函數(shù)畫下面相量圖 ua =1;ub cos(-2*pi/3)+sin(-2*pi/3)*i uc=cos(2*pi/3)+sin(2*pi/3)*i

14、; compass(ua,ub,uc,ua-ub,ub-uc,uc-ua) ua=1; ob=cos(-2*pi/3) ob=cos(-2*pi/3)+sin( -2*pi/3)*i; ub=cos(-2*pi/3)+sin( -2*pi/3)*i; uc=cos(2*pi/3)+sin(2*pi/3)*i; compass(ua,ub,uc,ua -ub,ub -uc,uc-ua) 20、三維空間曲線繪制 z=0:0.1:4*pi; x=cos(z); y=sin(z); plot3(x,y,z) z=0:0.1:4*pi; x=cos(z); y=sin(z); plot3(x,y,z)

15、21、用mesh或surf函數(shù)，繪制下面方程所 表示的三維空間曲面， x 和 y 的取值范 圍設(shè)為 -3 ， 3 。 x=-3:0.01:3; y=-3:0.01:3; x,y=meshgrid(x,y); z=-x.*x/10+y.*y/10; mesh(x,y,z) z=x.*x/10+y.*y/10; mesh(x,y,z) 式。 對分法的基本思想是： 若 f(x1)*f(x2)0 內(nèi)有實數(shù)解。 xm=(x1+x2)/2 ，判定 f(x1) 和 f(x2) 者 中 哪 一 個 與 f(xm) 異 號 ， f(x1)*f(xm) x0=0; x1=1; y=1; i=0; while y=

16、0 i=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*y00 x0=x3; elseif y*y1 a=3,4,7,2,9,12; y=roots(a) plot(y,*) 5 求代數(shù)方程 x 1 0 的 5 個根，并將其用 星號（*）標(biāo)記在復(fù)平面圖上。（用roots 和plot函數(shù)） 程序： a仁1,0,0,00-1; y1=roots(a1) plot(y1,*) 24、求下面函數(shù)在0.5,4區(qū)間內(nèi)的過零點。 （用 fzero 函） x=fsolve(f

17、u n,0, op timset(Dis play,off) y3=5*10(-11)*ex卩(1.6*對(1.38*10人(-4)* 27)-1) 函數(shù)： fun cti on y=fun(x) y=5*10A( -11)*exp(1.6*x1.38*10A( -4)*27)-1) +x-4; p 0.0 1.1 2.1 2.8 4.2 5.0 6.1 u 10 11 13 14 17 18 22 26、實驗數(shù)據(jù)處理：已知某壓力傳感器的測 試數(shù)據(jù)如下表 p為壓力值，u為電壓值，試用多項式 3 2 1 f (x) x 2x sin(X) 5x cos(x) x x=0.5:0.01:4; y2=x.A3-2*x.A2.*si n(x)+5*x.*cos(x)+1./x; plot(x,y2) 函數(shù)： fun ctio n y2=fz(x) y2=x.A3-2*x.A2.*si n(x)+5*x.*cos(x)+1./x; 分析： 根據(jù)圖像可知零點在 x=1和x=3附近 所以 u(p) ap3 bp2 cp d 來擬合其特性函數(shù)，求出a,b,c,d ，并 把擬合曲線和各個測試數(shù)據(jù)點畫在同 一幅圖上。 程序： fzero(fz,1) ans = 1.