离散线性系统模型（Discrete linear system model）.docVIP

离散线性系统模型（Discrete linear system model） 1. time domain model: 1) filter function: uses recursive or non recursive filters to filter data Y=filter (B, a, x) a system consisting of vectors B and a to filter input x, and the output of the system is y [y, zf]=filter (B, a, x, Zi) the initial state of the zi: input signal, the zf: function, the returned system, the final state vector 2) impz function: implements the unit impulse response (B, a) of discrete systems (y=impz) Eg. when the input and output system of differential equation: Y (n) -0.8*y (n-1) -0.5*y (n-2) =0.7*x (n) +0.3*x (n-1), and the input of the system for the unit impulse response function, respectively using filter and impz to get the system unit impulse response %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. frequency domain model: 1. transfer function: freqz function: to achieve the transfer function response, using the base 2 FFT algorithm [h, f]=freqz (B, a, N, FS, B, a) FS parameters for the discrete system for sampling frequency, if fs=1 frequency axis gives the normalized frequency n [0 fs/2] frequency point selection, often take a power of 2, f returned to the frequency axis vector Eg.1 implements the input output function of the system: 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));