控制系统数字仿真与CAD-第四次实验-附完整代码

我的微信公众号名称：深度学习与先进智能决策
微信公众号ID：MultiAgent1024
公众号介绍：主要研究强化学习、计算机视觉、深度学习、机器学习等相关内容，分享学习过程中的学习笔记和心得！期待您的关注，欢迎一起学习交流进步！

完整代码。

exc4.m

clear;
clear all;
clear close;
sa=2;sb=20;l=sb-sa;
e=1e-5;
x1=sa+(1-0.618)*l;
x2=sa+0.618*l;
q1=obj_asr(x1);
q2=obj_asr(x2);
count=0;
x_dot=[x1,x2];
while((sb-sa)/l>e)if(q1>q2)sa=x1;x1=x2;q1=q2;x2=sa+0.618*(sb-sa);x_dot=[x_dot,x2];q2=obj_asr(x2);elsesb=x2;x2=x1;q2=q1;x1=sa+(1-0.618)*(sb- sa);x_dot=[x_dot,x1];q1=obj_asr(x1);end
end
Kn_opt=(x1+x2)/2;
q_min=obj_asr(Kn_opt);
opt_asr_plot(Kn_opt);

fun4.m

function [result]=fun4(Kn)
tn=0.03;
Ki=0.269;
ti=0.067;
Ks=76;
Ts=0.00167;
R=6.58;
Tl=0.018;
Tm=0.25;
Ce=0.031;
Alpha=0.00337;
Beta=0.4;
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;
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 ]';
block_num=5;
h=0.001;
t_end=0.5;
t=0:h:t_end;
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) %若为积分比例，修正fai dFID(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
end
Y0=[0 0 0 0 0]';
n=length(t);
Y=Y0;
X=zeros(block_num,1);
result=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;X=FI'.*X+FIM'.*Uk+FIJ'.*Udot;Y=FIC'.*X+FID'.*Uf;Y(1)=satur(Y(1),c);result=[result,Y];
endend

obj_asr.m

function q=obj_asr(Kn)
result = fun4(Kn);
y_inf = 296.7350;
q=100*(max(result(5,:))-y_inf)/y_inf;
end

opt_asr_plot.m

function q=opt_asr_plot(Kn)
result = fun4(Kn);
h=0.001;
t_end=0.5;
t=0:h:t_end;
figure
subplot(4,1,1),plot(t,result(4,:)),grid,title('Id')
subplot(4,1,2),plot(t,result(3,:)),grid,title('Ud')
subplot(4,1,3),plot(t,result(2,:)),grid,title('ACR')
subplot(4,1,4),plot(t,result(1,:)),grid,title('ASR')
figure;
plot(t,result(5,:),'r');grid,title('n')
end

satur.m

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

控制系统数字仿真与CAD-第四次实验-附完整代码相关推荐

html页面数字滚动,Vue.js大屏可视化数据数字滚动翻转跳转效果（通俗易懂附完整代码）...
原文:Vue.js大屏可视化数据数字滚动翻转跳转效果(通俗易懂附完整代码) 大屏数字滚动翻转效果来源于最近工作中element后台管理页面一张大屏的UI图,该UI图上有一个模块需要有数字往上翻动的效果 ...
【控制系统数字仿真与CAD——实验报告】实验四：黄金分割法最优化PI调节器参数（文末附完整代码 + 实验结果）
一.实验目的掌握控制系统计算机辅助设计的原理与方法: 掌握黄金分割法的基本原理: 掌握黄金分割法在参数优化中的应用: 掌握MATLAB的基本编程,实现闭环系统的数字仿真,并增加参数优化过程: 掌握M ...
【控制系统数字仿真与CAD——实验报告】实验三：离散相似法数字仿真（文末附完整代码 + 实验结果）
一.实验目的 1. 了解离散相似法的基本原理 2. 掌握离散相似法仿真的基本过程 3. 应用离散相似法仿真非线性系统 4. MATLAB实现离散相似法的非线性系统仿真 5. 掌握SIMULINK仿真方 ...
[渝粤教育] 中国地质大学控制系统数字仿真复习题
控制系统数字仿真模拟题一.单选题 1.运行下列命令后A1=[123;456;789];A2=A1';A3=cat(1A1A2)系统输出结果为________. A. B. C. D. 2.可将两个状 ...
matlab zp2,matlab用于控制系统数字仿真
matlab用于控制系统数字仿真传递函数 n >>num=[6 12 6 10];den=[1 2 3 1 1]; n >>printsys(num,den) n 执行结果为 ...
【仿真】Carla之收集数据快速教程 (附完整代码)
收集过程可视化展示,随后进入正文: 参考与前言看到仿真群对这类任务下(用carla收集数据然后再做训练等) 需求量大,顺手马上写一个好了,首先收集数据需要考虑清楚: 收集什么数据,需要什么样的数据格 ...
基于 MATLAB 、Simulink 仿真和 TrueTime 工具箱的网络控制系统的仿真方法附完整代码
通过动手实践了解网络控制系统的架构及运行控制原理熟悉使用基于 MATLAB 平台.Simulink 仿真工具和 TrueTime 工具箱的网络控制系统的仿真方法考察不同丢包率下网络控制系统状态.输 ...
扶稳！四大步“上手”超参数调优教程，就等你出马了 | 附完整代码
作者 | Matthew Stewart 译者 | Monanfei 责编 | Jane 出品 | AI科技大本营(ID: rgznai100) [导读]在本文中,我们将为大家介绍如何对神经网络的超参 ...
【雷达】线性调频脉冲雷达仿真实验附matlab代码
1 内容介绍在MATLAB上实现线性调频信号(LMCW)的脉冲压缩(PC).动目标显示 (MTI)和动目标检测(MTD).解速度度模糊以及恒虚警率检测(CFAR)等信号处理算法.(本实验原内容为D ...
OpenCV SVM支持向量机和KNearest数字识别的实例(附完整代码)
OpenCV SVM支持向量机和KNearest数字识别的实例 OpenCV SVM支持向量机和KNearest数字识别的实例 OpenCV SVM支持向量机和KNearest数字识别的实例 #inc ...

控制系统数字仿真与CAD-第四次实验-附完整代码

控制系统数字仿真与CAD-第四次实验-附完整代码相关推荐

最新文章

热门文章