Matlab考題題整理帶答案

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

2、or 循環(huán)語句！A=123;456;789;rc=size(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 網站上獲得微軟公司股票的 2008 年 9 月的每日收盤價。6. 編寫 M文件，從 Yahoo網站批量讀取 60000.SH 至 600005.SH 在 2008 年 9 月份的每日收盤價(提 示：使用字符串函數)。7. 將金牛股份( 000937) 2

3、005 年 12 月 14 日至 2006 年 1 月 10 日的交易記錄保存到 Excel 中， 編寫程序將數據讀入 MATLAB中，進一步將數據讀入 Access 數據庫文件。8. 已知資產每日回報率為 0.0025 ，標準差為 0.0208 ，資產現在價值為 0.8 億，求 5 水平下資產的 10 天在險價值 (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).試用 MATLAB中最簡單的方法計算 b ，注意最簡單哦。1、求下列聯立方程的解3x+4y-7z-12w=45x-7y+4z+2w=-3x+8

4、z-5w=9-6x+5y-2z+10w=-8 求系數矩陣的秩； 求出方程組的解。解：( 1)>>a=34-7-12;5-742;108-5;-65-210;c=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.24292、設y=cos0.5+(3sinx)/(1+x2)把 x=02 間分為 101y 為縱坐標的曲線；解： >>x=linspace(0,2*pi,101);y=cos(0.5+3.*sin(x

5、)./(1+x.*x); plot(x,y)3、設 f(x)=x5-4x4+3x2-2x+6( 1)取 x=-2,8 之間函數的值(取 100 個點)，畫出曲線， (提示：用 polyval 函數)解： >>p=1-43-26;x=linspace(-2,8,100);y=polyval(p,x);plot(x,y);點，畫出以 x 為橫坐標，看它有幾個零點axis(-2,8,-200,2300);為了便于觀察，在 y=0 處畫直線，圖如下所示：與 y=0 直線交點有兩個，有兩個實根。（ 2）用 roots 函數求此多項式的根>>a=roots（p）a=3.0000,1

6、.6956,-0.3478+1.0289i,-0.3478-1.0289i4、在 -10 ， 10； -10 ， 10 范圍內畫出函數的三維圖形解： >>X,Y=meshgrid(-10:0.5:10);a=sqrt(X.2+Y.2)+eps;Z=sin(a)./a;mesh(X,Y,Z);matlab 試卷，求答案一、選擇或填空（每空 2 分，共 20 分） 1、標點符號（）可以使命令行不顯示運算結果，（）用來表示該行為注釋行。2 、下列變量名中（）是合法的。 （A）char_1;（B）x*y;（C）xy;（D）end3 、為，步長為的向量，使用命令（）創(chuàng)建。4、輸入矩陣，使用全

7、下標方式用（）取出元素“，”使用單下標方式用（）取出元素 5 、符號表達式中獨立的符號變量為（）。6、 M腳本文件和 M函數文件的主要區(qū)別是（）和（）。7 、在循環(huán)結構中跳出循環(huán)，但繼續(xù)下次循環(huán)的命令為（）。 （A）return;（B）break;（C）continue;（D）keyboad二、（本題 12 分）利用 MATLAB數值運算，求解線性方程組 （ 將程序保存為 test02.m 文件 ）三、（本題 20 分）利用 MATALAB符號運算完成（將程序保存為 test03.m 文件）：（ 1 ）創(chuàng)建符號函數（ 2 ）求該符號函數對的微分；（ 3 ）對趨向于求該符號函數的極限；（ 4 ）

8、求該符號函數在區(qū)間上對的定積分；（ 5 ）求符號方程的解。四、（本題 20 分）編寫 MATALAB程序，完成下列任務（將程序保存為test04.m 文件）：（ 1 ）在區(qū)間上均勻地取 20 個點構成向量；（ 2 ）分別計算函數與在向量處的函數值；（ 3 ）在同一圖形窗口繪制曲線與，要求曲線為黑色點畫線，曲線為紅色虛線圓圈；并在圖中恰當位 置標注兩條曲線的圖例；給圖形加上標題“ y1andy2 ?！蔽?、（本題 15 分）編寫 M函數文件，利用 for 循環(huán)或 while 循環(huán)完成計算函數的任務，并利用該函 數計算時的和（將總程序保存為 test05.m 文件）。六、（本題 13 分）已知求解線

9、性規(guī)劃模型：的 MATLAB命令為x=linprog （ c,A,b,Aeq,beq,VLB,VUB ）試編寫 MATLAB程序，求解如下線性規(guī)劃問題（將程序保存為test06.m 文件）：問題補充：卷子的地址看不見符號 , 能做就做了一些 .1、標點符號（ ; ）可以使命令行不顯示運算結果，（ %）用來表示該行為注釋行。2 、下列變量名中（ A）是合法的。（A）char_1;（B）x*y;（C）xy;（D）end3 、為，步長為的向量，使用命令（本題題意不清）創(chuàng)建。4、輸入矩陣，使用全下標方式用（本題題意不清）取出元素“，”使用單下標方式用（本題題意不清） 取出元素“ ?！? 、符號表達式中

10、獨立的符號變量為（）。6、 M腳本文件和 M函數文件的主要區(qū)別是（變量生存期和可見性）和（函數返回值）。7 、在循環(huán)結構中跳出循環(huán)，但繼續(xù)下次循環(huán)的命令為（C）。（A）return;（B）break;（C）continue;（D）keyboad二、（本題 12 分）利用 MATLAB數值運算，求解線性方程組 （ 將程序保存為 test02.m 文件 ）三、（本題 20 分）利用 MATALAB符號運算完成（將程序保存為 test03.m 文件）： （ 1 ）創(chuàng)建符號函數 symsx（ 2 ）求該符號函數對的微分；（ 3 ）對趨向于求該符號函數的極限；（ 4 ）求該符號函數在區(qū)間上對的定積分；（

11、 5 ）求符號方程的解。四、（本題 20 分）編寫 MATALAB程序，完成下列任務（將程序保存為test04.m 文件）：（ 1 ）在區(qū)間上均勻地取 20 個點構成向量；（ 2 ）分別計算函數與在向量處的函數值；（ 3 ）在同一圖形窗口繪制曲線與，要求曲線為黑色點畫線，曲線為紅色虛線圓圈；并在圖中恰當位 置標注兩條曲線的圖例；給圖形加上標題“ y1andy2 ?！蔽?、（本題 15 分）編寫 M函數文件，利用 for 循環(huán)或 while 循環(huán)完成計算函數的任務，并利用該函 數計算時的和（將總程序保存為 test05.m 文件）。六、（本題 13 分）已知求解線性規(guī)劃模型：的 MATLAB命令為

12、x=linprog （ c,A,b,Aeq,beq,VLB,VUB ）試編寫 MATLAB程序，求解如下線性規(guī)劃問題（將程序保存為test06.m 文件）：例 2.1 已知 SISO系統(tǒng)的狀態(tài)空間表達式為 （2-3） 式，求系統(tǒng)的傳遞函數。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)的傳遞函數 （2-4） 式求狀態(tài)空間表達式。num=153;den=1234;A,B,C,D=tf2ss(num,den) 例 2.3 對上述結果進行驗證編程。%將例 2

13、.2 上述結果賦值給 A、B、C、D陣；A=-2-3-4;100;010； B=1;0;0 ； C=153 ； D=0；專業(yè)資料整理num,den=ss2tf(A， B，C， D,1)例 2.4 給定系統(tǒng)32s2ss3G(s)，求系32統(tǒng) 的零極點增益模型和狀態(tài)空間模型，并求其s0.5s2s1單位脈沖響應及單位階躍響應。解：num=1213;den=10.521;sys=tf(num,den)% 系統(tǒng)的傳遞函數模型Transferfunction:s3+2s2+s+3s3+0.5s2+2s+1sys1=tf2zp(num,den)% 系統(tǒng)的零極點增益模型 sys1=形式sys2=tf2ss(s

14、ys)% 系統(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)定性，并繪制出時間響

15、應曲線來驗證上述判斷。 解：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)境工程學院

16、2008 級碩士研究生 MatLab 及其應用試題 注意，每題的格式均須包含 3 個部分a. 程序(含程序名及完整程序)：b. 運行過程：c. 運行結果：(1) 求解線性規(guī)劃問題：minZ 4x x 7x1 2 3s.t. x x x 51 2 33x x x 41 2 3x x 4x 71 2 3x, x 01 2問各 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

17、(f,A,b,Aeq,beq,lb,ub);xz=f'*x;fprintf('MINz=%fn',z); 運行結果：線性規(guī)劃問題求解Optimizationterminatedsuccessfully.x=2.25006.75004.0MINz=25.7500000.5x ，x2(2) 編寫一個函數，使其能夠產生如下的分段函數： f(x)1.50.25x2x60.56xx 2)曲線。(10 分)并調用此函數，繪制在 x0 ， 2 范圍的 f(x)f(WORD格式答 ： functiony=f(x) ifx<=2y=0.5*x;elseifx>6y=0.5;e

18、lsey=1.5-0.25*x;endend 運行結果 x=2 f(x)=1 x=0:0.05:2; y=diag(A2(x)'*A2(x+2); plot(x,y); xlabel('bfx'); ylabel('bfy');專業(yè)資料整理10(3) 將一個屏幕分 4 幅，選擇合適的步長在右上幅與左下幅繪制出下列函數的圖形 分)xy cos(x ， x，(曲線圖)； (-2x2,-4y4) )f(x,y)(曲面圖)。22 ；22答： >>subplot(2,2,2) >>ezplot('(cos(x)(1/2)',

19、-pi/2pi/2)>>ylabel('y')>>subplot(2,2,3)>>x=-2:0.5:2;>>y=-4:1:4;>>ezsurfc('x2/22+y2/42')10 分）（4）A 是一個維度 m×n 的矩陣 . 寫一段程序 , 算出 A中有多少個零元素（答： >>A=input（' 請輸入一個矩陣 '） m,n=size(A);sig=0;fori=1:mforj=1:nifA(i,j)=0sig=sig+1;endendend請輸入一個矩陣 012;

20、102;000A=012 102000>>sigsig=5(5) 向量 Aa, a, a 11n. 寫一段程序 , 找出 A 中的最小元素( 10 分)答： A=input（' 請輸入一個向量 '）m,n=sizeAmin=A(1,n); fori=1:n ifA(1,i)<minmin=A(1,i)endend請輸入一個向量 123-520A=123-520min=-5B.應用題( 50 分)根據專業(yè)方向特色和相關科研工作需求，經過與導師商量后，結合一個課題具 體任務，編寫一份 Matlab 應用工作報告。報告由： a 課題任務要求， b 技術路 線， c

21、程序， d 運行結果， e 總結、等部分構成，完成的報告經導師給出簡單評 語并簽字后繳來。a, 課題任務：研究了一種生物質，油菜秸稈對水溶液中金屬離子銅的吸附行為，分別從 pH，用量，溫度幾個方面考察秸稈的吸附性，并對分析的最佳條件進行了探 討。同時從吸附熱力學和吸附動力學角度探討吸附機理。結果表明，100ml 溶液 pH=5.30 ，秸稈用量 0.75g 時，秸稈對銅的吸附量可達到 6mg/g 左右。b, 技術路線：通過實驗，獲得一系列的數據，然后通過 Matlab 來做各種關 系圖。從圖中找到 g 各種關系式。c, 程序： x=2.202.723.444.135.38y=2.393.836

22、.076.396.84plot(x,y);xlabel('pH');ylabel(' 吸附量 ')圖1x=0.50.751.01.251.5y=6.056.195.334.694.02 plot（x,y）;xlabel（' 秸稈用量 g'）; ylabel（' 吸附量 '）圖2通過數據圖，得到比較理想的實驗條件pH 和秸稈用量，接下來做動力學和等溫線。>x=0.1670.5123458y=0.0332010.0860590.1697790.3220610.4807690.6441220.80906110. plot(x,y)

23、;xlabel(' 時間 t');ylabel(' 時間 / 吸附量 ')圖3x=0.23630.154960.136190.129060.133730.13315y=0.252180.047070.020140.012670.008810.00706 plot(x,y);xlabel('1/ 吸附量 ');ylabel('1/ 平衡濃度 ')圖4x=0.626540.809770.865850.88920.873770.87564y=0.598291.32731.695891.897372.055032.15149 plot(x

24、,y);吸附量 '); 平衡濃度 ')xlabel('Lgylabel('Lg圖5d, 總結：從圖 1 和圖 2，分析看可以得到比較理想的對于本次實驗的pH和秸稈用量。后面實驗是在前面的基礎上得到的。圖3 是吸附動力學反應速率圖，從圖中可以看到線性擬合程度很好，符合二級反應速率方程。圖4 和圖 5 是吸附等溫線作圖，看以看出圖 4 的線性擬合較圖 5 的好，說明符合 Langmuir 吸附等 溫模型。 例 2.1 已知 SISO 系統(tǒng)的狀態(tài)空間表達式為 (2-3) 式，求系統(tǒng)的傳遞函數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)的傳遞函數 (2-4) 式求狀態(tài)空間表達式。num=153;den=1234;A,B,C,D=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