matlab 微分方程边值问题.docVIP

我知道微分方程边值问题，把微分变成差分，然后解方程组可解。但是这个问题边界值不是函数的边界值，而是函数导数的边界值，请问怎么处理？看了2本书，也没答案 比如y+y+y=f(x) a1=x=a2 y(a1)=b1,y(a2)=b2 2011-5-14 08:23 最佳答案 这个方程最简单的方法是用shooting method.可wiki 之。 主要就假设一个y(a1), 然后看，y(a2)是否等于b2。 可以用matlab写两个小函数，一个用在ode45（fun1），然后y(a1)的值可以用fzero(fun2). fun1就是把二阶微分化成一阶方程组，fun2就是调用ode45并返回y(a2)-b2. 这样会自动求出初值的。 别外这个方程可能有解析解。可以利用格林函数来解，得到一个积分形式的解。 给你看一个我以前写的例子吧。适当改下就行了。 function out=fdif(eta,f) out=zeros(size(f)); out(1)=f(2); out(2)=f(3); out(3)=f(2)^2-1-f(1)*f(3); function result=soldif(x) [etaout,fout]=ode45(fdif,[0,10],[0,0,x]); result=fout(end,2)-1; 追问 您好，最近看了一下shooting method，发现我这个问题似乎不能应用此办法，我的方程齐次项是y+9/4y=0，shooting method是把原方程分解为y=u+c*v，c是常数，v就是齐次方程。x的区间是[-pi,pi]，这就导致了一个问题，因为shooting method要求v‘（-pi）=0，所以v’（pi）=0，c无法求解。我看的shooting method例子边界条件是y=。。。而不是y‘=。。。请问我该如何处理，两边求导？ 回答 你对shooting method的理解好像不太对(或者说和我说的不一样)。对于齐次方程，比如你说的y+9/4y=f(x)。应该像我之前说的先变成一次微分方程组。y1=y2; (y2)+9/4y1=f(x)所以下面函数改成 function out=fdif(eta,f) out=zeros(size(f)); out(1)=f(2); out(2)=f(x)-9/4*f(1); //注意这个f(x)改成相应的，如f(x)=x^2; then 改成eta^2 这样就可以调用ode45. 把下面代码改成 function result=soldif(x) [etaout,fout]=ode45(fdif,[-pi,pi],[x,0]); // [x,0]分别为假设初如值及其一阶导数 result=fout(end,2)-0; // 这里的0是最后一点导数。 最后利用fzero就可以求解，试试吧。Good luck 追问 您好，解二阶初值问题微分方程的方法我知道，也就是ode45之前我都知道。但是我看书上和wiki百科，关于线性问题，y=y1(t)+(y1-y1(t))/2*y2(t)，（篇幅太长，您可以自己上wiki上看看）。但是我这个问题y2(t)=0，所以我就不知道如何处理了。wiki百科上第一部分说的和您一样，是要解方程的根，但是我不知道是解什么方程的跟，他上边写的是:F(a)=y(t1;a)-y1,其中y(t0)=y0,y(t1)=y1,假设初值是y(t0)=y0,y(t0)=a。 回答 wiki那个问题是已经知道两个端点的值。所以先假设一个端点的值y0,一阶导数值为a,这样求解程后，会得到另一个端点的值，和已知的端点值y1比较，直到差为0.我上面替你改的程序可以直接用的。只要改下f(x)。为什么这样我大体说了。fzero就是求那个F(a)的。 追问 您还是没明白我的意思，我取了个特殊情况，能求出解析解，发现解里确实有个不确定常数，所以shooting method里的(y1-y1(t))/2*y2(t)确实求不出来，换成常数就能和解析解对应上。THX，anyway