Python 和matlab 关于DMD（动态模态分解）的实现和对比 21/06/08

参考资料：寇家庆，张伟伟. 动力学模态分解及其在流体动力学中的应用. 空气动力学报. 2018

 

通过numpy和matlab对于该论文中的典型算例1进行复现，并基于实数算例进行对比分析.

1 典型算例DMD方法代码和对比：

matlab:

X = [0.5,0.55-0.25i,0.48-0.55i,0.253-0.845i;1,1+0.1i,0.99+0.2i,0.97+0.299i;0.8,0.16-0.008i,0.0319-0.0032i,0.0064-0.001i];
Y = [0.55-0.25i,0.48-0.55i,0.253-0.845i,-0.1442-1.056i;1+0.1i,0.99+0.2i,0.97+0.299i,0.9401+0.396i;0.16-0.008i,0.0319-0.0032i,0.0064-0.001i,0.0013-0.0003i];
%X = [4,4,4,4; 3,9,27,81; 1,6,36,216];
%Y = [4,4,4,4; 9,27,81,243; 6,36,216,1476];

[u,s,v] = svd(X,'econ');
a1 = u' * Y *v * s^-1;
[V,D]= eig(a1);

其中，econ项将分解后的vt去除一列，从4*4矩阵变为4*3矩阵

python:

# python中特征值sigma是以列表形式出现，因此需要将其转化为对角阵

#numpy中linalg的SVD项目前似乎尚不支持econ项，因此需要手动进行去值。

# 两种工具进行SVD分解后的Vt项并不是互为转置，而是互为共轭转置。 因此，在某些情况下，U项