利用MATLAB绘制信号时域波形和信号的频谱

[x,Fs] = wavread('test1');%读取音频数据

x = x(:,1);

x = x';

N = length(x);%求取抽样点数

t = (0:N-1)/Fs;%显示实际时间

y = fft(x);%对信号进行傅里叶变换

f = Fs/N*(0:round(N/2)-1);%显示实际频点的一半

subplot(211);

plot(t,x,'g');%绘制时域波形

axis([0 max(t) -1 1]);

xlabel('Time / (s)');ylabel('Amplitude');

title('信号的波形');

grid;

subplot(212);

plot(f,abs(y(1:round(N/2))));

xlabel('Frequency / (s)');ylabel('Amplitude');

title('信号的频谱');

grid;

%% 以实际频点显示或者显示频率的一半

f = Fs/N*(0:N-1);

f = Fs/N*(0:round(N/2)-1);

%%以实际频点显示，显示的是Fs的一半

[h,w] = freqz(a,b,512,Fs);

%% w介于0--pi之间，除以pi，归一化到0--1，乘以奈奎斯特频率，实现显示频率的一半

[h,w] = freqz(a,b,512);

f = Fs*w/(2*pi);

%% 显示实际的时间，抽样点数除以抽样频率

t = (0:N-1)/Fs;

%% 信号频率，采样频率，关系

T/deltat = Fs/fs = T * Fs = N %抽样点数等于信号周期T乘以抽样频率，即每秒抽样多少次，乘以时间