1、MATLAB考試試題(1)產(chǎn)生一個1x10的隨機矩陣，大小位于(-55),并且按照從大到小的順序排列好！(注：要程序和運行結果的截屏)答案：a=10*rand(1,10)-5;b=sort(a,'descend')1 .請產(chǎn)生一個100*5的矩陣，矩陣的每一行都是123452 .已知變量：A='ilovematlab；'B='matlab'請找出：(A) B在A中的位置。(B) 把B放在A后面，形成C='ilovematlabmatlab'3 .請修改下面的程序，讓他們沒有for循環(huán)語句！A=123;456;789;rc=size

2、(A);fori=1:1:rforj=1:1:cif(A(i,j)>8|A(i,j)<2)A(i,j)=0;endendend4 .請把變量A=123;456;789寫到文件里(output.xls)，寫完后文件看起來是這樣的1234567895 .試從Yahoo網(wǎng)站上獲得微軟公司股票的2008年9月的每日收盤價。6 .編寫M文件，從Yahoo網(wǎng)站批量讀取60000.SH至600005.SH在2008年9月份的每日收盤價(提示：使用字符串函數(shù))。7 .將金牛股份(000937)2005年12月14日至2006年1月10日的交易記錄保存到Excel中,編寫程序將數(shù)據(jù)讀入MATLAB中

3、，進一步將數(shù)據(jù)讀入Access數(shù)據(jù)庫文件。8 .已知資產(chǎn)每日回報率為0.0025,標準差為0.0208,資產(chǎn)現(xiàn)在價值為0.8億，求5%水平下資產(chǎn)的試用MATLAB10天在險價值(Var)。9 .a=12345,b=a(1)*a(5)+a(2)*a(4)+a(3)*a(3)+a(4)*a(2)+a(5)*a(1).中最簡單的方法計算b,注意最簡單哦。1、求下列聯(lián)立方程的解3x+4y-7z-12w=45x-7y+4z+2w=-3x+8z-5w=9-6x+5y-2z+10w=-8求系數(shù)矩陣的秩；求出方程組的解。解：(1)>>a=34-7-12;5-742;108-5;-65-210;c=

4、4;-3;9;-8;b=rank(a)b=4(2) >>d=acd=-1.4841,-0.6816,0.5337,-1.2429即:x=-1.4841;y=-0.6816;z=0.5337;w=-1.2429為橫2、設y=cos0.5+(3sinx)/(1+xA2)把x=02兀間分為101點,畫出以x坐標，y為縱坐標的曲線；解：>>x=linspace(0,2*pi,101);y=cos(0.5+3.*sin(x)./(1+x.*x);plot(x,y)3、設f(x)=xA5-4xA4+3xA2-2x+6(1) Mx=-2,8之間函數(shù)的值(取100個點)，畫出曲線，看它

5、有幾個零點(提示：用polyval函數(shù))解：>>p=1-43-26;x=linspace(-2,8,100);y=polyval(p,x);plot(x,y);axisQ-2,8,-200,2300）;為了便于觀察，在y=0處畫直線，圖如下所示:與y=0直線交點有兩個，有兩個實根。（2）用roots函數(shù)求此多項式的根>>a=roots（p）a=3.0000,1.6956,-0.3478+1.0289i1.0289i4、在卜10,10;-10,10范圍內畫出函數(shù)的三維圖形,-0.3478-sin解:>>X,Y=meshgrid(-10:0.5:10);a=sq

6、rt(X.A2+Y.A2)+eps;Z=sin(a)./a;mesh(X,Y,Z);）用來表示該行為注釋行。）創(chuàng)建。matlab試卷，求答案一、選擇或填空（每空2分，共20分）1、標點符號（）可以使命令行不顯示運算結果,2、下列變量名中（）是合法的。(A)char_1;(B)x*y;(C)xy;(D)end3、為，步長為的向量，使用命令(4、輸入矩陣，使用全下標方式用（）取出元素“，”使用單下標方式用（）取出元素一5、符號表達式中獨立的符號變量為（）。6、M腳本文件和M函數(shù)文件的主要區(qū)別是（）和（）°7、在循環(huán)結構中跳出循環(huán)，但繼續(xù)下次循環(huán)的命令為（）。（A）return;（B）br

7、eak;（C）continue;（D）keyboad二、（本題12分）利用MATLAB數(shù)值運算，求解線性方程組（將程序保存為test02.m文件）三、（本題20分）利用MATALAB符號運算完成（將程序保存為test03.m文件）：（1）創(chuàng)建符號函數(shù)（2）求該符號函數(shù)對的微分；（3）對趨向于求該符號函數(shù)的極限；（4）求該符號函數(shù)在區(qū)間上對的定積分；（5）求符號方程的解。四、（本題20分）編寫MATALAB程序，完成下列任務（將程序保存為test04.m文件）：（1）在區(qū)間上均勻地取20個點構成向量；（2）分別計算函數(shù)與在向量處的函數(shù)值；（3）在同一圖形窗口繪制曲線與，要求曲線為黑色點畫線，曲線

8、為紅色虛線圓圈；并在圖中恰當位置標注兩條曲線的圖例；給圖形加上標題“y1andy2/五、（本題15分）編寫M函數(shù)文件，利用for循環(huán)或while循環(huán)完成計算函數(shù)的任務，并利用該函數(shù)計算時的和（將總程序保存為test05.m文件）。六、（本題13分）已知求解線性規(guī)劃模型：的MATLAB命令為x=linprog（c,A,b,Aeq,beq,VLB,VUB）試編寫MATLAB程序，求解如下線性規(guī)劃問題（將程序保存為test06.m文件）：問題補充：卷子的地址看不見符號，能做就做了一些.1、標點符號（；）可以使命令行不顯示運算結果，（）用來表示該行為注釋行。2、下列變量名中（A）是合法的。（A）cha

9、r_1;（B）x*y;（C）xy;（D）end3、為，步長為的向量，使用命令（本題題意不清）創(chuàng)建。4、輸入矩陣，使用全下標方式用（本題題意不清）取出元素“，”使用單下標方式用（本題題意不清）取出元素“?！?、符號表達式中獨立的符號變量為（）。6、 M腳本文件和M函數(shù)文件的主要區(qū)別是（變量生存期和可見性）和（函數(shù)返回值）。7、在循環(huán)結構中跳出循環(huán)，但繼續(xù)下次循環(huán)的命令為（C）。（A）return;（B）break;（C）continue;（D）keyboad二、(本題12分)利用MATLAB數(shù)值運算，求解線性方程組(將程序保存為test02.m文件)三、(本題20分)利用MATALAB符號運算完

10、成(將程序保存為test03.m文件)：(1)創(chuàng)建符號函數(shù)symsx(2)求該符號函數(shù)對的微分；(3)對趨向于求該符號函數(shù)的極限；(4)求該符號函數(shù)在區(qū)間上對的定積分；(5)求符號方程的解。四、(本題20分)編寫MATALAB程序，完成下列任務(將程序保存為test04.m文件)：(1)在區(qū)間上均勻地取20個點構成向量；(2)分別計算函數(shù)與在向量處的函數(shù)值；(3)在同一圖形窗口繪制曲線與，要求曲線為黑色點畫線，曲線為紅色虛線圓圈；并在圖中恰當位置標注兩條曲線的圖例；給圖形加上標題“y1andy2/五、(本題15分)編寫M函數(shù)文件，利用for循環(huán)或while循環(huán)完成計算函數(shù)的任務，并利用該函數(shù)計

11、算時的和(將總程序保存為test05.m文件)。六、(本題13分)已知求解線性規(guī)劃模型：的MATLAB命令為x=linprog(c,A,b,Aeq,beq,VLB,VUB)試編寫MATLAB程序，求解如下線性規(guī)劃問題(將程序保存為test06.m文件)：例2.1已知SISO系統(tǒng)的狀態(tài)空間表達式為(2-3)式，求系統(tǒng)的傳遞函數(shù)。A=010;001;-4-3-2;B=1;3;-6;C=100;D=0;num,den=ss2tf(a,b,c,d,u)num,den=ss2tf(A,B,C,D,1)例2.2從系統(tǒng)的傳遞函數(shù)(2-4)式求狀態(tài)空間表達式。num=153;den=1234;A,B,C,D=

12、tf2ss(num,den)例2.3對上述結果進行驗證編程。%鈿例2.2上述結果賦值給A、B、C、D陣；A=-2-3-4;100;010;B=1;0;0;C=153;D=0;num,den=ss2tf(A,B,C,D,1)s32s2s3例2.4給定系統(tǒng)G(s)=-2,求系統(tǒng)的零極點增益模型和狀態(tài)空間模型，并求其s30.5s22s1單位脈沖響應及單位階躍響應。解：num=1213;den=10.521;sys=tf(num,den)%系統(tǒng)的傳遞函數(shù)模型Transferfunction:sA3+2sA2+s+3sA3+0.5sA2+2s+1sys1=tf2zp(num,den)%系統(tǒng)的零極點增益模

13、型sys1=sys2=tf2ss(sys)%系統(tǒng)的狀態(tài)空間模型模型；或用a,b,c,d=tf2ss(num,den)形式impulse(sys2)%系統(tǒng)的單位脈沖響應step(sys2)%系統(tǒng)的單位階躍響應例3.1對下面系統(tǒng)進行可控性、可觀性分析。解：a=-1-22;0-11;10-1;b=201,;c=120Qc=ctrb(a,b)%生成能控性判別矩陣rank(Qc)%求矩陣Qc的秩ans=3%蔭秩，故系統(tǒng)能控Qo=obsv(a,c)%生成能觀測性判別矩陣rank(Qo)%求矩陣Qo的秩ans=3%蔭秩，故系統(tǒng)能觀測例3.2已知系統(tǒng)狀態(tài)空間方程描述如下：試判定其穩(wěn)定性，并繪制出時間響應曲線來

14、驗證上述判斷。解：A=-10-35-50-24;1000;0100;0010;B=1;0;0;0;C=172424;D=0;z,p,k=ss2zp(A,B,C,D,1);Flagz=0;n=length(A);fori=1:nifreal(p(i)>0Flagz=1;endenddisp('系統(tǒng)的零極點模型為');z,p,k系統(tǒng)的零極點模型為ifFlagz=1disp('系統(tǒng)不穩(wěn)定,);elsedisp('系統(tǒng)是穩(wěn)定的,);end運行結果為：系統(tǒng)是穩(wěn)定的step(A,B,C,D)%系統(tǒng)的階躍響應資源與環(huán)境工程學院2008級碩士研究生MatLab及其應用試題

15、注意，每題的格式均須包含3個部分a.程序(含程序名及完整程序)：b.運行過程：c.運行結果：(1)求解線性規(guī)劃問題：minZ=-4x1x27x3st.xix2-x3=53x1-x2x3£4x1x2-4x3-7x1,x2-0問各Xi分別取何值時，Z有何極小值。(10分)答：fprintf('線性規(guī)劃問題求解n');f=-4;1;7;A=3,-1,1;1,1,-4;b=4,-7'Aeq=1,1,-1;beq=5'lb=0,0,;ub=;x=linprog(f,A,b,Aeq,beq,lb,ub);xz=f*x;fprintf('MINz=%fn&#

16、39;,z);運行結果：線性規(guī)劃問題求解Optimizationterminatedsuccessfully.x=2.25006.75004.0000MINz=25.7500000.5x,x<2(2)編寫一個函數(shù)，使其能夠產(chǎn)生如下的分段函數(shù)：f(x)=<1.5-0.25x,2cxM6,0.5,6<x并調用此函數(shù)，繪制在乂=0,+2范圍的£他)“乂+2)曲線。(10分)答：functiony=f(x)ifx<=2y=0.5*x;elseifx>6y=0.5;elsey=1.5-0.25*x;endend運行結果x=2f(x)=1x=0:0.05:2;y=d

17、iag(A2(x)'*A2(x+2);plot(x,y);xlabel('bfx');ylabel('bfy');(3)將一個屏幕分4幅，選擇合適的步長在右上幅與左下幅繪制出下列函數(shù)的圖形(10分)_2Jcos(x),x=(曲線圖)；f(x,y)=222224；(-2<x<2,-4<y<44)(曲面m)o答:>>subplot(2,2,2)>>ezplot('(cos(x)A(1/2)',-pi/2pi/2)>>ylabel('y')>>subplot

18、(2,2,3)>>x=-2:0.5:2;>>y=-4:1:4;>>ezsurfc('xA2/2A2+yA2/4A2')(4) A是一個雉度亦n的矩陣.寫一段程序,算出A中有多少個零元素(10分)答：>>A=input(E青輸入一個矩陣)m,n=size(A);sig=0;fori=1:mforj=1:nifA(i,j)=0sig=sig+1;endendend請輸入一個矩陣012;102;000A=0121 02000>>sigsig=5(5)向量A=a1，a1，*.寫一段程序，找出A中的最小元素(10分)答：A=in

19、put('請輸入一個向量')m,n=sizeAmin=A(1,n);fori=1:nifA(1,i)<minmin=A(1,i)endend請輸入一個向量123-520A=123-520min=-5B.應用題(50分根據(jù)專業(yè)方向特色和相關科研工作需求，經(jīng)過與導師商量后，結合一個課題具體任務，編寫一份Matlab應用工作報告。報告由：a課題任務要求，b技術路線，c程序，d運行結果，e總結、等部分構成，完成的報告經(jīng)導師給出簡單評語并簽字后繳來。a,課題任務：研究了一種生物質，油菜秸稈對水溶液中金屬離子銅的吸附行為，分別從pH,用量，溫度幾個方面考察秸稈的吸附性，并對分析的最佳

20、條件進行了探討。同時從吸附熱力學和吸附動力學角度探討吸附機理。結果表明，100ml溶液pH=5.30,秸稈用量0.75g時，秸稈對銅的吸附量可達到6mg/g左右。b,技術路線：通過實驗，獲得一系列的數(shù)據(jù)，然后通過Matlab來做各種關系圖。從圖中找到g各種關系式。c,程序：x=2.202.723.444.135.38y=2.393.836.076.396.84Plot(x,y);xlabel('pH');ylabel('吸附量')x=0.50.751.01.251.5y=6.056.195.334.694.02plot（x,y）;xlabel（'秸稈用量

21、g'）;ylabel（'吸附量'）Gj5力百0£070.6口911.11.21.31.41.5秸稈用圖2通過數(shù)據(jù)圖，得到比較理想的實驗條件pH和秸稈用量，接下來做動力學和等溫線。>x=0,1670.5123458y=0.0332010.0860590.1697790.3220610.4807690.6441220.8090611.269841plot（x,y）;xlabel（'時間t'）;ylabel（'時間/吸附量'）x=0.23630.154960.136190.129060.133730.13315y=0.25218

22、0.047070.020140.012670.008810.00706plot（x,y）;xlabel（'1/吸附量'）;ylabel（'1/平衡濃度'）D.14Q.10.1S0.21碟附O.BQ,1I,1D.8.120.220.x=0.626540.809770.865850.88920.873770.87564y=0.598291.32731.695891.897372.055032.15149plot(x,y);xlabel('Lg吸附量');ylabel('Lg平衡濃度')00D.B口.昭a7O75QB口酣Q.3O.WLf

23、l破陽422.工圖5d,總結：從圖1和圖2,分析看可以得到比較理想的對于本次實驗的pH和秸稈用量。后面實驗是在前面的基礎上得到的。圖3是吸附動力學反應速率圖，從圖中可以看到線性擬合程度很好，符合二級反應速率方程。圖4和圖5是吸附等溫線作圖，看以看出圖4的線性擬合較圖5的好，說明符合Langmuir吸附等溫模型。例2.1已知SISO系統(tǒng)的狀態(tài)空間表達式為（2-3）式，求系統(tǒng)的傳遞函數(shù)。A=010;001;-4-3-2;B=1;3;-6;C=100;D=0;num,den=ss2tf(a,b,c,d,u)num,den=ss2tf(A,B,C,D,1)例2.2從系統(tǒng)的傳遞函數(shù)(2-4)式求狀態(tài)空間

24、表達式。num=153;den=1234;A,B,C,D=tf2ss(num,den)例2.3對上述結果進行驗證編程。%鈿例2.2上述結果賦值給AB、C、D陣；A=-2-3-4;100;010;B=1;0;0;C=153;D=0;num,den=ss2tf(A,B,C,D,1)s32s2s3例2.4給定系統(tǒng)G(s)2,求系統(tǒng)的零極點增益模型和狀態(tài)空間模s3+0.5s2+2s+1型，并求其單位脈沖響應及單位階躍響應。解：num=1213;den=10.521;sys=tf(num,den)Transferfunction:sA3+2sA2+s+3%系統(tǒng)的傳遞函數(shù)模型sA3+0.5sA2+2s+1

25、sys1=tf2zp(num,den)%系統(tǒng)的零極點增益模型sysl=sys2=tf2ss(sys)%系統(tǒng)的狀態(tài)空間模型模型；或用a,b,c,d=tf2ss(num,den)形式impulse(sys2)%系統(tǒng)的單位脈沖響應step(sys2)%系統(tǒng)的單位階躍響應例3.1對下面系統(tǒng)進行可控性、可觀性分析。解：a=-1-22;0-11;10-1;b=201,;c=120Qc=ctrb(a,b)%fe成能控性判別矩陣rank(Qc)麻矩陣Qc的秩ans=3%蔭秩，故系統(tǒng)能控Qo=obsv(a,c)%fe成能觀測性判別矩陣rank(Qo)麻矩陣Qo的秩ans=3%蔭秩，故系統(tǒng)能觀測例3.2已知系統(tǒng)狀

26、態(tài)空間方程描述如下：試判定其穩(wěn)定性，并繪制出時間響應曲線來驗證上述判斷。解：A=-10-35-50-24;1000;0100;0010;B=1;0;0;0;C=172424;D=0;z,p,k=ss2zp(A,B,C,D,1);Flagz=0;n=length(A);fori=1:nifreal(p(i)>0Flagz=1;endenddisp('系統(tǒng)的零極點模型為');z,p,k系統(tǒng)的零極點模型為ifFlagz=1disp('系統(tǒng)不穩(wěn)定,);elsedisp（'系統(tǒng)是穩(wěn)定的,）;end運行結果為：系統(tǒng)是穩(wěn)定的step（A,B,C,D）%系統(tǒng)的階躍響應。1

27、、使用下列哪一個函數(shù)可以產(chǎn)生單位矩陣（B）A.zerosB.eyeC.randD.diag2、下列哪一個函數(shù)是求模函數(shù)（D）A.remB.signC.fixD.mod3、使用下列哪一個函數(shù)可以交換矩陣左右對稱位置上的元素（A）A.fliplrB.flipdimC.flipudD.find4、使用下列哪一個函數(shù)可以比較字符串，且比較時忽略字符的大小寫（D）A.strncmpB.strcmpC.strncmpiD.strcmpi5、要利用圖形方式顯示元胞數(shù)組，則應該使用下列哪一個函數(shù)（B）A.cellfunB.cellplotC.celldispD.cell2mat6、下列哪一個函數(shù)可以獲取結構字

28、段的數(shù)據(jù)（B）A.fieldnamesB.getfieldC.setfieldD.rmfield7、執(zhí)行下列哪一條命令后，圖形窗體的軸將顯示坐標網(wǎng)格線（A）A.gridonB.holdonC.gridoffD.holdoff8、進行格式化繪圖時，使用哪一個函數(shù)可以添加圖例（B）A.titleB.legendC.labelD.text9、使用下列哪一條指令可以將圖形窗體分割成二行三列，并且將第一行第二列的繪圖區(qū)域設置為當前的繪圖區(qū)域（B）A.subplot（2,3,1）B.subplot（2,3,2）C.subplot（2,3,4）D.subplot（2,3,5）10、使用下列哪一個函數(shù)可以繪制

29、三維網(wǎng)線圖（C）A.surfB.plotC.meshD.plot31、>>A=123;456;B=25;83（2分）B=25832、假設向量A=9876543210>>A（1:3:5）（2分）ans=9 6>>A（135）（2分）ans=9753、>>A=ones（2,2）;A（:）=1:4;>>A*A'（2分）ans=10 141420>>B=A.*A'(2分)B=166164、使用三元組法，將下列滿陣轉變?yōu)橄∈杈仃?500220-150113000S=000-60091000000002800解：>

30、; >ir=14221351;> >jc=11234446;> >data=159111322-628-15;> >s=sparse(ir,jc,data,5,6)s=(1,1)15(4,1)91(2,(2) 11(2,(3) 3(2,(4) 22(3,4)-6(5,4)28(1,6)-155、>>A=reshape(1:24,4,6);>>A(:,234)=A=1 17212 18223 19234 20246、使用函數(shù)struct創(chuàng)建一個結構。此結構名為Student;有三個字段，分另1J為nameage、grade;有兩條

31、記錄，分別為'Way、23、3和'Deni'、21、1解：>>Student=struct('name','Way','Deni','age',23,21,'grade',3,1)7、繪出下幅圖>>plot(14285)8765432111,522.533.544.55三、寫出使以下這段文字成為字符串的MATLAB旨令。注意保持這段文字的格式。(10分)Areyouateacher?No,I'mastudent/解：>>a='"Areyouateacher?'''>>b='''No,I''ma_student.'''>>c=strvcat(a