Matlab_考题题整理_带答案

1、MATLAB考试试题产生一个1x10的随机矩阵，大小位于(-5 5 )，并且按照从大到小的顺序排列好！(注：要程序和运行结果的截屏)答案：a=10*ra nd(1,10)-5;b=sort(a,'desce nd')1 .请产生一个100*5的矩阵，矩阵的每一行都是1 2 3 4 52. 已知变量： A=' ilovematlab ' B=' matlab '请找岀：(A ) B在A中的位置。(B) 把 B 放在 A 后面，形成 C= ' ilovematlabmatlab '3. 请修改下面的程序，让他们没有for循环语句！A=

2、1 2 3; 4 5 6; 7 8 9;r c=size(A);for i=1:1:rfor j=1:1:cif (A(i,j)>8 | A(i,j)<2)A(i,j)=0;endendend4. 请把变量 A=1 2 3; 4 5 6; 7 8 9写到文件里(output.xls)，写完后文件看起来是这样的1 2 3 4 5 6 7 8 95. 试从Yahoo网站上获得微软公司股票的2008年9月的每日收盘价。6. 编写M文件，从Yahoo网站批量读取 60000.SH 至600005.SH 在2008年9月份的每日收盘价(提 示：使用字符串函数)。7. 将金牛股份(000937

3、 ) 2005年12月14日至2006 年1月10日的交易记录保存到Excel中，编写程序将数据读入MATLAB中，进一步将数据读入Access数据库文件。8. 已知资产每日回报率为0.0025 ，标准差为0.0208 ，资产现在价值为0.8亿，求5%水平下资产的10天在险价值(Var)。试用MATLAB 中9. a=1 2 3 4 5,b=a(1)*a(5)+a (2) *a (4)+a (3)*a (3)+a *a (2)+a(5)*a(1). 最简单的方法计算b，注意最简单哦。1、求下列联立方程的解 3x+4y-7z-12w=45x-7y+4z+ 2w=-3x+8z- 5w=9-6x+5

4、y-2z+10w=-8 求系数矩阵的秩； 求出方程组的解。解：（ 1）>> a=34-7 -12;5 -742 ;1 08-5;-6 5-210;c=4; -3; 9;-8;b=ra nk(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+xA2)把 x=02n 间分为 101 点，画出以 x 为横坐标，y为纵坐标的曲线；解:>> x=linspace(0,2*pi,101

5、);y=cos(0.5+3.*si n(x)./(1+x.*x);plot(x,y)3、设f(x)=xA5-4xA4+3xA2-2x+6(1)取x=-2,8之间函数的值(取100个点)，画出曲线，看它有几个零点(提示：用polyval函数) 解:>>p=1 -4 3 -2 6; x=li nspace(-2,8,100); y=polyval(p,x);Plot(x,y);axis(-2,8, -200,2300);为了便于观察，在y=0处画直线，图如下所示:与y=0直线交点有两个，有两个实根。(2)用roots函数求此多项式的根>> a=roots(p)a =3.00

6、00,1.6956, -0.3478 + 1.0289i, -0.3478 - 1.0289i4、在-10，10； -10，10范围内画出函数的三维图形。解:>>X,Y=meshgrid(-10 : 0.5 :10); a=sqrt(X.A2+Y.A2)+eps;Z=si n( a)./a;mesh(X,Y,Z);IO -1 matlab试卷，求答案、选择或填空(每空2分，共20分)1、标点符号()可以使命令行不显示运算结果，()用来表示该行为注释行。2、下列变量名中()是合法的。(A)char_1 ;(B)x*y ;(C) xy ;(D) end3、为，步长为的向量，使用命令()

7、创建。出元素5、符号表达式中独立的符号变量为（）6、 M脚本文件和M函数文件的主要区别是（）和（）。7、在循环结构中跳出循环，但继续下次循环的命令为（）（A） return; （B） break ;（C） continue ;（D） keyboad二、（本题12分）利用MATLAB数值运算，求解线性方程组（将程序保存为test02.m文件）三、（本题20分）利用MATALAB符号运算完成（将程序保存为test03.m文件）：（1）创建符号函数（2 ）求该符号函数对的微分；（3）对趋向于求该符号函数的极限；（4） 求该符号函数在区间上对 的定积分；（5）求符号方程的解。四、（本题20分）编写MA

8、TALAB程序，完成下列任务（将程序保存为test04.m文件）：（1）在区间 上均匀地取20个点构成向量；（2）分别计算函数与在向量处的函数值；（3） 在同一图形窗口绘制曲线与，要求 曲线为黑色点画线，曲线为红色虚线圆圈；并在图中恰当位置标注两条曲线的图例；给图形加上标题“y1 and y2。”五、（本题15分）编写M函数文件，利用for循环或while循环完成计算函数的任务，并利用该函数计算 时的和（将总程序保存为 test05.m文件）。六、（本题13分）已知求解线性规划模型： 的MATLAB命令为x=li nprog （ c,A,b,Aeq,beq,VLB,VUB ）试编写MATLAB

9、程序，求解如下线性规划问题（将程序保存为test06.m文件）问题补充：卷子的地址看不见符号，能做就做了一些.1、标点符号 （；）可以使命令行不显示运算结果,（%） 用来表示该行为注释行。(A)char_1;(B)x*y ;(C) xy ;(D) end3、为，步长为的向量,使用命令（本题题意不清）创建。2、下列变量名中（A ）是合法的。）取出元素“，使用单下标方式用4、 输入矩阵，使用全下标方式用（本题题意不清题题意不清）取出兀素“。"5、 符号表达式中独立的符号变量为（）变量生存期和可见性6、 M脚本文件和M函数文件的主要区别是（和（函数返回值）。7、在循环结构中跳出循环，但继续

10、下次循环的命令为（C ）(D) keyboad(A) return; (B) break ;(C) continue ;（本题12分）利用MATLAB数值运算，求解线性方程组（将程序保存为test02.m文件）三、（本题20分）利用MATALAB符号运算完成（将程序保存为test03.m文件）：（1）创建符号函数syms x（2 ）求该符号函数对的微分；（3）对趋向于求该符号函数的极限；（4） 求该符号函数在区间上对 的定积分；（5）求符号方程的解。四、（本题20分）编写MATALAB程序，完成下列任务（将程序保存为testO4.m文件）：（1） 在区间 上均匀地取20个点构成向量；（2）分别

11、计算函数与在向量处的函数值；（3） 在同一图形窗口绘制曲线 与，要求 曲线为黑色点画线，曲线为红色虚线圆圈；并在图中恰当位置标注两条曲线的图例；给图形加上标题“y1 and y2。”五、（本题15分）编写M函数文件，利用for循环或while循环完成计算函数的任务，并利用该函数计算 时的和（将总程序保存为 test05.m文件）。六、（本题13分）已知求解线性规划模型：的MATLAB命令为x=li nprog （ c,A,b,Aeq,beq,VLB,VUB ）试编写MATLAB程序，求解如下线性规划问题（将程序保存为test06.m文件）：例2.1已知SISO系统的状态空间表达式为(2-3)式

12、，求系统的传递函数。A=0 1 0;0 0 1;-4 -3 -2;B=1;3;-6;C=1 0 0;D=0;nu m,de n=ss2tf(a,b,c,d,u)n um,de n=ss2tf(A,B,C,D,1)例2.2从系统的传递函数(2-4)式求状态空间表达式。num =1 5 3;den =1 2 3 4;A,B,C,D=tf2ss( num,de n)例2.3对上述结果进行验证编程。%各例2.2上述结果赋值给A B、C、D 阵;A =-2 -3 -4;1 0 0; 0 1 0；B =1;0;0； C =1 5 3； D=0;num,den=ss2tf(A， B，C,D,1)例2.4给定

13、系统G（s）s3 2s2 s 3s30.5s2 2s !，求系统的零极点增益模型和状态空间模型，并求其单位脉冲响应及单位阶跃响应。解：num=1 2 1 3;den=1 0.5 2 1; sys=tf(num,den)%系统的传递函数模型Transfer function:sA3 + 2 sA2 + s + 3 sA3 + 0.5 sA2 + 2 s + 1sys1=tf2zp(num,den)%系统的零极点增益模型sys1 =sys2=tf2ss(sys)%系统的状态空间模型模型；或用a,b,c,d=tf2ss(num,den) 形式impulse(sys2)%系统的单位脉冲响应step(s

14、ys2)%系统的单位阶跃响应 例 3.1 对下面系统进行可控性、可观性分析。 解：a=-1 -2 2;0 -1 1;1 0 -1;b=2 0 1'c=1 2 0Qc=ctrb(a,b)%生成能控性判别矩阵rank(Qc)%求矩阵Qc的秩ans = 3%满秩，故系统能控Qo=obsv(a,c)%生成能观测性判别矩阵rank(Qo)%求矩阵Qo的秩ans = 3%满秩，故系统能观测 例 3.2 已知系统状态空间方程描述如下： 试判定其稳定性，并绘制出时间响应曲线来验证上述判断。 解：A=-10 -35 -50 -24;1 0 0 0;0 1 0 0;0 0 1 0;B=1;0;0;0;C=

15、1 7 24 24;D=0;z,p,k=ss2zp(A,B,C,D,1);Flagz=0; n=length(A);for i=1:nif real(p(i)>0Flagz=1;endenddisp(' 系统的零极点模型为 ');z,p,k 系统的零极点模型为if Flagz=1disp(' 系统不稳定 ');else disp（'系统是稳定的'）;end运行结果为：系统是稳定的step(A,B,C,D)%系统的阶跃响应资源与环境工程学院2008级硕士研究生MatLab及其应用试题注意，每题的格式均须包含3个部分a. 程序(含程序名及完整程

16、序)：b. 运行过程：c. 运行结果：求解线性规划问题：minZ 二-4x1 x2 7x3st. X1 X2 - X3 = 53X1 -X2X3 冬 4X1X2 -4X3 - -7X1,X2 _ 0问各Xi分别取何值时，Z有何极小值。（10 分） 答:fprintf（'线性规划问题求解n'）;f =卜4；1；7;A = 3,-1,1;1,1,-4;b = 4,-7'Aeq = 1,1,-1;beq ='；lb = 0,0,;ub =;x = lin prog（f,A,b,Aeq,beq,lb,ub）;xz = f * x;fprintf（'MIN z =

17、 %f n' , z）;运行结果：线性规划问题求解Optimization terminated successfully. x =2.25006.75004.0000MIN z = 25.750000编写一个函数，使其能够产生如下的分段函数:0.5x,x<2f(x)=<1.50.25x, 2<x 兰 6 ,0.5,6 兰 x并调用此函数，绘制 在x二02范围的f(xf(x 2)曲线。(10分)答：function y=f(x)if x<=2y=0.5*x;else if x>6y=0.5;else y =1.5-0.25*x;endend运行结果x=2f

18、(x)=1x = 0:0.05:2;y = diag(A2(x)'*A2(x+2);plot(x,y);xlabel('bfx');ylabel('bfy');将一个屏幕分4 幅,选择合适的步长在右上幅与左下幅绘制出下列函数的图形 分）2 2f（x,y）歩朱）（曲. COS(X),-，(曲线图)；面图)。答：>> subplot(2,2,2)>> ezplot('(cos(x)F(1/2)', -pi/2 pi/2)>> ylabel('y')>>subplot(2,2,3)

19、>> x=-2:0.5:2;>> y=-4:1:4;>>ezsurfc('xA2/2A2+yA2/4A2')(4) A是一个維度mKn的矩阵.写一段程序，算出A中有多少个零元素(10分) 答：>> A= in put (请输入一个矩阵')m, n= size(A);sig=0;for i=1:mfor j=1: nif A(i,j)=0sig = sig+1;endendend请输入一个矩阵A =0 1 2;1 0 2; 0 0 0012102000>> sigsig =向量A "“印，an.写一段程

20、序，找出A中的最小元素(10 分)答：A= in put ('请输入一个向量')m, n=sizeAmin =A(1, n);for i=1: nif A(1,i)<minmi n=A(1,i)endend请输入一个向量1 2 3 -5 2 0A =123-520min =B应用题（50分）根据专业方向特色和相关科研工作需求，经过与导师商量后，结合一个课题具 体任务，编写一份 Matlab应用工作报告。报告由：a课题任务要求，b技术路 线，c程序，d运行结果，e总结、等部分构成，完成的报告经导师给出简单评 语并签字后缴来。a, 课题任务：研究了一种生物质，油菜秸秆对水溶液

21、中金属离子铜的吸附行 为，分别从pH,用量，温度几个方面考察秸秆的吸附性，并对分析的最佳条件 进行了探讨。同时从吸附热力学和吸附动力学角度探讨吸附机理。结果表明， 100ml溶液pH=5.30,秸秆用量0.75g时，秸秆对铜的吸附量可达到6mg/g左右。b, 技术路线：通过实验，获得一系列的数据，然后通过Matlab来做各种关系图。从图中找到g各种关系式。c, 程序:x=2.20 2.72 3.44 4.13 5.38y=2.39 3.83 6.07 6.39 6.84Plot(x,y);xlabel('pH');ylabel('吸附量')x=0.5 0.75

22、1.0 1.25 1.5 y=6.05 6.19 5.33 4.69 4.02 plot（x,y）;xlabel（'秸秆用量g'）;ylabel（'吸附量'）Miff.B03 I口 E Q?5 6fi-通过数据图，得到比较理想的实验条件pH和秸秆用量，接下来做动力学和等温线。> x=0.167 0.5 1 2 3 4 5 8 y=0.033201 0.086059 0.169779 0.322061 0.480769 0.644122 0.809061 1.269841Plot（x,y）;xlabel（'时间 t'）;ylabel（时间/吸

23、附量'）IF I.8.&E D.mD.2DD.x=0.2363 0.15496 0.13619 0.12906 0.13373 0.13315 y=0.25218 0.04707 0.02014 0.01267 0.00881 0.00706 plot（x,y）;xlabel（'1/吸附量'）;ylabel（'1/平衡浓度'）匚14口 TE IB0.2>Z213.245 3 5 2 LH 1 -5 "匚 2匚0.1匸M0舞證i口 12图4x=0.62654 0.80977 0.86585 0.8892 0.87377 0.8756

24、4 y=0.59829 1.3273 1.69589 1.89737 2.05503 2.15149 plot（x,y）;xlabel（'Lg 吸附量'）;ylabel（'Lg平衡浓度'）t±请告0.0Q.h 6507口巧030 95Q.PU 5EL楝附8图5d, 总结：从图1和图2,分析看可以得到比较理想的对于本次实验的pH和秸秆用量。后面实验是在前面的基础上得到的。图3是吸附动力学反应速率图，从图中可以看到线性拟合程度很好，符合二级反应速率方程。图4和图5是吸附 等温线作图，看以看出图4的线性拟合较图5的好，说明符合Langmuir吸附等 温模型。

25、例2.1已知SISO系统的状态空间表达式为(2-3)式，求系统的传递函数。A=0 1 0;0 0 1;-4 -3 -2;B=1;3;-6;C=1 0 0;D=0;nu m,de n=ss2tf(a,b,c,d,u) n um,de n=ss2tf(A,B,C,D,1)例2.2从系统的传递函数(2-4)式求状态空间表达式。num =1 5 3;den =1 2 3 4;A,B,C,D=tf2ss( nu m,de n) 例2.3对上述结果进行验证编程。%各例2.2上述结果赋值给 A B、C、D阵;A =-2 -3 -4;1 0 0; 0 1 0；B =1;0;0； C =1 5 3； D=0;num,den=ss2tf(A， B，C，D,1)例2.4给定系统G(s)，求系统的零