time for parameter estimation using ode 45

1 次查看（过去 30 天）

nado 2022-10-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1830513-time-for-parameter-estimation-using-ode-45

回答： Rajiv Singh 2022-11-7

Hello everyone,

I used lsqnonlin to find unknown parameters of an ODE. However, when i run the script, it is taking a long time. Until now, it tooks 25 min and it's the first time i encounter this !

The system i'm solvin consists of 3x3 matrix.

Do you think the time for this is reasonable or there is something going wrong ???

I searched a lot on how to know whether the script is stuck somewhere but i didn't succeed.

Based on your experience, is there any method to know where the code is stuck ?

11 个评论
显示 9更早的评论隐藏 9更早的评论

nado 2022-10-19

[Ca1,resnorm,residual]=lsqnonlin(@(Ca1)odefit_calib_damping(Ca1),1e+5)

%% lsqnonlin is the function used to optimize the error between experiment and numerical model

function [ err ] = odefit_calib_damping(Ca1 )

time=0:0.1:630;

[t,y]=ode45(@(t,y)calib_damping(t,y,Ca1),time,[10;0;0;0;0;0]);

err=abs(y_exp_surge-y(:,1)).^2; %% y_exp is a vector containing data from experiment

end

function [dydt] = calib_damping( t,y,Ca1)

% this is the function where the system of coupled ODEs are implemented

M(1,1)=8.3163e+7;

M(1,2)=-2.6982e+9;

M(1,3)=1.8107e+6;

M(2,2)=5.6542e+11;

M(2,3)=3.4131e+8;

M(3,3)=1.4249e+6;

M(2,1)=M(1,2);

M(3,1)=M(1,3);

M(3,2)=M(2,3);

K(1,1)=687024.52188410972265443647316524;

K(1,2)=-4.861e+7;

K(1,3)=0;

K(2,1)=-4.861e+7;

K(2,2)=1.4604e+10-1467627741.2664321859450536919855;

K(2,3)=0;

K(3,1)=0;

K(3,2)=0;

K(3,3)=1.4249e+6;

C(1,1)=1.8448e+5 + Ca1;

C(2,2)=1.8286e+10;

C(3,3)=4.0008e+10;

C(1,2)=5.8081e+7;

C(1,3)=-8.5911e+7;

C(2,3)=-2.7048e+10;

C(2,1)=C(1,2);

C(3,1)=C(1,3);

C(3,2)=C(2,3);

M_K=-(M^(-1))*K;

M_C=-(M^(-1))*C;

dydt(1)= y(4);

dydt(2)=y(5);

dydt(3)=y(6);

dydt(4)= M_K(1,1)*y(1)+M_K(1,2)*y(2)+M_K(1,3)*y(3)+M_C(1,1)*y(4)+M_C(1,2)*y(5)+M_C(1,3)*y(6);

dydt(5)= M_K(2,1)*y(1)+M_K(2,2)*y(2)+M_K(2,3)*y(3)+M_C(2,1)*y(4)+M_C(2,2)*y(5)+M_C(2,3)*y(6);

dydt(6)= M_K(3,1)*y(1)+M_K(3,2)*y(2)+M_K(3,3)*y(3)+M_C(3,1)*y(4)+M_C(3,2)*y(5)+M_C(3,3)*y(6);

dydt=[dydt(1);dydt(2);dydt(3);dydt(4);dydt(5);dydt(6)];

end

nado 2022-10-19

Thank you @Sam Chak

nado 2022-10-19

@Torsten

I tried ode15s and it works !!!

Thanks a lot !!!

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Rajiv Singh 2022-11-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1830513-time-for-parameter-estimation-using-ode-45#answer_1093193

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

类别

MATLAB Programming

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by