%%%%%%傅里叶变换/逆变换/短时傅里叶变换%%%%%%
%[y,Fs]=wavread('C:\Users\HSF\Desktop\AgilentTechnologies.wav'); %读出信号,采样率和采样位数。
[y,Fs]=wavread('C:\Users\HSF\Desktop\sound\beiyou.wav');%读出信号,采样率和采样位数。
sound(y,Fs);
y=y(:,1);%我这里假设你的声音是双声道,我只取单声道作分析,如果你想分析另外一个声道,请改成y=y(:,2)
figure(1)
sigLength=length(y);
t=(0:length(y)-1)/Fs; %计算时间轴
subplot(3,1,1);plot(t,y);xlabel('Time(s)');%在第一个窗口画波形
grid on;
n=0:sigLength-1;
Y= fft(y,length(y));%在sigLength这个有限区间内做快速傅立叶变换
mag=abs(Y);
f=Fs*n/sigLength;
subplot(3,1,2);plot(f,abs(Y));
xlabel('Frequency(Hz)');
%在第一个窗口画率谱,
grid on;
subplot(3,1,3);
xifft=ifft(Y);
sound(xifft,Fs);
plot(t,xifft);xlabel('Time(s)')
grid on;
figure(2)
Nw=20; %窗函数长window length
L=Nw/2; %窗函数每次移动的样点数
Ts=round((sigLength-Nw)/L)+1; %计算把数据y共分成多少段
nfft=512; % FFT的长度
TF=zeros(Ts,nfft); %将存放三维图谱,先清零
for i=1:Ts
xw=y((i-1)*L+1:i*L+L); %取一段数据
temp=fft(xw,nfft); %FFT变换
%temp=fftshift(temp); %频谱以0频为中心
TF(i,:)=temp; %把谱图存放在TF中代表TF矩阵中的第i行
end
fn=(1:nfft)*Fs/nfft;
tn=(1:Ts)*Nw/2/Fs;
[T,F]=meshgrid(tn,fn);
mesh(T,F,abs(TF.')); %三维绘图
axis tight;
title('三维时频图');
ylabel('Frequency(Hz)');
xlabel('Time(s)');
grid on;
仿真结果图如下:?
本站仅提供存储服务,所有内容均由用户发布,如发现有害或侵权内容,请
点击举报。