微分方程数值解法,知道导数的边界值,而不是函数的边界值,我知道微分方程边值问题,把微分变成差分,然后解方程组可解。但是这个问题边界值不是函数的边界值,而是函数导数的边界值
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/23 08:15:12
![微分方程数值解法,知道导数的边界值,而不是函数的边界值,我知道微分方程边值问题,把微分变成差分,然后解方程组可解。但是这个问题边界值不是函数的边界值,而是函数导数的边界值](/uploads/image/z/10426085-53-5.jpg?t=%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E6%95%B0%E5%80%BC%E8%A7%A3%E6%B3%95%2C%E7%9F%A5%E9%81%93%E5%AF%BC%E6%95%B0%E7%9A%84%E8%BE%B9%E7%95%8C%E5%80%BC%2C%E8%80%8C%E4%B8%8D%E6%98%AF%E5%87%BD%E6%95%B0%E7%9A%84%E8%BE%B9%E7%95%8C%E5%80%BC%2C%E6%88%91%E7%9F%A5%E9%81%93%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E8%BE%B9%E5%80%BC%E9%97%AE%E9%A2%98%EF%BC%8C%E6%8A%8A%E5%BE%AE%E5%88%86%E5%8F%98%E6%88%90%E5%B7%AE%E5%88%86%EF%BC%8C%E7%84%B6%E5%90%8E%E8%A7%A3%E6%96%B9%E7%A8%8B%E7%BB%84%E5%8F%AF%E8%A7%A3%E3%80%82%E4%BD%86%E6%98%AF%E8%BF%99%E4%B8%AA%E9%97%AE%E9%A2%98%E8%BE%B9%E7%95%8C%E5%80%BC%E4%B8%8D%E6%98%AF%E5%87%BD%E6%95%B0%E7%9A%84%E8%BE%B9%E7%95%8C%E5%80%BC%EF%BC%8C%E8%80%8C%E6%98%AF%E5%87%BD%E6%95%B0%E5%AF%BC%E6%95%B0%E7%9A%84%E8%BE%B9%E7%95%8C%E5%80%BC)
微分方程数值解法,知道导数的边界值,而不是函数的边界值,我知道微分方程边值问题,把微分变成差分,然后解方程组可解。但是这个问题边界值不是函数的边界值,而是函数导数的边界值
微分方程数值解法,知道导数的边界值,而不是函数的边界值,
我知道微分方程边值问题,把微分变成差分,然后解方程组可解。但是这个问题边界值不是函数的边界值,而是函数导数的边界值,请问怎么处理?看了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;
使用matlab用龙格-库塔解法求数值解,可以调用函数,也可以利用simulink做。
简单的用matlab的符号运算可以得到解析解。
有限元中一般使用伽辽金法(余量加权法的一种),简单的例子可参考http://19se.com.cn/a/gongju/2011/0511/515.html
解微分方程组对工程问题来说十分重要,正不断学习中。您好,龙格库塔是求解初值问题的,而...
全部展开
使用matlab用龙格-库塔解法求数值解,可以调用函数,也可以利用simulink做。
简单的用matlab的符号运算可以得到解析解。
有限元中一般使用伽辽金法(余量加权法的一种),简单的例子可参考http://19se.com.cn/a/gongju/2011/0511/515.html
解微分方程组对工程问题来说十分重要,正不断学习中。
收起
这个是线性非其次方程,解出来的通解一定是 齐次式通解加特解形式。
特解,可以用你的数值解法来找一个。
这样就能得到原方程的通解。通解中有两个待定常数。
把y'(a1)=b1,y'(a2)=b2代入进去,得到两个方程。可以把两个常数解出来大哥,您搞笑吗,那个要是有解析解,谁求数值解。我的问题是怎么求数值解,谢谢y''+y'+y=0这主用一般方法可以解的吧 y''+y'+y=...
全部展开
这个是线性非其次方程,解出来的通解一定是 齐次式通解加特解形式。
特解,可以用你的数值解法来找一个。
这样就能得到原方程的通解。通解中有两个待定常数。
把y'(a1)=b1,y'(a2)=b2代入进去,得到两个方程。可以把两个常数解出来
收起