second order system analysis 自动控制原理 二阶系统的matlab仿真分析

本文主要是介绍second order system analysis 自动控制原理 二阶系统的matlab仿真分析,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

二阶系统的matlab仿真分析





二阶系统的matlab仿真分析如上图。

根据二阶函数对阶跃函数的响应函数,我们对参数epsilon进行分析讨论


由于临界阻尼和无阻尼的情况在现实生活中比较难出现,二阶方程的根几乎不可能恰好,实部为0,或者两个实部相同且虚部为0. 于是,并为对以上两种较特殊的情况进行讨论。
选择欠阻尼和过阻尼两种情况进行分析讨论。

可以看出,当epsilon比较小的时候,响应时间短,且伴随有明显的超调。
随着epsilon的增大,超调明显降低,epsilon在0.7(恰巧工程上的最佳阻尼系数是0.7!)之后就没有超调了。随着epsilon的增大响应时间变得越来越长。



Wn = (1/T);
二阶方程的根
实部:X = -epsilon_0*Wn
虚部:Y = j*(Wn ).*sqrt(1-(epsilon_0).^2)

下图是上面十条曲线对应的根分布



可以看出图中所有的点均位于Y轴左侧,说明最终响应都将收敛
虚部为0的点很好的对应了过阻尼状态(同样颜色的是一对实根),图中关于X轴对称的点,系统处于欠阻尼状态

下图是上升时间和阻尼系数epsilon之间的关系(欠阻尼状态下)



可以看出,随着阻尼系数的增大,上升时间变长!

matlab和本文相关代码:

%%*************************************************************
% code writer: EOF
% code date:2014.03.18
% e-mail: jasonleaster@gmail.com
% code purpose : 
%           I just want to share with someone who is interesting
% in adaptive control. This code is to help people to understand
% second order system.
%%**************************************************************
clear all
clc
syms s f t m;
K01 = 1;
K02 = 1;
K0 = (K01*K02)./(1+K01*K02);
hold on;
figure(1);
T0 = 1;
for epsilon_0 = 0.1:0.2:2T = T0./(1+K0);
%     epsilon_0 = 0.5*(1/(K01*K02*T0));epsilon = epsilon_0./(1+K0);K = K0/(1+K0);f = (K./((T.^2).*(s.^2)+2*epsilon_0.*T.*s+K0)).*(1./s);m = ilaplace(f);ezplot(m,[0,120]);axis([0 60 0 1.2]);
end
legend('0.1','0.3','0.5','0.7','0.9','1.1','1.3','1.5','1.7','1.9');
hold off;
figure(2);
hold on;
for epsilon_0 = 0.1:0.2:2if epsilon_0 <1plot(-epsilon_0.*(1./T),(1./T).*sqrt(1-(epsilon_0).^2),'*');plot(-epsilon_0.*(1./T),-(1./T).*sqrt(1-(epsilon_0).^2),'*');elseif abs(epsilon_0-1.1) < 0.1  plot(-epsilon_0.*(1./T)+(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','r');plot(-epsilon_0.*(1./T)-(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','r');elseif abs(epsilon_0-1.3) < 0.1 plot(-epsilon_0.*(1./T)+(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','g');plot(-epsilon_0.*(1./T)-(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','g');elseif abs(epsilon_0-1.5) < 0.1 plot(-epsilon_0.*(1./T)+(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','b');plot(-epsilon_0.*(1./T)-(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','b');elseif abs(epsilon_0-1.7) < 0.1 plot(-epsilon_0.*(1./T)+(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','y');plot(-epsilon_0.*(1./T)-(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','y');elseif abs(epsilon_0-1.9) < 0.1 plot(-epsilon_0.*(1./T)+(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','k');plot(-epsilon_0.*(1./T)-(1./T).*sqrt((epsilon_0).^2-1),0,'*','Color','k');end
end
axis([-6 1 -3 3]);
hold off;
figure(3);
hold on;
Wn = (1./T);
for epsilon_0 = 0.1:0.2:2 Wd = Wn.*(1-(epsilon_0).^2);if epsilon_0 < 1belta = acos(epsilon_0);tr = (pi-belta)./Wd;plot(epsilon_0*10,tr,'*');end
end
legend('epsilon = 0.1','epsilon = 0.3','epsilon = 0.5','epsilon = 0.7','epsilon = 0.9');


这篇关于second order system analysis 自动控制原理 二阶系统的matlab仿真分析的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1035827

相关文章

Java程序进程起来了但是不打印日志的原因分析

《Java程序进程起来了但是不打印日志的原因分析》:本文主要介绍Java程序进程起来了但是不打印日志的原因分析,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录Java程序进程起来了但是不打印日志的原因1、日志配置问题2、日志文件权限问题3、日志文件路径问题4、程序

Java字符串操作技巧之语法、示例与应用场景分析

《Java字符串操作技巧之语法、示例与应用场景分析》在Java算法题和日常开发中,字符串处理是必备的核心技能,本文全面梳理Java中字符串的常用操作语法,结合代码示例、应用场景和避坑指南,可快速掌握字... 目录引言1. 基础操作1.1 创建字符串1.2 获取长度1.3 访问字符2. 字符串处理2.1 子字

Spring Boot循环依赖原理、解决方案与最佳实践(全解析)

《SpringBoot循环依赖原理、解决方案与最佳实践(全解析)》循环依赖指两个或多个Bean相互直接或间接引用,形成闭环依赖关系,:本文主要介绍SpringBoot循环依赖原理、解决方案与最... 目录一、循环依赖的本质与危害1.1 什么是循环依赖?1.2 核心危害二、Spring的三级缓存机制2.1 三

C#中async await异步关键字用法和异步的底层原理全解析

《C#中asyncawait异步关键字用法和异步的底层原理全解析》:本文主要介绍C#中asyncawait异步关键字用法和异步的底层原理全解析,本文给大家介绍的非常详细,对大家的学习或工作具有一... 目录C#异步编程一、异步编程基础二、异步方法的工作原理三、代码示例四、编译后的底层实现五、总结C#异步编程

Python 迭代器和生成器概念及场景分析

《Python迭代器和生成器概念及场景分析》yield是Python中实现惰性计算和协程的核心工具,结合send()、throw()、close()等方法,能够构建高效、灵活的数据流和控制流模型,这... 目录迭代器的介绍自定义迭代器省略的迭代器生产器的介绍yield的普通用法yield的高级用法yidle

Go 语言中的select语句详解及工作原理

《Go语言中的select语句详解及工作原理》在Go语言中,select语句是用于处理多个通道(channel)操作的一种控制结构,它类似于switch语句,本文给大家介绍Go语言中的select语... 目录Go 语言中的 select 是做什么的基本功能语法工作原理示例示例 1:监听多个通道示例 2:带

鸿蒙中@State的原理使用详解(HarmonyOS 5)

《鸿蒙中@State的原理使用详解(HarmonyOS5)》@State是HarmonyOSArkTS框架中用于管理组件状态的核心装饰器,其核心作用是实现数据驱动UI的响应式编程模式,本文给大家介绍... 目录一、@State在鸿蒙中是做什么的?二、@Spythontate的基本原理1. 依赖关系的收集2.

利用Python快速搭建Markdown笔记发布系统

《利用Python快速搭建Markdown笔记发布系统》这篇文章主要为大家详细介绍了使用Python生态的成熟工具,在30分钟内搭建一个支持Markdown渲染、分类标签、全文搜索的私有化知识发布系统... 目录引言:为什么要自建知识博客一、技术选型:极简主义开发栈二、系统架构设计三、核心代码实现(分步解析

C++ Sort函数使用场景分析

《C++Sort函数使用场景分析》sort函数是algorithm库下的一个函数,sort函数是不稳定的,即大小相同的元素在排序后相对顺序可能发生改变,如果某些场景需要保持相同元素间的相对顺序,可使... 目录C++ Sort函数详解一、sort函数调用的两种方式二、sort函数使用场景三、sort函数排序

kotlin中const 和val的区别及使用场景分析

《kotlin中const和val的区别及使用场景分析》在Kotlin中,const和val都是用来声明常量的,但它们的使用场景和功能有所不同,下面给大家介绍kotlin中const和val的区别,... 目录kotlin中const 和val的区别1. val:2. const:二 代码示例1 Java