matlab与控制系统仿真期末复习文档格式.docx
《matlab与控制系统仿真期末复习文档格式.docx》由会员分享,可在线阅读,更多相关《matlab与控制系统仿真期末复习文档格式.docx(11页珍藏版)》请在冰豆网上搜索。
(2)而乘*,这个和我们在数学中学到的矩阵相乘是一个意义。
a*c
ans=
a*b
mtimes
Innermatrixdimensionsmustagree.
由于a的列不等于b的行,所以两者不可相乘。
Shift+enter:
一行多命令输入
…:
连接符可以把断开的命令连起来,但是…后的不在被识别,是字符要单独用‘隔开
;
:
表示不显示结果
A=5;
A=5,B=6
\矩阵左除
/矩阵右除
x=a/b即x=a*inv(b)a乘以b的逆矩阵。
x=a\b,即x=inv(a)*ba的逆矩阵乘以b矩阵。
AB是标量A/B等于B\A
AB是矩阵A/B不等于A\B
标量:
单个数据
a=2;
向量:
一行或一列数据
a=[3,5,7,8]或a=[3578]行向量,在方括号内的一行数据
%:
注释
closeall关闭图形窗口
clear清除工作空间
clc清屏
who:
who命令只列出当前工作空间的变量名
whos:
以列表的形式,给当前工作空间的变量名以及变量的其它信息(大小,字节,类等)
figure(n):
figure是建立图形的意思,系统自动从1,2,3,4...来建立图形,数字代表第几幅图形,figure
(1),figure
(2)就是第一第二副图的意思
saveacb.matAB–append
将AB变量保存成文件名为acb.mat的文件。
loadacd.matAB装载当前路径下的数据文件
B
冒号表达式A=0:
pi/180:
2*pi范围(0,2*pi)步长为pi/180
B=0:
10
B=
012345678910
linespace(0,2*pi,181)
x=linspace(a,b,n);
a为下限,b为上限,n为将a到b等分为(n-1)份即分为n个点
若没有n,则默认为将a到b等分为99份
ones(m,n):
ones(m,n):
生成m*n的全1矩阵
eye(m,n):
eye(m,n):
生成m*n的矩阵E,对角线元素为1,其他为0
rand(m,n):
生成m*n的随机矩阵函数rand函数用来生成[0,1]之间均匀分布的随机函数
round\fix\floor\ceil
Matlab中的取整函数fix,floor,ceil与round
1.fix向0取整
fix([-0.5-0.40.40.5])==>
[0000]
2.floor向-inf取整
floor([-0.5-0.40.40.5])==>
[-1-100]
3.ceil向+inf取整
ceil([-0.5-0.40.40.5])==>
[0011]
4.round四舍五入
round([-0.5-0.40.40.5])==>
[-1001]
频率法超前校正零点小于极点(零点加在极点左侧),中频段,作用是提高动态性能,ts
滞后校正零点大于极点(零点加在极点右侧),低频段,作用是提高稳态性能
常用simulink模块来源:
1
Continuous(连续模块)
Integrator
输入信号积分
Derivative
输入信号微分
State-Space
线性状态空间系统模型
Transfer-Fcn
线性传递函数模型
Zero-Pole
以零极点表示的传递函数模型
Memory
存储上一时刻的状态值
Transport
Delay
输入信号延时一个固定时间再输出
Variable
输入信号延时一个可变时间再输出
2
Discrete
(离散模块)
Discrete-time
离散时间积分
Filter
IIR与FIR滤波器
离散状态空间系统模型
离散传递函数模型
以零极点表示的离散传递函数模型
First-Order
Hold
一阶采样和保持器
Unit
一个采样周期的延时
3
Function&
Tables(函数和表格模块)
Fcn
用自己定义的函数(表达式)进行运算
MATLAB
利用MATLAB的现有函数进行运算
S-Function
调用自编的S函数的程序进行运算
Look-Up
Table
建立输入信号的查询表(线性峰值匹配)
(2-D):
建立两个输入信号的查询表(线性峰值匹配)
4
Math
Operations(数学运算模块)
Sum
加减运算
Product
乘运算
Dot
点乘运算
Gain
比例运算
Function
包括指数函数、对数函数、求平方、开根号等常用数学函数
Trigonometric
三角函数,包括正弦、余弦、正切等
MinMax
最值运算
Abs
取绝对值
Sign
符号函数
Logical
Operator
逻辑运算
Relational
关系运算
Complex
to
Magnitude-Angle
由复数输入转为幅值和相角输出
由幅值和相角输入合成复数输出
Real-Imag
由复数输入转为实部和虚部输出
由实部和虚部输入合成复数输出
5
Signals&
Systems
(信号和系统模块)
In1:
输入端
Out1:
输出端
Mux:
将多个单一输入转化为一个复合输出
Demux:
将一个复合输入转化为多个单一输出
Ground:
连接到没有连接到的输入端
Terminator:
连接到没有连接到得输出端
SubSystem:
建立新的封装(Mask)功能模块
6
Sinks
(接收器模块)
Scope:
示波器
XY
Graph:
显示二维图形
To
Workspace:
将输入写入MATLAB的工作空间
File(.mat):
将输出写入数据文件
7
Sources
(输入源模块)
Constant
常数信号
Clock:
时钟信号
From
来自MATLAB的工作空间
来自数据文件
Pulse
Generator:
脉冲发生器
Repeating
Sequence:
重复信号
Signal
信号发生器,可以产生正弦、方波、锯齿波及随意波
Sine
Wave
正弦波信号
Step:
阶跃波信号
根轨迹:
开环传递函数G=11.9/s(s+10)
直接利用rlocus()函数可绘制根轨迹
具体程序如下:
s=tf('
s'
);
%定义传递函数算子
G=11.9/(s*(s+10));
rlocus(G)
二、选择题(20分)
三、问答题(60分)
1.方程组的求解P39
AX=BB一般是列向量B=[0;
0;
1]X=inv(A)*B=A\B
2.矩阵的组合和分解、ones、eye
例一:
a=[12];
b=[34];
c=cat(1,a,b);
结果:
c=
1
2
3
4
例二:
a=[12;
34];
b=[1122;
3344];
d=cat(2,a,b);
11
22
33
44
d
=
2
4
3.图形的绘制plot,subplot,holdon/off,text,title第五章
a.plot(x,y,’r’)图形修饰title,text
b.图形保持(holdon)c.子图
4.符号表达式的操作symsyms简化,分解P45
a.定义symsymsb.求解c.积分,微分
5.超前滞后校正(自动控制原理课后)
6.系统稳定性判定feedbackrootif
练习
1.在MATLAB的命令行窗口中按上方向键可以恢复已执行或已输入的上一行命令。
2.将已计算变量T1和T2存入已有数据的数据文件daa.mat中,数据文件的原数据不丢失,在命令行中输入的指令为savedaa.datT1T2,从daa.mat文件中导入T1,T2的命令为:
loaddaa.matT1T2
3.如果一个表达式或命令过长,在MATLAB命令行中的一行书写不下,可以使用连接符…将断开的命令连起来。
4.MATLAB中变量名由字母/数字/下划线等组成,区分字母的大小写;
变量中不能包含有标点符号。
5.MATLAB的预定义变量inf表示无穷大,pi表示圆周率,ans为默认变量名,虚数单位用i或者j表示。
6.为了调整数据的输出格式为15位定点数,需在命令行窗口执行formatlong命令。
7.在公式后加符号;
可以使命令行不显示运算结果,在表达式或文字前加符%可以将其转换为注释。
8.在使用MATLAB时,要查找名称中有plot的函数用lookforplot命令查找;
要查看plot的功能,使用helpplot
9.A矩阵和B矩阵进行点乘运算,两个矩阵的行列数应相等,A和B进行矩阵运算,其乘法A(m行n列),B(n行p列)才符合。
10.对于除法运算,左除和右除的区别:
、有标量A和B,及矩阵C和D,除法运算A/B等于B\A,C/D不等于D\C
11.用冒号表达式建立从数值A到B,间隔为0.1的线性行向量C的命令为C=A:
0.1:
B。
用linspace产生1到10,间隔为1的向量的表达式为linspace(1,10,11)
12.查找出矩阵A大于4的元素的位置的matlab命令为find(A>
4)
13.已知矩阵A=[123;
246;
345];
,其逆矩阵计算命令为inv(A),算其特征矢量c和特征值d的命令为[c,d]=eig(a)
14.对于一次输入多行命令语句的情况,可以在每条命令输入完后,按组合键shift+enter后输入下一条命令,直至按enter键后一次性执行。
15.在MATLAB的工作空间中查看所有已存在变量的名称、sizes可以用以下whos命令。
16.在MATLAB的主窗口中,可以在历史命令面板中利用已执行的命令建立m文件。
17.Clc为清屏,clear为清除变量工作空间,clf为关闭图形窗口
18.ones(A)为全1矩阵,eye(m,n)为对角线为1的矩阵,rand(n)为随机矩阵
19.已知数组a,则运行mean(a)命令是计算a每列的平均值,det(a)为求矩阵a的行列式,inv(a)为其逆运算,a‘为其转置运算
20.已知x向量,其对数函数有以e为底的log(x),以10为底的log10(x),以2为底的log2(x),
21.在取整函数中:
ceil:
向正无穷方向舍入;
fix:
向零方向舍入;
floor:
向负无穷方向舍入;
round为四舍五入。
22.a=0:
4,b=1:
5,a和b可以做加法和减法,可以求转置,但不可以做矩阵的乘法和除法,点乘和点除可以。
23.利用conv函数计算三个多项式s1、s2和s3的乘积,则算式为conv(conv(s1,s2),s3)
24.已知方程
,通过符号方程求解x
symsabcx
s=a*x^2+b*x+c;
r=solve(s)
已知方程
,求解方程的根
r=roots(a)
对下式进行部分分式展开:
b=[32546];
[r,s,k]=residue(b,a)
求矩阵
的行列式值、逆
AD=det(A)
AI=inv(A)
计用符号方法求下列极限或导数