utopai 发表于 2006-4-18 19:57

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

<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 08:44

回复:(utopai)[求助]这种形式的微分方程组该怎么解...

这个怎么看着像是常微分方程?另外你写的两个边界条件有什么不同吗
[此贴子已经被作者于2006-4-19 8:45:12编辑过]

utopai 发表于 2006-4-19 09:56

回复:(suffer)

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

happy 发表于 2006-4-20 19:05

回复:(utopai)[求助]这种形式的微分方程组该怎么解...

<P>用ode45好像求解不了<BR>试试用差分法求解</P>
页: [1]
查看完整版本: [求助]这种形式的微分方程组该怎么解算啊?