信号与系统MATLAB课稿.pptVIP

例2 已知某连续系统的状态方程和输出方程为 其初始状态和输入分别为 求该系统的系统函数矩阵H(s) 和输出。 计算系统函数矩阵H(s) A=[2 3;0 -1];B=[0 1; 1 0]; C=[1 1; 0 -1];D=[1 0; 1 0]; [num1,den1]=ss2tf(A,B,C,D,1) [num2,den2]=ss2tf(A,B,C,D,2) 运行结果 num1 =1 0 -1 1 -2 0 den1 =1 -1 -2 num2 =0 1 1 0 0 0 den2 =1 -1 -2 所以系统函数矩阵H(s)为 计算输出 %Program 9_1 A=[2 3;0 -1];B=[0 1; 1 0]; C=[1 1; 0 -1];D=[1 0; 1 0]; q0=[2 -1]; dt=0.01;t=0:dt:2; x(:,1)=ones(length(t),1); x(:,2)=exp(-3*t); sys=ss(A,B,C,D); y=lsim(sys,x,t,q0); subplot(2,1,1);plot(t,y(:,1),r); ylabel(y1(t));xlabel(t); subplot(2,1,2);plot(t,y(:,2)); ylabel(y2(t));xlabel(t); 运行结果 例3 已知某离散系统的状态方程和输出方程为 其初始状态和输入分别为 求该系统的输出。 计算输出 %Program 9_2 A=[0 1; -2 3];B=[0;1]; C=[1 1; 2 -1];D=zeros(2,1); q0=[1; -1]; N=10;x=ones(1,N); y=dlsim(A,B,C,D,x,q0); subplot(2,1,1);y1=y(:,1); stem((0:N-1),y1); xlabel(k);ylabel(y1); subplot(2,1,2);y2=y(:,2); stem((0:N-1),y2); xlabel(k);ylabel(y2); 运行结果 例3 试用数值方法近似计算三角波信号的频谱 function y=sf1(t,w);y=(t=-1 t=1).*(1-abs(t)).*exp(-j*w*t); w=linspace(-6*pi,6*pi,512); N=length(w);X=zeros(1,N); for k=1:N X(k)=quad8(sf1,-1,1,[],[],w(k)); end figure(1); plot(w,real(X));title() xlabel(\omega);ylabel(‘X(j\omega)); figure(2); plot(w,real(X)-sinc(w/2/pi).^2); xlabel(\omega);title(计算误差); 例3 试用数值方法近似计算三角波信号的频谱 运行结果 利用MATLAB进行系统频域分析 连续系统频率响应的计算 周期信号通过系统的响应 离散系统频率响应的计算 一、连续系统频率响应的计算 b 分子多项式系数 a 分母多项式系数 w 需计算的H(jw)的抽样点 (数组w中少需包含两个w的抽样点)。 计算频响的MATLAB函数 H=freqs(b,a,w) 一、连续系统频响特性的计算 例1 三阶归一化的Butterworth低通滤波器的系统函数为 w=linspace(0,5,200); b=[1];a=[1 2 2 1]; h=freqs(b,a,w); subplot(2,1,1); plot(w,abs(h)); subplot(2,1,2); plot(w,angle(h)); 试画出|H(jw)| 和?(w)。 一、连续系统频响特性的计算 三阶Butterworth低通滤波器的幅度响应和相位响应 二、周期信号通过系统的响应 例2 周期方波通过RC系统的响应。 二、周期信号通过系统的响应 例2 周期方波通过RC系统的响应。 %p5_2 Periodic signal pass LTI system T=4;w0=2*pi/T;RC=0.1; t= -6:0.01:6;N=51; c0=0.5;xN=c0*ones(1,length(t)); %dc for n=1:2:N % even harmonics are zero H=abs(1/(1+j*RC*w0*n)); phi=angle(1/(1+j*RC*w0*n)); xN=xN