このような以下のプロ​グラムを実行しました​.図は得られたdとL​chをFFTしたもの​なのですが,dはxを​単純遅延させただけな​ので本来形は同じにな​るはずなのですが,半​円を何度も描いている​ようなグラフになって​います.特にプログラ​ムに問題はなさそうな​のですが,どうしてで​しょうか.

2 次查看（过去 30 天）
Kaede
Kaede 2020-5-27
评论： Kaede 2020-6-3
clear;
[data,fs]=audioread('2000rpm.wav');
a=[0,1,0,0,0,0,0,0];
L=8;
Lch=data(:,1);
N=numel(Lch);
d=zeros(1,N);
for n=L:N
d(n)=a*Lch(n-L+1:n);
end
f=fs/1024:fs/1024:fs;
fft_Lch=fft(Lch,1024);%fftは同様にしてdも行いました
dB_Lch=mag2db(abs(fft_Lch));
plot(f,dB_Lch)
  4 个评论
显示 2更早的评论
Kenta
Kenta 2020-5-27
なるほど、理解しました。ありがとうございます。ちなみに、2000rpm.wavを添付していただくことは可能ですか？
Kaede
Kaede 2020-5-27
お返事ありがとうございます.wavファイルを添付できないみたいなのですが,解決策はございますでしょうか？

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Shoumei
Shoumei 2020-5-29
単純に0データを挿入して遅延させると、0データから、有効データの境界で不連続点が発生し、それが高調波成分となります。それがFFTの結果に現れているものと推測します。
予想しない結果が得られることはよくあります。今後も沢山そのような場面に遭遇すると思います。
そういうときは、まず単純なデータ（例えばSin波）などで、同じようなことをやってみて、現象を理解することをお勧めします。
  5 个评论
显示 3更早的评论
Kaede
Kaede 2020-6-3
いつもコメントいただきありがとうございます
Kaede
Kaede 2020-6-3
これはeが発散しているのが原因でした

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Fast Fourier Transforms 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!