[求助]这种形式的微分方程组该怎么解算啊？

utopai

<P>请问这样的微分方程组这么解算啊？</P>
<P>
<P><FONT face="Times New Roman">  </FONT></P><FONT face="Times New Roman">dy1/dx</FONT>＝<FONT face="Times New Roman">f (x, y1, y2)</FONT>；<BR><FONT face="Times New Roman">dy2/dx</FONT>＝<FONT face="Times New Roman">g (x, y1, y2)</FONT>；
<br>
<P>    区间为（<FONT face="Times New Roman">0, L</FONT>）；</P>
<P>    边界条件为：<BR><FONT face="Times New Roman">  x</FONT>＝<FONT face="Times New Roman">0</FONT>时，<FONT face="Times New Roman">y1</FONT>＝<FONT face="Times New Roman">a </FONT>；<BR><FONT face="Times New Roman">  x</FONT>＝<FONT face="Times New Roman">L</FONT>时，<FONT face="Times New Roman">y2</FONT>＝<FONT face="Times New Roman">b </FONT>；</P>
<P>
<P><FONT face="Times New Roman"></FONT></P>
<p>
<P>我查了很多的资料，一般资料上的算例所给的边界条件为都是（<FONT face="Times New Roman">x</FONT>＝<FONT face="Times New Roman">0</FONT>时，<FONT face="Times New Roman">y1</FONT>＝<FONT face="Times New Roman">a </FONT>；<FONT face="Times New Roman">x</FONT>＝<FONT face="Times New Roman">0</FONT>时，<FONT face="Times New Roman">y2</FONT>＝<FONT face="Times New Roman">b </FONT>）。而对于边界条件为（<FONT face="Times New Roman">x</FONT>＝<FONT face="Times New Roman">0</FONT>时，<FONT face="Times New Roman">y1</FONT>＝<FONT face="Times New Roman">a </FONT>；<FONT face="Times New Roman">x</FONT>＝<FONT face="Times New Roman">L</FONT>时，<FONT face="Times New Roman">y2</FONT>＝<FONT face="Times New Roman">b </FONT>）时的基本没有。希望那位知道的朋友能够指点一下，先谢谢了啊！</P>

suffer

这个怎么看着像是常微分方程？另外你写的两个边界条件有什么不同吗

[此贴子已经被作者于2006-4-19 8:45:12编辑过]

utopai

是常微分方程组。两个边界条件是不一样的啊，一个是给的都是x＝0的边界值，另一个则分别给的x＝0和x＝L时的边界值，是不一样的。我看了一些资料说这种常微分方程组可以用ode45解算。但是用ode45的形式只有这样：<BR>    y0=[a,b];<BR>    [x,y]=ode45(fun,[0:0.02L:L],y0);<BR>上边的形式只能解算都是x＝0的边界值的情况，却不能解算x＝0和x＝L时的边界值的情况。我现在就是不知道怎么解算当边界条件分别为x＝0和x＝L时的常微分方程组。

happy

<P>用ode45好像求解不了<BR>试试用差分法求解</P>