tanai
2009-03-21, 19:16
初次遇到这样的方程，困扰很久，不知道如何求解，希望各位高手赐教，谢谢！：）
da/dt=da^2/dx^2+f(a,u,j) (1) 其中f(a,u,j)=-a+exp(a)-exp(-j+a)+12.5*u
dj/dt=-j-10*(j+exp(a)-1) (2)
du/dt=-u-2*(u-a) (3)
边界条件：
da(0,t)/dx=da(L,t)/dx=0 变量x的取值范围[0,L]
我现在只会用pdepe函数解方程(1), 但加了 (2)， (3)两个约束方程后，就不知道如何编程了。
da/dt=da^2/dx^2+f(a,u,j) (1) 其中f(a,u,j)=-a+exp(a)-exp(-j+a)+12.5*u
dj/dt=-j-10*(j+exp(a)-1) (2)
du/dt=-u-2*(u-a) (3)
边界条件：
da(0,t)/dx=da(L,t)/dx=0 变量x的取值范围[0,L]
我现在只会用pdepe函数解方程(1), 但加了 (2)， (3)两个约束方程后，就不知道如何编程了。