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 循环语句！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 、为，步长为的向量，使用命令（）创建。4、输入矩阵，使用全

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 ）

8、求该符号函数在区间上对的定积分；（ 5 ）求符号方程的解。四、（本题 20 分）编写 MATALAB程序，完成下列任务（将程序保存为test04.m 文件）：（ 1 ）在区间上均匀地取 20 个点构成向量；（ 2 ）分别计算函数与在向量处的函数值；（ 3 ）在同一图形窗口绘制曲线与，要求曲线为黑色点画线，曲线为红色虚线圆圈；并在图中恰当位 置标注两条曲线的图例；给图形加上标题“ y1andy2 。”五、（本题 15 分）编写 M函数文件，利用 for 循环或 while 循环完成计算函数的任务，并利用该函 数计算时的和（将总程序保存为 test05.m 文件）。六、（本题 13 分）已知求解线

9、性规划模型：的 MATLAB命令为x=linprog （ c,A,b,Aeq,beq,VLB,VUB ）试编写 MATLAB程序，求解如下线性规划问题（将程序保存为test06.m 文件）：问题补充：卷子的地址看不见符号 , 能做就做了一些 .1、标点符号（ ; ）可以使命令行不显示运算结果，（ %）用来表示该行为注释行。2 、下列变量名中（ A）是合法的。（A）char_1;（B）x*y;（C）xy;（D）end3 、为，步长为的向量，使用命令（本题题意不清）创建。4、输入矩阵，使用全下标方式用（本题题意不清）取出元素“，”使用单下标方式用（本题题意不清） 取出元素“ 。”5 、符号表达式中

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

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

12、x=linprog （ c,A,b,Aeq,beq,VLB,VUB ）试编写 MATLAB程序，求解如下线性规划问题（将程序保存为test06.m 文件）：例 2.1 已知 SISO系统的状态空间表达式为 （2-3） 式，求系统的传递函数。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 从系统的传递函数 （2-4） 式求状态空间表达式。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；专业资料整理num,den=ss2tf(A， B，C， D,1)例 2.4 给定系统32s2ss3G(s)，求系32统 的零极点增益模型和状态空间模型，并求其s0.5s2s1单位脉冲响应及单位阶跃响应。解：num=1213;den=10.521;sys=tf(num,den)% 系统的传递函数模型Transferfunction:s3+2s2+s+3s3+0.5s2+2s+1sys1=tf2zp(num,den)% 系统的零极点增益模型 sys1=形式sys2=tf2ss(s

14、ys)% 系统的状态空间模型模型；或用 a,b,c,d=tf2ss(num,den) impulse(sys2)% 系统的单位脉冲响应step(sys2)% 系统的单位阶跃响应 例 3.1 对下面系统进行可控性、可观性分析。 解：a=-1-22;0-11;10-1;b=201'c=120Qc=ctrb(a,b)% 生成能控性判别矩阵 rank(Qc)%求矩阵 Qc 的秩 ans=3%满秩，故系统能控Qo=obsv(a,c)%生成能观测性判别矩阵 rank(Qo)%求矩阵 Qo 的秩ans=3%满秩，故系统能观测 例 3.2 已知系统状态空间方程描述如下： 试判定其稳定性，并绘制出时间响

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(' 系统的零极点模型为 ');z,p,k系统的零极点模型为ifFlagz=1disp(' 系统不稳定 ');elsedisp(' 系统是稳定的 ');end运行结果为 :系统是稳定的step(A,B,C,D) 系统的阶跃响应资源与环境工程学院

16、2008 级硕士研究生 MatLab 及其应用试题 注意，每题的格式均须包含 3 个部分a. 程序(含程序名及完整程序)：b. 运行过程：c. 运行结果：(1) 求解线性规划问题：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(' 线性规划问题求解 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); 运行结果：线性规划问题求解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');专业资料整理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 分)根据专业方向特色和相关科研工作需求，经过与导师商量后，结合一个课题具 体任务，编写一份 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 系统的状态空间表达式为 (2-3) 式，求系统的传递函数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 从系统的传递函数 (2-4) 式求状态空间表达式。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