




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、MATLAB 考试试题 (1)产生一个 1x10 的随机矩阵,大小位于( -5 5 ),并且按照从大到小的顺序排列好! (注:要程序和运行结 果的截屏)答案:a=10*rand(1,10)-5;b=sort(a,descend)1. 请产生一个 100*5 的矩阵,矩阵的每一行都是 1 2 3 4 52. 已知变量: A=ilovematlab ; B=matlab,请找出:(A ) B 在 A 中的位置。(B ) 把 B 放在 A 后面,形成 C= ilovematlabmatlab3. 请修改下面的程序,让他们没有 for 循环语句!A=1 2 3; 4 5 6; 7 8 9;r c=si
2、ze(A);for i=1:1:rfor j=1:1:cif (A(i,j)8 | A(i,j) a=3 4 -7 -12;5 -7 4 2 ;108-5;-65-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);y=cos(0.5+3
3、.*s in (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),-0.3478 -, -0.3478 + 1.0289i ,1.6
4、956 3.0000 a =1.0289i解:a=sqrt(X.A2+Y.A2)的三维图形。范围内画出函数-10,10、在4-10 , 10;X,Y=meshgrid(-10 : 0.5 :10);+eps; Z=sin(a)./a;mesh(X,Y,Z);matlab试卷,求答案一、 选择或填空(每空1、 标点符号( 行。2、下列变量名中2分,共20分)可以使命令行不显示运算结果,)用来表示该行为注释)是合法的。(A) char_1 ;)创建。4、输入矩阵取出元素”5、符号表达式(B)x*y ;(C) xy ;(D) end(的向量,使用命令)取出元素“使用全下标方式用(,步长为为、3 .”
5、,使用单下标方式用中独立的符号变量为(6、M脚本文件和 M函数文件的主要区别是)和(7、在循环结构中跳岀循环,但继续下次循环的命令为(A) return; (B) break ;(C) continue ;二、(本题12分)利用(B) break ;MATLAB数值运算,求解线性方程组(D)keyboad(将程序保存为test02.m 文件)三、(本题20分)利用 MATALAB符号运算完成(将程序保存为test03.m 文件):(1 )创建符号函数(2 )求该符号函数对 的微分;(3 )对趋向于求该符号函数的极限;(4 )求该符号函数在区间上对 的定积分;(5 )求符号方程的解。test04
6、.m 文件):四、(本题20分)编写 MATALAB程序,完成下列任务(将程序保存为 (1 )在区间上均匀地取20个点构成向量 ;(2 )分别计算函数 与 在向量 处的函数值;(3 )在同一图形窗口绘制曲线与 ,要求 曲线为黑色点画线, 曲线为红色虚线圆圈;并在图中恰当位置标注两条曲线的图例;给图形加上标题“y1 and y2 ”。五、(本题 15 分)编写 M 函数文件,利用 for 循环或 while 循环完成计算函数 的任务,并利用 该函数计算 时的和(将总程序保存为 test05.m 文件)。六、(本题 13 分)已知求解线性规划模型:的 MATLAB 命令为x=linprog ( c
7、,A,b,Aeq,beq,VLB,VUB )试编写 MATLAB 程序,求解如下线性规划问题(将程序保存为 test06.m 文件): 问题补充:卷子的地址2 、下列变量名中 ( A ) 是合法的。(A)char_1 ; (B)x*y ;(C) xy ;(D) end3、为,步长为的向量,使用命令 (本题题意不清)创建。)取出元素“ ”,使用单下标方式看不见符号 ,能做就做了一些 .1 、标点符号 (; )可以使命令行不显示运算结果,(% ) 用来表示该行为注释行。4 、输入矩阵 ,使用全下标方式用 ( 本题题意不清 用 ( 本题题意不清 )取出元素“ ”。变量生存期和可见)。C ) 。( 将
8、程序保存为 test02.m 文件 )5 、符号表达式 中独立的符号变量为 ( )6、 M 脚本文件和 M 函数文件的主要区别是 ( 性 ) 和( 函数返回值 7、在循环结构中跳出循环,但继续下次循环的命令为((D) keyboad(A) return; (B) break ; (C) continue ;本题 12 分)利用 MATLAB 数值运算,求解线性方程组三、(本题 20 分)利用 MATALAB 符号运算完成(将程序保存为 test03.m 文件):(1 )创建符号函数syms x(2 )求该符号函数对 的微分;(3 )对 趋向于 求该符号函数的极限;(4 )求该符号函数在区间 上
9、对 的定积分;(5 )求符号方程 的解。四、(本题 20 分)编写 MATALAB 程序,完成下列任务(将程序保存为 test04.m 文件):(1 )在区间 上均匀地取 20 个点构成向量 ;(2 )分别计算函数 与 在向量 处的函数值;(3 )在同一图形窗口绘制曲线与 ,要求 曲线为黑色点画线, 曲线为红色虚线圆圈;并在图中恰当位置标注两条曲线的图例;给图形加上标题“ y1 and y2 ”。五、(本题 15 分)编写 M 函数文件,利用 for 循环或 while 循环完成计算函数 的任务,并利用 该函数计算 时的和(将总程序保存为 test05.m 文件)。六、(本题 13 分)已知求
10、解线性规划模型:的MATLAB命令为x=li nprog (c,A,b,Aeq,beq,VLB,VUB )试编写MATLAB程序,求解如下线性规划问题(将程序保存为 test06.m文件):例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;n um,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
11、 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)32 s s3 2sG(s),求系统的零极点增益模型和状态空间模型,并求其例2.4给定系统 23 2s s. s051单位脉冲响应及单位阶跃响应。解:n um=1 2 1 3;de n=1 0.5 2 1;sys=tf( num,de n)%系统的传递函数模型Tran sfer fun cti on:sA3 + 2 sA2 + s + 3 sA3 + 0.5
12、 sA2 + 2 s + 1%系统的零极点增益模型sys1=tf2zp( num,de n) sys1 =sys2=tf2ss(sys) %系统的状态空间模型模型;或用 a,b,c,d=tf2ss(num,den)形式 %系统的单位脉冲响应impulse(sys2)%step(sys2)系统的单位阶跃响应例3.1对下面系统进行可控性、可观性分析。解: a=-1 -2 2;0 -1 1;1 0 -1;b=2 0 1;c=1 2 0Qc=ctrb(a,b)生成能控性判别矩阵% Qc的秩rank(Qc) % 求矩阵ans = 3满秩,故系统能控%Qo=obsv(a,c) %生成能观测性判别矩阵ran
13、k(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=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=1);系统不稳定 disp(else disp( 系统是稳定的
14、);end运行结果为 : 系统是稳定的 step(A,B,C,D) 系统的阶跃响应资源与环境工程学院 2008 级硕士研究生 MatLab 及其应用试 题注意,每题的格式均须包含 3 个部分a. 程序(含程序名及完整程序) :b. 运行过程:c. 运行结果:(1)求解线性规划问题:minZ 4x x 7x 3125 3213x x x 4321x x 4x 7321x,x 021问各 x 分别取何值 时,Z有何极小值。(10分)i答:fprintf(线性规划问题求解n);f = -4;1;7;A = 3,-1,1;1,1,-4;b = 4,-7;Aeq = 1,1,-1;beq = 5;lb
15、= 0,0,;ub = ;x = linprog(f,A,b,Aeq,beq,lb,ub);xz = f * x;fprintf(MIN z = %f n , z);运行结果:线性规划问题求解Optimizati on termi nated successfully.x =2.25006.75004.0000MIN z = 25.750000x 20.5x, f(x) 1.5 0.25x, 2 x 6,编写一个函数,使其能够产生如下的分段 函数: 0.5, 6 x 2范围的f(x) f(x 2)曲线在x 0,( 10并调用此函数,绘制 分) 0 答:function y=f(x)if x6y
16、=0.5;else y =1.5-0.25*x;endend运行结果x=2f(x)=1x = 0:0.05:2;y = diag(A2(x)*A2(x+2);Plot(x,y);xlabel(bfx);ylabel(bfy);(3)将一个屏幕分4幅,选择合适的步长在右上幅与左下幅绘制出下列函数的图 形。分)10 (.22 yx; ,xy) f(y(-2 x 2,-4 4)(曲面,x)x cos(;(曲线图). 224222图)。答: subplot(2,2,2) ezplot(cos(x)F(1/2) ,-pi/2 pi/2) ylabel(y)subplot(2,2,3) x=-2:0.5:
17、2; y=-4:1:4;ezsurfc(xA2/2A2+yA2/4A2).(4) A是一个維度mx n的矩阵.写一段程序,算出A中有多少个零元素(10分) 答: A= in put (请输入一个矩阵)m, n= size(A);sig=0;for i=1:mfor j=1:n if A(i,j)=0 sig = sig+1;endendend请输入一个矩阵A =0 i 2;i 0 2; 0 0 00i2i02000sigsig5Aa,a,a?.写一段程序 向量,找出A中的最小元素(10分)nii答:A= in put(请输入一个向量 )m,n=sizeAmin =A(i,n);for i=i:
18、nif A(i,i) 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.8090611.269841Plot(x,y);xlabel(时间 t);ylabel(时间/吸附量).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/平衡浓度) x=0.62654 0.80977
19、 0.86585 0.8892 0.87377 0.87564 y=0.59829 1.3273 1.69589 1.89737 2.05503 2.15149Plot(x,y);xlabel(Lg 吸附量);ylabel(Lg平衡浓度)图5d,总结:从图1和图2,分析看可以得到比较理想的对于本次实验的pH和秸秆用量。后面实验是在前面的基础上得到的。图3是吸附动力学反应速率图,从图 中可以看到线性拟合程度很好,符合二级反应速率方程。图4和图5是吸附等温 线作图,看以看出图4的线性拟合较图5的好,说明符合Langmuir吸附等温模 型。例2.1已知SISO系统的状态空间表达式为(2-3)式,求系
20、统的传递函数。A=0 1 0;0 0 1;-4 -3 -2;B=1;3;-6;C=1 0 0;D=0;n um,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)32 s s
21、3 2sG(s),求系统的零极点增益模型和状态空间模给定系统例 2.4. 23S 0.5s 2s 1型,并求其单位脉冲响应及单位阶跃响应。解:n um=1 2 1 3;de n=1 0.5 2 1;sys=tf( num,de n)%系统的传递函数模型Tran sfer fun cti on:sA3 + 2 sA2 + s + 3 sA3 + 0.5 sA2 + 2 s + 1%系统的零极点增益模型sys1=tf2zp( num,de n)sys1 =sys2=tf2ss(sys) %系统的状态空间模型模型;或用a,b,c,d=tf2ss(num,den)形式%系统的单位脉冲响应impulse
22、(sys2)%step(sys2)系统的单位阶跃响应例3.1对下面系统进行可控性、可观性分析。解: a=-1 -2 2;0 -1 1;1 0 -1;b=2 0 1;c=1 2 0Qc=ctrb(a,b) % 生成能控性判别矩阵%求矩阵Qc的秩ran k(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;
23、B=1;0;0;0;C=1 7 24 24;D=0;z,p,k=ss2zp(A,B,C,D,1);Flagz=0; n=length(A); for i=1:n if real(p(i)0Flagz=1;end end disp( 系统的零极点模型为 );z,p,k 系统的零极点模型为 if Flagz=1); 系统不稳定 disp( else disp( 系统是稳定的 ); end 运行结果为 : 系统是稳定的 step(A,B,C,D) 系统的阶跃响应。 1、使用下列哪一个函数可以产生单位矩阵( B )A.zeros B.eye C.rand D.diag2、下列哪一个函数是求模函数 (
24、D ) A.rem B.sign C.fix D.mod3、使用下列哪一个函数可以交换矩阵左右对称位置上的元素( A )A.fliplr B.flipdim C.flipud D.find4、使用下列哪一个函数可以比较字符串,且比较时忽略字符的大小写( D )A.strncmp B.strcmp C.strncmpi D.strcmpi5、要利用图形方式显示元胞数组,则应该使用下列哪一个函数( B )A.cellfun B.cellplot C.celldisp D.cell2mat6、下列哪一个函数可以获取结构字段的数据( B )A.fieldnames B.getfield C.setfi
25、eld D.rmfield7、执行下列哪一条命令后,图形窗体的轴将显示坐标网格线( A )A.grid on B.hold on C.grid off D.hold off8、进行格式化绘图时,使用哪一个函数可以添加图例( B )A.title B.legend C.label D.text9、使用下列哪一条指令可以将图形窗体分割成二行三列,并且将第一行第二列 的绘图区域设置为当前的绘图区域 ( B )A.subplot(2,3,1) B.subplot(2,3,2)C.subplot(2,3,4) D.subplot(2,3,5)10、使用下列哪一个函数可以绘制三维网线图( C )A.sur
26、f B.plot C.mesh D.plot31、 A=1 2 3;4 5 6;B=2 5;8 3(2 分)B =2 58 32、假设向量 A=9 8 7 6 5 4 3 2 1 0 A(1:3:5)(2 分)ans =9 6 A(1 3 5)(2 分)ans =9 7 53、 A=ones(2,2);A(:)=1:4; A*A (2 分)ans =10 1414 20 B=A.*A (2分 )B =1 66 164、使用三元组法,将下列满阵转变为稀疏矩阵1500220-150113000S= 000-60091000000002800解: ir=1 4 2 2 1 3 5 1; jc=1 1 2 3 4 4 4 6; data=15 91 11 3 22 -6 28 -15; s=sparse(ir,jc,data,5,6)s =(1,1) 15(4,1) 91(2,2) 11(2,3) 3(1,4) 22(3,4) -6(5,4) 28(1,6) -155、 A=reshape(1:24,4,6); A(:,2 3 4)=117212182231923420246、使用函数struct创建一个结构此结构名为Student ;有三个字段,分别为name age、grade ;有两条记录,分别为 Way、23、3和Deni、21、1 解: Stude nt
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 职业道德与要求试题及答案
- 基于强化学习的复杂系统建模与仿真实验-洞察及研究
- 康复训练技术创新应用-洞察及研究
- 社保缴费基数监管-洞察及研究
- 2025个体工商户经营资金短缺借款合同
- 物业冲突法律解决策略-洞察及研究
- 2025年学校食堂餐饮服务承包合同协议
- 冲床安全生产教育培训课件
- 冲压模具车间安全培训课件
- 2025合同违约是否需要承担责任
- GB/T 19839-2005工业燃油燃气燃烧器通用技术条件
- GB/T 19478-2018畜禽屠宰操作规程鸡
- (完整版)人工智能介绍课件
- 陶瓷材料的制备课件
- 中职统计基础知识课件
- 预防校园欺凌-共创和谐校园-模拟法庭剧本
- 《人间词话》十则公开课
- 质量管理学课件第1章
- 磁刺激仪技术参数
- 通用机场建设审批程序
- 城市雕塑工程工程量清单计价定额
评论
0/150
提交评论