一、实验目的

1. 了解离散相似法的基本原理
　　2. 掌握离散相似法仿真的基本过程
　　3. 应用离散相似法仿真非线性系统
　　4. MATLAB实现离散相似法的非线性系统仿真
　　5. 掌握SIMULINK仿真方法，应用于非线性系统的仿真，并对实验结果进行分析比较

二、实验原理

　　在ASR的输出增加限幅装置（饱和非线性，饱和界为c=8 ）。 Ce=0.031，其它参数不变。输入为单位阶跃，用离散相似法求系统各环节的输出。

　　要求：采用零阶保持器和一阶保持器离散化系统，分别完成本实验。

1、各环节的参数：

　　由5个典型环节组成：

　　A=[0 0 1 1 0];

　　B=[tn ti Ts Tl Tm*Ce];
　　C=[Kn Ki Ks 1/R R];
　　D=[Kn*tn Ki*ti 0 0 0];

　　还有一饱和非线性环节：c=8;

2、各环节的离散化系数矩阵

3、各环节的输入作用

（1）u(n)可通过联接矩阵直接求得：

u(n)=[u1(n), u2(n),…,un(n)]为各环节的输入量， n为环节数。 Y(n)=[Y1(n), Y2(n),…,Yn(n)]为各环节的输出量， r为外中参考输入量。

（2）U(n)由近似表达式求得：
　　　　　　

（3）u(n+1)用折线法近似求得：

4、状态和输出计算

（1）一阶保持器
　　　　　　X=FI'.*X+FIM'.*Uk+FIJ'.*Udot;
　　　　　　Y=FIC'.*X+FID'.*Uf;

（2）零阶保持器
　　　　　　X=FI'.*X+FIM'.*Uk;
　　　　　　Y=FIC'.*X+FID'.*Uf;

5、饱和非线性环节

看作环节1（ASR）的一部分。建立satur.m文件：

function [uo]=satur(ui,c)if (abs(ui)<=c)uo=ui;elseif ( ui > c )uo = c;elseuo=-c;end
end

三、实验过程

1、新建脚本文件，命名为satur.m

function [uo]=satur(ui,c)if (abs(ui)<=c)uo=ui;elseif ( ui > c )uo = c;elseuo=-c;end
end

2、新建脚本文件，命名为test3.m

　　完整代码：

clc;
clear;
% ******  各环节参数  ****** %
Kn=26.7;
tn=0.03;
Ki=0.269;
ti=0.067;
Ks=76;
Ts=0.00167;
R=6.58;
T1=0.018;
Tm=0.25;
Ce=0.031;
Alpha=0.00337;
Beta=0.4;
A=[0 0 1 1 0];
B=[tn ti Ts T1 Tm*Ce];
C=[Kn Ki Ks 1/R R];
D=[Kn*tn Ki*ti 0 0 0];
c=8;
r=1;
W=[0 0 0 0 -Alpha;1 0 0 -Beta 0;0 1 0 0 0;0 0 1 0 -Ce;0 0 0 1 0];
W0=[1 0 0 0 0]';
h=0.001;
t_end=0.5;
t=0:h:t_end;
n=length(t);
% ******  各环节离散化系数  ****** %
block_num=5;
for k=1:block_numif(A(k)==0)FI(k)=1;FIM(k)=h*C(k)/B(k);FIJ(k)=h*h*C(k)/B(k)/2;FIC(k)=1;FID(k)=0;if(D(k)~=0)FID(k)=D(k)/B(k);endelseFI(k)=exp(-h*A(k)/B(k));FIM(k)=(1-FI(k))*C(k)/A(k);FIJ(k)=h*C(k)/A(k)-FIM(k)*B(k)/A(k);FIC(k)=1;FID(k)=0;if(D(k)~=0)FIC(k)=C(k)/D(k)-A(k)/B(k);FID(k)=D(k)/B(k);endend
endY0=[0 0 0 0 0]';
Y=Y0;
X=zeros(block_num,1);
result1=Y;
Uk=zeros(block_num,1);
Ub=Uk;for m=1:(n-1)Ub=Uk;Uk=W*Y+W0*r;Uf=2*Uk-Ub;Udot=(Uk-Ub)/h;%******  零阶保持器  ******%X=FI'.*X+FIM'.*Uk;Y=FIC'.*X+FID'.*Uf;Y(1)=satur(Y(1),c);result1=[result1,Y];
endY0=[0 0 0 0 0]';
Y=Y0;
X=zeros(block_num,1);
result2=Y;
Uk=zeros(block_num,1);
Ub=Uk;
for m=1:(n-1)Ub=Uk;Uk=W*Y+W0*r;Uf=2*Uk-Ub;Udot=(Uk-Ub)/h;%******  一阶保持器  ******%X=FI'.*X+FIM'.*Uk + FIJ'.*Udot;Y=FIC'.*X + FID'.*Uf;Y(1)=satur(Y(1),c);result2=[result2,Y];
endplot(t,result1(5,:),'-.',t,result2(5,:),'--',t,ScopeData.signals.values,'k');
legend('零阶保持器','一阶保持器','Simulink');

3、在Simulink中绘制仿真图

注意：Simulink中的变量名和工作区变量关联方法请点击：这里

四、实验结果

【控制系统数字仿真与CAD——实验报告】实验三：离散相似法数字仿真（文末附完整代码 + 实验结果）相关推荐

【控制系统数字仿真与CAD——实验报告】实验四：黄金分割法最优化PI调节器参数（文末附完整代码 + 实验结果）
一.实验目的掌握控制系统计算机辅助设计的原理与方法: 掌握黄金分割法的基本原理: 掌握黄金分割法在参数优化中的应用: 掌握MATLAB的基本编程,实现闭环系统的数字仿真,并增加参数优化过程: 掌握M ...
分布式数据库九大发展趋势｜文末附完整报告下载
随着数据成为第五大生产要素,新兴软硬件技术快速发展,各行各业的数据呈指数式增长,传统集中式数据库面对海量数据存储.高并发等场景的局限性愈发凸显.分布式数据库凭借高可用.高可扩展性.高性价比等优势,已经 ...
数字2DPSK频带传输系统的设计（文内附完整代码）
使用MATLAB实现二进制数字频带系统2DPSK,完成基带信号产生.差分变换.调制解调.抽样判决.码反变换等功能. 2DPSK -- 二进制差分相移键控码率:24Kbps 信道:高斯信道解调采用相 ...
计算机控制系统AD实验报告,实验一AD转换实验报告
实验一AD转换实验报告 (6页) 本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦! 9.9 积分计算机控制系统实验报告实验一实验名称:A/D转换实验二零一六年三 ...
计算机控制系统AD实验报告,实验一AD转换实验报告.docx
文档介绍: 实验一AD转换实验报告.docx计算机控制系统实验报告实验一实验名称:A/D转换实验二零一六年三月实验_A/D转换实验.实验要求利用ADC做A/D转换器,将模拟量转换成二进制数字量,用延时 ...
计算机仿真技术应用报告,计算机仿真技术实验报告-实验二
计算机仿真技术实验报告-实验二 (7页) 本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦! 9.9 积分 <仿真技术与应用>实验报告计算机仿真技术实验 ...
c语言cad实验报告,1.电子线路CAD实验报告实验
电子线路CAD实验报告实验电子线路CAD实验报告实验序号:实验一实验名称:Altium Designer 基本操作实验日期:2014.2.18 专业班级:电子信息工程姓名:陈蕾学号:122 ...
dsp 实验c语言,dsp实验报告实验1 ccs入门实验2（c语言的使用）
1. 1 DSP 实验报告实验报告院 ( 系 ) 名称自动化科学与电气工程学院学生姓名学号任课老师吴冠 2014 年 6 月 DSP 实验报告 1 ...
实验c语言程序数据类型,C语言程序设计实验报告——实验2--数据类型、运算符和表达式.doc...
C语言程序设计实验报告--实验2--数据类型.运算符和表达式.doc 下载提示(请认真阅读)1.请仔细阅读文档,确保文档完整性,对于不预览.不比对内容而直接下载带来的问题本站不予受理. 2.下载的文档 ...

【控制系统数字仿真与CAD——实验报告】实验三：离散相似法数字仿真（文末附完整代码 + 实验结果）