MATLAB使用的说明 .ppt

由 [A,B,C,D]=tf2ss([1],[1 5 10]) 可得 所以系统的状态方程为 C=[0 1] D=0 例2 已知某连续系统的状态方程和输出方程为 其初始状态和输入分别为 求该系统的系统函数矩阵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]; x0=[2 -1]; dt=0.01;t=0:dt:2; f(:,1)=ones(length(t),1); f(:,2)=exp(-3*t); sys=ss(A,B,C,D); y=lsim(sys,f,t,x0); 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); x0=[1; -1]; N=10;f=ones(1,N); sys=ss(A,B,C,D,[]); y=lsim(sys,f,[],x0); 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); F(jw)= Sa2(w / 2) 解: 图示三角波可表示为 三角波信号频谱的理论值为 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);F=zeros(1,N); for k=1:N F(k)=quad8(sf1,-1,1,[],[],w(k)); end figure(1); plot(w,real(F));title() xlabel(\omega);ylabel(F(j\omega)); figure(2); plot(w,real(F)-sinc(w/2/pi).^2); xlabel(\omega);title(计算误差); 运行结果 连续系统频率响应的计算 周期信号通过系统的响应 离散系统频率响应的计算 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+H*cos(w0*n*t+phi)*sinc(n*0.5); end plot(t,x