No.6ベストアンサー
- 回答日時:
書き方にちょっと誤解があったかもしれませんので。
>やはりFFTが正統なのでしょうね
ガウス-ザイデルの方法は、FFTは使いません。
直接デコンボリューション(内部ではコンボリューションもしている)を行う手法です。
> #3にも書きましたとおり,補完してFFTでやってみました.
> 今のところ,期待した結果とはなっておりません.
単純にFFTして、わり算して求めても、測定値が実測の場合はノイズの影響でうまく復元しないでしょう。
実測データに対してやらなければならない手順は、
1.データ補間
No.5の方の言われる窓関数を使った方法はかなり有効です。
詳しくはNumerical Recipes in C などの本を参照して下さい。
2.ノイズフィルタで可能な限りノイズの除去
窓関数を使う方法をとれば、同時にフィルタの役目もしますので、これは必要なくなります。
保管方法によってはこれがないとまともな結果が得られません。
3.デコンボリューション
もっとも数学的にわかりやすいのはガウスザイデルの方法(ガウスの消去法で連立方程式を解く)ですね。
これ以外にも沢山の方法が考案されています。
h(t)が推定とのことですが、それだけではそれらしい結果が得られない理由にはなりません。
うまくデコンボリューション出来ない理由の大半はノイズによる影響です。
このデコンボリューションというのはかなり難しく、色々なやり方が研究されています。
デコンボリューション(deconvolution)、ガウス、ウィナーフィルタなどでwebを検索すると色々出てくると思いますよ。
では、がんばって下さい。
No.5
- 回答日時:
yの帯域に上限があることがはっきりしている場合には
直線補間より高度な
sin(x)/xを使う方法があります
これだと無限個のサンプリング値を必要とするので
窓関数w(x)をかけて
w(x)・sin(x)/x
を補間に使うと良いでしょう
するとサンプリング値を補間点近辺に限定できます
サンプリング間隔が大きいと単なる直線補間は余りにも荒すぎます
No.4
- 回答日時:
>逆畳み込み積分は応答関数の複素共役と畳み込み積分をするということで良いのでしょうか?
どのようにイメージされているかが、今一つ理解できませんでしたが、
y(t) = f(t) * x(t)
の操作をしてy(t)を得ているという点はよろしいわけですね。
原理的には、No.1の方が書いたように、フーリエ変換で表すと、x(t)をy(t)とx(t)から逆算出来ることになります。
(これが最も簡単で早い解法です)
が、もし実測波形の場合は、
y(t) = f(t)*x(t) + n(t)
とノイズの項n(t)が入ってしまうために、最もx(t)に近い波形を推定しなければ求めることが出来ません。
ノイズn(t)も実測できれば良いのでしょうけど、事実上無理ですよね?
また、わずかなノイズでもx(t)に大きく影響してしまうためにうまく復元できません。
なお、どの手法を選ぶにしてもデータは等間隔でないと、扱いがとても面倒(私にはその扱い方は理解できませんでした)なので、普通は補間などの方法で、等間隔にマッピングし直して処理します。
>単に相関を取ったのではなぜだめなのでしょうか?
デコンボリューションは、応答関数f(t)によってゆがめられた信号から、元の波形を推定することになりますので、相関(コリレーションの意味ですか?)では解けないと思いますが。
それとも相関という意味を、最小2乗法でfittingするようなイメージで使われていますか?
ならば、ガウス-ザイデル法がそれに該当します。
y1(t) = f(t) * x1(t)
x1(t)を初期値として、y(t)-y1(t)の残差から、次のx2(t)を推定して、また残差をみて、、、と繰り返して、残差が小さくなるようなxn(t)を求めることで、元波形を推定します。
では。
またまた回答ありがとうございました.経過報告です.
本屋にて信号処理関係を漁りましたが,直接デコンボルーションをする方法については全く見つかりませんでした.やはりFFTが正統なのでしょうね.
#3にも書きましたとおり,補完してFFTでやってみました.
今のところ,期待した結果とはなっておりません.
元データが2万5000点を,補完して2万点にしてしまったので,この時点でおかしいかも知れません.もう少しやってみます.
又,f(t)またはf(ω)が正確に判っているわけでは無く,判っているのは時定数だけなのです.mHz程度の応答であるため,手持ちの測定器では調べることもままならずというところです.
このセンサ自体にヒステリシスも含まれているようなので,計測方法自体も見直すことにします.
No.3
- 回答日時:
出力信号を直線補間をして等間隔にあるデータを求めfftを行えばいいと思います
もし精度を上げたかったら高次補間をすればいいと思います
xが未知なのだからfは予測値で我慢するしかありません
補完処理->FFTしてみました.
f(t)は一次遅延で時定数のみが判っており,そのまま計算したのでは期待した結果は得られませんでした.
きっと周波数特性が一致していないのだと思います.他にもヒステリシスが含まれていそうで,これらを元信号からできる限り除去しないと,信号処理だけでは厳しいようでした.
助言ありがとうございました.
No.2
- 回答日時:
その操作は、一般に「逆畳み込み積分」または「デコンボリューション」と呼ばれます。
やり方は実に様々あるので、お探しになってみて下さい。一番簡単なのはガウス-ザイデルの方法です。
(連立一次方程式を解く形となります)
データが等時間間隔でない場合は、扱いが大変面倒なので、一度等時間間隔の空間へ補間などを用いてマッピングし直してから計算すると良いでしょう。
観測波形にノイズなどがあると、強調されてしまいますので、事前にフィルタリングが必要になったりします。
(単純なFFTによる復元では、このノイズのために実用にはなりません。)
では。
回答ありがとうございます.
まだよく判っていないのですが,逆畳み込み積分は応答関数の複素共役と畳み込み積分をするということで良いのでしょうか?
単に相関を取ったのではなぜだめなのでしょうか?
No.1
- 回答日時:
Y(ω)=∫dt・exp(-i・ω・t)・y(t)
X(ω)=∫dt・exp(-i・ω・t)・x(t)
F(ω)=∫dt・exp(-i・ω・t)・f(t)
とする
Y(ω)=F(ω)・X(ω)
であるから
X(ω)=Y(ω)/F(ω)
である
これより
x(t)=1/(2・π)・∫dω・exp(i・ω・t)・X(ω)
=1/(2・π)・∫dω・exp(i・ω・t)・Y(ω)/F(ω)
もしy(t)が
y(t)=Σ(n)・x[n]・δ(t-t[n])
と表されるならば
Y(ω)=Σ(n)・x[n]・exp(-i・ω・t[n])
である
ただしiは虚数単位である(jにしたほうがいいと思いますが)
x(t),y(t),f(t)のうち何が時系列データになるのですか?
adc、dacがどこかには行っているのでしょうか?
もっと詳しい構成を教えてください
回答ありがとうございます.
得られているのはtとy(t)だけで,f(t)というのは部品カタログからの予測値です.
y(t)の測定が等時間間隔で無い(ある条件から外れたときはデータを取り込まないようになっているため),fftは使えないのです.
tが判っているので,フーリエ積分をコツコツやればいいのかと思ったのですが,やはり等時間間隔で無いでない為,積分値に誤差が生じてしまいそうです.
そのためフーリエ積分を用いないでやればいいのかと思ったのでした.
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 工学 制御工学の問題について 1 2022/10/07 20:25
- 工学 制御工学の問題について 2 2022/10/30 20:20
- 工学 制御工学の問題について 1 2022/11/01 23:45
- 工学 制御工学の問題について 1 2022/11/01 09:12
- エアコン・クーラー・冷暖房機 20年使ったエアコンを買換えても電気代は安くならない? 11 2023/06/29 14:54
- アルバイト・パート 超人手不足のパートの辞め方【アイデア求む!】 現在フランチャイズのリラクゼーションサロンで、リラクゼ 1 2023/04/22 21:59
- 工学 1つのタクトスイッチで複数の並列回路を閉じ、アースにつなげることができるスイッチはありますか。 2 2022/08/06 10:38
- その他(コンピューター・テクノロジー) 量子コンピュータの動作原理がわかりません。同じビットが、1でも0でも有って良いだろうか? 3 2023/02/04 03:20
- 工学 エミッタ接地増幅回路 入出力特性、周波数特性、位相特性それぞれの特性3種に与える回路定数の影響を答え 1 2023/05/29 22:43
- 加湿器・除湿機 加湿器の水量? 1 2023/04/26 03:59
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
力学の問題です
-
RL-C並列回路のインピーダ...
-
電気回路-EとIが同相なるとき
-
回転するコイルの誘導起電力に...
-
この交流回路の電源周波数をf1...
-
1.027の求め方について教えて下...
-
最高角加速度の求め方
-
力学
-
物理の単振動について
-
1質点系の運動方程式について
-
物理;角運動量
-
写真の図の問題で、t<0では、ス...
-
線形非同次微分方程式に関して
-
摩擦力の働く単振動について
-
複素誘電率と電磁波の損失エネ...
-
伝達関数とゲインについて
-
半径がr[m]のタイヤが角速度ω[r...
-
RL直列回路の電流ベクトルの...
-
この電気回路の時間関数が分か...
-
微分方程式 重ね合わせの原理
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
遠心力について。
-
大学の物理が難しすぎることに...
-
RL-C並列回路のインピーダ...
-
減衰係数の単位換算
-
単極誘導の説明の間違い
-
連続系一次遅れの離散系バージ...
-
1.027の求め方について教えて下...
-
半径がr[m]のタイヤが角速度ω[r...
-
複素振幅ってなんですか?
-
フーリエスペクトルの振幅について
-
周波数スペクトル図の、マイナ...
-
2物体の単振動
-
電荷qの荷電粒子が角速度ω、半...
-
物理の回路の問題です (2)の一...
-
太陽の赤緯の求め方教えてくだ...
-
物理の問題がどうしても解けま...
-
ローパスフィルタのカットオフ...
-
回転するコイルの誘導起電力に...
-
復元力の運動方程式F=-KxのK=mω...
-
2自由度系の固有振動数
おすすめ情報