0
1

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?

ノイズに鈍感なDMD

Last updated at Posted at 2025-04-30

初めに

DMDの中心的な考え方は,大雑把に理解したつもりである(こちら).

DMDはノイズに敏感だと聞く.本当だろうか.どう改良するのか.

MAIN

Data

$\mathrm{Re} = 100$の流れに,Gaussianなノイズを加える.適当に$\sim \mathcal{N} (0, 1)$とした.なお,前回は$0.01 \ [\mathrm{s}]$毎のスナップショットを取ったが,今回は$0.1 \ [\mathrm{s}]$毎とした.

$\text{Data}$
$\text{w/o noise}$ wo_noise.gif
$\text{w/ noise}$ w_noise.gif

POD: Proper Orthogonal Decomposition

w/o noise

スナップショットの撮影間隔が前回と違うので,念の為,確認しておく.

singular_value.png
POD_mode_0.png
POD_mode_1.png
POD_mode_2.png
POD_mode_3.png
POD_mode_4.png
POD_mode_5.png
POD_mode_6.png
POD_mode_7.png
POD_mode_8.png
POD_mode_9.png
POD_mode_10.png
POD_mode_11.png
POD_mode_12.png
POD_mode_13.png
POD_mode_14.png
POD_mode_15.png
POD_mode_16.png
POD_mode_17.png
POD_mode_18.png
POD_mode_19.png

時間モードの解像能が違うが,殆ど同じような挙動をしている.

w/ noise

ノイズが入ると,かなりのモードが影響を受ける.本当に少しだけのleading modeが存在し,その後ひどくplateauしている.

singular_value.png
POD_mode_0.png
POD_mode_1.png
POD_mode_2.png
POD_mode_3.png
POD_mode_4.png
POD_mode_5.png
POD_mode_6.png
POD_mode_7.png
POD_mode_8.png
POD_mode_9.png
POD_mode_10.png
POD_mode_11.png
POD_mode_12.png
POD_mode_13.png
POD_mode_14.png
POD_mode_15.png
POD_mode_16.png
POD_mode_17.png
POD_mode_18.png
POD_mode_19.png

平均場とほんの幾つかのleading modeを除けば,殆どがノイズに支配されている.ノイズのことをエネルギーと解釈してしまっている.PODはノイズに敏感と言えるだろう.

singular_value.png
singular_value.png
singular_value.png
singular_value.png

DMD: Dynamic Mode Decomposition

w/o noise

前回も実施したが,是非に何度でも.

eigenvals.png
DMD_mode_0.png
DMD_mode_1.png
DMD_mode_2.png
DMD_mode_3.png
DMD_mode_4.png
DMD_mode_5.png
DMD_mode_6.png
DMD_mode_7.png
DMD_mode_8.png

mode 7, 8を除けば,殆どがunityであり,エネルギーを保存している.

w/ noise

$\sim \mathcal{N} (0, 1)$だけのノイズを加えるとどうなるか?

eigenvals.png
DMD_mode_0.png
DMD_mode_1.png
DMD_mode_2.png
DMD_mode_3.png
DMD_mode_4.png
DMD_mode_5.png
DMD_mode_6.png
DMD_mode_7.png
DMD_mode_8.png

固有値はunit circle上に留まるはずだが(そうあってほしいが),内側に押し込まれている.エネルギーが減衰している.

piDMD: physics-informed DMD

Baddoo2023による.

w/o noise

先ずは,ノイズを入れずに.

eigenvals.png
DMD_mode_0.png
DMD_mode_1.png
DMD_mode_2.png
DMD_mode_3.png
DMD_mode_4.png
DMD_mode_5.png
DMD_mode_6.png
DMD_mode_7.png
DMD_mode_8.png

素朴なDMDと殆ど同じだが,特にmode 7-8は,より保存性が良い.

w/ noise

ノイズを入れるとどれくらい劣化するか?

eigenvals.png
DMD_mode_0.png
DMD_mode_1.png
DMD_mode_2.png
DMD_mode_3.png
DMD_mode_4.png
DMD_mode_5.png
DMD_mode_6.png
DMD_mode_7.png
DMD_mode_8.png

驚くべき結果である.

w/ stronger noise

何処まで耐えるか.

eigenvals.png
eigenvals.png
eigenvals.png
eigenvals.png

なお驚くべき結果である.但し,固有ベクトルには多大なノイズが残ったことを強調しておく.

終わりに

探索するMatrix Manifoldを積極的に制約している.

Misc

こういうのに適用するとどうなるだろうか?なんて思ったが,尾が長いだけだろう.

45_degrees.gif

0
1
0

Register as a new user and use Qiita more conveniently

  1. You get articles that match your needs
  2. You can efficiently read back useful information
  3. You can use dark theme
What you can do with signing up
0
1

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?