stft1.m

基于matlab的功率谱分析以及短时傅立叶变换

[x1,fs1,bits1]=wavread('C:\Documents and Settings\hp\My Documents\MATLAB\ak47-1.wav');
%????
plot(real(x1));
hold on;
plot(imag(x1),'r*--');
legend('??','??');
xlabel('?? t');
ylabel('?? A');
%?????????
[s,f,t]=spectrogram(x1);
%??
figure(2)
plot(t,abs(s));
xlabel('?? t');
ylabel('?? A');
%??
figure(3)
plot(t,angle(s));
xlabel('?? t');
ylabel('?? phi');
%????
figure(4)
contour(t,f,abs(s));
xlabel('?? t');
ylabel('?? f');