Clear Filters
Clear Filters

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

17 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
Naoya
Naoya on 11 Aug 2020
fft関数は2つめの引数を省略する場合、信号の点数でfftをしますので、上記で良さそうです。
Kaede
Kaede on 12 Aug 2020
確認いただきありがとうございます.

Sign in to comment.

Answers (0)

Products

Community Treasure Hunt

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

Start Hunting!