- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
4-2设典型闭环结构控制系统如图4-47所示,当阶跃输入幅值 时,用sp4_1.m求取输出的响应。
解:用sp4_1.m求解过程如下:
在MATLAB语言环境下,输入以下命令语句
a=[0.016 0.864 3.27 3.42 1];
b=[30 25];
X0=[0 0 0 0]; %系统状态向量初值为零
V=2; %反馈系数
n=4;
T0=0;Tf=10;
h=0.01;R=20 ; %仿真步长h=0.01,阶跃输入幅值
sp4_1 %调用sp4_1.m函数
plot(t,y)
运行结果为:
附:sp4_1.m函数为
b=b/a(1);a=a/a(1);A=a(2:n+1);
A=[rot90(rot90(eye(n-1,n)));-fliplr(A)];
B=[zeros(1,n-1),1];
m1=length(b);
C=[fliplr(b),zeros(1,n-m1)];
Ab=A-B*C*V;
X=X0;
y=0;t=T0;
N=round((Tf-T0)/h);
for i=1:N
K1=Ab*X+B*R;
K2=Ab*(X+h*K1/2)+B*R;
K3=Ab*(X+h*K2/2)+B*R;
K4=Ab*(X+h*K3)+B*R;
X=X+h*(K1+2*K2+2*K3+K4)/6;
y=[y,C*X];
t=[t,t(i)+h];
end
4-4系统结构图如图4-48,写出该系统的联结矩阵和,并写出联结矩阵非零元素阵。
解:根据图4-48中,拓扑连结关系,可写出每个环节输入受哪些环节输出的影响,
现列如入下:
把环节之间的关系和环节与参考输入的关系分别用矩阵表示出来,
即=,=,
4-6若系统为图4-5b 双输入-双输出结构,试写出该系统的联接矩阵,,说明应注意什么?
解:根据图4-5b中,拓扑连结关系,可列写如下关系式:
转换成矩阵形式为
所以联接矩阵=,=
此时应注意输入联接矩阵变为型。
4-8求图4-49非线性系统的输出响应y(t),并与无非线性环节情况进行比较。
解:(1)不考虑非线性环节影响时,求解过程如下:
先将环节编号标入图中。
2) 在MATLAB命令窗口下,按编号依次将环节参数输入P阵;
P=[0.1 1 0.5 1;0 1 20 0;2 1 1 0;10 1 1 0];
3) 按各环节相对位置和联接关系,有联接矩阵如下:
, ,所以非零元素矩阵
WIJ=[1 0 1;1 4 -1;2 1 1;3 2 1;4 3 1] ;
4)由于不考虑非线性影响,则非线性标志向量和参数向量均应赋零值;
Z=[0 0 0 0];S=[0 0 0 0];
5)输入运行参数:开环截至频率约为1,故计算步长h取经验公式值,即
,取h=0.01;每0.25秒输出一点。故取=25。
h=0.01;
L1=25;
n=4;
T0=0
Tf=20;
nout=4;
Y0=10;
sp4_4;
plot(t,y,r)
hold on
运行结果如图中红色实线所示。
(2)考虑非线性环节N影响时,只需将非线性标志向量Z和参数向量S的相应分量正确输入即可。
在MATLAB命令窗口中输入下列语句:
Z=[4 0 0 0];S=[5 0 0 0]; %第一个线性环节后有饱和非线性,参数值为5。
sp4_4;
plot(t,y,--)
运行结果如图中蓝色虚线所示。
从图中可以清楚的地看出,饱和非线性环节对线性系统输出响应的影响。
附:sp4_4函数为:
A=P(:,1);B=P(:,2);
C=P(:,3);D=P(:,4);
m=length(WIJ(:,1));
W0=zeros(n,1);W=zeros(n,n);
for k=1:m
if (WIJ(k,2)==0); W0(WIJ(k,1))=WIJ(k,3);
else W(WIJ(k,1),WIJ(k,2))=WIJ(k,3);
end;
end;
for i=1:n
if(A(i)==0);
FI(i)=1;
FIM(i)=h*C(i)/B(i);
FIJ(i)=h*h*(C(i)/B(i))/2;
FIC(i)=1;FID(i)=0;
您可能关注的文档
- 冷冲压工艺与模具设计 教学课件 作者 匡余华 习 题 及答案.doc
- 客户服务导论与呼 叫中心实务 第3版 教学课件 作者 978 7 302 21568 4 FLA.doc
- 客户服务导论与呼 叫中心实务 第3版 教学课件 作者 978 7 302 21568 4 FLB.doc
- 客户服务导论与呼 叫中心实务 第3版 教学课件 作者 978 7 302 21568 4 FLC.doc
- 客户服务导论与呼 叫中心实务 第3版 教学课件 作者 978 7 302 21568 4 FLD.doc
- 客户服务导论与呼 叫中心实务 第3版 教学课件 作者 978 7 302 21568 4 FLE.doc
- 客户服务导论与呼 叫中心实务 第3版 教学课件 作者 978 7 302 21568 4 FLF.doc
- 客户服务导论与呼 叫中心实务 第3版 教学课件 作者 978 7 302 21568 4 FLG.doc
- 客户服务导论与呼 叫中心实务 第3版 教学课件 作者 978 7 302 21568 4 ML.doc
- 客户服务实务 教案 第八周 第二次课.doc
- 控制系统数字仿真与CAD 第3版 教学课件 作者 张晓华 主编 第五章习题答案.doc
- 控制系统数字仿真与CAD 第3版 教学课件 作者 张晓华 主编 第一章习题答案.doc
- 控制系统数字仿真与CAD 第3版 教学课件 作者 张晓华 主编 吊车防摆技术的研究.pdf
- 控制系统数字仿真与CAD 第3版 教学课件 作者 张晓华 主编 非连续系统的Simulink仿真方法研究.doc
- 控制系统数字仿真与CAD 第3版 教学课件 作者 张晓华 主编 概率统计课程中的数学思想方法研究.pdf
- 控制系统数字仿真与CAD 第3版 教学课件 作者 张晓华 主编 回收塑料_木材纤维复合材料的工艺及性能.PDF
- 控制系统数字仿真与CAD 第3版 教学课件 作者 张晓华 主编 机器人运动控制仿真系统研究与实现.PDF
- 控制系统数字仿真与CAD 第3版 教学课件 作者 张晓华 主编 基于滑模方法的桥式吊车系统的抗摆控制.PDF
- 控制系统数字仿真与CAD 第3版 教学课件 作者 张晓华 主编 集成电路设计中的实物仿真技术.PDF
- 控制系统数字仿真与CAD 第3版 教学课件 作者 张晓华 主编 局内电梯调度问题与竞争算法.PDF
最近下载
- 中小学书法培训中级班课件04横折折撇-附教案.pptx VIP
- 一种治疗阿尔茨海默症的中药组合物及其制备方法与应用.pdf VIP
- 新标准大学英语综合教程4__课后答案解析.docx VIP
- 广东河源电厂×600MW机组工程主要危险源辨识评价与控制措施.doc
- 康明斯KTA50发动机零件图册.pdf VIP
- 均值不等式练习题.docx VIP
- ISO 19880-3-2018 气态氢 - 加气站-第3部分:阀门.pdf
- 通风管道技术规程T 141-2017.doc VIP
- 《系统工程 第6版》 课件 第3、4章 系统模型与模型化、 系统仿真及系统动力学方法.pptx
- 新版FMEA表单模板(DFMEA和PFMEA)案例.xlsx VIP
文档评论(0)