搜搜考试网微分方程数值解法,知道导数的边界值,而不是函数的边界值,
来源：学生作业帮 编辑：搜搜考试网作业帮 分类：数学作业 时间：2024/05/02 07:16:15
微分方程数值解法,知道导数的边界值,而不是函数的边界值,
我知道微分方程边值问题，把微分变成差分，然后解方程组可解。但是这个问题边界值不是函数的边界值，而是函数导数的边界值，请问怎么处理？看了2本书，
比如y''+y'+y=f(x) a1
我知道微分方程边值问题，把微分变成差分，然后解方程组可解。但是这个问题边界值不是函数的边界值，而是函数导数的边界值，请问怎么处理？看了2本书，
比如y''+y'+y=f(x) a1
这个方程最简单的方法是用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
主要就假设一个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