How to get the output for the following code into the workspace?

Question

Urvi 2012-10-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/51293-how-to-get-the-output-for-the-following-code-into-the-workspace

My code is as follows:

function H = newmain 
global b1 b2 b3;
b1=0;
b2=2.0;
b3=0.33;
options=odeset('InitialStep', 0.01, 'MaxStep', 0.01, 'RelTol', 10., 'AbsTol',10.);
[t2,y2]=ode45(@equation,[0:0.05:0.1],[1 2 0], options);
H = [t2 y2]
end
function dy=equation(t,y)
global b1 b2 b3 
dy=zeros(3,1);
b1=b1+1-exp(-b3)
b2=b1-y(1)
b3=b1+b2+y(2)
dy(1)=-b1*y(1);
dy(2)=b3*y(1)+b2*y(2);
dy(3)=sqrt(b1)+y(1)+y(3);
end

How can directly get the values of b1 b2 b3 y1 y2 y3 at the same time instant in the workspace?

Thanks!

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

请先登录，再进行评论。

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

Answer 1

Azzi Abdelmalek 2012-10-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/51293-how-to-get-the-output-for-the-following-code-into-the-workspace#answer_62533

在 MATLAB Online 中打开

function dy=equation(t,y)
global c b1 b2 b3 y1 y2 y3
y1=y(1);y2=y(2);y3=y(3);
dy=zeros(3,1);
b1=b1+1-exp(-b3)
b2=b1-y(1)
b3=b1+b2+y(2)
dy(1)=-b1*y(1);
dy(2)=b3*y(1)+b2*y(2);
dy(3)=sqrt(b1)+y(1)+y(3);
c=[c; b1 b2 b3 y1 y2 y3]

add in your main program

global c b1 b2 b3 y1 y2 y3
y1=0;y2=0;y3=0;c=[];

2 个评论
显示无隐藏无

Urvi 2012-10-19

why have you given y1=0 y2=0 y3=0? I need them from my odes.

Azzi Abdelmalek 2012-10-19

it's just an initialization, not effect

请先登录，再进行评论。

How to get the output for the following code into the workspace?

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

回答（1 个）

2 个评论
显示无隐藏无

另请参阅

类别

标签

Community Treasure Hunt

How to get the output for the following code into the workspace?

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

回答（1 个）

2 个评论 显示 无隐藏 无

另请参阅

类别

标签

Community Treasure Hunt

WeChat

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

2 个评论
显示无隐藏无