matlab求解迭代方程例子

胥峰莎3994用matlab自行编程实现:采用迭代法的思想求方程 x3 - 3.2x2+1.9x+0.2=0的在根区间[1,4.2]内的根,要求准确到小数点后第4位.初值x0选1. -
赫采初19775677932 ______ x0=1; y=x0^3-3.2*x0^2+1.9*x0+0.2; while (y<0) x0=x0+0.0001; y=x0^3-3.2*x0^2+1.9*x0+0.2; end

胥峰莎3994matlab简单迭代公式编程 -
赫采初19775677932 ______ N=1000; x=[-.8,zeros(1,N)]; for n=0:N-1 x(n+1+1)=x(n+1)^2-2;%matlab从1开始编号,所以再+1 end plot(0:N,x)

胥峰莎3994求高手!用牛顿迭代法解三元三次方程的MATLAB程序! -
赫采初19775677932 ______ 用fsolve可解出来: 先构造函数: function output=solveproblem(X) c=X(1); m=X(2); y=X(3); output(1)=(1-c)*(1-y)*(1-m)*10.61+y*(1-c)*(1-m)*30.56+c*(1-y)*(1-m)*60.67+m*(1-c)*(1-y)*0.95+y*m*(1-c)*28.45+y*c*(1-m) *8.19+c*m*(1-y) *2.51...

胥峰莎3994MATLAB中用迭代法求x=sqrt(a) -
赫采初19775677932 ______ xn=1;xn1=2;a=16; while abs(xn1-xn)>1e-3 xn1=(xn+a/xn)/2; xn=(xn1+a/xn1)/2; end x=xn

胥峰莎3994用matlab解三元三次方程组,牛顿迭代法.求程序,及结果 -
赫采初19775677932 ______ syms X Y Z eq1=X+Y+4*(Z-40)-4410; eq2=Z*X^3-(Z-40)*(X-80)^3+80*Z^3-105304178*10^4; eq3=Z*Y^3-(Z-40)*(Y-80)^3+80*Z^3-306400178*10^4; d=solve(eq1,eq2,eq3) d.X d.Y d.Z

胥峰莎3994MATLAB求解方程系数,重复迭代
赫采初19775677932 ______ 试试. 建立一个m文件mysolve.m如下: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function chanshu=mysolve(index) % index是被抽中的四个数的下标组成的向量.可能为[3,1,4,6]...

胥峰莎3994用MATLAB求解SOR迭代程序
赫采初19775677932 ______ 以a=[6,2,1,-1;2,4,1,0;1,1,4,-1;-1,0,-1,3];b=[6;1;5;-5];为例 function[y,r,n]=sor(a,b,x0,w,e,N) D=diag(diag(a));U=-triu(a,1);L=-tril(a,-1); G=(D-w*L)\((1-w)*D+w*U);f=w*((D-w*L)\b); y=G*x0+f;n=1; while(norm(y-x0)>=e&n<N) x0=y;y=G*x0+f;n=n+1; end n;y;...

胥峰莎3994MATLAB用牛顿迭代求解非线性方程的程序 -
赫采初19775677932 ______ 给你一个完整版: % 牛顿法求解非线性方程组 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(...

胥峰莎3994matlab求助 怎样用循环迭代求解 -
赫采初19775677932 ______ x=0:1000; x(1)=0.3; //因为在MATLAB的数组中没有0下标,所以从1开始即X0=X1,所求的X1000=X1001 //循环语句中默认步长为1 for n=1:1000 x(n+1)=(x(n)+2.5)*x(n); end x(1001)

胥峰莎3994Matlab 牛顿加速迭代法 急…………Matlab 牛顿加速迭代法 急…………记住是 牛顿“加速”迭代法解方程是x=e^( - x)在x=0.5附近的根 -
赫采初19775677932 ______[答案] clc;clear f='x-exp(-x)' df=diff(f) x=0.5; ac=[0.7;0.9;1];%加速因子 ac2=0.9; for k=1:10 x=x-ac.*subs(f/df); X(k,:)=x; end format long g X plot(X,'o-') 结果: X = 0.546417702238053 0.559679902877496 0.566311003197218 0.560870862876879 0....