三阶常微分方程边值问题 用有限差分法和打靶法的matlab 程序有么?或者思路
t^2x'''+2tx''+t^2x'+2tx=6t,[0.001,3]区间,x(0.001)=5,x'(0.001)=-2,x(3)=1.373这个问题。。ode45可...
t^2x'''+2tx''+t^2x'+2tx=6t,[0.001,3]区间,x(0.001)=5,x'(0.001)=-2,x(3)=1.373
这个问题。。ode45可以做了,有限差分法和打靶法的程序有么? 展开
这个问题。。ode45可以做了,有限差分法和打靶法的程序有么? 展开
2个回答
展开全部
ode45只能解初值问题,不能解边值问题
x(0.001)=5,x'(0.001)=-2,x(3)=1.373
这是个边值问题
你可以用matlab bvp4c解或者打靶法,
具体的就是猜一个x=0.001的二阶导数值x''(0.001)=a,用ode45求解,求出x(3)的值与1.373比较,如果误差不允许,就在迭代变换a值,直到解出的x(3)的值与1.373误差可接受为止
x(0.001)=5,x'(0.001)=-2,x(3)=1.373
这是个边值问题
你可以用matlab bvp4c解或者打靶法,
具体的就是猜一个x=0.001的二阶导数值x''(0.001)=a,用ode45求解,求出x(3)的值与1.373比较,如果误差不允许,就在迭代变换a值,直到解出的x(3)的值与1.373误差可接受为止
更多追问追答
追问
我的那道题是这样的:1.用打靶法求解2.用ode45求解3.用有限差分法求解。。并比较三者结果。
打靶法和有限差分法求解三阶的微分方程 完全找不到思路。。。网上都是二阶的。诶。
追答
这个题目有问题,ode45只能解初值问题,你若不行用ode45解试试,保证和你给条件对不上。
打靶法我已经跟你说思路了
有限差分法很简单,主要是三阶差分表示形式问题,和边界条件的处理
三阶差分形式,采用了迎风格式和五点格式,这个你可以查下,或者由二阶的推导
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询