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

完整代码。

实验内容

1. 分别用欧拉法和四阶龙格库塔法求解如下：

（1）R=4K，C=1F，L=1H，仿真时间：0.03s

（2）R=280，C=8F，L=2H，仿真时间：0.1s

两组参数下电路的单位阶跃响应。（初始条件均为零）。

熟悉并理解掌握MATLAB模型转换以及参数提取指令（针对上述系统写出它们的运行结果）

G=tf (num, den)

G1=tf（G） G=zpk (z, p, k)

G1=zpk (G) G=ss (A, B, C, D)

G3=ss (G)

[num, den] = tfdata (G)

[Z, P,K] = zpkdata (G)

[A, B, C, D] = ssdata (G)

1.参数提取：

% -----参数提取----------------
clear;
clc;
L = 2;
C = 0.000008;
R = 280;
num = [1];
den = [L*C, R*C, 1];
[z,p,k] = tf2zp(num,den)
[A,B,C,D] = tf2ss(num,den)

2.欧拉法求解（1）R=4K，C=1F，L=1H，仿真时间：0.03s ：

function [z]=fun(X)
R =4000;
C =0.000001;
L = 1;
A = [0, 1/C;-1/L,-R/L];
B = [0;1/L];
z = A*X+B;
end

clear;
clc;
% ---------------欧拉法-------------
h = 0.0001;
y = [0;0];
for i = 0:1:300t = i*h;x(i+1) = t;k1 = fun(y);k2 = fun(y+h*k1);y = y + (k1 +k2)*h/2;output1(i+1) = y(1,1);
end
plot(x,output1)
%------------四阶-龙格-库塔-----------
% h = 0.0001;
% y = [0;0];
% for i = 0:1:1000
%     t = i*h;
%     x(i+1) = t;
%     k1 = fun(y);
%     k2 = fun(y+h*k1/2);
%     k3 = fun(y+h*k2/2);
%     k4 = fun(y+h*k3);
%     y = y + (k1 + 2*k2 +2*k3 + k4)*h/6;
%     output1(i+1) = y(1,1);
% end
% plot(x,output1)

运行结果如下：

四阶龙格库塔（2）R=280，C=8F，L=2H，仿真时间：0.1s：

function [z]=fun(X)
R =280;
C =0.000008;
L = 2;
A = [0, 1/C;-1/L,-R/L];
B = [0;1/L];
z = A*X+B;
end

clear;
clc;
% ---------------欧拉法-------------
% h = 0.0001;
% y = [0;0];
% for i = 0:1:300
%     t = i*h;
%     x(i+1) = t;
%     k1 = fun(y);
%     k2 = fun(y+h*k1);
%     y = y + (k1 +k2)*h/2;
%     output1(i+1) = y(1,1);
% end
% plot(x,output1)
%------------四阶-龙格-库塔-----------
h = 0.0001;
y = [0;0];
for i = 0:1:1000t = i*h;x(i+1) = t;k1 = fun(y);k2 = fun(y+h*k1/2);k3 = fun(y+h*k2/2);k4 = fun(y+h*k3);y = y + (k1 + 2*k2 +2*k3 + k4)*h/6;output1(i+1) = y(1,1);
end
plot(x,output1)

运行结果如下：

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

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

【控制系统数字仿真与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 执行结果为 ...
JavaScript实现切换数字的符号switchSign算法（附完整源码）
JavaScript实现切换数字的符号switchSign算法(附完整源码) switchSign.js完整源代码 switchSign.js完整源代码 /*** Switch the sign of ...
【仿真】Carla之收集数据快速教程 (附完整代码)
收集过程可视化展示,随后进入正文: 参考与前言看到仿真群对这类任务下(用carla收集数据然后再做训练等) 需求量大,顺手马上写一个好了,首先收集数据需要考虑清楚: 收集什么数据,需要什么样的数据格 ...
html页面数字滚动,Vue.js大屏可视化数据数字滚动翻转跳转效果（通俗易懂附完整代码）...
原文:Vue.js大屏可视化数据数字滚动翻转跳转效果(通俗易懂附完整代码) 大屏数字滚动翻转效果来源于最近工作中element后台管理页面一张大屏的UI图,该UI图上有一个模块需要有数字往上翻动的效果 ...
Three.js实例详解___旋转的精灵女孩(附完整代码和资源)（一）
Three.js实例详解___旋转的精灵女孩(附完整代码和资源)(一) 本文目录: 一.[旋转的精灵女孩]案例运行效果二.Three.js简介三.Three.js代码正常运行显示条件 (1)不载入 ...
对“科大讯飞2021丨广告点击率预估挑战赛 Top1方案(附完整代码)_Jack_Yang的博客-CSDN博客”的补充。
这篇文章的初衷是针对科大讯飞2021丨广告点击率预估挑战赛 Top1方案(附完整代码)_Jack_Yang的博客-CSDN博客进行补充. 博客的信息量很少,对任务背景的介绍也不太对,说实话令人费解.我 ...
Three.js实例详解___旋转的精灵女孩(附完整代码和资源)（二）
Three.js实例详解___旋转的精灵女孩(附完整代码和资源)(二) 本篇目录: 五.实例中所使用的代码语法详细解释 (1).构建一个三维空间场景 (2).选择一个透视投影相机作为观察点 (a).创 ...

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

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

最新文章

热门文章