matlab有意思程序,求助一段matlab程序的意思

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

[y,Fs,bits]=wavread('file name');

y=y(:,2); sigLength=length(y);

Y =

fft(y,sigLength);

Pyy = Y.* conj(Y) /

sigLength;

halflength=floor(sigLength/2);

f=Fs*(0:halflength)/sigLength;

figure;plot(f,Pyy(1:halflength+1));xlabel('Frequency(Hz)');

t=(0:sigLength-1)/Fs;

figure;plot(t,y);xlabel('Time(s)');

第二段：

[k,Fs,bits]=wavread('file

name');

k=k(:,1);

y_temp=k(1:90000);

dfactor=3;

y=decimate(y_temp,dfactor);

[b,a] =

butter(10,200/(Fs/(dfactor*2)),'high');

y=filter(b,a,y);

y=interp(y,dfactor);

sigLength=length(y);

Y = fft(y,sigLength);

Pyy = Y.* conj(Y) / sigLength;

halflength=floor(sigLength/2);

f=Fs*(0:halflength)/sigLength;

figure;plot(f,Pyy(1:halflength+1));xlabel('Frequency(Hz)');

sigLength=length(y_temp);

Y = fft(y_temp,sigLength);

Pyy = Y.*

conj(Y) / sigLength;

halflength=floor(sigLength/2);

f=Fs*(0:halflength)/sigLength;

figure;plot(f,Pyy(1:halflength+1));xlabel('Frequency(Hz)');

t=(0:sigLength-1)/Fs;

figure;plot(t,y,t,y_temp);xlabel('Time(s)');

wavplay(y,Fs);

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！