%上面是一段利用FFT变换检测噪音的程序,x是随机加噪的 t=0:0.001:1;    %采样频率1000Hzx=sin(2*pi*100*t)+sin(2*pi*200*t)+rand(size(t));%x长度1000  subplot(2,1,1)plot(x(1:50));axis on;xlabel('Variable t');       % 加X轴说明ylabel('Variable x(t)'); y=fft(x,512);   %程序中只是选取前面512个点做FFT，其实点数可以任意选取，一般不小于数据长度，% 甚至FFT长度可以不选取2的幂次。如果选取2的幂次意味着可以使用碟型展开计算FFT，容易硬件实现。仿真时可不拘泥于2的幂次。 f=1000*(0:256)/512; %采样频率1000Hz，数据窗长度512，频率分辨率1000/512，所以横坐标为：(0:256)×1000/512。 subplot(2,1,2)plot(f,y(1:257));axis on;xlabel('Variable f');       % 加X轴说明ylabel('Variable A(f)');

评论