example2_13.m

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

n=-5:10;
x=sin(pi*n/2);
k=-100:100;w=pi/100*k;
X=x*(exp(-j*pi/100)).^(n'*k);
[xe,xo,m]=evenodd(x,n);
XE=xe*(exp(-j*pi/100)).^(m'*k);
XO=xo*(exp(-j*pi/100)).^(m'*k);
XR=real(X);
error1=max(abs(XE-XR));
XI=imag(X);
error2=max(abs(XO-j*XI));
subplot(2,2,1);plot(w/pi,XR);grid;axis([-1,1,-2,2]);
xlabel('frequency');ylabel('Re(X)');
subplot(2,2,2);plot(w/pi,XI);grid;axis([-1,1,-10,10]);
xlabel('frequency');ylabel('Im(X)');
subplot(2,2,3);plot(w/pi,real(XE));grid;axis([-1,1,-2,2]);
xlabel('frequency');ylabel('XE');
subplot(2,2,4);plot(w/pi,imag(XO));grid;axis([-1,1,-10,10]);
xlabel('frequency');ylabel('XO');
disp(error1);
disp(error2);