微分要素のインパルス応答の周波数領域こんにちはゲストさん。先日MATLABで次のような波形を作りインパルス応答を調べました。出力y(t)=微分フィルタのよって出力された波形(cos波)X=fftshift(fft(x)),Y=fftshift(fft(y))としてH=Y./Xとしたところ、添付画像のような結果が得られました。時間領域でのインパルス応答h(t)は以前同じような質問がなされていたので理解できましたが、この周波数領域(ラプラス変換をするならばS領域)でこのような結果が得られる理由をご教授いただきたく思います。※t=((1:N)-1)*(T/N)omega=((1*N)-N/2)*(2*pi/T)と設定しました。今日は、インパルス応答は入力信号にユニット・インパルス信号を入力したときの出力の応答でそれが系の伝達関数になります。MATLABは経験がありませんのでoboro7411さんの書かれた式X=fftshift(fft(x)),Y=fftshift(fft(y))は良く理解できませんが、少なくとも入力信号X(t)はsin波ではなくユニット・インパルス信号を使用しなければなりません。H(s)になるのでインパルス信号を入力したときの出力はその系の伝達関数を表しているというのがインパルス応答の特徴です。したがって、入力信号にユニット・インパルス以外の信号を入力してもその系の伝達関数は得られませんので計算結果はH(s)の周波数特性を表してはいないということになります。返事が遅れて申し訳ありません。質問内容に微分フィルタと書かせていただきましたが実際のプログラムを見るとy(t)=(x(t+⊿t)-x(t))/⊿tというような式よって微分フィルタを表現されていたということを記入し忘れておりました。これにより理論通りの微分フィルタにはならないためより深く考察すれば図のような結果が得られそうだという結論に達しました。回答していただいたのに申し訳ない限りです。ありがとうございました。返事が遅れて申し訳ありません。いまは事情により手元にプログラムがないため再実験はできません。回答ありがとうございます。出力はその微分つまりy(t)=cosω0tになるかと思います。iωX(ω)となると思うのです。matlabで勉強中なんですが、以下の操作をするとどうなるのか教えて頂きたいです。ゆらぎについて研究していています。matlabで信号処理を学んでいます。Reserved.掲載情報の著作権は提供元企業等に帰属します。