MATLAB在计算机控制系统中的应用第二课.ppt

文档介绍：

试绘制单位阶跃响应曲线。
解:%MATLAB程序如下:
num=10*[1 1];
den=[1 8 7 10];
step(num,den)

式中ωn为无阻尼振荡频率;ζ为阻尼比。
当ωn=6, ζ=,,,,1,2时,试绘制单位阶跃响应曲线;
当ζ=, ωn=2,4,6,8,10,12时,试绘制单位阶跃响应曲线。
解:%MATLAB程序如下:
wn=6;
zeta=[::,,];
figure(1)
hold on
for zet=zeta
num=wn.^2;
den=[1,2*zet*wn,wn.^2];
step(num,den)
end
hold off
②解:%MATLAB程序如下:
w=2:2:12;
zeta=;
figure(1)
hold on
for wn=w
num=wn.^2;
den=[1,2*zeta*wn,wn.^2];
step(num,den)
end
hold off


当ωn=5,ζ=,试绘制单位冲击响应曲线。
解:%MATLAB程序如下:
wn=6;
zeta=;
num=wn.^2;
den=[1,2*zeta*wn,wn];
impulse(num,den)

:
已知a=1,b=3,c=2,x=0,初始条件:y(1)=0,y(0)=0,求响应y(n)。
解(1)%对差分方程两边进行Z变换,求对应的变换方程的MATLAB程序如下:
syms a b c n ;
a=1;b=3;c=2;
l=a*sym('y(n+2)')+b*sym('y(n+1)')+c*sym('y(n)');
L=factor(ztrans(l));
R=0;
S=L-R
(2)%解变换方程,求出未知函数的象函数表达式,并进行Z反变换,程的MATLAB程序如下:
syms z Y y0 y1;
y0=0;y1=1;
S=z^2*Y-y0*z^2-y1*z+3*z*Y-3*y0*z+2*Y;
Y=solve(S,Y);
y=iztrans(Y)

=1s,脉冲传递函数为:
输入采样信号的Z变换为
试求系统输出信号y(n)。
解:%求系统输出采样信号Z变换C(z)的MATLAB
%语句如下:
syms z;
G=(+*z^(-1))/(1-*z^(-1));
R=z/(z-1);
C=R*G;
y=factor(iztrans(C))