数值分析 非线性方程组的数值方法.pptVIP

局部收敛性原理原理的局限性：（1）收敛域很难找（2）对非线性问题这是一个充分性原理，不是充分必要的，只有对线性问题，才是充分必要条件.如：第29页，共62页，星期日，2025年，2月5日第30页，共62页，星期日，2025年，2月5日P=1，C1为线性收敛，P=2为平方收敛。收敛速度

第31页，共62页，星期日，2025年，2月5日第四节非线性方程组的Newton型算法

一、Newton-Raphson方法的迭代格式

第32页，共62页，星期日，2025年，2月5日第33页，共62页，星期日，2025年，2月5日第34页，共62页，星期日，2025年，2月5日第35页，共62页，星期日，2025年，2月5日第36页，共62页，星期日，2025年，2月5日二、同伦算法第37页，共62页，星期日，2025年，2月5日第38页，共62页，星期日，2025年，2月5日第39页，共62页，星期日，2025年，2月5日第40页，共62页，星期日，2025年，2月5日第41页，共62页，星期日，2025年，2月5日第42页，共62页，星期日，2025年，2月5日求解非线性方程组的同伦算法%文件名：Homotopy.mfunctionroot=Homotopy(funcF,funcH,x,N,tol,Nmax)%功能:求解非线性方程组的同伦算法%输入:%funcF-----原始方程函数句柄%funcH-----同伦方程函数句柄%x----初始迭代点%N(可选)----同伦参数t的划分数（默认是10）%tol(可选)---精度要求（默认是1e-4）%Nmax(可选)---最大迭代次数（默认100次）%输出:%root----解向量第43页，共62页，星期日，2025年，2月5日ifsize(x,1)==1;x=x;endifnargin6;Nmax=100;endifnargin5|isempty(tol);tol=1e-4;endifnargin4|isempty(N);N=10;endk=0;x0=x;n=length(x);dx=tol+1;f0=tol+1;form=1:Nt=m/N;while(norm(dx)tolnorm(f0)tol)kNmax[jac,f0]=JacobianH(x,funcF,funcH,x0,t);dx=jac\(-f0);x=x+dx;k=k+1;Enddx=tol+1;f0=tol+1;root=x;k=0;end第44页，共62页，星期日，2025年，2月5日数值分析数值分析非线性方程组的数值方法第1页，共62页，星期日，2025年，2月5日非线性方程组解的复杂性第2页，共62页，星期日，2025年，2月5日clear,clfx1=-2:.2:2;y2=-2:.2:2;y1=f1(x1);x2=f2(y2);plot(x1,y1,r:,x2,y2,b)xlabel(x),ylabel(y)(3)a=0(4)a=-1(1)a=1(2)a=1/4第3页，共62页，星期日，2025年，2月5日几类典型非线性问题第4页，共62页，星期日，2025年，2月5日第5页，共62页，星期日，2025年，2月5日例：半线性椭圆型边值问题解:(1)剖分求解域.YN+1N:210012….NN+1X第6页，共62页，星期日，2025年，2月5日(2)对微分算子进行离散.在每个点(xi,yj)上的有限差分方程为在边界上第7页，共62页，星期日，2025年，2月5日对非边界点进行编号:顺序为-----从下往上,从左往右相应的解向量和右端向量分别为第8页，共62页，星期日，2025年，2月5日第9页，共62页，星期日，2025年，2月5日多元向量值函数的导数第10页，共62页，星期日，2025年，2月5日第11页，共62页，星期日，2025年，2月5日第12页，共62页，星期日，2025年，2月5日多元向量值函数的导数第13页，共62页，星期日，2025年，2月5日第14页，共62页，星期日，2025年，2月5日第15页，共62页，星期日，2025年，2月5日多元实函数的高阶导数第16页，共62页，星期日，2025年，2月5日