湖北汽车工业学院MATLAB单输入单输出系统习题答案

合集下载

Matlab编程与应用习题和一些参考答案

Matlab编程与应用习题和一些参考答案

Matlab编程与应用习题和一些参考答案Matlab 上机实验一、二3.求下列联立方程的解⎪⎪⎩⎪⎪⎨⎧=+-+-=-+=++-=--+41025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10];>> b=[4;4;9;4];>> c=a\b4.设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=81272956313841A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=793183262345B ,求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。

>> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8];>> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7];>> C1=A*B'>> C2=A'*B>> C3=A.*B>> inv(C1)>> inv(C2)>> inv(C3)5.设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y ,把x=0~2π间分为101点,画出以x 为横坐标,y 为纵坐标的曲线。

>> x=linspace(0,2*pi,101);>> y=cos(x)*(0.5+(1+x.^2)\3*sin(x));>> plot(x,y,'r')6.产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。

并求该矩阵全体数的平均值和均方差。

(mean var )a=randn(8,6)mean(a)var(a)k=mean(a)k1=mean(k)i=ones(8,6)i1=i*k1i2=a-i1i3=i2.*i2g=mean(i3)g2=mean(g)10.利用帮助查找limit 函数的用法,并自己编写,验证几个函数极限的例子。

MATLAB作业1 参考答案

MATLAB作业1 参考答案

MATLAB 作业1参考答案1 、在你的机器上安装MATLAB 语言环境,并键入demo 命令,由给出的菜单系统和对话框原型演示程序,领略MATLAB 语言在求解数学问题方面的能力与方法。

【求解】略.2、启动MATLAB 环境,并给出语句tic, A=rand(500); B=inv(A); norm(A*B-eye(500)),toc ,试运行该语句,观察得出的结果,并利用help 命令对你不熟悉的语句进行帮助信息查询,逐条给出上述程序段与结果的解释。

【求解】在MATLAB 环境中感触如下语句,则可以看出,求解500500⨯随机矩阵的逆,并求出得出的逆矩阵与原矩阵的乘积,得出和单位矩阵的差,得出范数。

一般来说,这样得出的逆矩阵精度可以达到1210-。

>> tic, A=rand(500); B=inv(A); norm(A*B-eye(500)), tocans =1.2333e-012Elapsed time is 1. seconds.3、试用符号元素工具箱支持的方式表达多项式5432()34236f x x x x x x =+++++,并令11s x s -=+,将f(x) 替换成s 的函数。

【求解】可以先定义出f 函数,则由subs() 函数将x 替换成s 的函数>> syms s x;f=x^5+3*x^4+4*x^3+2*x^2+3*x+6;F=subs(f,x,(s-1)/(s+1)) F =(s-1)^5/(s+1)^5+3*(s-1)^4/(s+1)^4+4*(s-1)^3/(s+1)^3+2*(s-1)^2/(s+1)^2+3*(s-1)/(s+1)+6syms x sx=(s-1)./(s+1);y=x.^5+3*x.^4+4*x.^3+2*x.^2+3*x+6;ysimple(y)y=19 - (72*s^4 + 120*s^3 + 136*s^2 + 72*s + 16)/(s + 1)^54、用MATLAB 语句输入矩阵A 和B123414233241432141322314,234123324114324132234114j j j j j j j j A B j j j j j j j j ++++⎡⎤⎡⎤⎢⎥⎢⎥++++⎢⎥⎢⎥==⎢⎥⎢⎥++++⎢⎥⎢⎥++++⎣⎦⎣⎦前面给出的是4 ×4 矩阵,如果给出A(5; 6) = 5 命令将得出什么结果?【求解】用课程介绍的方法可以直接输入这两个矩阵>> A=[1 2 3 4; 4 3 2 1; 2 3 4 1; 3 2 4 1]A =1 2 3 44 3 2 12 3 4 13 24 1若给出A(5,6)=5 命令,虽然这时的行和列数均大于B 矩阵当前的维数,但仍然可以执行该语句,得出>> A(5,6)=5A =1 2 3 4 0 04 3 2 1 0 02 3 4 1 0 03 24 1 0 00 0 0 0 0 5复数矩阵也可以用直观的语句输入>> B=[1+4i 2+3i 3+2i 4+1i; 4+1i 3+2i 2+3i 1+4i;2+3i 3+2i 4+1i 1+4i; 3+2i 2+3i 4+1i 1+4i];B =1.0000 + 4.0000i2.0000 +3.0000i 3.0000 + 2.0000i4.0000 + 1.0000i4.0000 + 1.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 1.0000 + 4.0000i2.0000 +3.0000i 3.0000 + 2.0000i4.0000 + 1.0000i 1.0000 + 4.0000i3.0000 + 2.0000i 2.0000 + 3.0000i4.0000 + 1.0000i 1.0000 + 4.0000i5、假设已知矩阵A ,试给出相应的MATLAB 命令,将其全部偶数行提取出来,赋给B 矩阵,用A =magic(8) 命令生成A 矩阵,用上述的命令检验一下结果是不是正确。

matlab基础测习题(带答案)

matlab基础测习题(带答案)

精心整理Matlab基础练习题常量、变量、表达式1、MATLAB中,下面哪些变量名是合法的?()2、3、4、5、)6、小数7、8、在MATLAB中,a=1,b=i,则a占_8__个字节,b占_16_个字节,c占________字节。

9、在MATLAB中,inf的含义是__无穷大__,nan的含义是__非数(结果不定)___。

数组1、在MATLAB中,X是一个一维数值数组,现在要把数组X中的所有元素按原来次序的逆序排列输出,应该使用下面的()指令。

(A)X[end:1] (B)X[end:-1:1] (C)X(end:-1:1)(D)X(end:1) 2、在MATLAB中,A是一个字二维数组,要获取A的行数和列数,应该使用的MATLAB的命令是()。

(A)class(A) (B)sizeof(A) (C)size(A)(D)isa(A)3、4、这5、6、___linspace(0,5,50)___7、在MATLAB中,A=[0:1/2:2]*pi,那么sin(A)=___[010-10]_____。

8、在MATLAB中,A=[1,2,3;4,5,6;7,8,0],B=[2,1,6;8,5,2;14,2,1]。

写出下面MATLAB语句执行的结果:(为节省篇幅,把矩阵写成mat2str的形式)(1)A==B [000;010;000](2)A.*B [2218;322512;98160](3)A(:)' [147258360](4)A(1,:)*B(:,3) 139、在MATLAB中,写出下面MATLAB语句执行的结果:(1)c lear,A=ones(2,6)A=(5)B=A(:,end:-1:1)B=211713951231915117310、请编写一段matlab程序,完成以下功能:(1) 生成一个100行,200列的二维随机数组;A=rand(100,200);(2) 找出数组A 中所有大于0.49且小于0.51的元素的单下标;Idx=find(A(:)>0.49&A(:)<0.51)(3) 数组A 中满足(2)中的条件的元素有多少个?11、12、 22cos 0.7sin (1)x y x x ⎡⎤=+⎢⎥+⎣⎦,编写一段matlab 程序,要求如下(1) 在[0,2]π区间,每隔0.01取一x 数值,计算出相应的y 的函数值;x=0:0.01:2*pi;y=(0.7+2*cos(x)./(1+x.^2)).*sin(x)(2)根据MATLAB计算出的数据,找出在[0,2]π内该函数的极小值的坐标。

湖北汽车工业学院MATLAB单输入单输出系统习题答案

湖北汽车工业学院MATLAB单输入单输出系统习题答案

湖北汽车⼯业学院MATLAB单输⼊单输出系统习题答案考察下列单输⼊单输出系统:x’’+3x’+2x=ux(0)=1 ,x’(0)=0输⼊:u=u(t), 输出:x=x(t), 0编程做以下四件事:1、u=1时,⽤符号法求输出(dslove)2、u=1时,⽤数值法求输出(ode45)3、⽤传递函数求输出阶跃响应(tf,step)4、⽤Simulink法求输出阶跃响应(mdl)给出源程序代码或模块,以及运⾏结果分析。

题⽬解答如下:1、运⾏代码:x=dsolve('D2x+3*Dx+2*x=1','x(0)=1,Dx(0)=0','t')>> ezplot(x,[0,10])输出结果:x =exp(-t) - exp(-2*t)/2 + 1/22、在脚本⽂件中建⽴函数:function xdot=DyDt(t,x)mu=3;xdot=[x(2);1-2*x(1)-mu*x(2)];在命令窗⼝中运⾏指令:tspan=[0,10];x0=[1,0];[tt,xx]=ode45(@DyDt,tspan,x0);plot(tt,xx(:,1));xlabel('t'),title('x(t)')由于函数ode45采⽤的数值积分法,故tt,xx的结果为⼀系列插值点,结果不易处理,在此只输出图形。

3、在命令窗⼝中运⾏指令:clear all;G1=tf([1],[1,3,2])step(G1,'m*');legend('*G1')输出结果:G1 =1-------------s^2 + 3 s + 24、建⽴siso.mdl⽂件,在⽂件中建⽴模型如下:。

自动控制原理习题全解及MATLAB实验 第6章习题解答

自动控制原理习题全解及MATLAB实验 第6章习题解答

系统开环传递函数为 G0 s
s0.1s
K
10.2s
1
,要求:
(1)系统响应斜坡信号 r(t)=t 时,稳态误差 ess 0.01 ;
(2) 系统相位裕量 ' 40 。
试用分析法设计一个串联滞后-超前校正装置。
解:(1)系统为Ⅰ型系统,在单位斜坡信号下
分稳态误差为
essr
1 k
令 essr
稳态性能与动态性能? 答:PID 兼有 PI、PD 控制的特点,它相当于提供了一个积分环节与两个一阶微分环节。
积分环节改善稳态性能,两个一阶微分环节改善动态性能。 试分别叙述利用比例负反馈和微分负反馈包围振荡环节所起到的作用。
答:二阶振荡环节的频率特性为
1
T 2S 2 2 S 1
用比例负反馈 H(s)=h
0.2s 1 0.0143s 1
(5) Gc (s)
s 1 14s 1
0.2s 1 0.0143s
(6)
G
k
(s)
s(14s
100(s 1) 1)(0.1s 1)(0.0143s
1)
' 180 [90 arctan 7 arctan(14 7) arctan(0.1 7) arctan(0.0143 7)] 41.9 40
10lg( 12)
6dB
,
最后得出 c' m 4.47rad/s>4.4rad/s
(4) 确定校正装置的转折频率
1 m
2.2rad/s ,2 m
8.8rad/s ,T 1 0.45s , 1
G(s)=
s
2.2 s
1 1
0.45s 0.11s
1 1

(完整word版)含答案《MATLAB实用教程》

(完整word版)含答案《MATLAB实用教程》

第二章 MATLAB 语言及应用实验项目实验一 MATLAB 数值计算三、实验内容与步骤1.创建矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a(1(2)用(3)用(42.矩阵的运算(1)利用矩阵除法解线性方程组。

⎪⎪⎩⎪⎪⎨⎧=+++=-+-=+++=+-12224732258232432143214321421x x x x x x x x x x x x x x x 将方程表示为AX=B ,计算X=A\B 。

(2)利用矩阵的基本运算求解矩阵方程。

已知矩阵A 和B 满足关系式A -1BA=6A+BA ,计算矩阵B 。

其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=7/10004/10003/1A ,Ps: format rata=[1/3 0 0;0 1/4 0;0 0 1/7];b=inv(a)*inv(inv(a)-eye(3))*6*a(3)计算矩阵的特征值和特征向量。

已知矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=1104152021X ,计算其特征值和特征向量。

(4)Page:322利用数学函数进行矩阵运算。

已知传递函数G(s)=1/(2s+1),计算幅频特性Lw=-20lg(1)2(2w )和相频特性Fw=-arctan(2w),w 的范围为[0.01,10],按对数均匀分布。

3.多项式的运算(1)多项式的运算。

已知表达式G(x)=(x-4)(x+5)(x 2-6x+9),展开多项式形式,并计算当x 在[0,20]内变化时G(x)的值,计算出G(x)=0的根。

Page 324(2)多项式的拟合与插值。

将多项式G(x)=x 4-5x 3-17x 2+129x-180,当x 在[0,20]多项式的值上下加上随机数的偏差构成y1,对y1进行拟合。

对G(x)和y1分别进行插值,计算在5.5处的值。

Page 325 四、思考练习题1.使用logspace 函数创建0~4π的行向量,有20个元素,查看其元素分布情况。

Ps: logspace(log10(0),log10(4*pi),20) (2) sort(c,2) %顺序排列 3.1多项式1)f(x)=2x 2+3x+5x+8用向量表示该多项式,并计算f(10)值. 2)根据多项式的根[-0.5 -3+4i -3-4i]创建多项式。

MATLAB复习题(答案仅供参考)1

MATLAB复习题(答案仅供参考)1

一、填空题1. MA TLAB命令窗口中可用clear 命令清除工作区中的变量;用clc 命令清除命令窗口中的内容。

2. MA TLAB中的运算符包括数学运算符、矩阵运算符、数组运算符。

3 MA TLAB的含义为矩阵运算。

4. MA TLAB中的函数类型有嵌套函数、主函数、子函数、私有函数和匿名函数。

5. MA TLAB中用系统预定义的特殊变量nargin 表示函数输入参数个数。

6. size命令的功能确定数组大小。

7. MA TLAB中用linspace 符号定义单元数组。

8. MA TLAB中用hold on 命令打开图形保持功能。

9 MA TLAB中如果一个变量名为123,该变量名是否合法否(填是或否)。

10. MA TLAB中用inv 函数完成矩阵的求逆运算,用det 函数求矩阵的行列式。

11.已知C={2+5i,’Hi’;1:4,[1 2]} 输入s=C(3) 即s=__'_Hi__'___.12.先输入P=struct(‘name’,’ Jack‘,’num’,43) ,再输入s= 即s=___Jack___.13.当x=[ 1 2 3; 4 5 6] , y=[ 5 8; 5 8 ] 时,__不能__(能﹨不能)使用plot 命令绘图。

14.sym 或syms _能_______(能﹨不能)创建符号矩阵。

15.已知x=3:2:24,则x有__11_个元素。

16.if结构的开始是“if”命令,结束是_end____命令。

17.已知x为一个向量,计算ln(x)的MA TLAB命令是计算_log(x)___。

18.命令:>>subplot(3,3,x)中有__9_个子图。

19.已知s=‘显示“hey”’,则s的元素个数是___9__。

20.function y=f()global WW=2y=f1(5)function y=f()global Wy=w+W创建以上函数文件,在命令窗口中运行“y=f”命令则显示__出错___。

湖北汽车工业学院MATLAB试题样本

湖北汽车工业学院MATLAB试题样本

☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆☆ ☆☆☆☆☆☆☆☆ ☆☆☆☆☆☆☆☆ ☆☆☆☆☆☆☆☆ ☆☆☆☆☆☆☆☆湖北汽车工业学院 《系统辨识与计算机仿真》考试试卷(题)课程编号 020300 考试形式考试日期题号 一 二 三 四 五 六 七 八 九 十 总分 得分 阅卷一、 简答题(每题6分,共48分)1.在MATLAB 计算机仿真过程中,符号运算与数值运算的区别是什么?2. 在MATLAB 计算机仿真过程中,插值运算的目的是什么?3.请叙述用示波器、workspace I/O 页和To workspace 模块把计算机仿真结果的数据送入到MATLAB 工作空间,其变量名各为什么。

4.MATLAB 计算机仿真过程中,在装机过程中没有勾选自动控制系统工具箱,将会产生什么后果运行?5. 控制系统的仿真过程是什么?6. 在MATLAB 命名规则中,变量名、函数名是对大小写是否敏感?命名规则是什么?7. 请叙诉数组运算与矩阵运算的区别?8. 怎么提高SIMULINK 系统集成度?二、 根据如下命令画出图形(10分)x=0:0.05:7; title(‘cos(x)’) y1=cos(x); subplot(1 2 2) y2=2*sin(x); plot(x,y2); subplot(1 2 1) title(‘sin(x)’) plot(x,y1);考生班号 学号 姓名 (座位号 )装订线内不得答题三、 某正弦信号为y=0.8+sin(5πt),试用函数dtrend 对信号进行处理,使之均值为零,输出处理前与处理后的信号曲线。

MATLAB系统仿真试题复习资料,有答案

MATLAB系统仿真试题复习资料,有答案

PART I 教材复习要点第一章 MATLAB 入门测试1.11.MATLAB 的命令窗口的作用是什么?编辑/调试窗口?图象窗口? MATLAB 命令窗口是MATLAB 起动时第一个看到的窗口,用户可以在命令窗口提示符">>"后面输入命令,这些命令会被立即执行。

命令窗口也可以用来执行M 文件。

编辑/调试窗口是用来新建,修改或调试M 文件的。

图像窗口用来显示MATLAB 的图形输出。

2.列出几种不同的得到MATLAB 帮助的方法。

(1)在命令窗口中输入help <command_name>,本命令将会在命令窗口中显示关于些命令的有信息。

(2)在命令窗口中输入lookfor <keyword>,本命令将会在命令窗口中显示所有的在第一注释行中包含该关键字的命令和函数。

(3)通过在命令窗口输入helpwin 或helpdesk 启动帮助浏览器,或者是在启动板中选择"Help"。

帮助浏览器包含了基于超文本的MATLAB 所有特性的描述,HTML 或PDF 格式的在线手册,这是MATLAB 最全面的帮助资源。

3.什么是工作区?在同一工作区内,你如何决定它里面存储了什么?工作区是命令、M 文件或函数执行时被MATLAB 使用的变量或数组的收集器,所有命令都在命令窗口(所有的脚本文件也是从命令窗口执行)共享公共工作区,因此它们也共享所有变量,工作区的内容可以通过whos 命令来查看,或者通过工作区浏览器来图形化地查看。

4.你怎样清空MATLAB 工作区内的内容?要清除工作区的内容,只需在命令窗口中输入clear 或clear variables 即可。

第二章MATLAB 基础测试2.11.数组,矩阵,向量有什么区别?数组是在内存中被组织成行和列的数据集合,只有一个名称,数据要通过在数组名后面圆括号里加上表示数据所在行和列的数字来访问。

大一matlab考试题及答案

大一matlab考试题及答案

大一matlab考试题及答案一、选择题(每题2分,共20分)1. MATLAB中用于创建一个3x3的单位矩阵的命令是:A. eye(3)B. ones(3)C. zeros(3)D. identity(3)答案:A2. 下列哪个函数用于计算矩阵的行列式?A. detB. invC. eigD. rank答案:A3. MATLAB中,用于绘制函数f(x)=x^2的图形的命令是:A. plot(x, x^2)B. plot(x^2)C. plot(x, f(x))D. plot(f(x))答案:A4. 以下哪个命令用于计算矩阵的特征值?A. eigB. invC. detD. rank5. MATLAB中,用于求解线性方程组Ax=b的命令是:A. solve(A, b)B. A\bC. b/AD. A/b答案:B6. MATLAB中,用于创建一个1到10的向量的命令是:A. 1:10B. linspace(1, 10)C. ones(1, 10)D. zeros(1, 10)答案:A7. 下列哪个函数用于计算矩阵的逆?A. detB. invC. eigD. rank答案:B8. MATLAB中,用于计算向量元素的总和的命令是:A. sumB. meanC. maxD. min答案:A9. MATLAB中,用于计算两个矩阵的点积的命令是:B. *C. .*D. cross答案:A10. MATLAB中,用于生成一个随机数矩阵的命令是:A. randB. randnC. randiD. randperm答案:A二、填空题(每题3分,共15分)1. MATLAB中,用于创建一个5x5的全零矩阵的命令是____。

答案:zeros(5)2. MATLAB中,用于计算矩阵的秩的命令是____。

答案:rank3. MATLAB中,用于计算向量的范数的命令是____。

答案:norm4. MATLAB中,用于计算矩阵的迹的命令是____。

matlab基础练习题及答案

matlab基础练习题及答案

matlab基础练习题及答案第1章MATLAB基础1.4MATLAB操作桌面有几个窗口?如何使某个窗口脱离桌面成为独立窗口?又如何将脱离出去的窗口重新放置到桌面上?答:在MATLAB操作桌面上有五个窗口。

在每个窗口的右上角有两个小按钮,一个是关闭窗口的Cloe按钮,一个是可以使窗口成为独立窗口的Undock按钮,点击Undock按钮就可以使该窗口脱离桌面成为独立窗口。

在独立窗口的Dektop菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。

1.5如何启动M文件编辑/调试器?答:在操作桌面上选择“建立新文件”或“打开文件”操作时,M文件编辑/调试器将被启动。

在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。

1.6存储在工作空间中的数组能编辑吗?如何操作?答:存储在工作空间的数组可以通过数组编辑器进行编辑:在工作空间浏览器中双击要编辑的数组名打开数组编辑器,再选中要修改的数据单元,输入修改内容即可。

1.7命令历史窗口除了可以观察前面键入的命令外,还有什么用途?答:命令历史窗口除了用于查询以前键入的命令外,还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。

1.8如何设置当前目录和搜索路径,在当前目录上的文件和在搜索路径上的文件有什么区别?答:当前目录可以在当前目录浏览器窗口左上方的输入栏中设置,搜索路径可以通过选择操作桌面的file菜单中的SetPath菜单项来完成。

在没有特别说明的情况下,只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用,如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件,如果没有特别说明,数据文件将存储在当前目录上。

1.9在MATLAB中有几种获得帮助的途径?答:(1)帮助浏览器:选择view菜单中的Help菜单项或选择Help 菜单中的MATLABHelp菜单项可以打开帮助浏览器。

(完整word版)Matlab考试题库+答案

(完整word版)Matlab考试题库+答案

填空题1、标点符号; %—用来表示该行为注释行。

可以使命令行不显示运算结果,2、x为0~4pi,步长为0.1pi的向量,使用命令 x=0:0.1*pi:4*pi创建。

3、输入矩阵A=,使用全下标方式用A(2,2)取出元素“-5”,使用单下标方式用A(5)取出元素“-5”。

4、符号表达式sin(2*a+t)+m中独立的符号变量为t。

5、M脚本文件和M函数文件的主要区别是M脚本文件没有函数定义和M函数文件有函数定义_______。

6. 设x是一维数组,x的倒数第3个元素表示为x(_end-2_)设y为二维数组,要删除y的第34行和48列,可使用命令y(34,:)=[] ;y(:,48)=[];7. 将变量x以Ascii文本格式存储到文件fname.txt,应使用命令save _x ;8. 在while 表达式, 语句体, End 循环语句中,表达式的值非零时表示循环条件为真,语句体将被执行,否则跳出该循环语句;9.要从键盘读入一个字符串并赋值给变量x,且给出提示“Who is she?”,应使用命令x=input(‘Who is she?’,’s’);10.设A=和B=和C=均为m*n矩阵,且存在于WorkSpace中,要产生矩阵D=,可用命令D=(A-C)/B.^C,计算可用命令det(inv(A’*B)11. 在MATLAB命令窗口中的“>>”标志为MATLAB的命令行提示符,“│”标志为输入提示符。

12.已知A=[1 2 3;4 5 0;7 8 9];B=[1 0 3;1 5 0;0 1 2];写出下列各指令运行的结果。

A+B;A.*B;A==B ;ans= [2,2,6;5,10,0;7,9,11]ans= [1,0,9;4,25,0;0,8,18]ans= [1,0,1;0,1,1;0,0,0]13.已知A是矩阵,求A的对角矩阵函数是 diag ;求A的下三角矩阵函数是 tril。

matlab课后习题答案1到6章

matlab课后习题答案1到6章

欢迎共阅习题二1.如何理解“矩阵是MATLAB 最基本的数据对象”? 答:因为向量可以看成是仅有一行或一列的矩阵,单个数据(标量)可以看成是仅含一个元素的矩阵,故向量和单个数据都可以作为矩阵的特例来处理。

(2) 删除矩阵A 的第7号元素。

答:A(7)=[](3) 将矩阵A 的每个元素值加30。

答:A=A+30;(4) 求矩阵A 的大小和维数。

答:size(A);ndims(A);(5) 将向量 t 的0元素用机器零来代替。

答:t(find(t==0))=eps; (6) 将含有12个元素的向量 x 转换成34⨯矩阵。

答:reshape(x,3,4); (7) 求一个字符串的ASCII码。

答:abs(‘123’); 或E 。

答:B=A(1:3,:); C=A(:,1:2); D=A(2:4,3:4);E=B*C;(2) 分别求E<D 、E&D 、E|D 、~E|~D 和find(A>=10&A<25)。

答:E<D=010001⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,E&D=110111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,11⎡⎤⎢⎥答:student(1).id='0001';student(1).name='Tom'; student(1).major='computer';student(1).grade=[89,78,67,90,86,85]; 8.建立单元矩阵B 并回答有关问题。

B{1,1}=1;B{1,2}='Brenden';B{2,1}=reshape(1:9,3,3); B{2,2}={12,34,2;54,21,3;4 ,23,67};(1)size(B)和ndims(B)的值分别是多少?答:size(B) 的值为2, 2。

ndims(B) 的值为2。

(2)B(2)和B(4)的值分别是(2)建立5×6随机矩阵A,其元素为[100,200]范围内的随机整数。

MATLAB仿真课后习题答案

MATLAB仿真课后习题答案

MATLAB仿真课后习题答案MATLAB仿真课后习题答案第一章1.要求在闭区间[0,2Π]上产生具有10个等间距采样点的一维数组。

试用两种不同的指令实现。

解答:方法1:a=0:2*pi/9:2*pi方法2:a1=linspace(0,2*pi,10)4.任意建立矩阵A,然后找出在[10,20]区间的元素的位置。

解答:A=[4,15,-45,10,6;56,0,17,-45,0]A =4 15 -45 10 656 0 17 -45 0>> find(A>=10&A<=20)ans =367第二章M文件的2种形式:命令文件和函数文件。

命令文件没有输入输出参数,执行时只需在命令窗口中键入文件名回车即可;而函数文件是一条以function语句作为引导。

即文件的第一行为function[返回参数1,返回参数2,…]=函数名(输入参数1,输入参数2,…)。

这一行的有无是区分命令文件与函数文件的重要标志。

函数文件可以接受输入变量,还可以返回输出变量,执行时需在命令窗口中以固定格式调用函数方可。

习题22-1.编写程序,建立向量N=[1,2,3,4,5],然后利用向量N产生下列向量:(1) 2,4,6,8,10(2) 1/2, 1, 3/2, 2, 5/2(3) 1, 1/2, 1/3, 1/4, 1/5(4) 1, 1/4, 1/9, 1/16, 1/25解答:>> N1=2*N >> N2=N/2 >> N3=1./N >> N4=1./N.^2 2-2从键盘输入一个3位数的整数,将它反向输出。

如输入639,输出936。

输入一个百分制成绩,要求输出成绩等级A、B、C、D、E。

其中90~100分为A,80~89分为B,70~79分为C,60~69分为D,60分以下为E。

要求:(1)分别用if语句和switch语句实现。

(2)输入百分制成绩后要判断该成绩的合理性,对不合理的成绩应输出出错信息。

MATLAB 语言及其用 实验(答案)

MATLAB 语言及其用 实验(答案)

《MATLAB 语言及其用》实验指导书目录实验一Matlab 使用方法和程序设计........................实验二控制系统的模型及其转换.............................实验三控制系统的时域、频域和根轨迹分析...........实验四动态仿真集成环境-Simulink.........................实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法;2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容:1、帮助命令使用help命令,查找 sqrt(开方)函数的使用方法;在 CommandWindowL里输入help,接在在search里输入sqr即可。

sqrtSquare rootSyntaxB = sqrt(X)DescriptionB = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results.TipsSee sqrtm for the matrix square root. Examplessqrt((-2:2)')ans =0 + 1.4142i0 + 1.0000i1.00001.4142See Alsonthroot | realsqrt | sqrtm2、矩阵运算(1)矩阵的乘法已知A=[1 2;3 4]; B=[5 5;7 8];求A^2*BA=[1 2;3 4];B=[5 5;7 8];C=A^2*B>> format compactC =105 115229 251(2)矩阵除法已知 A=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];A\B,A/BA=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];C=A\B,D=A/BC =1.0e+016 *0.3152 -1.2609 0.9457-0.6304 2.5218 -1.89130.3152 -1.2609 0.9457D =1.0000 1.0000 1.00004.0000 2.5000 2.00007.0000 4.0000 3.0000(3)矩阵的转置及共轭转置已知A=[5+i,2-i,1;6*i,4,9-i];求A.', A'A=[5+i,2-i,1;6*i,4,9-i];B=A.', C=A'B =5.0000 + 1.0000i 0 +6.0000i2.0000 - 1.0000i 4.00001.0000 9.0000 - 1.0000iC =5.0000 - 1.0000i 0 -6.0000i2.0000 + 1.0000i 4.00001.0000 9.0000 + 1.0000i(4)使用冒号表达式选出指定元素已知: A=[1 2 3;4 5 6;7 8 9];求A中第3列前2个元素;A中所有列第2,3行的元素;A=[1 2 3;4 5 6;7 8 9];B1=A([1,2],[3])B2=A([2,3],:)B1 =36B2 =4 5 67 8 9方括号[]用magic函数生成一个4阶魔术矩阵,删除该矩阵的第四列A=magic(4)B=A(:,[1,2,3])或A=magic(4)A(:,4)=[]A =16 2 3 135 11 10 89 7 6 124 14 15 1B =16 2 35 11 109 7 64 14 153、多项式(1)求多项式4xxp的根=x)2(3--Y=[1 0 -2 -4];S=roots(Y)S =2.0000-1.0000 + 1.0000i-1.0000 - 1.0000i(2)已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ,求矩阵A的特征多项式;把矩阵A作为未知数代入到多项式中;A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4]P=poly(A)polyval(P,A)A =1.2000 3.0000 5.0000 0.90005.0000 1.7000 5.00006.00003.0000 9.0000 0 1.00001.00002.00003.00004.0000P =1.0000 -6.9000 -77.2600 -86.1300 604.5500ans =1.0e+003 *0.3801 -0.4545 -1.9951 0.4601-1.9951 0.2093 -1.9951 -2.8880-0.4545 -4.8978 0.6046 0.43530.4353 0.0840 -0.4545 -1.16174、基本绘图命令(1)绘制余弦曲线 y=cos(t),t∈[0,2π](2)在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5),t∈[0,2π](1)t=[0:0.05:2*pi];y=cos(t);plot(t,y)-1-0.8-0.6-0.4-0.20.20.40.60.81(2)t=[0:0.05:2*pi];y1=cos(t-0.25);y2=sin(t-0.5);plot(t,y1)hold onplot(t,y2)5、基本绘图控制绘制[0,4π]区间上的x1=10sint曲线,并要求:(1)线形为点划线、颜色为红色、数据点标记为加号;(2)坐标轴控制:显示范围、刻度线、比例、网络线(3)标注控制:坐标轴名称、标题、相应文本;t=[0:0.1:4*pi];x1=10*sin(t);plot(t,x1,'r-.+'); %画图,显示红色、点划线、标记加号;axis([0,15,-10,10]); %定义显示范围,横轴为[0,15],纵轴为[-10,10];title('曲线x1=10sint'); %显示标题;xlabel('T轴');ylabel('X1轴'); %显示坐标轴名称;set(gca,'xminortick','on');set(gca,'yminortick','on'); %显示刻度线;grid on %显示网络线T 轴X 1轴6、基本程序设计(1)编写命令文件:计算1+2+…+n<2000 时的最大n 值; (2)编写函数文件:分别用for 和while 循环结构编写程序,求2的0到n 次幂的和。

MATLAB平时练习题(附答案)

MATLAB平时练习题(附答案)

一、填空题1、MATLAB常用操作界面包括、工作空间窗口、、、内存数组编辑器、M文件编辑/调试器、帮助导航/浏览器、图形窗口等。

2、MATLAB中Inf或inf表示、NaN或nan表示、nargout表示。

3、MATLAB中逗号主要用作;用作输入量与输入量之间的分隔符;。

4、工作空间浏览器主要用于内存变量的、和。

5、MATLAB实现将全下标转换为单下标的指令为、据单下标换算出全下标的指令为。

6、二维数组的标识有、、“逻辑1”标识。

二、简述下列命令在MATLAB中的作用1.clf %清除图对象clear %清除工作空间内的所有变量clc %清除当前屏幕上显示的所有内容,但不清除工作空间中的数据2.ceil %沿+∞方向取整factor %符号计算的因式分解3.sph2cart %球坐标变为直角坐标pow2 %2的幂4.digits(d) %设置当前的可变算术精度的位数为整数 d 位vpa(pi,100) %设置当前pi的小数位数为1005.box on %打开框状坐标轴开grid off %网格关6.logspace %对数分度向量cat %串接成高维数组7.sym2poly %符号多项式转变为双精度多项式系数向量poly2sym %双精度多项式系数转变为向量符号多项式8.plot3 %三维线图poly2str %以习惯方式显示多项式9.bar %二维直方图pie %二维饼图10.zoom on %打开图形缩放模式 edit %M文件编辑11.whos %对当前工作空间变量的信息进行列表figure %生成图形窗口12.cart2sph %直角坐标变为球坐标pol2cart % 极或柱坐标变为直角坐标13.diff %数值差分、符号微分14.ezplot3 %画三维曲线的简捷指令fix %向零取整factor %符号计算的因式分解三简答题1.叙述MATLAB的主要窗口包括哪些?其作用是什么?命令窗口(Commend Window)——供输入命令,实现计算或绘图功能。

【免费下载】自动控制系统计算机仿真习题参考答案

【免费下载】自动控制系统计算机仿真习题参考答案
3 1 2
0 A1= -0.2
0.4
-0.6667 1
-0.3333
0.3333
-0.2
0.0667

以及特征值
2 4 3 求出 B 阵的秩为 3、转置为 B= 1 1 3 、行列式为 56 、逆矩阵
8 4 2
B
1
=
-0.1786

0.0714
y 1.3 2
3X
1.3(s 0.7692 1.31i)(s 0.7692 1.31i)
(s 0.1153 1.1642i)(s 0.1153 1.1642i)(s 0.7307)
对全部高中资料试卷电气设备,在安装过程中以及安装结束后进行高中资料试卷调整试验;通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关,系电,力根通保据过护生管高产线中工敷资艺设料高技试中术卷资,配料不置试仅技卷可术要以是求解指,决机对吊组电顶在气层进设配行备置继进不电行规保空范护载高与中带资负料荷试下卷高问总中题体资,配料而置试且时卷可,调保需控障要试各在验类最;管大对路限设习度备题内进到来行位确调。保整在机使管组其路高在敷中正设资常过料工程试况中卷下,安与要全过加,度强并工看且作护尽下关可都于能可管地以路缩正高小常中故工资障作料高;试中对卷资于连料继接试电管卷保口破护处坏进理范行高围整中,核资或对料者定试对值卷某,弯些审扁异核度常与固高校定中对盒资图位料纸置试,.卷保编工护写况层复进防杂行腐设自跨备动接与处地装理线置,弯高尤曲中其半资要径料避标试免高卷错等调误,试高要方中求案资技,料术编试交写5、卷底重电保。要气护管设设装线备备置敷4高、调动设中电试作技资气高,术料课中并3中试、件资且包卷管中料拒含试路调试绝线验敷试卷动槽方设技作、案技术,管以术来架及避等系免多统不项启必方动要式方高,案中为;资解对料决整试高套卷中启突语动然文过停电程机气中。课高因件中此中资,管料电壁试力薄卷高、电中接气资口设料不备试严进卷等行保问调护题试装,工置合作调理并试利且技用进术管行,线过要敷关求设运电技行力术高保。中护线资装缆料置敷试做设卷到原技准则术确:指灵在导活分。。线对对盒于于处调差,试动当过保不程护同中装电高置压中高回资中路料资交试料叉卷试时技卷,术调应问试采题技用,术金作是属为指隔调发板试电进人机行员一隔,变开需压处要器理在组;事在同前发一掌生线握内槽图部内 纸故,资障强料时电、,回设需路备要须制进同造行时厂外切家部断出电习具源题高高电中中源资资,料料线试试缆卷卷敷试切设验除完报从毕告而,与采要相用进关高行技中检术资查资料和料试检,卷测并主处且要理了保。解护现装场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况,然后根据规范与规程规定,制定设备调试高中资料试卷方案。

MATLAB及应用参考答案

MATLAB及应用参考答案

《MATLAB及应用》上机作业学院名称:机械工程学院专业班级:测控1201学生姓名:学生学号:201 年 4 月《MATLAB及应用》上机作业要求及规范一、作业提交方式:word文档打印后提交。

二、作业要求:1.封面:按要求填写学院、班级、姓名、学号,不要改变封面原有字体及大小。

2.内容:只需解答过程(结果为图形输出的可加上图形输出结果),不需原题目;为便于批阅,题与题之间应空出一行;每题答案只需直接将调试正确后的M文件内容复制到word 中(不要更改字体及大小),如下所示:%作业1_1clcA=[1 2 3 4;2 3 5 7;1 3 5 7;3 2 3 9;1 8 9 4];B=[1+4*i 4 3 6 7 7;2 3 3 5 5 4+2*i;2 6+7*i 5 3 4 2;1 8 9 5 4 3];C=A*BD=C(4:5,4:6)三、大作业评分标准:1.提交的打印文档是否符合要求;2.作业题的解答过程是否完整和正确;3.答辩过程中阐述是否清楚,问题是否回答正确;4.作业应独立完成,严禁直接拷贝别人的电子文档,发现雷同者都以无成绩论处。

作业11、用MATLAB 可以识别的格式输入下面两个矩阵12342357135732391894A ⎛⎫⎪⎪ ⎪= ⎪⎪ ⎪⎝⎭,144367723355422675342189543i i B i +⎛⎫⎪+⎪= ⎪+ ⎪⎪⎝⎭再求出它们的乘积矩阵C ,并将C 矩阵的右下角23⨯子矩阵赋给D 矩阵。

赋值完成后,调用相应的命令查看MATLAB 工作空间的占有情况。

解:A=[1 2 3 4;2 3 5 7;1 3 5 7 ;3 2 3 9 ;1 8 9 4;]B=[1+4i 4 3 6 7 7;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;] B=[1+4i 4 3 6 7 7;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;] C=A*B D=C(4:5,4:6);2、设矩阵16231351110897612414152⎛⎫⎪⎪ ⎪ ⎪⎝⎭,求A ,1A -,3A ,12A A -+,1'3A A --,并求矩阵A 的特征值和特征向量。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

考察下列单输入单输出系统:
x’’+3x’+2x=u
x(0)=1 ,x’(0)=0
输入:u=u(t), 输出:x=x(t), 0<t<10
编程做以下四件事:
1、u=1时,用符号法求输出(dslove)
2、u=1时,用数值法求输出(ode45)
3、用传递函数求输出阶跃响应(tf,step)
4、用Simulink法求输出阶跃响应(mdl)
给出源程序代码或模块,以及运行结果分析。

题目解答如下:
1、运行代码:x=dsolve('D2x+3*Dx+2*x=1','x(0)=1,Dx(0)=0','t')
>> ezplot(x,[0,10])
输出结果:
x =
exp(-t) - exp(-2*t)/2 + 1/2
2、在脚本文件中建立函数:function xdot=DyDt(t,x)
mu=3;
xdot=[x(2);1-2*x(1)-mu*x(2)];
在命令窗口中运行指令:
tspan=[0,10];
x0=[1,0];
[tt,xx]=ode45(@DyDt,tspan,x0);
plot(tt,xx(:,1));
xlabel('t'),title('x(t)')
由于函数ode45采用的数值积分法,故tt,xx的结果为一系列插值点,结果不易处理,在此只输出图形。

3、
在命令窗口中运行指令:
clear all;
G1=tf([1],[1,3,2])
step(G1,'m*');
legend('*G1')
输出结果:
G1 =
1
-------------
s^2 + 3 s + 2
4、建立siso.mdl文件,在文件中建立模型如下:。

相关文档
最新文档