版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

辛迪熊啦

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 54.8
帖子: 6
在线: 3小时
虫号: 2645279
注册: 2013-09-10
性别: MM
专业: 资源环境政策与管理

[求助] 牛顿迭代法接非线性方程组求解已有1人参与

新手，菜鸟
有一组比较复杂的方程组需要求解，含有15个未知数，15个方程，找到一个案例照着编程了一下
运行结果是Warning: Matrix is close to singular or badly scaled. Results may be inaccurate.
RCOND =  1.322942e-30
请问这种现象是由于方程本身有问题，还是求解算法有问题，还是初值给的有问题？
如果是算法有问题，还有什么其他算法可以推荐呢
(还有我是新手所以金币不多如果不满意金币数量我可以后期再多赚金币补偿^-^)

以下是用到的编程案例，我是照着这个案例编程的，把我这15个方程套进去，最后得出了上面的结论
matlab实现牛顿迭代法求解非线性方程组
已知非线性方程组如下
3*x1-cos(x2*x3)-1/2=0
x1^2-81*(x2+0.1)^2+sin(x3)+1.06=0
exp(-x1*x2)+20*x3+(10*pi-3)/3=0
求解要求精度达到0.00001
---------------------------------------------------------分--割--线---------------------------------------------------------
首先建立函数fun
储存方程组编程如下将fun.m保存到工作路径中:
function f=fun(x);
%定义非线性方程组如下
%变量x1 x2 x3
%函数f1 f2 f3
syms x1 x2 x3
f1=3*x1-cos(x2*x3)-1/2;
f2=x1^2-81*(x2+0.1)^2+sin(x3)+1.06;
f3=exp(-x1*x2)+20*x3+(10*pi-3)/3;
f=[f1 f2 f3];
---------------------------------------------------------分--割--线---------------------------------------------------------
建立函数dfun
用来求方程组的雅克比矩阵将dfun.m保存到工作路径中:
function df=dfun(x);
%用来求解方程组的雅克比矩阵储存在dfun中
f=fun(x);
df=[diff(f,'x1');diff(f,'x2');diff(f,'x3')];
df=conj(df');
---------------------------------------------------------分--割--线---------------------------------------------------------
编程牛顿法求解非线性方程组将newton.m保存到工作路径中:
function x=newton(x0,eps,N);
con=0;
%其中x0为迭代初值eps为精度要求N为最大迭代步数con用来记录结果是否收敛
for i=1:N;
f=subs(fun(x0),{'x1' 'x2' 'x3'},{x0(1) x0(2) x0(3)});
df=subs(dfun(x0),{'x1' 'x2' 'x3'},{x0(1) x0(2) x0(3)});
x=x0-f/df;
for j=1:length(x0);
      il(i,j)=x(j);
end
if norm(x-x0)<eps
      con=1;
      break;
end
x0=x;
end

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2015-02-24 17:25:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

匿名

用户注销 (文坛精英)

应助: 164 (高中生)
金币: 9550.6
散金: 34344
红花: 361
沙发: 31
帖子: 20580
在线: 3633.3小时
虫号: 0
注册: 2011-12-15
专业: 建筑环境与结构工程

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
辛迪熊啦: 金币+10, ★有帮助, 谢谢，我再去查找看看pinv怎么使用:) 2015-02-24 22:42:31

本帖仅楼主可见

赞一下

2楼2015-02-24 21:44:35

已阅申请仿真EPI 回复此楼编辑查看我的主页

相关版块跳转我要订阅楼主辛迪熊啦的主题更新