matlab上機(jī)練習(xí)附答案

1、MATLAB上機(jī)練習(xí)1 .以下兩種說法對嗎？(1) MATLAB進(jìn)行數(shù)值的表達(dá)精度與其指令窗中的數(shù)據(jù)顯示精度相同。(2) MATLAB指令窗中顯示的數(shù)據(jù)有效位數(shù)不超過七位。2 .歷史指令窗所記錄的內(nèi)容與diary指令所產(chǎn)生的“日志”內(nèi)容有什么不同？DIARYfilenamecausesacopyofallsubsequentcommandwindowinputandmostoftheresultingcommandwindowoutputtobeappendedtothenamedfile.Ifnofileisspecified,thefile'diary'isused.DIA

2、RYOFFsuspendsit.DIARYONturnsitbackon.DIARY,byitself,togglesthediarystate.UsethefunctionalformofDIARY,suchasDIARY('file'),whenthefilenameisstoredinastring.3 .如何把用戶自己的“工作目錄”永久地設(shè)置在MATLAB的搜索路徑上？“位于搜索路徑上的目錄”與“當(dāng)前目錄”在MATLAB中的功用相同嗎？4 .如何向MATLAB工作空間輸入一個含有100個左右元素的一維或二維數(shù)值數(shù)組？用直接鍵入法？用數(shù)組編輯器？用M文件編輯器？5 .運(yùn)用

3、數(shù)組算術(shù)運(yùn)算符去掉下面程序里的for/end循環(huán)：x=11:15fork=1:length(x)z(k)=x(k42+2.3*x(k)A0.5;endx=11:15z1=x.A2+2.3*x.A0.56 .不使用數(shù)組算術(shù)運(yùn)算符，重寫下面的程序代碼：x=214z=1./(1+x.A2)x=2;k=1;whilei<=4,z2(k)=1/(1+iA2);i=i+1;x=x+1;end7 .某公司銷售電腦打印機(jī)的價格方案如下：()如果顧客只買一臺打印機(jī)，則一臺的基本價格為$150。()如果顧客購買兩臺以上打印機(jī)，則第二臺價格為$120。()第三臺以后，每臺$110。寫一段程序分別計(jì)算出購買1-

4、10臺打印機(jī)所需的錢數(shù)。打印機(jī)臺數(shù)可以在程序開始處指定，或通過input命令讀入。運(yùn)行程序，計(jì)算出購買10臺打印機(jī)的總價格。寫出程序，生成分別購買1-10臺打印機(jī)所需價格的圖表(使用fprintf命令輸出圖表,不允許手算)。x=1:10;ifx<=1y(x)=150*x;elseifx<=2y(x)=150+120*(x-1);elsey(x)=150+120+110*(x-2);endy(10)plot(x,y(x),'r*-)'8 .編寫一個函數(shù)M文件fun_es(x)，計(jì)算如下函數(shù)y=0.心/3一x2sinxfunctiony=es(x)y=0.5*exp(x

5、/3)-x.A2.*sin(x)其中參數(shù)可以為標(biāo)量，也可以為向量。在MATLAB里鍵入如下命令檢驗(yàn)此函數(shù)：fun_es(3)fun_es(123)1259 .定義一個矩陣m=312413_1m=125;312;413;如何用語句input將m讀入MATLAB?在MATLAB上驗(yàn)證此過程是否正確。在MATLAB上sum(m),max(m)和min(m)的結(jié)果分別是什么？1010.寫一段程序計(jì)算S=£ndn1(a)使用一個for/end循環(huán)，但不能使用數(shù)組運(yùn)算符和sum。(b)使用數(shù)組運(yùn)算符和sum，但不能使用任何for/end循環(huán)，S的值是多少？%第1哂aS=0;forn=1:10S=

6、S+n/(n+1);end%第1哂bn=1:10;a=n./(n+1);sum(a)11 .(a)寫一個函數(shù)M文件以計(jì)算方程ax2+bxc0的解。規(guī)定此函數(shù)形式為quad_rt(a,b,c),其中a,b和c可以為向量。quad_rt(a,b,c)x1=(-b+sqrt(b.A2-4*a.*c)./(2*a),x2=(-b-sqrt(b.A2-4*a.*c)./(2*a),(b)用a=3,b=1,c=1檢驗(yàn)此函數(shù)。(c)用a=312,b=1-49,c=13同檢驗(yàn)此函數(shù)。12 .編寫一個名為fun_xb的M文件以計(jì)算下面級數(shù)的和：232xxn1xf(x)=x-工-III(-1)23nx和n的值由i

7、nput給出。與手算的x=1,n=4時的結(jié)果相比，檢驗(yàn)此函數(shù)。上述級數(shù)是log(1+x)的截?cái)帑溈藙诹终归_，當(dāng)-1<x<1時，此級數(shù)收斂。利用該知識選定一些x值檢驗(yàn)此函數(shù)。例如：可取x=-0.5,0.5,n=1,2,3,5,10,20,50,將結(jié)果與log(1+x)進(jìn)行比較(當(dāng)x接近-1或1時，收斂將變得越來越困難)。13 .寫出下面MATLAB程序所描述的數(shù)學(xué)表達(dá)式。(1)forj=1:n(2)forj=1:nfori=1:my=x(j)*A(:,j)+yy(i)=A(i,j)*x(j)+y(i)endendend14 .寫出下面MATLAB程序所描述的數(shù)學(xué)表達(dá)式。(1)fori

8、=1:m(2)forj=1:nforj=1:mA(i,j)=A(i,j)+x(i)*y(j)A(:,j)=A(:,j)+y(endendend15 .用一條plot命令，在同一幅圖中繪制下面兩個函數(shù)的圖形:yJx-1)(x-2)(x-4)(x-5),0x(3-1)(3-2)(3-4)(3-5)yJx-2)(x-3)(x-4)(x-5),0，ix,6(1-2)(1-3)(1-4)(1-5)使用兩次plot命令和holdon重復(fù)上述任務(wù)。x=0:0.1:6;y1=(x-1).*(x-2).*(x-4).*(x-5)/(3-1)*(3-2)*(3-4)*(3-5);y2=(x-2).*(x-3).*

9、(x-4)*(x-5)/(1-2)*(1-3)*(1-4)*(1-5);plot(x,y1,'r-',x,y2,'g')plot(x,y1,'r-'),holdonplot(x,y2,'g')，holdoff2216 .繪制如下方程給出的橢圓：(x2)+2(y-x-1)=3,把y看做x的函數(shù)解方程。這樣做將得到兩個方程，一個是橢圓的上半部分，另一個是橢圓的下半部分。繪制這兩個關(guān)于x的方程，得到橢圓。x=-4:0.01:4;y1=-sqrt(3-(x-2).A2)/2)+x+1;y2=sqrt(3-(x-2).A2)/2)+x+1;

10、plot(x,y1,'g-',x,y2,'r-')17 .將兩個設(shè)計(jì)參數(shù)x和y定為0Vx<5和0<y<5。產(chǎn)品的成本為：一2.2f=x-8xy-6y-0.1xy50用mesh命令繪圖，近似地找出使成本最低的最優(yōu)參數(shù)，并寫出最低成本值。x,y=meshgrid(0:0.01:5,0:0.01:5);f=x.A2-8*x+y.A2-6*y-0.1*x.*y+50;mesh(f),m,n=find(f=min(min(f);x(m,n),y(m,n)18 .用contour命令重做上題。P98-1編寫程序Gauss-Siedel程序functiony=gs(A,b,x0)D=diag(diag(A);U=-triu(A,1);L=-tril(A,1);G=(D-L)U;f=(D-L)b;y=G*x0+f;n=1;whilenorm(y-x0)>=10e-6x0=y;y=G*x0+f;n=n+1;endynA=10,-1,0;-1,10,-2;0,-2,10;b=9;7;6;gs(A