exe0410_03.m

来自「Matlab程序为清华大学胡广书版信号处理的上机习题答案」· M 代码 · 共 28 行

clear all；

% 构造三个不同频率的正弦信号的叠加作为试验信号
N=64;
f1=10.8;f2=11.75;f3=12.55;fs=40;
n=0:N-1;
t=2*pi*n/fs;
x=sin(f1*t)+sin(f2*t)+sin(f3*t);

% 第一组参数 
M=50;
f0=9;
DELf=0.2;
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(211);plot(n2,abs(Y3));grid on;

% 第二组参数 
M=60;
f0=8;
DELf=0.12;
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(212);plot(n2,abs(Y3));grid on;