如何解二维常微方程组？

-0.0002t+0.0037=(0.0891685780439393*2*x+0.938619232671706*y)*(dx/dt)
-0.0004t+0.0089=(0.0650571169521575*2*y+0.938619232671706*x)*(dy/dt)
如何解出x和y关于t的表达式，已知x(0)=0.3327,y(0)=0.0497;望大虾们不吝赐教！

可以用matlab解。

求大侠具体算法！
>> syms x y t
>> [x,y]=dsolve('-0.0002t+0.0037=(0.0891685780439393*2*x+0.938619232671706*y)*(dx/dt)','-0.0004t+0.0089=(0.0650571169521575*2*y+0.938619232671706*x)*(dy/dt)','x(0)=0.3327','y(0)=0.0497')
显示出错，望大侠指教！

我也不知道怎么解。但我用系统动力学软件却能得到x和y随时间的图形

4# nkunku
麻烦楼上的大侠，你那个图是怎么画出来的？既然能画出来，x和y肯定能表示成t的表达式吧？能写出来吗？

李宏海 发表于 2011-3-23 09:48
4# nkunku
麻烦楼上的大侠，你那个图是怎么画出来的？既然能画出来，x和y肯定能表示成t的表达式吧？能写出来吗？

虽然能画出图像，但未必就能得到解析解。我觉得它是利用数值解画出来的。

另外，我用Mathematica和Maple也没能得到解析解，但得到数值解并画出图像，确实是都没问题的。

哎，我觉得自己还是水平有限。不好意思。

“另外，我用Mathematica和Maple也没能得到解析解，但得到数值解并画出图像，确实是都没问题的。”{:3_46:}
麻烦大侠把两个程序都给我看看吧，太谦虚了，呵呵！

In[6]:= DSolve[{x'[t] ==(0.0037 - 0.0002 t)/(0.178337 x[t] + 0.938619 y[t]),y'[t] ==(0.0089 - 0.0004 t)/(0.938619 x[t] + 0.130114 y[t])},{x[0]==0.3327,y[0]==0.049},{x[t],y[t]},t]
楼上的大侠，这是我用Mathematic编写的算法，但是解不出来呀，望大侠不吝赐教！

见以下图形

感谢了，呵呵{:3_59:}