zchirpz.m

用MATLAB实现的DSP函数

%Chirp Z变换
close all;clear all;clc;
h=fir1(50,125/300,boxcar(51));
fs=1000;
f1=100;
f2=200;
m=1024;
w=exp(-j*2*pi*(f2-1)/(m*fs));
a=exp(j*2*pi*f1/fs);
y=fft(h,1024);
z=czt(h,m,w,a);
fy=((0:length(y)-1)'*1000/length(y));
fz=((0:length(z)-1)'*(f2-f1)/length(z)+f1);
figure, subplot(1,2,1),plot(fy(1:512),abs(y(1:512))),title('FFT的频率响应');
subplot(1,2,2),plot(fz,abs(z)),title('Chirp z的频率响应');