控制系统数字仿真

matlab学习笔记3（控制系统仿真）

1、典型二阶系统

,绘制出

,

分别为0.1、0.2、0.3、···、1.0和2.0时的单位阶跃响应。

clear
clc
Wn=0.6;
kesai=[0.1:0.1:1,2];
hold on
for k=kesai
    num=Wn^2;
    den=[1 2*k*Wn Wn^2];
    step(num,den);
end
title('step response');           

2、对1中的典型二阶系统，绘制出当

，

时的单位阶跃响应。

clear
clc
W=2:2:12;
kesai=0.7;
hold on
for wn=W
    num=wn^2;
    den=[1 2*kesai*wn wn^2];
    step(num,den);
end
title('step response');
hold off           

3、求三阶系统

的单位阶跃响应。

clear
clc
num=[5 25 30];
den=[1 6 10 8];
step(num,den);
title('step response');           

4、求典型二阶系统

,当

clear
clc
Wn=6;
kesai=0.7;
num=Wn^2;
den=[1 2*kesai*Wn Wn^2];
step(num,den);
title('step response');           

5、有高阶系统$\left\lbrace \begin{array}{l}

\dot{x} =\left\lbrack \begin{array}{cccc}

-1.6 & -0.9 & 0 & 0\\

0.9 & 0 & 0 & 0\\

0.4 & 0.5 & -5.0 & -2.45\\

0 & 0 & 2.45 & 0

\end{array}\right\rbrack x+\left\lbrack \begin{array}{l}

1\

0\

\end{array}\right\rbrack u\

y=\left\lbrack \begin{array}{llll}

1 & 1 & 1 & 1

\end{array}\right\rbrack x

\end{array}\right.

x_0 ={\left\lbrack 1\;,1,\;1,-1\right\rbrack }^T$。

a0=[-1.6,-0.9,0,0;0.9,0,0,0;0.4,0.5,-5.0,-2.45;0,0,2.45,0];
b0=[1 0 1 0]';
c0=[1 1 1 1];
d0=[0];
t=0.5;   %采样周期
[a,b,c,d]=c2dm(a0,b0,c0,d0,t,'tustin');%连续系统离散化
subplot(2,2,1);
dstep(a,b,c,d);%离散系统的单位阶跃响应
title('discrete step response');
subplot(2,2,2);
dimpulse(a,b,c,d);%离散系统的单位冲激响应
title('discrete impulse response');
subplot(2,2,3);
x0=[1 1 1 -1]';%离散系统初始状态
dinitial(a,b,c,d,x0);
axis([0 6 -0.5 2.5]);
title('discrete initial response');
subplot(2,2,4);
[z,p,k]=ss2zp(a,b,c,d,1);%状态方程转化为零极点增益形式
zplane(z,p);
title('discrete ploe-zero map');           

6、多输入多输出系统$\left\lbrace \begin{array}{l}

2.25 & -5 & -1.25 & -0.5\\

2.25 & -4.25 & -1.25 & -0.25\\

0.25 & -0.5 & -1.25 & -1\\

1.25 & -1.75 & -0.25 & -0.75

\end{array}\right\rbrack x+\left\lbrack \begin{array}{cc}

4 & 6\

2 & 4\

2 & 2\

0 & 2

y=\left\lbrack \begin{array}{cccc}

0 & 0 & 0 & 1\\

0 & 2 & 0 & 2

\end{array}\right.$，求单位阶跃响应和单位冲激响应。

clear
clc
a=[2.25 -5 -1.25 -0.5;
   2.25 -4.25 -1.25 -0.25;
   0.25 -0.5 -1.25 -1;
   1.25 -1.75 -0.25 -0.75];
b=[4,6;2,4;2,2;0,2];
c=[0 0 0 1;0 2 0 2];
d=zeros(2,2);
subplot(1,2,1);
step(a,b,c,d);
title('step response');
subplot(1,2,2);
impulse(a,b,c,d);
title('impule response');           

a0=[-1.6,-0.9,0,0;0.9,0,0,0;0.4,0.5,-5.0,-2.45;0,0,2.45,0];

b0=[1 0 1 0]';

c0=[1 1 1 1];

d0=[0];

t=0.5; %采样周期

[a,b,c,d]=c2dm(a0,b0,c0,d0,t,'tustin');%连续系统离散化

subplot(2,2,1);

dstep(a,b,c,d);%离散系统的单位阶跃响应

title('discrete step response');

subplot(2,2,2);

dimpulse(a,b,c,d);%离散系统的单位冲激响应

title('discrete impulse response');

subplot(2,2,3);

x0=[1 1 1 -1]';%离散系统初始状态

dinitial(a,b,c,d,x0);

axis([0 6 -0.5 2.5]);

title('discrete initial response');

subplot(2,2,4);

[z,p,k]=ss2zp(a,b,c,d,1);%状态方程转化为零极点增益形式

zplane(z,p);

title('discrete ploe-zero map');

matlab学习笔记4（控制系统仿真）

1、将5中的连续系统，以t=0.5 取样周期，采用双线性变换算法转换成离散系统，然后求出离散系统的单位阶跃响应、单位冲激响应及零输入响应，初始状态$x_0 ={\left\lbrack 1\;1\;1-1\right\rbrack }^T$ 。

clear
clc
a=[-1.6 -0.9 0 0;0.9 0 0 0;0.4 0.5 -5 -2.45;0 0 2.45 0];
b=[1,0,1,0]';
c=[1 1 1 1];
d=[0];
H=ss(a,b,c,d);
t=0.5;%采样周期
Gd=c2d(H,t,'tustin');%连续时间系统转化为离散时间系统
[ad,bd,cd,dd]=ssdata(Gd);
subplot(2,2,1);
dstep(ad,bd,cd,dd);
title('discrete step response');
subplot(2,2,2);
dimpulse(ad,bd,cd,dd);
title('discrete impulse response');
subplot(2,2,3);
x0=[1 1 1 -1]';
dinitial(ad,bd,cd,dd,x0);
axis([0 6 -0.5 2.5]);
title('discrete initial response');
subplot(2,2,4);
[z,p,k]=ss2zp(ad,bd,cd,dd,1);
zplane(z,p);
title('discrete ploe-zero map');           

2、设开环传递函数$H(s)=\frac{k(3s+1)}{s(2s+1)}$,绘制出通过单位负反馈构成的闭环系统的根轨迹.

clear
clc
num=[3,1];
den=[2,1,0];
rlocus(num,den);
sgrid
title('根轨迹');           

3、设开环系统$H\left(S\right)=\frac{k\left(S+5\right)}{S\left(S+2\right)\left(S+3\right)}$，绘制出闭环系统的根轨迹，并确定交点处的增益K。

clear
clc
num=[1 5];
den=[1 5 6 0];
rlocus(num,den);
sgrid
title('根轨迹');
[k,p]=rlocfind(num,den);           

Select a point in the graphics window

selected_point = - 1.9896 - 0.09931

10、已知开环传递函数为$H\left(S\right)=\frac{k}{S^4 +{16S}^3 +36S^2 +80S}$，绘制处闭环系统根轨迹，并分析其稳定性。

clearclcnum=[1];den=[1 16 36 80 0];rlocus(num,den);title('根轨迹');           

Select a point in the graphics window

select_point = 0.0152 + 3.13651

figure(2);k1=55;num1=k1*num;den1=conv(den,den);G1=tf(num1,den1);[num_F1,den_F1]=feedback(G1,1,-1);%单位负反馈G1_0=tf(num_F1,den_F1);impulse(G1_0);title('impulse response(k=55)');           

figure(3);
k2=56;
num2=k2*num;
den2=conv(den,den);
G2=tf(num2,den2);
[num_F2,den_F2]=feedback(G2,1,-1);
G2_0=tf(num_F2,den_F2);
impulse(G2_0);
title('impulse response(k=56)');