声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1141|回复: 7

[编程技巧] 关于方程组估计解的问题

[复制链接]
发表于 2007-7-22 19:32 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
众所周知,用fsolve解方程组时要估计解,这个解的估计有技巧吗?不解 出方程如何知道解大概是多少?若解出来了也就不用MATLAB解了,请高手指点
回复
分享到:

使用道具 举报

发表于 2007-7-22 20:16 | 显示全部楼层
可以使用隐函数绘图,看曲线的大致走势。然后在过零点附近取个值作为解的估计值
 楼主| 发表于 2007-7-23 10:46 | 显示全部楼层
Thank you!虽然不大懂,只能怪我功力不够了
发表于 2007-7-23 10:55 | 显示全部楼层
花如月的方法适合三维以下的非线性方程组,如果维数多了,就得具体情况具体分析了。或者使用1stOpt来求解。

评分

1

查看全部评分

发表于 2007-7-23 10:59 | 显示全部楼层

回复 #3 wwg.no.19 的帖子

举个简单的例子
要求sin(3*x)=0的解,可以先用fplot(@(x) sin(3*x),[0.5 6])画出图形可以看到在1,2,3,4左右都有解。
>> x = fsolve(@(x) sin(3*x),1,optimset('Display','off'))

x =

    1.0472

>> x = fsolve(@(x) sin(3*x),2,optimset('Display','off'))

x =

    2.0944

>> x = fsolve(@(x) sin(3*x),3,optimset('Display','off'))

x =

    3.1416

>> x = fsolve(@(x) sin(3*x),4,optimset('Display','off'))

x =

    4.1888

>>  x = fsolve(@(x) sin(3*x),[1 2 3 4],optimset('Display','off'))

x =

    1.0472    2.0944    3.1416    4.1888
发表于 2007-7-23 11:07 | 显示全部楼层
是啊,多维非线性方程组的情形就复杂了。本身不是学数学的,解估计的技巧也就没资格谈了。能想到的这个方法应用确实很受限制。路过的高手也谈谈自己的经验吧

[ 本帖最后由 花如月 于 2007-7-23 11:10 编辑 ]
发表于 2007-7-23 13:19 | 显示全部楼层
原帖由 花如月 于 2007-7-23 11:07 发表
是啊,多维非线性方程组的情形就复杂了。本身不是学数学的,解估计的技巧也就没资格谈了。能想到的这个方法应用确实很受限制。路过的高手也谈谈自己的经验吧



关于fsolve和1stOpt的问题好像版面之前就有讨论:非线性方程组-1stOpt与fsolve的比较!

[ 本帖最后由 eight 于 2007-7-23 13:21 编辑 ]
发表于 2008-3-27 18:12 | 显示全部楼层

easy

try 1stOpt!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-18 08:35 , Processed in 0.054455 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表