example5_2.m

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

%例5-2程序，example5_2
%信号产生
t = 0:0.025:3.1750; % 采样频率 200Hz . 128个样本点
x =sin(2*pi*7*t)+sin(2*pi*7.5*t)+sin(2*pi*8*t)+sin(2*pi*8.5*t);% 信号频率由 8Hz 8.22Hz和 9Hz的正弦信号组成

figure(1)
subplot(2,2,1);
plot(1:128,x)
title('x(n)')
xlabel('t')
grid on

%FFT变换
Y_FFT = fft(x,128);
m1 = abs(Y_FFT);% p1 = unwrap(angle(Y_FFT));  幅度和相位
f1 = (0:length(Y_FFT)-1)*40/length(Y_FFT); % Frequency vector
subplot(2,2,2);
plot(f1(1:64),m1(1:64));
title('FFT');
xlabel('Hz')
grid on


% CZT_DFT变换
Y_CZT_DFT=czt(x);
m2 = abs(Y_CZT_DFT);%  p2 = unwrap(angle(Y_CZT_DFT)); 幅度和相位
f2 = (0:length(Y_CZT_DFT)-1)*40/length(Y_CZT_DFT); % Frequency vector
subplot(2,2,3);
plot(f2(1:64),m2(1:64));
title('用CZT做DFT');
xlabel('Hz')
grid on

% CZT变换
fs =40; f1 = 6; f2 = 10;  % in hertz
M = 50;
w = exp(-j*2*pi*(f2-f1)/(M*fs));
a = exp(j*2*pi*f1/fs);
z = czt(x,M,w,a);
subplot(2,2,4);
fz = ((0:length(z)-1)'*(f2-f1)/length(z)) + f1;
plot(fz,abs(z)); 
title('CZT');
xlabel('Hz')
grid on