[讨论]利用牛顿迭代求解非线性方程组问题，麻烦帮忙改一下吧。 — 编程爱好者论坛
主题：[讨论]利用牛顿迭代求解非线性方程组问题，麻烦帮忙改一下吧。
[专家分：0]
我的思路是这样的：输入样本D为一个n乘3的矩阵，其中第一二三列分别表示q、r、v，然后代入方程组;&&f1=sum((v(i)-k*q(i)^a*r(i)^b)*q(i)^a*r(i)^b);&f2=sum((v(i)-k*q(i)^a*r(i)^b)*q(i)^a*r(i)^b*log(q(i)));f3=sum((v(i)-k*q(i)^a*r(i)^b)*q(i)^a*r(i)^b*log(r(i)));中求解k、a、b并且使(x-x0)/x0误差最小。下面是编写的程序，请帮忙改正一下吧，谢了。clear&clcglobal&x&v&k&q&a&b&r&e1;D=[5.6,15.62,1.431;5.6,22.36,1.067;5.6,29.73,0.582;5.6,36.4,0.516;5.6,43.17,0.348;9.4,20.59,1.433;9.4,27.86,0.964;12.8,27.31,1.204];%数据样本q=D(:,1)';r=D(:,2)';v=D(:,3)';%计算输入样本的组数n=length(q);x0=[100,0.5,-0.5];%初始值%编写循环程序for&i=1:n&&&&&&if&abs(e1)&1e-4&(w&=500)&&&&&&f1=sum((v(i)-k(w)*q(i)^a(w)*r(i)^b(w))*q(i)^a(w)*r(i)^b(w));&&&&&&f2=sum((v(i)-k(w)*q(i)^a(w)*r(i)^b(w))*q(i)^a(w)*r(i)^b(w)*log(q(i)));&&&&&&f3=sum((v(i)-k(w)*q(i)^a(w)*r(i)^b(w))*q(i)^a(w)*r(i)^b(w)*log(r(i)));&&&&&&&&&&&f=[f1;f2;f3];x=[k,a,b];&&&&&&df=jacobian(f,x);%利用jacobian函数求微分&&&&&&x=newton(f,df,x0,e,N);&&&&&&x=x0-feval(f,x0)/feval(df,x0);&&&&&%利用norm函数求精度&&&&&&e1=norm((x-x0)./x,inf);%取dx的无穷范数&&&&&&&&&&&&&&&w=w+1;&&&&&&&&x0=x;&&&&&&&&disp(x);&&&&else&if&w&N&&&&warning('已达迭代次数上限');&&&&end&&&i=i+1;&&&&&&endenddisp(x)
最后更新于： 08:45:00
本帖地址：&
回复列表 （共7个回复）
[专家分：430]
请将仿真参数“&x&v&k&q&a&b&r&w&e1&e&N”给一下，以便我验证一下
[专家分：0]
我在过程中都给出了，e是精度1e-4，N为500，其他的都有了。麻烦了，谢谢！
[专家分：430]
我按以下验算，没有错误，但没有结果。&e=1e-4;N=500;global&x&v&k&q&a&b&r&w&e1&e&N;D=[5.6,15.62,1.431;5.6,22.36,1.067;5.6,29.73,0.582;5.6,36.4,0.516;5.6,43.17,0.348;9.4,20.59,1.433;9.4,27.86,0.964;12.8,27.31,1.204];%数据样本q=D(:,1)';r=D(:,2)';v=D(:,3)';%计算输入样本的组数n=length(q);x0=[100,0.5,-0.5];%初始值%编写循环程序for&i=1:n&&&&&&if&abs(e1)&1e-4&(w&=500)&&&&&&f1=sum((v(i)-k(w)*q(i)^a(w)*r(i)^b(w))*q(i)^a(w)*r(i)^b(w));&&&&&&f2=sum((v(i)-k(w)*q(i)^a(w)*r(i)^b(w))*q(i)^a(w)*r(i)^b(w)*log(q(i)));&&&&&&f3=sum((v(i)-k(w)*q(i)^a(w)*r(i)^b(w))*q(i)^a(w)*r(i)^b(w)*log(r(i)));&&&&&&&&&&&f=[f1;f2;f3];&&&&&&x=[k,a,b];&&&&&&df=jacobian(f,x);%利用jacobian函数求微分&&&&&&x=newton(f,df,x0,e,N);&&&&&&x=x0-feval(f,x0)/feval(df,x0);&&&&&%利用norm函数求精度&&&&&&e1=norm((x-x0)./x,inf);%取dx的无穷范数&&&&&&&&&&&&&&&w=w+1;&&&&&&&&x0=x;&&&&&&&&disp(x);&&&&else&if&w&N&&&&warning('已达迭代次数上限');&&&&end&&&i=i+1;&&&&&&endendx%disp(x)结果：&Warning:&The&value&of&local&variables&may&have&been&changed&to&match&the&&&&&&&&&globals.&&Future&versions&of&MATLAB&will&require&that&you&declare&&&&&&&&&a&variable&to&be&global&before&you&use&that&variable.&&In&fenxianxing&at&5Warning:&The&value&of&local&variables&may&have&been&changed&to&match&the&&&&&&&&&globals.&&Future&versions&of&MATLAB&will&require&that&you&declare&&&&&&&&&a&variable&to&be&global&before&you&use&that&variable.&&In&fenxianxing&at&5x&=&&&&&[]
[专家分：430]
全局变量没有给我。global&x&v&k&q&a&b&r&w&e1
[专家分：0]
你好，其实我是个初学者，程序里面肯定有很多的错误，至于那里面的全局变量，我是这样假定的：刚开始输入一个数据样本，其中此样本的第一二三列分别为q、r、v，需要求算的是k、a、b，程序里的w为求解次数。麻烦了。
[专家分：430]
我的机子没有newton()库函数，但是没有报错，只有警告。建议你看一本书，高教社，冯烟利翻译《数值分析》，里面有你类似方程的解例子。但是看了，你的程序结构，应该没有错误。我无能为力，另请高明解答。
[专家分：0]
您尚未登录，请登录后再回复。错误提示- _星空游戏网
你的位置：&&&&
作者：游戏玩家inone &&&更新日期：
错误提示……
0X000000该内存不能为read的解决方法 出现这个现象有方面的,一是硬件,即内存方面有问题,二...进不了某个网页 提示错误代码500……
500系列错误代码是对方网站服务器问题,以下是500常见错误 HTTP 500 - 内部服务器错误 ...win7系统上的防火墙打不开?点击打开提示错误代码0x……
按windows开始键-》运行-》调出运行窗口,输入services.msc,找到Windows F...word2007打开时总是出现“向程序发送命令时出现问题”的错误提示……
跟我的情况一模一样啊 以前装的office2007序列号可能无法验证,于是就在注册表里将原来的sn信...在Excel中出现这样的错误提示各代表什么意思?……
“####”:单元格长度不够显示 #VALUE!值计算错误,用非数值参与计算 #DIV/0!:被除数...错误提示691……
提示错误691真正意义上来讲可能是以下原因: 1:域名出现错误,(用户名或密码输入错误)。 2:电话...Altium design中PCB错误提示?……
楼主新手吧,这个绿色是DRC错误,楼主用AD打开你的PCB,点键盘L,出现对话框,将系统颜色下面DR...我的eclipse突然没有错误提示怎么办啊?……
Preferences -& General -& Editors -& Text Editors ...IE脚本错误提示……
确认活动脚本、ActiveX 和 Java 未被阻止 确认您计算机上的 Internet Explo...暴风影音提示错误代码40601……
1.卸载暴风影音; 2.在控制面板中添加删除不必要的软件及插件; 3.检查电脑是否安装多个防病毒软件...
Copyright & 2014
新 版权所有
文字内容来自网友交流，仅供参考，不表示本站同意或赞成其观点，对其准确性和真实性不作任何的担保。新生咨询，帮我一下吧，谢谢～_洛阳师范学院吧_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：79,530贴子：
新生咨询，帮我一下吧，谢谢～收藏
我统考考了214，文化大概430左右吧，能上你们学校艺术设计吗？我洛阳的。谢谢各位了～
1楼 05:29&|
应该差不多肯定管上！是美术专业么
2楼 05:29&|
应该差不多肯定管上！
3楼 05:29&|
是美术啊，最好能上装潢吧，别的也行～谢谢你
4楼 05:29&|
这个，美术，今年应该可以的，装潢、数码跟环艺的分数线是分开的，一般装潢跟数码的高，你加起来640多，可以报装潢，不过数码的学长建议你报数码，，，写生机会多……
5楼 05:29&|
师院的装潢其实算平面，数码是动画，……所以……平面现在市场饱和……你懂的哦……
6楼 05:29&|
师院的装潢其实算平面，数码是动画，……所以……平面现在市场饱和……你懂的哦……
7楼 05:29&|
师院的装潢其实算平面，数码是动画，……所以……平面现在市场饱和……你懂的哦……
8楼 05:29&|
师院的装潢其实算平面，数码是动画，……所以……平面现在市场饱和……你懂的哦……
9楼 05:29&|
师院的装潢其实算平面，数码是动画，……所以……平面现在市场饱和……你懂的哦……
10楼 05:29&|
小白你这是卡了吗
11楼 05:29&|
卡卡卡卡卡……
12楼 05:29&|
谢谢LS了，真有爱。分数稍低点的都是什么专业呢？分数够的话会考虑数码的～
13楼 05:29&|
登录百度帐号我的游戏推荐游戏
后查看最近玩过的游戏
内&&容:使用签名档&&
想了解更多关于 ”洛阳师范学院“的信息,请&或精彩推荐：
微信扫一扫到手机
随时看帖,分享到朋友圈
手机阅读→
发表于： 13:44
用户被禁言,该主题自动屏蔽!
本帖评分记录
【常州宝达】售价不变配置更高——2013款迈腾正式上市楼主可以看一下这个 帖子，然后打电话过去咨询一下
外星的牛牛
九品芝麻官
新人来此。。。。。。。。。
～～ 婚期7.7
婚纱照：巴黎心情婚纱摄影
婚庆：巴黎爱妃 看贴回贴好习惯！谁删我贴跟谁急！
化龙巷网站·可爱MM奖
新人进步奖
化龙巷网站·新人进步奖
最近买一起大众的朋友多的嘛，听说钟楼一汽大众在化龙巷的口碑一直不错，你可以电话去咨询他们啊 。好像是
化龙巷网站·个人认证
全款买，省事。
回 楼主(神州4号) 的帖子
我有个朋友在武进汽车城一汽大众的，你如果需要可以联系我QQ
做那个50/50的贷款&&还可以 省事么直接全款
　　　前几天刚从迈腾4S店里购完车，非常认真的在买车前做了笔记哈哈。简单给LZ说说：　　　他们的金融方案一共五项，五选一。第一种是首付最低20%，1年0利率。但貌似能审批下来的都没有那么低，个人认为30%左右的首付审批是完全木有问题滴。第二种是首付50%，一年后结清另外50%，月供最低288。也可以选择信用卡分期，有建行和招行可以选择，具体利率还是咨询下4S店。如果LZ有二手车，也可以用其他品牌的二手车置换迈腾，补贴八千。 第五种就是全款购买迈腾，赠送价值五千的三大主险。我选的是置换，当时交款时候直接补贴了。基本就这些了，希望能帮到你！楼主也可以去他们官网看看。
&&&&二手车置换，补贴8千太亏了吧？&&&&&&&&&&&&&&&&
&&&&有钱的话，当然是一次付款最好了。&&&&
大众车性能很好，节能省油，比较不错&&&&&&&&&&&&
&&&&&&&&城市代步的话，迈腾不错了，好车~~&&&&&&&&
最后的感谢
二手车闲杂越来越不值钱，能不贴八千不错了
大家一起来聊聊 哪种方式最划算呢？
使用（可批量传图、插入视频等）
[完成后可按 Ctrl+Enter 发布]
广告热线：400 970 0519 转 8888 投诉受理：400 970 0519 转 9999 法律顾问：江苏日月泰律师事务所 费小强 律师
信息产业部备案/许可证编号：苏ICP备号&经营性ICP：苏B2-号&电子公告服务许可文件号:苏通[号}