| 查看: 2409 | 回复: 5 | ||
[求助]
Matlab数值求解非线性常微分方程 已有1人参与
|
|
现在遇到一个问题:求解非线性微分方程两点边值问题的数值解。如何用Matlab求解？（对matlab不熟悉），文中说采用差分方法将其离散化，然后用Newton迭代求解。可是如何用差分方法离散？如何迭代求解？请大侠指教.... 方程为: d/dr(dy/dr)=d/dr(secy)*2/A-B*d/dr(tany/r) (其中，A和B都是常数) 如下图所示： 边值条件是:当r=0时，y=0 当r=5时，y=0.2  20130831103501088.jpg |
回复此楼» 猜你喜欢
filecode
已经有10人回复
-
听说今天filecode变了
已经有3人回复
-
filecode
已经有12人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有5人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有11人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有30人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教一个matlab求解非线性方程组的问题
已经有9人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- Matlab数值积分结果是虚数
已经有6人回复
- matlab求解常微分方程
已经有6人回复
- 偏微分方程数值计算书籍推荐
已经有10人回复
- 二阶非线性微分方程求解
已经有14人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 求助一阶非线性微分方程的解
已经有7人回复
- 各位朋友，如何用matlab编写这个方程，并求解？
已经有15人回复
- 一阶微分方程求解，并优化参数值
已经有14人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- 隐式欧拉法求解一阶常微分方程
已经有7人回复
- matlab求解非线性方程组
已经有16人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- MATLAB解方程出现??? Subscripted assignment dimension mismatch.
已经有3人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】求解隐函数
已经有6人回复
- 【求助】matlab 解矩阵方程 重金悬赏解决方案
已经有10人回复
- 【讨论】CFD软件是求解PDE的，matlab中也可以求解PDE
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】求教这种微分方程的详细求解过程
已经有14人回复
2楼2013-09-01 15:10:59
z65668396
金虫 (正式写手)
- 应助: 14 (小学生)
- 金币: 1120.8
- 散金: 1682
- 红花: 3
- 帖子: 469
- 在线: 369.6小时
- 虫号: 1774509
- 注册: 2012-04-23
- 性别: GG
- 专业: 机械动力学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
1592203609: 金币+2, 谢谢回帖 2013-09-05 16:00:29
815292578: 金币+10, ★有帮助, 感谢你的回复！ 2013-09-05 16:43:55
1592203609: 金币+2, 谢谢回帖 2013-09-05 16:00:29
815292578: 金币+10, ★有帮助, 感谢你的回复！ 2013-09-05 16:43:55
|
我最近也在用牛顿迭代解一个非线性方程组 Newton迭代有程序 我拷给你你 编写Newton迭代法的M文件NewtonIterate.如下： function s=NewtonIterate(x,eps) %Newton迭代法求非线性方程组 %x为迭代初值，eps为允许误差值 if nargin==1 eps=1.0e-6; elseif nargin<1 error return end x1=fx1(x);%非线性方程组 x2=-dfx1(x);%非线性方程组导数 x3=inv(x2); x0=x3*x1'; while norm(x0)>=eps %循环迭代 x=x0'+x; x1=fx1(x); x2=-dfx1(x); x3=inv(x2); x0=x3*x1'; end s=x0'+x; return 例： x1.^2-10x1+x2.^2+b=0 x1*x2.^2+x1-10x2+8=0 用MATLAB实现，编写非线性方程组的M文件fx1.m如下所示： function y=fx1(x) y(1)=x(1)*x(1)-10*x(1)+x(2)*x(2)+8; y(2)=x(1)*x(2)*x(2)+x(1)-10*x(2)+8; y=[y(1) y(2)]; 编写非线性方程组导数的M文件dfx1.m function y=dfx1(x) y(1)=2*x(1)-10; y(2)=2*x(2); y(3)=x(2)*x(2)+1; y(4)=2*x(1)*x(2)-10; y=[y(1) y(2);y(3) y(4)]; |
3楼2013-09-05 11:34:17
4楼2013-09-05 16:55:26
z65668396
金虫 (正式写手)
- 应助: 14 (小学生)
- 金币: 1120.8
- 散金: 1682
- 红花: 3
- 帖子: 469
- 在线: 369.6小时
- 虫号: 1774509
- 注册: 2012-04-23
- 性别: GG
- 专业: 机械动力学
5楼2013-09-05 20:41:43
CSU_HUSTer
新虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2705.1
- 红花: 3
- 帖子: 713
- 在线: 169.4小时
- 虫号: 2836351
- 注册: 2013-11-29
- 性别: GG
- 专业: 工程热物理相关交叉领域
6楼2014-03-10 18:19:27
50