|
Matlab:
- Fs = 1024; % 采样率
- N = 1024; % 采样点数
- n = 0:N-1; % 采样序列
- t = 0:1/Fs:1-1/Fs; % 时间序列
- f = n * Fs / N; %真实的频率
- %波形是由直流分量,50Hz正弦波正弦波组成
- x = 1024 + 1024*sin(2*pi*50*t + pi/3) ;
- y = fft(x, N); %对原始信号做FFT变换
- Mag = abs(y);
- subplot(2,1,1);
- plot(f, Mag);
- title('幅值谱');
- xlabel('频率');
- ylabel('幅值');
- subplot(2,1,2);
- plot(f, 20*log10(Mag));
- title('对数谱');
- xlabel('频率');
- ylabel('幅值/db')
复制代码
效果:
|
|