全データのフーリエ変換について

7 次查看（过去 30 天）

显示更早的评论

Kaede 2020-8-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/577168-

评论： Kaede 2020-8-12

在 MATLAB Online 中打开

このように2つの.wavファイルをフーリエ変換することによって周波数特性を見たいのですが,Nを1024というデータの特定部分だけをみるようなfftしかわかりません.

全データに対してフーリエ変換する方法を教えていただきたいです.

久しぶりにMATLABを触り忘れてしまっている部分がありますので初歩的な質問で申し訳ございません.

clear;
[p1,fs]=audioread('40kph_RN.wav');
[p2,fs]=audioread('2000rpm_Knock.wav');
fs=44100;
N=1024;
p2(448512,[1 2]) = [0,0];
smix=p1+p2;
%soundsc(smix,fs)
p1_fft=fft(p1,N);
p2_fft=fft(p2,N);
A=zeros(1,N/2+1);
B=zeros(1,N/2+1);
frequency=zeros(1,N/2+1);
for k=1:N/2+1;
    A(k)=abs(p1_fft(k));
    B(k)=abs(p2_fft(k));
    frequency(k)=(k-1)*fs/N;
end
dB_A=mag2db(A);
dB_B=mag2db(B);
plot(frequency,dB_A);
hold on
plot(frequency,dB_B);
hold off
legend('roadnoise','2000rpm.knock');