example2_5.m

原名：数字信号处理基础及MATLAB实现,解压缩密码：[e-mail]

%例2-5 程序:example2_5
n=-5:10;x=sin(pi*n/2);    
k=-100:100;w=(pi/100)*k;              %-pi到pi之间的频率
X=x*(exp(-j*pi/100)).^(n'*k);         %x的FT
%信号分解
[xe,xo,m]=evenodd(x,n);               %利用evenodd分解奇偶部
XE=xe*(exp(-j*pi/100)).^(m'*k);       %xe的FT
XO=xo*(exp(-j*pi/100)).^(m'*k);       %xo的FT
%校验
XR=real(X);                           %X的实部
error1=max(abs(XE-XR))                %XE与XR是否相等
XI=imag(X);                           %X的虚部
error2=max(abs(XO-j*XI))              %XO与XI是否相等
%用图形验证
subplot(1,1,1)
subplot(2,2,1);plot(w/pi,XR);grid
axis([-1,1,-2,2])
xlabel('频率（以pi为单位）');ylabel('Re(X)')
title('X的实部')
subplot(2,2,2);plot(w/pi,XI);grid
axis([-1,1,-10,10])
xlabel('频率（以pi为单位）');ylabel('Im(X)')
title('X的虚部')
subplot(2,2,3);plot(w/pi,real(XE));grid
axis([-1,1,-2,2])
xlabel('频率（以pi为单位）');ylabel('XE')
title('x偶部的FT')
subplot(2,2,4);plot(w/pi,imag(XO));grid
axis([-1,1,-10,10])
xlabel('频率（以pi为单位）');ylabel('XO')
title('x奇部的FT')