信号的频谱分析
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
信号的频谱分析
music1.wav的处理
[x, b]=audioread('music1.wav');
% 原始
x = x(:,1);
x = x';
N = length(x);%求取抽样点数
t = (0:N-1)/b;%显示实际时间
y = fft(x);%对信号进行傅里叶变换
f = b/N*(0:round(N/2)-1);%显示实际频点的一半subplot(3,2,1);
plot(t,x,'g');%绘制时域波形
xlabel('Time / (s)');ylabel('Amplitude');
title('信号的波形');
grid;
subplot(3,2,2);
plot(f,abs(y(1:round(N/2))));
xlabel('Frequency / (s)');ylabel('Amplitude'); title('信号的频谱');
grid;
sound(x ,b);
% 变低
x1=resample(x,10,8);%变低
N = length(x1);%求取抽样点数
t = (0:N-1)/b;%显示实际时间
y = fft(x1);%对信号进行傅里叶变换
f = b/N*(0:round(N/2)-1);%显示实际频点的一半subplot(3,2,3);
plot(t,x1,'g');%绘制时域波形
xlabel('Time / (s)');ylabel('Amplitude');
title('信号的波形');
grid;
subplot(3,2,4);
plot(f,abs(y(1:round(N/2))));
xlabel('Frequency / (s)');ylabel('Amplitude'); title('信号的频谱');
grid;
sound(x1,b);
% 变高
x2=resample(x,10,13);%变高
N = length(x2);%求取抽样点数
t = (0:N-1)/b;%显示实际时间
y = fft(x2);%对信号进行傅里叶变换
f = b/N*(0:round(N/2)-1);%显示实际频点的一半subplot(3,2,5);
plot(t,x2,'g');%绘制时域波形
xlabel('Time / (s)');ylabel('Amplitude');
title('信号的波形');
grid;
subplot(3,2,6);
plot(f,abs(y(1:round(N/2))));
xlabel('Frequency / (s)');ylabel('Amplitude'); title('信号的频谱');
grid;
sound(x2,b);
music2.wav的处理
[x, b]=audioread('music2.wav');
% 原始
x = x(:,1);
x = x';
N = length(x);%求取抽样点数
t = (0:N-1)/b;%显示实际时间
y = fft(x);%对信号进行傅里叶变换
f = b/N*(0:round(N/2)-1);%显示实际频点的一半subplot(3,2,1);
plot(t,x,'g');%绘制时域波形
xlabel('Time / (s)');ylabel('Amplitude');
title('信号的波形');
grid;
subplot(3,2,2);
plot(f,abs(y(1:round(N/2))));
xlabel('Frequency / (s)');ylabel('Amplitude'); title('信号的频谱');
grid;
sound(x ,b);
% 儿童
x1=resample(x,10,13);%变低
N = length(x1);%求取抽样点数
t = (0:N-1)/b;%显示实际时间
y = fft(x1);%对信号进行傅里叶变换
f = b/N*(0:round(N/2)-1);%显示实际频点的一半subplot(3,2,3);
plot(t,x1,'g');%绘制时域波形
xlabel('Time / (s)');ylabel('Amplitude');
title('信号的波形');
grid;
subplot(3,2,4);
plot(f,abs(y(1:round(N/2))));
xlabel('Frequency / (s)');ylabel('Amplitude'); title('信号的频谱');
grid;
sound(x1,b);
% 男人
x2=resample(x,10,8);%变高
N = length(x2);%求取抽样点数
t = (0:N-1)/b;%显示实际时间
y = fft(x2);%对信号进行傅里叶变换
f = b/N*(0:round(N/2)-1);%显示实际频点的一半subplot(3,2,5);
plot(t,x2,'g');%绘制时域波形
xlabel('Time / (s)');ylabel('Amplitude');
title('信号的波形');
grid;
subplot(3,2,6);
plot(f,abs(y(1:round(N/2))));
xlabel('Frequency / (s)');ylabel('Amplitude'); title('信号的频谱');
grid;
sound(x2,b);