fs=2048;%设定采样频率
N=1024;%设定数据长度
i=0:N-1;
t=i/fs;
f=100;%设定正弦信号频率
%生成正弦信号
x=cos(3*pi*f*t);
x1=x.*exp(-2*j*pi*f0/fs*t);
bw=fs/(2*50);b=fir1(32,2*bw/fs);
y=filter(b,1,x);
yd=zeros(1,1024);
for p=1:length(yd)
yd(p)=y(p*50);
end
Yd=fft(yd,1024) ;
其实就是用个滤波器滤波 不知道怎么错了 按照论文抄的
