exp3_16.m

是Digital signal processing--using Matlab中的一些函数

a=[1 -1.76 1.1829 -0.2781];
b=[0.0181 0.0543 0.0543 0.0181];
m=0:length(a)-1;
n=0:length(b)-1;
k=0:1:500;
w=pi*k/500;
x=b*exp(-j*m'*w)./(a*exp(-j*n'*w));%系统函数
%直接写出有理分式，没有调用其他函数
realx=real(x);
imagx=imag(x);
magx=abs(x);
angx=angle(x);
subplot(2,2,1),plot(w/pi,realx),title('realx');grid on;
subplot(2,2,2),plot(w/pi,imagx),title('imagx');grid on;
subplot(2,2,3),plot(w/pi,magx),title('magx');grid on;
subplot(2,2,4),plot(w/pi,angx),title('angx');grid on;