# Global min optimisation for aorta 3-element windkessel model

8 ビュー (過去 30 日間)
Hussam 2024 年 3 月 28 日
コメント済み: William Rose 2024 年 4 月 1 日
So, I'm currently trying to optimise a 3-element windkessel model for the aorta using lsqcurvefit (only optimises 1 parameter) and fminsearch (okay, but very dependent on the initial parameters and does not match the pressure or flow exactly or very close). So I am thinking of implementing a global min optimiser, such as the Monte Carlo method but without the large computational cost. Any recommendations?

サインインしてコメントする。

### 採用された回答

William Rose 2024 年 3 月 28 日
Use fmincon(). I have used it to optimize a model of the circulation with up to 9 parameters, in order to match presure and flow traces recorded from the radial, carotid, and femoral arteries.
A windkessel model, by itself, cannot predict both pressure and flow. It can predict pressure from flow, or flow from pressure, but not both. To get both pressure and flow, you also must have a model of the heart itself, such as a time-varying elastance model. See, for example, this recent article by me and colleagues:
##### 9 件のコメント7 件の古いコメントを表示7 件の古いコメントを非表示
Hussam 2024 年 3 月 30 日
Hi @William Rose, thank you very much for your feedback. With number of points, I meant the number of points for multistart. Sometimes fmincon() gets stuck with certain points.
Regarding C I did put down 1e-10 m^3/Pa (where your suggested calculcation is 6e-9, not very far off).
Regarding my pressure decay, would you recommend me to use a higher order model (4 element for example), or should I try to have a pressure waveform that decays to zero. What do you think? Thank you!
William Rose 2024 年 4 月 1 日
You're welcome.
"Would you recommend me to use a higher order model (4 element for example), or should I try to have a pressure waveform that decays to zero. What do you think?"
To answer your question about how to proceed (for example, use a 4-4lement Windkessel), I would like to understand where the P and F data came from*, and what your goals are. Please email me securely by clicking on the envelope icon in the pop-up window that appears when you click on the "WR" next to my posts. I may not answer for a day or two, due to other obligations.
* I see "CFD/FSI" in a comment in your code. Did your data come from a CFD simulation involving fluid-structure interaction?

サインインしてコメントする。

### カテゴリ

Help Center および File ExchangeProgramming についてさらに検索

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by