作者fafasisi (change)
看板MATLAB
标题不同initial condition ODE 求解
时间Sun Apr 14 02:54:00 2013
大家好 我是matlab超级新手
对它还很不熟悉
最近需要解一些微分方程式所以开始接触
请大家帮帮忙 :)
我想求不同initial condition下
二微ODE的解
写了个for loop
vi=initial condition
i=1:50
每个initial condition下照理说都会产生一组解
但每次算完 都只有出现i=50时的解
实在不知如何是好阿~~~
code 如下 请大家替我解惑!!
ode function:
function dydt=fun(t,y)
Pe=5;
r=25;
Da=12;
B=0.5;
dydt(1)=y(2);
dydt(2)=Pe*(y(2)-Da*(1+B-y(1))*exp(r*(1-1/y(1))));
dydt=dydt';
command window:
clear all
clc
vi=linspace(1,1.5,50);
for i=1:50
[t,y]=ode45(@fun,[1 0],[vi(i) 0]);
end
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 98.223.161.93
※ 编辑: fafasisi 来自: 98.223.161.93 (04/14 02:55)
1F：→ rex0707:t(i),y(i) 04/14 11:24