example2_6.m

来自「西交大刘树棠写的《数字信号处理》那本书的所有源代码」· M 代码 · 共 14 行

w=[0:1:500]*pi/500;
X=exp(j*w)./(exp(j*w)-0.5*ones(1,501));
magX=abs(X);angX=angle(X);
realX=real(X);imagX=imag(X);
subplot(2,2,1);plot(w/pi,magX);grid;
xlabel('frequency');ylabel('amplitude');title('amplitude');
subplot(2,2,2);plot(w/pi,angX);grid;
xlabel('frequency');ylabel('phase');title('phase');
subplot(2,2,3);plot(w/pi,realX);grid;
xlabel('frequency');ylabel('real');title('real');
subplot(2,2,4);plot(w/pi,imagX);grid;
xlabel('frequency');ylabel('image');title('image');