1、Matlab程序設計上機報0.982924(3) y 2a 1 -O.8333,其中 o=3.32 ,蘆-7.9.a=3.32;b=-7.9;y=2*pi*aA2*(1-pi/4)*b-(0.8333-pi/4)*a)1. 求以下表達式的值.(1) w .、2 (1 0.34245 10 6)>> w=sqrt(2)*(1+0.34245*10A(-6)1.4142(4)z - e2t In t . 1 t2,其中221 3it=.50.65-128.4271(2) xbe 2eabctan (b c) a,其中>> t=2,1-3i;5,-0.65;z=0.5*exp

2、(2*t)*log(t+sqrt(1+t*t)a=3.6,b=5,c=-9.8>> a=3.5;b=5;c=-9.8;x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2)/(ta n( b+c) +a)1.0e+040.0057 - 0.0007i0.0049 - 0.0027i1.9884 - 0.3696i1.7706 - 1.0539i2.154831A078 ,B2533617320x =123726-15737>> AA2-B+eye(3)ans =求以下表達式的值(1) A+6B 和 A2-B+l(2) A*B、A.*B 和 B*A(3) A/

3、B 及 B/A(4) A,B和A(1,3,:);BA2>> A=-1,5,-4;0,7,8;3,61,7;>> B=8,3,-1;2,5,3;-3,2,0; >> A+6*B-18 -217172253310921867526>> A*Bans =ans =4723-10141416-105121125328180>> A.*Bans =-815403524-91220>> A.*B23100.77804145655A32503269.54543.14完成以下操作：ans =-815403524-91220(1)輸出 A

4、在10,25范圍內的全部元素.(2) 取出A前3行構成矩陣B,前兩列構成 矩陣C,右下角3X2只矩陣構成矩陣D, B 與 C 的乘積構成矩陣 E.(3) 分別求表達式 E<D、E&D、E|D和E|D 的值.>> find(A>=10&A<=25)ans =>> A/Bans =151.2234-0.92552.9787>>B=A(1:3,:)-0.94682.3511-0.95744.61703.872313.8936B =>> B/A23.000010.0000-0.7780 0ans =41.0000-45.

5、000065.0000 5.000032.00005.00000-5.3173-3.56610.894232.0000-1.2768-0.47430.2411>>C=A(:,1:2)2.12291.3173-0.2924>> A,BC =ans =23.000010.000041.0000-45.0000-15-483-132.00005.00000782 536.0000-9.54003617-320>> “1,3,:)砂2>>D=A(2:4,3:4)ans =D =-15-465.00005.00003617032.00007337154.0

6、0003.1400173713>> E=B*C-20193.E =1.0e+03 *>> H=hilb(5);P=pascal(5);0.9141-0.2239Hh=det(H)1.20802.71231.1330-0.2103Hh =E<D3.7493e-12ans =>> Hp=det(P)01Hp =00011>> E<D>> Th=cond(H) ans =Th =0 10 0 4.7661e+05 0 1>> Tp=cond(P) >> E&DTp = ans =8.5175e+0

7、31 1Tp 比 Th 好,更接近于 1 .0 15. ：1 129 618A 20 512>> E|D885ans =求A的特征值及特征向量,并分析其數學意義.1 1>> A=-29,6,18;20,5,12;-8,8,5;1 1V,D=eig(A)1 1>> E|DV =0.71300.28030.2733-0.6084-0.78670.87250.34870.55010.4050D =-25.3169000-10.51820ans =0 01 00 04產生5階希爾伯特矩陣H和5階帕斯卡矩 陣P,且求其行列式的值Hh和Hp以及它們 的條件數Th和Tp,

8、判斷哪個矩陣性能更好, 為什么？0016.8351設A為n階矩陣,假設存在常數A及非零的n維向 量x,使得Ax=入x,那么稱 尿矩陣A的特征值,x是A屬于特征值入 的特征向量1.從鍵盤輸入一個3位整數,將它反向輸出.>> a=input'輸入一個3位整數:','s'fliplra輸入一個3位整數:345ans =5432. 輸入一個百分之成績,要求輸出等級A、B、C、D、E.其中 90100 分為 A,8089分為B,70分79分為C,60分69分為D,60分以下 為E.要求：1分別用if語句和switch語句實現.2輸入百分制成績后要判斷該成績的合

9、 理性,對不合理的成績應輸出出錯信息.>> m=input'請輸入一個百分制的成績 ：'；if m>=0 & m<=100switch fixm/10case 9,10disp'這個成績所屬的成績等級:A'case 8disp'這個成績所屬的成績等級:B'case 7disp'這個成績所屬的成績等級：C'case 6disp'這個成績所屬的成績等級:D'otherwisedisp'這個成績所屬的成績等級:E'endelsedisp'輸入有誤'enddi

10、sp('這個成績所屬的成績等級:A');elseif m>=80&m<=89disp('這個成績所屬的成績等級:B');elseif m>=70&m <=79disp('這個成績所屬的成績等級:C');elseif m>=60&m<=69disp('這個成績所屬的成績等級：C');selseif m>=0&m <=59disp('這個成績所屬的成績等級:E');elsedisp('錯誤:輸入的成績不合理');end3.

11、輸入20個數,求其中最大數和最小數. 要求分別 用循環結構和調用 matlab的max函數、min函數來 實>>for i=1:20s=strcat('輸入第',num2str(i),'個數字');a(i)=i nput(s);endb=a(1);c=a(1);for i=1:20if a(i)>bb=a(i);endif a(i)<cc=a(i);endendbc>> max(a)>> mi n(a)0.3a0.3ae e0.3 a5. ysin(a 0.3) ln,當 a取2 2-3.0、-2.9、-2.8&g

12、t; -2.8> 2.8、2.9、3.0 時,求各點的函數值.要求分別用順序結構和循環結構實現.>> a=-3:0.1:3;b=(exp(0.3.*a)-exp(-0.3.*a)./2);c=si n( a+0.3);y=b.*c+log(0.3+a)./2)>> m=input'請輸入一個百分制的成績：'if m>=90&m<=100>> a=-3:0.1:3;for i=1:61y=(exp(0.3*a(i)-exp(-0.3*a(i)/2*si n(a(i)+0.3)+log(0.3+a(i)/2)End5當n

13、分別取100、1000、10000時,求以下格式的值：要求分別用順序結構和向量運算(使用sum函數)來實現.(1)12213"n2(2-)6>> A=0,0,0;n=100,1000,10000; for i=1:3for j=1: n(i)A(i)=1/(jA2)+A(i); end1.56691.57041.57086.建立5X 6矩陣,要求輸出矩陣第n行元素.當n值超過矩陣的行數時,自動轉為輸出矩陣最后一行 元素,并給出出錯信息.>> A=ra nd(5,6)n=i nput('輸入第幾行？)；if n>5n=5;disp('超出范

14、圍,輸出最后一行的元素)；A(n,:)elseA(n,:)endendA1.繪制以下曲線1.63501.64391.6448>> A=0,0,0;n=100,1000,10000;for i=1:3j=1:1: n( i);A(i)=sum(1./j.A2);endA(1) y3x3!>> x=0:pi/100:2*pi;>> y=x-x.A3/factorial(3);>> plot(x,y)1.63501.64391.6448(2)2 24 466(2n)(2 n)133 557 ."(2n 1)(2n 1)2>> A=

15、1,1,1;n=100,1000,10000;for i=1:3for j=1: n(i)A(i)=A(i)*(2*j)A2)/(2*j-1)*(2*j+1);endend>> A(2) x2 2y2 64>> ezplot('x.A2+2*y.A2-64',-8,8,-6,6)眇*2= Q理-_R II、III B-64-202& B B胃2,設y亠,t ,1 e在同一圖形窗口采用子圖的形式繪制條形 圖、階梯圖、桿圖、和對數坐標圖等不同圖 形,并對不同圖形加標注說明>> t=-pi:pi/100:pi;>> y=1./(

16、1+exp(-t);>> subplot(2,2,1), stairs(t,y),title '階梯圖'>> t=-pi:pi/10:pi;>> y=1./(1+exp(-t);>> subplot(2,2,1), stairs(t,y),title '階梯圖'>> subplot(2,2,2), bar(t,y),title '條形圖'>> subplot(2,2,3), stem(t,y),title '桿圖'>> subplot(2,2,4)

17、, semilogx(y),grid on,title '對 數坐標圖(2)-JJ -2 C 73繪制以下及坐標圖15 cos 4>> x=-2*pi:pi/20:pi;>> y=5.*cos(x)+4;>> polar(x,y)5si n2cos ' 3>> x=-pi/3:pi/30:pi/3;>> y=5.*(si n(x)42)./cos(x);>> polar(x,y)4.繪制以下三維圖形tx e 20 costt(1) y e20sint,0 t 2.>> t=0:pi/100:2*

18、pi; >>x=exp(-t/20).*cos(t); >>y=exp(-t/20).*si n(t); >>z=t;>>plot3(x,y,z)co =21.建立一個圖形窗口,使之背景顏色為紅色, 并在窗口上保存原有的菜單項,而且在按下 鼠標左鍵后顯示出Left Button Pressed字樣.>> co=figure('Color','r','Wi ndowButton Dow nFcn ','axis off;text(0.1,0.5,"Left Button

19、Pressed","fontsize",30)')(2) z=5,|x0,|y|<5.要求應用插值著色處 理.>> x,y=meshgrid(-10:0.1:10);>>z=x;>>z(:,:)=5;>>shad ing in terp>> i=fin d(abs(x)>5&abs(y)>5); >>z1=z;>>z1(i)=NaN;>>surf(x,y,z1);>>shad ing in terp;5播放一個直徑不斷變化的

20、球體.>> n=30;x,y,z=sphere; m=moviei n(n);for j=1: nsurf(j*x,j*y,j*z);m(:,j)=getframe;>> end>> movie(m,2)2.利用圖形對象回執以下曲線,要求先利用默認屬 性繪制曲線,然后通過圖形句柄操作來改變曲線的 顏色、線性和線寬,并利用文本對象給曲線添加文 字標注.1 : 2(1) y In x . 1 x2>> x=0:0.1:20;y=Iog(x+sqrt(1+x.A2)./2;fi=plot(x,y);set(fi,'color','

21、;r','li nestyle',':','li newidth',1)1 fl -in.1 2 -H1-rg -J0i2 -2 I t 9 iu14 仙 18£(2)x t2y 5t2>> t=0:0.1:20;x=t.A2;y=5.*t.A2;fi=plot3(t,x,y);set(fi,'color','g','l in estyle','-.','li newidth',2)4.利用圖形對象繪制以下三維圖形,要求對圖形進 行光照處理.(1) z x y 5sin(xy)>> x,y=meshgrid(0:0.1:10); z=x.*2+y.*2-5*si n(x.*y);>> surf(x,y,z); light('style','local')(2)zy3>> x,y=meshgrid(0:0.1:10);>> z=y.A3;>> surf(x,y,z); light('posi'