牛顿迭代法matlab例题
作者&投稿:魏沾 (若有异议请与网页底部的电邮联系)
数值分析有什么步骤?
第四步,有代数方程以后,我们还需要求解它。这里面就有很多方法可供选择,比如直接法如高斯法和各种不同的迭代方法。一般除非我们是做求解器的,我们可以用很多现成的工具去求解这个代数方程,而不用再写一个求解器来求解。比如,在MATLAB中就可以使用一个再简单不过的X=A\\B命令来求解。第五步,求出...
米易县夏枯回答:[答案] 用牛顿迭代法,求导x=0.29644 >> x0=0;tol=1e-6;x1=newton(x0,tol) n = 6 x1 = 0.29644 >>syms x,ezplot(sin(4*x^2-4*x+1)-(3/4)*x-1/(10*x+3)+9/40),grid on >>hold on,plot(double(x1),0,'p'),text(0.5,1,'Zeros Point') 代码见附件 图形见下图
宠雁14730328901问: 用牛顿法求解x³+10x - 20=0在区间【1,2】内的一个根,要求迭代4次 - ?
米易县夏枯回答:[答案] 使用迭代表达式: x_{n+1} = x_n - f(x_n)/f'(x_n) 其中 f(x) = x^3 + 10 x - 20; f'(x) = 3x^2 + 10 为f(x)的导数; 初始值x_0 可以使用[1, 2]的中点, 1.5.
宠雁14730328901问: 牛顿迭代法求一个方程的解 MATLAB用MATLAB求方程y=x - c*sin(y+k)其中,x,c,k都为已知就是求y就好了 - ?
米易县夏枯回答:[答案] 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)>em=n;f0=feval(f,...
宠雁14730328901问: matlab中用牛顿法求方程的根用matlab求解以下问题:用牛顿法求解方程x=e^ - x在x0=0.5附近的根(ε=10^ - 5) - ?
米易县夏枯回答:[答案] 代码如下:function rtn=newton1(fx,dfx,x0,tol,N)% Newton Method% The first parameter fx is a external function with respect to viable x.% The second parameter dfx is the first order diffential function ...
宠雁14730328901问: 用牛顿迭代法三次求方程 f(x)=x^5 - x^2+x - 30=0,在区间[1,3]中的近似值 - ?
米易县夏枯回答:[答案] 明显f(1)=-29 f(3)=207 所以(1,3)之间一定有零点. 而且明显更靠近x=3 f(2)=-2 f(2.5)=63.9 f(2.2)=19.2 f(2.1)=8.53 f(2.01)=0.78 f(2.001)=0.08 f(2.0001)=0.008 应该满足要求了
宠雁14730328901问: Matlab编程问题:编写牛顿法程序来求解下面的方程组: (x - 2)^2+(y - 3+2x)^2=5 2(x - 3)^2+(y/3)^2=4 - ?
米易县夏枯回答: ^% 用牛顿迭代法解非线性方程组% F1=(x-2)^e79fa5e98193e78988e69d83313332646631612+(y-3+2x)^2-5=0% F2=2(x-3)^2+(y/3)^2-4=0% F=[F1;F2]=[5*x^2 + 4*x*y - 16*x + y^2 - 6*y + 8% 2*x^2 - 12*x + y^2/9 +14 ]%迭代初值设为:X0=[x,y]=[0,0]'...
宠雁14730328901问: 用牛顿迭代法求解x*x*x - x*x - 1=0 - ?
米易县夏枯回答:[答案] #include void main() { float x,x0,f,f1; x0=0.5; do { f=x0*x0*x0-x0*x0-1; f1=3*x0*x0-2*x0; x=x0-f/f1; x0=x; }while((fabs(x-x0)
宠雁14730328901问: 牛顿迭代法求一个方程的解 MATLAB - ?
米易县夏枯回答: 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)>em=n;f0=feval(f,m);f2=feval(f,m+h);f1=feval...
宠雁14730328901问: 用MATLAB的迭代法求解x^3 - x - 1=0在x0=1.5附近的一个根? - ?
米易县夏枯回答: %牛顿迭代法 解方程y=x.^3-x-1 x=1.5; format long; x1=x-func1_1(x)/func1_1_1(x); if(abs(x1)<1.5) delt=abs(x1-x); else delt=abs((x1-x)/x1); end while(delt>1e-6|abs(func1_1(x))>1e-6) x=x1; x1=x-func1_1(x)/func1_1_1(x); if(abs(x1)<1.5) ...
宠雁14730328901问: 牛顿迭代法和阻尼牛顿迭代法求极小值的matlab程序 - ?
米易县夏枯回答: 阻尼牛顿迭代法没听说过牛顿迭代法的程序 function [r,n]=mulNewton(x0,eps) if nargin==1eps=1.0e-4; endr=x0-myf(x0)/dmyf(x0); n=1; tol=1; while tol>epsx0=r;r=x0-myf(x0)/dmyf(x0);tol=norm(r-x0);n=n+1;if(n>100000)disp('迭代步数太多,可能不收敛!');return;end end