概述
[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乘以抽样频率,即每秒抽样多少次,乘以时间
最后
以上就是高兴咖啡为你收集整理的利用MATLAB绘制信号时域波形和信号的频谱的全部内容,希望文章能够帮你解决利用MATLAB绘制信号时域波形和信号的频谱所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
发表评论 取消回复