Trying to plot 3-d

Question

Avan Al-Saffar 2014-8-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/146817-trying-to-plot-3-d

评论： Yu Jiang 2014-9-3

My equation is

dxdt=N0*sin(omega*t)*x*(1-x/K);
N0=1:1:10;
omega=1:1:10;

I want to plot 3-D for N0,omega,Xmax and Xmean.

I am trying with the following code but it does not work.

function RunOsciliationsky3D
N0all= 1:1:10;
N=length(N0all);
omegaall= 1:1:10;
M=length(omegaall);
Pmax=zeros(N,M);
Pmean=zeros(N,M);
Pall=[Pmax,Pmean];
x=size(Pall);
for i=1:N
  for j =1:N
      [t,x]=ode45(@osciliation,[0 100],0.1,[],N0all(i),10,omegaall(j));
      Pall(i,j)=x;
  end
end
[N0x,omegay]=meshgrid(N0all,omegaall);
h=mesh(N0x,omegay,Pall);

1;

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Andy L 2014-8-8

What error are you getting? I have tried copying and pasting but the function in your optimisation function isn't included.

请先登录，再进行评论。

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

Answer 1

Yu Jiang 2014-8-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/146817-trying-to-plot-3-d#answer_147971

See the files I revised as attached. It can generate Pmean. Similarly you can plot Pmax. Notice that it is not clear what K is, so I let it be 0.15. I also force N0 to take negative values since positive values make the system unstable.

Hope it helps.

- Yu

2 个评论
显示无隐藏无

Avan Al-Saffar 2014-9-1

在 MATLAB Online 中打开

Dear Yu Jiang

Thank you for your previous help, but I need more help please if you have any advice?

I've chosen K=10 and N0 to get a positive values. I got the result which I am not happy with it and I got this error message:

Warning: Failure at t=8.083150e+01. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (2.273737e-13)at time t.

> In ode45 at 308

In Runanswertest at 13

Warning: Failure at t=8.701882e+01. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (2.273737e-13) at time t.

> In ode45 at 308

In Runanswertest at 13

Warning: Failure at t=3.011651e+01. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (5.684342e-14) at time t.

> In ode45 at 308

In Runanswertest at 13

Warning: Failure at t=5.698339e+00. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.421085e-14)at time t.

> In ode45 at 308

In Runanswertest at 13

Note: I thought to remove the initial transient period or just to reduce the time step but I am really confused about this code .

Yu Jiang 2014-9-3