Question

急急急！！！MATLAB求解微分方程组，需要得到数值解

MATLAB用ode45求解微分方程组，需要编写m文件，怎样直接在命令窗口直接编写代码求解微分方程组，而不编写m文件，或者怎样把m文件展开直接写在代码里面，前提是需要得到数值解，请给一个格式实例，谢谢
能解决问题请回答，非诚勿扰，谢谢

Answer 1

对于不是很复杂的微分方程（最好不包含逻辑判断之类的），可以使用匿名函数（7.0之前使用内联函数inline）定义微分方程，然后用ode*系列函数求解。

例如，作为ode*函数文档中的例子：
function dy = vdp1000(t,y)
dy = zeros(2,1); % a column vector
dy(1) = y(2);
dy(2) = 1000*(1 - y(1)^2)*y(2) - y(1);

可以用内联函数表示为：
vdp1000=inline('[y(2); 1000*(1 - y(1)^2)*y(2) - y(1)]','t','y');
[T,Y] = ode15s(vdp1000,[0 3000],[2 0]);
plot(T,Y(:,1),'-o')

追问

还是不对，哪里出错了，我在MATLAB7.0和MATLAB2010b都试了，不行

追答

上次忘说了，对于7.0之后的版本，最好使用匿名函数（由于我当时的电脑中正在用6.5，就随手使用了inline）。

匿名函数与inline函数的最大差别是，匿名函数可以使用当前工作区中已定义的变量，而inline由于是以字符串形式定义的，它不能直接使用工作区中的变量。也正因匿名函数相比之下更方便、更强大，较新的版本中计划删除inline函数，而建议用户使用匿名函数代替。

匿名函数的使用也很简单。例如，上面的内联函数如果用匿名函数定义如下：
vdp1000=@ (t,y)[y(2); 1000*(1 - y(1)^2)*y(2) - y(1)];
注意不要把方括号给丢了。

追问

您好 非常感谢您的回答 这个问题就解决了 我想问问MATLAB代码画图直接运行结果显示正常，做成m文件运行为啥显示的图有些地方的线就变成虚线了 我觉得应该是精度的问题，编写m函数文件的采用ode45求解时可以用里面可以用options 但是这种不知道怎么改 求解 谢谢您

追答

没看明白你的意思。
最好说具体一点，例如直接运行是怎样做的，编M文件又是怎样做的，结果分别是什么样。

就我所知道的而言，直接运行与写成M文件不应该有什么不同（相同的代码），而options在命令窗口中也是可以直接使用的，并不需要写成文件。

追问

直接运行MATLAB代码做图plot生成的线是实的，转换成m文件（function ***（参数）  。。。。end）然后在命令窗口输入***（给参数赋给数值）生成的图，线条显示不清楚

如图（第一个是m文件运行的结果，第二个是直接运行的结果）

追答

你的具体做法还是没说清楚。
请把代码通过私信发给我，或者传到网盘并告诉我链接。
不然，我不知道你是怎么做的，想帮你也无从下手。

Answer 2

我可以用Mathematica试试，这个软件和Matlab有很多共同点

追问

搞不懂你给我的这是什么回答，。。。。。。。不用换软件，就用MATLAB求解，不会咱就不用谈这些无关痛痒的东西，谢谢您