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

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

Melden Sie sich an, um zu kommentieren.

Antworten (0)

Produkte

Community Treasure Hunt

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

Start Hunting!