%上面是一段利用FFT变换检测噪音的程序,x是随机加噪的

t=0:0.001:1;    %采样频率1000Hz
x=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)');

评论