数字信号处理的思考
请看代码：
1 clf;fs=100; %采样频率 2 Ndata=128; %数据长度 3 N=256; %FFT的数据长度 4 n=0:Ndata-1;t=n/fs; %数据对应的时间序列 5 %x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); %时间域信号 6 x=0.5*sin(2*pi*20*t); 7 y=fft(x,N); %信号的Fourier变换 8 mag=abs(y); %求取振幅 9 f=(0:N-1)*fs/N; %真实频率 10 % subplot(2,2,1), 11 plot(f(1:N/2),mag(1:N/2)*2/N); %绘出Nyquist频率之前的振幅 12 xlabel('频率/Hz');ylabel('振幅'); 13 title('Ndata=128 Nfft=256');grid on;
请看图：
http://blog.163.com/fei_lai_feng/blog/static/9289962200971751114547/
