[問題] 一題關於BVP-ODE的問題已回收
我需要解四個連立的BVP-ODE方程式
dy(1)/dx= -x/2*(k1*y(1)-k2*y(3))
dy(2)/dx= x/2*(k1*y(3)-k2*y(4))
dy(3)/dx= -y(1)/D+y(3)/Ct*(y(1)+y(2))
dy(4)/dx= -y(2)/D+y(4)/Ct*(y(1)+y(2))
y1(0)=y2(0)=0 y3(1)=1 y4(1)=0
其中
Ct=1;
D=10^(-5);
k1=1;
k2=0.1;
每次執行完 就會出現
??? Input argument "yc" is undefined.
Error in ==> fourbc at 2
res=[ya(1)
Error in ==> bvp4c at 195
testBC = feval(bc,y(:,1),y(:,end),ExtraArgs{:});
請問是哪裡出現了問題
我的code如下
solinit=bvpinit(linspace(0,1,4),[1,0]);
sol=bvp4c(@fourode,@fourbc,solinit);
------------------------------------
function res=fourbc(ya,yb,yc,yd)
res=[ya(1)
yb(1)
yc(1)-1
yd(1)];
------------------------------------
function dydx=fourode(x,y)
Ct=1;
D=10^(-5);
k1=1;
k2=0.1;
dydx=[ -x/2*(k1*y(1)-k2*y(3))
x/2*(k1*y(3)-k2*y(4))
-y(1)/D+y(3)/Ct*(y(1)+y(2))
-y(2)/D+y(4)/Ct*(y(1)+y(2))];
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 216.47.139.238
推
03/16 00:54, , 1F
03/16 00:54, 1F
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
