sphm.net
当前位置:首页 >> 求mAtlAB编程 用二分法和牛顿迭代法求根号A >>

求mAtlAB编程 用二分法和牛顿迭代法求根号A

对于求平方根,变成方程模式为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; 牛顿迭代法: ...

牛顿迭代法(Newton's method)又称为牛顿-拉夫逊(拉弗森)方法(Newton-Raphson method),它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。 牛顿迭代法(Newton's method)又称为牛顿-拉夫逊(拉弗森)方法(Newton-Raphso...

function test clear clc %实验方程:3*x.^2+x+2*exp(x)=0 %原函数 f=@(x)3*x.^2+x-2*exp(x); %导函数 df=@(x)6*x+1-2*exp(x); %原函数在[-1 0]上的图像(有根范围) fplot(f,[-1 0]) hold on %牛顿切线法 [x1,n1]=fnewton(f,df,-0.5); disp(spr...

1. 创建一个函数 %牛顿法求立方根 function x=cube_newton(a) f=@(x)x^3-a; df=diff(sym('x^3-a')); if a==0; x1=a; else x0=a; x1=x0-f(x0)/subs(df,x0); while abs(x1-x0)>1e-6 x0=x1; x1=x0-f(x0)/subs(df,x0); end end x=x1; 2. 调用求解 >>...

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 x0=2; f=x^3-3*x-1; eps=1e-6; maxcnt=1000; fx=diff(f,x); x1=x0; cnt=1; while cnt

function [n,x]=sor22(A,b,X,nm,w,ww) %用超松弛迭代法求解方程组Ax=b %输入:A为方程组的系数矩阵,b为方程组右端的列向量,X为迭代初值构成的列向量,nm为最大迭代次数,w为误差精度,ww为松弛因子 %输出:x为求得的方程组的解构成的列向量,n...

m=0;%起始点 e=0.00001;%精度 h=0.000001;%步长 f=inline('1-y-2*sin(y+3)','y'); %x=1,c=2,k=3代入具体数值 t=0; f0=feval(f,m); f2=feval(f,m+h); f1=feval(f,m-h); n=m-2*h*f0/(f2-f1); while abs(1-m/n)>e m=n; f0=feval(f,m); f2=feval(f,m...

这个有例子 区间二分法: 与对分查找法相同 1 区间二分法求出的仅仅是方程的一个单根,如果方程有重根或者多个根时,在做区间二分法时就会出现分叉,这样方程有几个根,就会产生几个实数序列,每一个实数序列的极限便是方程的一个根 2 通常用区...

:牛顿迭代法的原理是, 实现过程如下: function y = fun_main(x) %定义函数 y =x^3+4*x^2-10; end function y = dfun_main(x) %定义导函数 y =x^3+4*x^2-10; y = 3*x^2 + 8*x; end syms x ezplot('y=x^3+4*x^2-10',[-5 4 -11 5]),grid on,hold o

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