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

Kaede

2020 8 9

0 个回答

7 次查看（30 天）

0 个投票

このように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');