连续/离散的控制系统阶跃测试(包括MATLAB里的step()函数)

阶跃测试

只要是连续时间系统,无论是传递函数还是连续状态空间形式的模型,直接可以用**step()**做阶跃测试;但是对于离散系统而言,不能用step()函数,可以自行编写代码,如下。

1、离散系统:x(k+1)=Ax(k)+Bu(k)


离散系统阶跃响应不能用step()函数,可以自行编写

复制代码
close all;
clear all;
A=[1.1 2;0 0.95];B=[0;0.079];

total = 100;
x = zeros(2,total+1);
u = zeros(1,total);
u(1) = 1;
t = zeros(1,total);

for i=1:total
    
    t(i) = i-1;
    x(:,i+1) = A*x(:,i) + B*u(:,i);
end

figure(1);
subplot(3,1,1);
plot(t,x(1,1:total));
ylabel('x1')
hold on;
grid on;

subplot(3,1,2);
plot(t,x(2,1:total));
ylabel('x2')
hold on;
grid on;

subplot(3,1,3);
plot(t,u(1:total));
ylabel('u')
hold on;
grid on;

可见系统自身不稳定,因为离散系统的特征值有在单位圆外,如下

2、连续系统:x'=Ax+Bu或y=G(s)u

一个是连续状态空间模型,一个是传递函数模型

2.1 传递函数形式:

MATLAB里写成传函形式,并做阶跃测试**step()**函数

2.2 化成连续状态空间形式

再做阶跃

不管是传函,还是连续状态空间,step()测试结果相同,说明系统稳定,因为状态矩阵特征值都<0;

eig(A)

ans =

-5.4178 + 0.0000i

-0.2911 + 1.3270i

-0.2911 - 1.3270i

2.3 连续状态空间形式化成离散状态空间形式

(1)对传递函数进行离散

得到的是离散时间传递函数,z变换
(2)对连续状态空间模型进行离散

得到的是离散状态空间模型
(3)对离散状态空间模型做阶跃测试

离散系统阶跃响应不能用step()函数,可以自行编写,就是本文第1节,下面继续写一遍:

复制代码
A = [0.5307 -0.2069 -0.1864;0.1491 0.978 -0.02058;0.01647 0.1985  0.9986];B = [0.1491;0.01647;0.001152];C = [0 2.5 1.25];%得到的离散状态空间模型

total = 100;
x = zeros(3,total+1);
u = zeros(1,total);
u(1) = 1;
y = zeros(1,total);
t = zeros(1,total);

for i=1:total
    
    t(i) = i-1;
    x(:,i+1) = A*x(:,i) + B*u(:,i);
    y(:,i+1) = C*x(:,i+1);
end

figure(1);
subplot(3,1,1);
plot(t,x(1,1:total));
ylabel('x1')
hold on;
grid on;

subplot(3,1,2);
plot(t,x(2,1:total));
ylabel('x2')
hold on;
grid on;

subplot(3,1,3);
plot(t,x(3,1:total));
ylabel('x3')
hold on;
grid on;

figure(2)
subplot(2,1,1);
plot(t,u(1:total));
ylabel('u')
hold on;
grid on;

subplot(2,1,2);
plot(t,y(1:total));
ylabel('y')
hold on;
grid on;

图形如下:

化成离散系统后,仍然稳定。

相关推荐
侃侃_天下2 天前
最终的信号类
开发语言·c++·算法
echoarts2 天前
Rayon Rust中的数据并行库入门教程
开发语言·其他·算法·rust
Aomnitrix2 天前
知识管理新范式——cpolar+Wiki.js打造企业级分布式知识库
开发语言·javascript·分布式
每天回答3个问题2 天前
UE5C++编译遇到MSB3073
开发语言·c++·ue5
伍哥的传说2 天前
Vite Plugin PWA – 零配置构建现代渐进式Web应用
开发语言·前端·javascript·web app·pwa·service worker·workbox
小莞尔2 天前
【51单片机】【protues仿真】 基于51单片机八路抢答器系统
c语言·开发语言·单片机·嵌入式硬件·51单片机
我是菜鸟0713号2 天前
Qt 中 OPC UA 通讯实战
开发语言·qt
JCBP_2 天前
QT(4)
开发语言·汇编·c++·qt·算法
Brookty2 天前
【JavaEE】线程安全-内存可见性、指令全排序
java·开发语言·后端·java-ee·线程安全·内存可见性·指令重排序
百锦再2 天前
[特殊字符] Python在CentOS系统执行深度指南
开发语言·python·plotly·django·centos·virtualenv·pygame