dfttest.m

自己编写的自适应notch fileter算法

n=0:1/3000:1; d=0.1*sin(2500*pi*n);
subplot(311)
%q=length(t);
plot(n,d);
%axis([0,t(q),-0.3,0.3]);
title('zhengxuan信号源曲线图');

%X=dft(d,1000);
%magX=abs(X(1:551));
%subplot(312);
k=0:99;
w=2*pi/100*k;
%plot(w/pi,magX);

subplot(313);
Y=fft(d,100);
magY=abs(Y(1:100));
plot(w/pi,magY,'r');
xlabel('以pi为单位的频率');
%axis([0,1,0,7]);
grid;