离散线性系统模型.docVIP

离散线性系统模型 1.时域模型： 1) filter function:利用递归滤波器或非递归滤波器对数据进行滤波 y=filter(b,a,x) 由向量b和a组成的系统对输入x进行滤波，系统的输出为y [y,zf]=filter(b,a,x,zi) zi:输入信号的初始状态 zf:函数返回的系统最终状态向量 2) impz function:实现离散系统的单位冲激响应 y=impz(b,a) eg.当系统的输入输出差分方程为：y(n)-0.8*y(n-1)-0.5*y(n-2)=0.7*x(n)+0.3*x(n-1)，并且系统的输入为单 位冲激响应函数时，分别利用filter和impz得到系统单位冲激响应 %filter and impz function realize unit pulse response clear all; pulse=[1 zeros(1,63)]; % generate unit sample sequence b=[0.3 0.2]; % system parameter a=[1 -0.4 -0.7]; % system parameter h1=filter(b,a,pulse); h2=impz(b,a,64); subplot(211) stem(h1) title(filter function) subplot(212) stem(h2) title(impz function) 2.频域模型： 1.传递函数：freqz function:实现传递函数响应，采用基2的FFT算法 [h,f]=freqz(b,a,n,fs) b,a为离散系统的参数 fs为采样频率,若fs=1频率轴给出归一化频率 n为[0 fs/2]内选取的频率点数,常取2的幂次方,f为返回频率轴坐标向量 eg.1 实现系统的输入输出函数：H(z)=( 0.3+0.2*z^(-1) )/( 1-0.4*z^(-1)-0.7*z^(-2) ) % transfer function repinse clear all; fs=1000; % sample frequency unit in Hz b=[0.3 0.2]; % system parameter a=[1 -0.4 -0.7]; % system parameter [h,f]=freqz(b,a,256,fs); mag=abs(h); % compute magnitude ph=angle(h); % compute phase ph=ph*180/pi; % unit in degree % ph=unwrap(ph); subplot(211) plot(f,mag) grid xlabel(frequency(Hz)); ylabel(magnitude); subplot(212) plot(f,ph) grid xlabel(frequency(Hz)); ylabel(phase); 2.零极点增益： eg.1 当系统的输入输出差分方程为：y(n)-0.8*y(n-1)-0.5*y(n-2)=0.7*x(n)+0.3*x(n-1)，并且系统的输入 为单位冲激响应函数时,求其零极点及零极点增益 % zero/pole gain clear all; b=[0.3 0.2]; a=[1 -0.4 -0.7]; zr=roots(b); % or g=b(1)/a(1); pk=roots(a); % zplane(b,a); g=b(1)/a(1); 3.状态空间： eg 系统的输入输出函数：H(z)=( 0.3+0.2*z^(-1) )/( 1-0.4*z^(-1)-0.7*z^(-2) ) 的状态空间表示 % state-space method clear all; b=[0.3 0.2]; a=[1 -0.4 -0.7]; [A,B,C,D]=tf2ss(b,a) 4.带佘数的部分分式展开： 5.二次分式： 6.Lattice结构： eg.1 对于全极点的IIR系统，其传递函数为： H(z)=1/( 1-0.4*z^(-1)-0.7*z^(-2) ) 求其Lattice结构 %Latttice structure IIR system clear all; a=[