File Exchange

image thumbnail

POD_DMD

version 1.2.0.1 (4.67 MB) by Jonathan Higham
A tool to connect POD and DMD

24 Downloads

Updated 18 Jan 2019

GitHub view license on GitHub

As featured in the Journal of Hydraulic Research Paper "Implications of the selection of a particular modal decomposition technique for the analysis of shallow flows". This algorithm connects the spatially orthogonal Proper Orthogonal Decomposition with the temporally orthogonal Dynamic Mode Decomposition. From flow visualisation images, PIV or CFD vector fields it is possible to extract coherent structures which are quasi-spatially and temporally orthogonal.
Two example cases are included with the script. Please cite:
Higham, J.E., Brevis, W & Keylock, C.J. (2018) - Implications of the selection of a particular modal decomposition technique for the analysis of shallow flows , Journal of Hydraulic Research, DOI:10.1080/00221686.2017.1419990
https://www.researchgate.net/publication/322722097_Implications_of_the_selection_of_a_particular_modal_decomposition_technique_for_the_analysis_of_shallow_flows

Cite As

Jonathan Higham (2020). POD_DMD (https://github.com/jonnyhigham/POD_DMD), GitHub. Retrieved .

Comments and Ratings (9)

Hello thank you for the great work. Can we apply this method on 1-D/time series signals?

Aryan Singh

I got it, Jonathan, I totally missed unzipping the data_cylinder zip file.

Thanks

Ali sameer

FC Wu

Hi Jonathan, Do you have any tutorial to work with this program? I appreciate you if you can prepare a tutorial on youtube or anywhere.

Best,
Omid

Dipanjan

Thanks Yonatan, I have changed the script as per your suggestion. There is no change in functionality. It was just me being lazy and not checking the manual for findpeaks better!

I think there is a bug in the script for Matlab 2013: X is the vector index of the peaks (not a frequency value), and x is a frequency value.

Original Modified
[Y, X]=findpeaks(px,fx); [Y, X]=findpeaks(px);
[v, l]=min(abs(X-x)); [v, l]=min(abs(fx(X)-x));
hold on
plot(X(l),Y(l),'ro') plot(fx(X(l)),Y(l),'ro')

% Find the matching frequency from POD
[~,l]=min(abs(X(l)-f)); [~,l]=min(abs(fx(X(l))-f));

Updates

1.2.0.1

Added to github

1.2.0.0

Check for 2018a - works fine.

1.2.0.0

Change to description and code edit. No change in functionality.

1.2.0.0

Small change to script. For this case no change in functionality.

1.1.0.0

Paper Link

1.1.0.0

Change to syntax as suggested by Yonatan Cadavid. There is no affect to the functionality.

1.0.0.0

Change to text

1.0.0.0

Added reference.

MATLAB Release Compatibility
Created with R2018b
Compatible with any release
Platform Compatibility
Windows macOS Linux