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

21 views (last 30 days)
Kaede
Kaede on 9 Aug 2020
Commented: Kaede on 12 Aug 2020
このように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');
  3 Comments
Kaede
Kaede on 12 Aug 2020
確認いただきありがとうございます.

Sign in to comment.

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!