sphm.net
当前位置:首页 >> mAtlAB 牛顿法 >>

mAtlAB 牛顿法

首先,建立原函数的自定义函数文件 function y = fun(x) y=x^4-3*x^3+5*cos(x)+8; 其次,建立导函数的自定义函数文件 fu...

这是两部分代码: function eq=exam1_5(x) eq=(x^7+2*x.^5+3*x.^3+5*x+7)/(7*x^6+10*x^4+9*x.^2+5); 这两行是函数部分,参数为x,返回值为eq。 下面的命令是在命令窗口输入的: x0=0.1; x1=x0-exam1_5(x0); while abs(x1-x0)>0.0001 x0=x1; x1=x...

function [ A ] = cal( a,b,v )%a,b表示区间,v是精度 i=1; x = (a+b)/2; A=[i x]; t = x-(x^3-x-1)/(3*x^2-1);%迭代函数 while(abs(t-x)>v) i=i+1; x = t; A = [A;i x]; t = x-(x^3-x-1)/(3*x^2-1);%迭代函数 end A = [A;i+1 t]; end 运行结果...

syms x f=x^x-10; df=diff(f,x); eps=1e-6; x0=10; cnt=0; MAXCNT=200; %最大循环次数 while cnt

对于求平方根,变成方程模式为f(x)=x^2-a,即求此方程的实根; 下面编写了两个function函数,可以直接调用。 二分法: function x=sqrt_bisect(a) f=@(x)x^2-a; if a0 xb=x; elseif f(xa)*f(x)>0 xa=x; else break end end end x; 牛顿迭代法: ...

定义函数 function y=nd(x) y= 0.036-((x/2090.7).^(1/0.1585))-x/182000 function y=nd0(x) y=- (1/0.1585)*(x/2090.7).^(1/0.1585-1)*(1/2090.7) - 1/182000 主程序 x=0.5; %迭代初值 i=0; %迭代次数计数 while i10^(-5); %收敛判断 x=y; else...

function [ A ] = cal( a,b,v )%a,b表示区间,v是精度 i=1; x = (a+b)/2; A=[i x]; t = x-(x^3-x-1)/(3*x^2-1);%迭代函数 while(abs(t-x)>v) i=i+1; x = t; A = [A;i x]; t = x-(x^3-x-1)/(3*x^2-1);%迭代函数 end A = [A;i+1 t]; end 运行结果...

这个牛顿法程序对吗?给你一个可用的:function main() clc; clear all; f = @(x)log(x+sin(x)); df = @(x)(1+cos(x))/(x+sin(x)); x0 = 0.1; x = TestNewton(f, df, x0)function x = TestNewton(fname, dfname, x0, e, N) % 用途:Newton迭代法...

Newton's method is the classic algorithm for finding roots of functions. The breif introduction is written below. The analytic solution is The code is clear allx = zeros(1,2);f = zeros(1,2);tol = 1E-7;invJ = zeros(2,2);x(1) = 0...

Newton-Raphson 求解非线性方程组matlab源程序 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+...

网站首页 | 网站地图
All rights reserved Powered by www.sphm.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com