不同initial condition ODE 求解
大家好 我是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)
→
04/14 11:24, , 1F
04/14 11:24, 1F