Newton-Raphson牛顿法求解非线性方程组matlab源程序

Newton-Raphson牛顿法求解非线性方程组matlab源程序 function hom [P,iter,err]=newton('f','JF',[7.8e-001;4.9e-001; 3.7e-001],0.01,0.001,1000); disp(P); disp(iter); disp(err); function Y=f(x,y,z) Y=[x^2+y^2+z^2-1; 2*x^2+y^2-4*z; 3*x^2-4*y+z^2]; function y=JF(x,y,z) f1='x^2+y^2+z^2-1'; f2='2*x^2+y...

function hom [P,iter,err]=newton('f','JF',[7.8e-001;4.9e-001; 3.7e-001],0.01,0.001,1000); disp(P); disp(iter); disp(err); function Y=f(x,y,z) Y=[x^2+y^2+z^2-1; 2*x^2+y^2-4*z; 3*x^2-4*y+z^2]; function y=JF(x,y,z) f1='x^2+y^2+z^2-1'; f2='2*x^2+y^2-4*z'; f3='3*x^2-4*y+z^2'; df1x=diff(sym(f1),'x'); df1y=diff(sym(f1),'y'); df1z=diff(sym(f1),'z'); df2x=diff(sym(f2),'x'); df2y=diff(sym(f2),'y'); df2z=diff(sym(f2),'z'); df3x=diff(sym(f3),'x'); df3y=diff(sym(f3),'y'); df3z=diff(sym(f3),'z'); j=[df1x,df1y,df1z;df2x,df2y,df2z;df3x,df3y,df3z]; y=(j); function [P,iter,err]=newton(F,JF,P,tolp,tolfp,max) %输入P为初始猜测值，输出P则为近似解 %JF为相应的Jacobian矩阵 %tolp为P的允许误差 %tolfp为f(P)的允许误差 %max：循环次数 Y=f(F,P(1),P(2),P(3)); for k=1:max J=f(JF,P(1),P(2),P(3)); Q=P-inv(J)*Y; Z=f(F,Q(1),Q(2),Q(3)); err=norm(Q-P); P=Q; Y=Z; iter=k; if (err function homework4 [P,iter,err]=newton('f','JF',[7.8e-001;4.9e-001; 3.7e-001],0.01,0.001,1000); disp(P); disp(iter); disp(err); function Y=f(x,y,z) Y=[x^2+y^2+z^2-1; 2*x^2+y^2-4*z; 3*x^2-4*y+z^2]; function y=JF(x,y,z) f1='x^2+y^2+z^2-1'; f2='2*x^2+y^2-4*z'; f3='3*x^2-4*y+z^2'; df1x=diff(sym(f1),'x'); df1y=diff(sym(f1),'y'); df1z=diff(sym(f1),'z'); df2x=diff(sym(f2),'x'); df2y=diff(sym(f2),'y'); df2z=diff(sym(f2),'z'); df3x=diff(sym(f3),'x'); df3y=diff(sym(f3),'y'); df3z=diff(sym(f3),'z'); j=[df1x,df1y,df1z;df2x,df2y,df2z;df3x,df3y,df3z]; y=(j); function [P,iter,err]=newton(F,JF,P,tolp,tolfp,max) %输入P为初始猜测值，输出P则为近似解 %JF为相应的Jacobian矩阵 %tolp为P的允许误差 %tolfp为f(P)的允许误差 %max：循环次数 Y=f(F,P(1),P(2),P(3)); for k=1:max J=f(JF,P(1),P(2),P(3)); Q=P-inv(J)*Y; Z=f(F,Q(1),Q(2),Q(3)); err=norm(Q-P); P=Q; Y=Z; iter=k; if (err

                    本文档为【Newton-Raphson牛顿法求解非线性方程组matlab源程序】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
                    图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。 
 该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。

                    [版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。

                    本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。

                    网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
                

下载需要：免费已有0 人下载

立即下载

你可能还喜欢

最新资料

资料动态

专题动态

is_899298

暂无简介~

格式：doc

大小：42KB

软件：Word

页数：3

分类：互联网

上传时间：2011-05-14

浏览量：219

热点搜索

江苏省建设工程施工现场主要管理人员变更申请表品质基本理论看性治疗师童嵩珍如何在北京电台《和100个人谈性》高中化学竞赛中级无机化学质子理论溶剂体系理论及其拉平效应和区分效应(共24张PPT) 顶管顶力技术计算来料加工合同,英文氟尿嘧啶类化疗药物PPT课件承接查验流程 8.用最小的代价获得最大的回报——什么是长期生产函数？初三作文那一刻我很幸福作文1000字 (华师版)初一数学最全最经典习题册湘夫人ppt课件 word快捷键大全(完整版) 发声亮剑表态发言稿150字江苏省建设工程施工现场主要管理人员变更申请表品质基本理论看性治疗师童嵩珍如何在北京电台《和100个人谈性》高中化学竞赛中级无机化学质子理论溶剂体系理论及其拉平效应和区分效应(共24张PPT) 顶管顶力技术计算来料加工合同,英文氟尿嘧啶类化疗药物PPT课件承接查验流程 8.用最小的代价获得最大的回报——什么是长期生产函数？初三作文那一刻我很幸福作文1000字 (华师版)初一数学最全最经典习题册湘夫人ppt课件 word快捷键大全(完整版) 发声亮剑表态发言稿150字