利⽤matlab对⾳频做FFT进⾏频谱分析利⽤matlab对⾳频的左右声道做频谱分析。 syms y fs;
syms in_put;
in_put= 'C:\Users\getao\Desktop\⽴体声.wav'
%in_put 为⾳频⽂件所在的位置,.wav⽂件为⽂件的名称。
[y,fs] = audioread(in_put);
info=audioinfo(in_put) ;
sound(y,fs);
T=1/fs;
t=(0:length(y)-1)*T;
f=(0:length(y)-1)*fs/length(y);
figure(1);
yz=y(:,1);
subplot(2,2,1);
plot(t,yz);
xlabel('时间');
ylabel('振幅');
yz2=y(:,2);
%提取⾳频信号的右声道数据
subplot(2,2,2);
plot(t,yz2);
title('原始信号时域(右声道)');
xlabel('时间');
ylabel('振幅');
subplot(2,2,3);
n=length(yz);
y1=fft(yz,n);
F=fs/length(yz);
plot(f,abs(y1));
title('原始信号频谱(左声道)');
xlabel('F(Hz)');
subplot(2,2,4);
n2=length(yz2);
y1=fft(yz2,n2);
F=fs/length(yz2);
支撑架
plot(f,abs(y1));
title('原始信号频谱(右声道)');
xlabel('F(Hz)');
ylabel('H(jw)');眼模
syms y fs;
syms in_put;
in_put= 'C:\Users\getao\Desktop\右声道.wav';
[y,fs] = audioread(in_put);
info=audioinfo(in_put) ;
stc2052sound(y,fs);
T=1/fs;
t=(0:length(y)-1)*T;
f=(0:length(y)-1)*fs/length(y);
figure(1);
音调控制电路
yz=y(:,1);
subplot(2,2,1);
plot(t,yz);
title('原始信号时域(左声道)');
xlabel('时间');
ylabel('振幅');
yz2=y(:,2);
subplot(2,2,2);
plot(t,yz2);
title('原始信号时域(右声道)');
xlabel('时间');
ylabel('振幅');
subplot(2,2,3);
n=length(yz);
y1=fft(yz,n);
F=fs/length(yz);
plot(f,abs(y1));
title('原始信号频谱(左声道)');
薄膜发电xlabel('F(Hz)');
ylabel('H(jw)');
subplot(2,2,4);
n2=length(yz2);
y1=fft(yz2,n2);
F=fs/length(yz2);
plot(f,abs(y1));
title('原始信号频谱(右声道)');
xlabel('F(Hz)');
ylabel('H(jw)');
下图分别为三段⾳频的分析,三段⾳频⽂件对同⼀⾸⾳乐的某⼀段进⾏录⾳得到。⽴体声左右声道都录制了声⾳,左右声道的⽂件则通过录⾳设备的两端进⾏分别录制。
⽴体声 ⾳频分析
左声道 ⾳频分析 右声道 ⾳频分析
从上图可以看出,该段声⾳的频率在0-10KHz范围内,在未录制声⾳的⾳道有杂⾳的⼲扰使其频谱有尖峰的存在。