MATLAB仿真及电子信息应用.ppt

例6-12中周期矩形波的脉冲宽度不变，将重复周期由原来的6变为10，重绘其幅度谱，与上图6-9相比较，可以看出，当脉冲宽度不变、周期增大时，谱线的幅度随之减小，谱线之间的间隔变小，谱线更密。函数功能调用格式tf2ss传递函数模型转换为状态空间模型[A,B,C,D]=tf2ss(b,a)tf2zp传递函数模型转换为零-极点增益模型[z,p,k]=tf2zp(b,a)ss2tf状态空间模型转换为传递函数模型[b,a]=ss2tf(A,B,C,D,iu)iu为输入量ss2zp状态空间模型转换为零-极点增益模型[z,p,k]=ss2zp(A,B,C,D,iu)iu为输入量zp2tf零-极点增益模型转换为传递函数模型[b,a]=zp2tf(z,p,k)zp2ss零-极点增益模型转换为状态空间模型[A,B,C,D]=zp2ss(z,p,k)表6-2系统模型转换函数b=[156];a=[2604];[A,B,C,D]=tf2ss(b,a)A=-30-2100010B=100C=0.50002.50003.0000D=0[z,p,k]=tf2zp(b,a)z=-3.0000-2.0000p=-3.19580.0979+0.7850i0.0979-0.7850ik=0.5000例6-4某连续系统的传递函数为求该系统的状态空间模型和零-极点增益模型。例6-5某系统的传递函数为求该系统的系数向量和状态空间模型。z=[-1-2];p=[-3-45];k=2;[b,a]=zp2tf(z,p,k)b=0264a=12-23-60[A,B,C,D]=zp2ss(z,p,k)A=5.0000001.0000-7.0000-3.464103.46410B=100C=2.0000-8.0000-5.7735D=0A=[-21;-30];B=[21];C=[10];D=[1];[b,a]=ss2tf(A,B,C,D)b=1.00004.00004.0000a=1.00002.00003.0000[z,p,k]=ss2zp(A,B,C,D)z=-2.0000+0.0000i-2.0000-0.0000ip=-1.0000+1.4142i-1.0000-1.4142ik=16.3连续时间系统的时域分析LTI连续时间系统在时域中以常系数微分方程来描述：可利用dsolve函数求解常系数微分方程来得到。解：先求零输入响应，再求零状态响应，然后两部分加在一起。MATLAB程序eq1=D2y+5*Dy+6*y=0;%设定零输入条件下的微分方程ic1=y(0)=1,Dy(0)=1;%输入初始状态yzi=dsolve(eq1,ic1);%求解微分方程，得到零输入响应yzi=simplify(yzi)%化简零输入响应eq2=D2y+5*Dy+6*y=exp(-t)*Heaviside(t);%设定给定输入条件下的微分方程ic2=y(-0.01)=0,Dy(-0.01)=0;%设定初始状态为0yzs=dsolve(eq2,ic2);%求解微分方程，得到零状态响应yzs=simplify(yzs)%化简零状态响应y=yzi+yzs%全响应运行程序，得到：yzi=4*exp(-2*t)-3*exp(-3*t)yzs=1/2*heaviside(t)*(exp(-t)-2*exp(-2*t)+exp(-3*t))y=4*exp(-2*