数字信号处理实验报告

x=sin(f1*t)+sin(f2*t)+sin(f3*t); M=N;

W=exp(-j*2*pi/M);

% A=1时的czt变换 A=1;

Y1=czt(x,M,W,A); subplot(311)

plot(n1,abs(Y1(1:N/2)));grid on; % DTFT Y2=abs(fft(x)); subplot(312)

plot(n1,abs(Y2(1:N/2)));grid on; % 详细构造A后的czt M=60; f0=7.2; DELf=0.05; A=exp(j*2*pi*f0/fs); W=exp(-j*2*pi*DELf/fs); Y3=czt(x,M,W,A);

n2=f0:DELf:f0+(M-1)*DELf; subplot(313);plot(n2,abs(Y3));grid on;

13

（三）、实验题目

设x(n)由三个实正弦组成，频率分别是8Hz,9Hz,10Hz,抽样频率为60 Hz,时域取256点，作CZT变换、IFFT变换和FFT变换，观察波形，更改参数，得出不同参数下的CZT变换波形。给出实验报告。

3.1.1实现CZT变换、IFFT变换和FFT变换任务的程序为： clear all；

% 构造三个不同频率的正弦信号的叠加作为试验信号 N=256;

f1=8;f2=9;f3=10;fs=60; a1=5;a2=3;a3=4; stepf=fs/N; n=0:N-1; t=2*pi*n/fs; n1=0:stepf:fs/2-stepf;

x=a1*sin(f1*t)+a2*sin(f2*t)+a3*sin(f3*t); M=N;

W=exp(-j*2*pi/M);% A=1时的czt变换 A=1;

Y1=czt(x,M,W,A); subplot(321)

14

plot(n1,abs(Y1(1:N/2))); title('A=1时的czt变换'); grid on; % DTFT Y2=abs(fft(x)); subplot(322)

plot(n1,abs(Y2(1:N/2))); title('DTFT');

grid on;% 详细构造A后的czt M=60; f0=7.2; DELf=0.05; A=exp(j*2*pi*f0/fs); W=exp(-j*2*pi*DELf/fs); Y3=czt(x,M,W,A);

n2=f0:DELf:f0+(M-1)*DELf; subplot(323); plot(n2,abs(Y3)); title('详细构造A后的czt'); grid on；

% 应用FFT 求频谱； subplot(324);

plot(x(1:N/4));%作图向量x的第一到第N/4个值；

15

f=-0.5:1/N:0.5-1/N;

X=fft(x); % 快速傅里叶变换； y=ifft(X);% 快速傅里叶逆变化 subplot(325); plot(f,fftshift(abs(X))); title('快速傅里叶变换'); grid on; subplot(326);

plot(real(y(1:N/4))); %作图向量y的第一到第N/4个值的实部； title('作图向量y的第一到第N/4个值的实部');

3.1.2实现CZT变换、IFFT变换和FFT变换任务的结果为：

16