您正在关注此问题

您将在关注的内容源中看到更新。
您可能会收到电子邮件，具体取决于您的通信预设项。

how to give initial condition as array in ode45

4 次查看（过去 30 天）

显示更早的评论

abhijit kulkarni 2014-4-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45

⋮

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45

评论： Jan 2017-2-20

采纳的回答： Star Strider

I want to enter initial condition of 6x6 matrix in ode45.

How to do that? Any idea.

Please help

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_430198

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_430198

Sir, how do i solve an ode with an array , using ode45? ODE is udY/dx=P, Y,P are arrays.

Jan 2017-2-20

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_430204

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_430204

@DIP: Please do not attach a new question to an existing thread. Thanks.

请先登录，再进行评论。

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

采纳的回答

Star Strider 2014-4-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#answer_134964

⋮

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#answer_134964

The ode45 initial conditions have to be a vector.

You have my interest, though. Why do you want them to be a matrix?

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210869

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210869

In my program, I have one differential equation of say : Hdot= something.

Now size od Hdot is 6x6

So i have to give initial condition of 6x6

abhijit kulkarni 2014-4-29

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210870

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210870

I want use it in Kalman Filter equation.

Star Strider 2014-4-29

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210871

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210871

编辑：Star Strider 2014-4-29

No, you don’t need a matrix of initial conditions.

You need a (6x1) vector of initial conditions, since you have a system of 6 differential equations, each an expression for the derivative of one of your variables.

The vector should work.

EDIT — I’ll have to go back to my Kalman filter books to understand what you might be doing (too tired tonight), but what I mentioned about the initial condition vector for ode45 is correct.

abhijit kulkarni 2014-4-29

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210873

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210873

If I send you the code is, it possible to rectify the error?

Please help me.

abhijit kulkarni 2014-4-29

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210874

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210874

Pdot = F * P + P * F' + Q - P * H' *inv® * H * P;

Here P dot is a 6x6 matrix. Yes it is diagonal matrix. Still I dont know how to enter the initial condition.

abhijit kulkarni 2014-4-29

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210875

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210875

Ok take rest.

Thanks for immediate reply. Good night

Star Strider 2014-4-29

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210876

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210876

Don’t send it to me, instead post it as an attachment to your original Question (use the ‘paper-clip’ icon in the toolbar above the window). That way, if I can’t answer it, there are others here on MATLAB Answers who probably can, and will be able to read and run your code. If you send it only to me, others will not have access to it.

Star Strider 2014-4-29

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210879

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210879

在 MATLAB Online 中打开

I don’t know exactly what the initial conditions for your system are, but a vector of zeros will likely work unless you know they will not. If you have a function file or anonymous function for Pdot, call ode45 this way:

% For a function file for |Pdot| named |Pdot.m|:
[T, P] = ode45(@Pdot, tspan, zeros(6,1));
% For an anonymous function such as: 
F = [matrix];
P = [matrix];
Q = [matrix];
R = [matrix];
H = [matrix];
Pdot = (t,p,F,P,Q,R,H) F * P + P * F' + Q - P * H' * R \ H * P;
[T, P] = ode45(@(t,p) Pdot(t,p,F,P,Q,R,H), tspan, zeros(6,1));

You will have to define F, P, Q, R, H. I cannot test this code because I do not have all your information. However, it should work. (Note that I replaced inv( R ) * H by R \ H. That is more efficient and will produce a more accurate result.)

abhijit kulkarni 2014-4-29

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210913

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210913

Respected Star Thank you very much!

I thought I took much time of your rest.

Well I have Pdot equation as one of the two more differential equations.

What I did is that as you said : I put initial condition of 6 variables (diagonal) of P initial matrix. And then inside program, I equated them in for loop with eye(y).

Any way your help and attitude to help is simply great.

In gratitude Abhijit

Star Strider 2014-4-29

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210942

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127554-how-to-give-initial-condition-as-array-in-ode45#comment_210942

My pleasure!

I am happy it worked.

请先登录，再进行评论。

另请参阅

类别

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

选择网站

选择网站以获取翻译的可用内容，以及查看当地活动和优惠。根据您的位置，我们建议您选择：。

(English)
(Deutsch)
(Français)

（简体中文）
(English)

您也可以从以下列表中选择网站：

如何获得最佳网站性能

选择中国网站（中文或英文）以获得最佳网站性能。其他 MathWorks 国家/地区网站并未针对您所在位置的访问进行优化。

美洲

América Latina (Español)
Canada (English)
United States (English)

欧洲

Belgium (English)
Denmark (English)
Deutschland (Deutsch)
España (Español)
Finland (English)
France (Français)
Ireland (English)
Italia (Italiano)
Luxembourg (English)

Netherlands (English)
Norway (English)
Österreich (Deutsch)
Portugal (English)
Sweden (English)
Switzerland
United Kingdom (English)

亚太

Australia (English)
India (English)
New Zealand (English)
中国
- 简体中文Chinese
- English
日本Japanese (日本語)
한국Korean (한국어)

联系您当地的办事处

how to give initial condition as array in ode45

2 个评论
显示无隐藏无

采纳的回答

10 个评论
显示 8更早的评论隐藏 8更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

how to give initial condition as array in ode45

2 个评论 显示 无隐藏 无

采纳的回答

10 个评论 显示 8更早的评论隐藏 8更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

WeChat

2 个评论
显示无隐藏无

10 个评论
显示 8更早的评论隐藏 8更早的评论