Matlab上機(jī)實驗答案

1、Matlab 上機(jī)實驗答案實驗一 MATLAB 運(yùn)算基礎(chǔ)1.先求下列表達(dá)式的值，然后顯示 MATLAB 工作空間的使用情況并保存全部變量z1=2*sin(85*pi/180)/(1+exp(2)z1=0.2375212i22)z2=-ln(x+x2),其中x=J22-0.455x=21+2i;-0.455;z2=1/2*log(x+sqrt(1+xA2)z2=0.7114-0.0253i0.8968+0.3658i0.2139+0.9343i1.1541-0.0044i(1)乙=2sin8501e2e0.3ao3a0.3a(3)4=2sin(a0.3)In2a=-3.0,-2.9,|,2.9,

2、3.0a=-3.0:0.1:3.0;z3=(exp(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=Columns1through50.7388+3.1416i0.7696+3.1416i0.7871+3.1416i0.7913+3.1416i0.7822+3.1416iColumns6through100.7602+3.1416i0.7254+3.1416i0.6784+3.1416i0

3、.6196+3.1416i0.5496+3.1416iColumns11through150.4688+3.1416i0.3780+3.1416i0.27753.1416i0.1680+3.1416i0.0497+3.1416iColumns16through20-0.0771+3.1416i-0.2124+3.1416i-0.5104+3.1416i-0.6752+3.1416iColumns21through25-0.8536+3.1416i-1.0497+3.1416i-1.5271+3.1416i-1.8436+3.1416iColumns26through30-2.2727+3.14

4、16i-2.9837+3.1416i-3.0017+0.0000i-2.3085+0.0000iColumns31through35-1.8971+0.0000i-1.5978+0.0000i-1.1531+0.0000i-0.9723+0.0000i-0.3566+3.1416i-1.2701+3.1416i-37.0245+0.0000i-1.3575+0.0000iColumns36through40-0.8083+0.0000i-0.6567+0.0000i-0.3819+0.0000i-0.2561+0.0000iColumns41through45-0.1374+0.0000i-0

5、.0255+0.0000i0.1766+0.0000i0.2663+0.0000iColumns46through50Columns51through55Columns56through60-0.5151+0.0000i0.0792+0.0000i0.3478+0.0000i0.4206+0.0000i0.4841+0.0000i0.5379+0.0000i0.5815+0.0000i0.6145+0.0000i0.6366+0.0000i0.6474+0.0000i0.6470+0.0000i0.6351+0.0000i0.4774+0.0000i0.4126+0.0000iColumn61

6、0.3388+0.0000ir2_t20t1(4)z4=Jt2_11Mt2,其中t=0:0.5:2.52一，一一t-2t+12tt=0:0.5:2.5;z4=(t=0&t=1&t=2&t A=1234-4;34787;3657; B=13-1;203;3-27; A+6*Bans=1852-10467105215349 I=eye(3); A-B+Ians=1231-3328840671A*Bans=684462309-72596154-5241 A.*Bans=1210246802619-13049 AA3ans=372262338244860424737014918

7、860076678688454142118820 A.A3ans=172839304-643930434365850327274625343 A/Bans=16.4000-13.60007.600035.8000-76.200050.200067.0000-134.000068.0000 BAans=109.4000-131.2000322.8000-53.000085.0000-171.0000-61.600089.8000-186.2000 A,Bans=1234-413-13478720336573-27 A(1,3,:);BA2ans=1234-436574511101920-5403

8、,設(shè)有矩陣 A 和 B(1)求它們的乘積 Co(2)將矩陣 C 的右下角 3X2 子矩陣賦給 D(3)查看 MATLAB 工作空間的使用情況。A=(reshape(1:25,5,5);B=3016;17-69;023-4;970;41311;C=A*BC=9315077258335237423520397588705557753890717一123451678910111213 1415161718 1920212223 2425130161176 9B=023-497041311-A=D=C(3:5,2:3)(1)求100,999之間能被 21 整除的數(shù)的個數(shù)。(2)建立一個字符串向量，刪除

9、其中的大寫字母n=100:999;l=find(rem(n,21)=0);length(l)ans=520397705557890717whosNameSizeBytesClassA5x5200doubleB5x3120doubleC5x3120doubleD3x248double4.完成下列操作：Attributes43ch=aegbBOIEG0je23RGnc;wz=find(ch=A&chch(wz)=ch=aegb0je23nc實驗二 MATLAB 矩陣分析與處理cERRS陣和對角陣，試通過數(shù)值1t算驗證A2=J9PS2_E=eye(3);R=rand(3,2);O=zeros

10、(2,3);S=diag(1:2);1.設(shè)有分塊矩陣 A=,EOR32S22,其中 E、R、O、S 分別為單位矩陣、隨機(jī)矩陣、零矩A=E,R;O,S1.0000000.81470.913401.000000.90580.6324001.00000.12700.09750001.0000000002.0000AA2ans=1.0000001.62942.740101.000001.81161.8971001.00000.25400.29260001.0000000004.0000E,R+R*S;O,SA2ans=1.0000001.62942.740101.000001.81161.897100

11、1.00000.25400.29260001.0000000004.0000AA2=E,R+R*S;O,SA2ans=1111111111111111111111111cERRS通過驗證，矩陣A=|PS212.產(chǎn)生 5 階希爾伯特矩陣 H 和 5 階帕斯卡矩陣 P,且求其行列式的值 Hh 和 Hp 以及它們的條件數(shù) Th 和 Tp,判斷哪個矩陣性能更好。為什么？H=hilb(5)1.00000.50000.33330.25000.20000.50000.33330.25000.20000.16670.33330.25000.20000.16670.14290.25000.20000.16670

12、.14290.12500.20000.16670.14290.12500.1111P=pascal(5)P=111111234513610151410203515153570成立。Hh=3.7493e-12Hp=det(P)Hp=1Th=cond(H)Th=4.7661e+05Tp=cond(P)Tp=8.5175e+03答：5 階帕斯卡矩陣 P 的性能好。矩陣的性能是由條件數(shù)決定的，條件數(shù)越接近于 1 其性能就越好。由上機(jī)操作求得 Th=4.7661e+005,Tp=8.5175e+003。Tp 的值更接近于 1 則其性能要好。 所以 5 階帕斯卡矩陣 P 的性能好。3 .建立一個 5X5

13、矩陣，求它的行列式值、跡、秩和范數(shù)。A=rand(5)0.27850.97060.42180.03570.74310.54690.95720.91570.84910.39220.95750.48540.79220.93400.65550.96490.80030.95950.67870.17120.15760.14190.65570.75770.7060det(A)ans=-0.1322trace(A)ans=3.4127ans=norm(A)ans=3.27594 .已知-29618A=20512I885求 A 的特征值及特征向量，并分析其數(shù)學(xué)意義A=-29,6,18;20,5,12;-8,8

14、,5-2961820512-885V,D=eig(A)0.71300.28030.2733-0.6084-0.78670.87250.34870.55010.4050D=-25.3169000-10.518200016.8351在數(shù)學(xué)上,特別是線性代數(shù)中,對于一個給定的線性變換，它的特征向量（本征向量或稱正規(guī)正交向量）是這樣一個斗零的向量v：當(dāng)v經(jīng)過這個線性變換的作用之后得到的新向量（長度也許改變）仍然與原來的v保持在同一條線上。一個特征向量的長度在該線性變換下縮放的比例稱為其特征值（本征值）。如果特征值為正，則表示v在經(jīng)過線性變換的作用后方向也不變；如果特征值為負(fù)，說明方向會反轉(zhuǎn)；如果特征值

15、為0,則是表示縮回零點。但無論怎樣，仍在同一條直線上。5.下面是一個線性方程組12131(1)求方程的解。(2)將方程右邊向量元素 b3改為 0.53 再求解，并比較 b3的變化和解的相對變化。(3)計算系數(shù)矩陣 A 的條件數(shù)并分析結(jié)論。formatrat%用分?jǐn)?shù)格式顯示A=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/61/21/31/41/31/41/51/41/51/6format%恢復(fù)默認(rèn)格式b=0.95;0.67;0.520.95000.67000.5200X=AbX=1.20000.60000.6000b2=0.95;0.67;0.53b2=0.95000.6

16、7000.5300X2=Ab2X2=3.0000-6.60006.6000D=cond(A)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)程序設(shè)計1 .求分段函數(shù)的值。x2+x6x0且x#-3y=x25x+60Mx5x02及xr3x2-x-1其他用 if 語句實現(xiàn)，分別輸出 x=-5.0,-3.0,1.0,2.0,2.5,3

17、.0,5.0 時的 y 值。程序設(shè)計：clear;x=-5.0-3.01.02.02.53.05.0;ifx=0&x=90&score=80&score=70&score=60&score=0&score=100%判斷成績合理性switchfix(score/10)case9,10grade=A;case8grade=B;case7grade=C;case6grade=D;otherwisegrade=E;endgradeelsedisp(error)%成績不合理時輸出出錯信息end3 .硅谷公司員工的工資計算方法如下：(1)工作時數(shù)超過 120

18、 小時者，超過部分加發(fā) 15%。(2)工作時數(shù)低于 60 小時者，扣發(fā) 700 元。(3)其余按每小時 84 元計發(fā)。試編程按輸入的工號和該號員工的工時數(shù)，計算應(yīng)發(fā)工資。程序設(shè)計：clear;num=input(請輸入員工工號：)；time=input(請輸入員工工時數(shù)：)；iftime120wages=84*120+(time-120)*84*(1+0.15);運(yùn)行結(jié)果：請輸入員工工號：disp=員工工號：4 應(yīng)發(fā)工資：22404 .設(shè)計程序，完成兩位數(shù)的加、減、乘、除四則運(yùn)算，即產(chǎn)生兩個兩位隨機(jī)整數(shù)一個運(yùn)算符號，做相應(yīng)的運(yùn)算，并顯示相應(yīng)的結(jié)果。解：程序設(shè)計：clear;a=floor(ra

19、nd(1)*90+10)b=floor(rand(1)*90+10)c=input(,請輸入運(yùn)算符號：,s);ifc=+elsewages=time*84;enddisp=(員工工號：num2str(num)應(yīng)發(fā)工資：num2str(wages)請輸入員工工時數(shù):35s=a+belseifc=-s=a-belseifc=*s=a*belseifc=s=abelsedisp(error)end運(yùn)行結(jié)果：3077請輸入運(yùn)算符號s=-475 .建立 5X6 矩陣，要求輸出矩陣第 n 行元素。當(dāng) n 值超過矩陣的行數(shù)時，自動轉(zhuǎn)為輸出矩陣最后一行元素，并給出出錯信息。程序設(shè)計：clear;A=rand(

20、5,6)n=input(n=?)ifn=1&n=5B=A(n,1:6);elsedisp(erro!)B=A(5,1:6);endB運(yùn)行結(jié)果：A=0.77920.01190.52850.68920.91330.07820.93400.33710.16560.74820.15240.44270.12990.16220.60200.45050.82580.10670.56880.79430.26300.08380.53830.96190.46940.31120.65410.22900.99610.0046n=?3實驗四循環(huán)結(jié)構(gòu)程序設(shè)計二1111M11.根據(jù)二=工十震十三十川二，求兀的近似值

21、。當(dāng) n 分別取 100、1000、6122232n20.12990.16220.60200.45050.82580.106710000時，結(jié)果是多少？要求：分別用循環(huán)結(jié)構(gòu)和向量運(yùn)算(使用 sum 函數(shù))來實現(xiàn)。程序設(shè)計(一)：clear;forn=100,1000,10000sum=0;forx=1:nsum=sum+1/(x.A2);endresult=sqrt(6*sum)end運(yùn)行結(jié)果：result=3.1321result=3.1406result=3.1415程序設(shè)計(二)：clear;forn=100,1000,10000 x=1:n;result=sqrt(6*sum(1./x

22、./x)end運(yùn)行結(jié)果：result=3.1321result=3.1406result=3.1415111.2 .根據(jù) y=1+為 Ij+-,求:3 52n-1(1) y3 時的最大 n 值。(2)與(1)的 n 值對應(yīng)的 y 值。程序設(shè)計：y=0;n=1;whiley10A(-5)&n0)negtivenum=length(find(ffushu2請輸入任意一個復(fù)數(shù) A：4i-0.6536-0.7568i1.3863+1.5708i0.0000+27.2899i27.30822.一物理系統(tǒng)可用下列方程組來表示fm1cos6-sin601p1-01nsin60cos80 a2二mg0

23、mb-sin60 N1一000cos81. .N2一_m2gl從鍵盤輸入 mi、m2和。的值，求 ai、a2、Ni和 N2的值。其中 g 取 9.8,輸入。時以角度為單位。要求：定義一個求解線性方程組 AX=B 的函數(shù)文件，然后在命令文件中調(diào)用該函數(shù)文functiona1,a2,N1,N2=xxfcz(m1,m2,t)t=t*pi/180;A=mi*cos(t)-mi-sin(t)0;m1*sin(t)0cos(t)0;0m2-sin(t)0;00-cos(t)1;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(請輸入角度t:);a1,a2,N1,N2=xxfcz(m1,m2,t);a1a2N1N2運(yùn)行結(jié)果：xxfcz2請輸入 m1：5請輸入 m2:3請輸入角度 t:90al=9.8000a2=3.7505e-16N1=1.1251e-15N2=29.40003.一個自然數(shù)是素數(shù)，且它的數(shù)字位置經(jīng)過任意對換后仍為素數(shù)。例如 13 是絕對素數(shù)。試求所