| 查看: 1086 | 回复: 5 | |||
chuzhenhui金虫 (小有名气)
|
[交流]
【求助】用fsolve求解非线性一元方程 已有2人参与
|
|
我用fsolve（）求解非线性方程，程序如下： >> t=@(D)(8/pi^2)*exp(-12000000*pi^2*D)-0.092034; >> [D,f,flag,out]=fsolve(t,0) 运行后结果如下： Optimization terminated: norm of relative change in X is less than max(options.TolX^2,eps) and sum-of-squares of function values is less than sqrt(options.TolFun). D = 1.8370e-008 f = -8.6133e-006 flag = 2 out = iterations: 22 funcCount: 36 algorithm: 'trust-region dogleg' firstorderopt: 44.0855 message: [1x167 char] 请问红色部分是什么意思？我的这个解好像不是很准确，怎么改进啊？ 要用optmiset设置参数吗？ 大侠指导下，感激不尽！！ |
回复此楼» 猜你喜欢
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何使用matlab求解非线性方程组的所有整数解？
已经有9人回复
- matlab求解非线性方程组
已经有16人回复
- 非线性方程组matlab拟合
已经有9人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- MATLAB非线性方程组求解~~~反应釜模拟
已经有9人回复
- 【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！
已经有5人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】MATLAB 如何解矩阵方程（达人救命）
已经有6人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
怎么成了非线性方程了啊?
2楼2010-05-13 19:34:56
chuzhenhui
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 600.7
- 散金: 35
- 帖子: 129
- 在线: 10.2小时
- 虫号: 763256
- 注册: 2009-05-05
- 性别: GG
- 专业: 食品&包装
★
zzuwangshilei:给个红包，欢迎积极关注自己的帖子 2010-05-14 08:11:52
nono2009(金币+1):红包补上。 2010-05-14 13:06:56
zzuwangshilei:给个红包，欢迎积极关注自己的帖子 2010-05-14 08:11:52
nono2009(金币+1):红包补上。 2010-05-14 13:06:56
|
谢谢，不好意思我粗心了，其实我想解的方程是有好多项加起来的，我图方便拿了其中的一项，具体的方程是 s=0 for n=0:inf s=s+(-8/pi^2)*(1/(2*n+1)^2)*exp(-(2*n+1)^2*pi^2*12000000*D) end 要解的方程是由上面得到的 s-0.907966=0 我拿了其中几项用同样的方法运行程序，得到上面的红色字的提示。 t=@(D)(8/pi^2)*(exp(-12000000*pi^2*D)+(1/9)*exp(-9*12000000*pi^2*D)+(1/25)*exp(-25*12000000*pi^2*D))-0.907966; [D,f,flag,out]=fsolve(t,0) Optimization terminated: norm of relative change in X is less than max(options.TolX^2,eps) and sum-of-squares of function values is less than sqrt(options.TolFun). D = 9.2703e-011 f = -1.7886e-005 flag = 2 我想知道的是红色字提示的是什么意思，我得到的解与精确解差的多不多，有什么提高的方法。希望得到指点。 [ Last edited by chuzhenhui on 2010-5-13 at 22:42 ] |
3楼2010-05-13 22:27:25
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+1):鼓励应助！ 2010-05-14 13:12:52
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+1):鼓励应助！ 2010-05-14 13:12:52
|
你可以将12000000*D作为新未知数D,结果最后除去12000000即可 >>t=@(D)(8/pi^2)*(exp(-pi^2*D)+(1/9)*exp(-9*pi^2*D)+(1/25)*exp(-25*pi^2*D))-0.907966; >>[D,f,flag,out]=fsolve(t,0) D = 0.001111597655874 f = 5.475275788313638e-011 flag = 1 out = iterations: 3 funcCount: 8 algorithm: 'trust-region dogleg' firstorderopt: 1.163031023008677e-009 message: [1x695 char] >> D/12000000 ans = 9.263313798950567e-011 [ Last edited by change0618 on 2010-5-14 at 12:47 ] |
4楼2010-05-14 12:46:40
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
5楼2010-05-14 14:24:40
chuzhenhui
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 600.7
- 散金: 35
- 帖子: 129
- 在线: 10.2小时
- 虫号: 763256
- 注册: 2009-05-05
- 性别: GG
- 专业: 食品&包装
6楼2010-05-16 12:54:50
