《MATLAB数值计算基础与实例教程 》课件（共十章）.pptVIP

10.1单变量函数的求解第10章【例10-1】使用二分法求解f(x)=x3cosx+2x2-2sinx在区间[1,4]内的根根据题意可知，先进行方程的绘制如图10-1（见下页），输入以下代码进行实现。fplot([x^3*cos(x)+2*x^2-2*sin(x)],[1,4])接着，创建M文件，文件名为fun.m，输入以下代码，进行保存。functionf=fun(x)f=x^3*cos(x)+2*x^2-2*sin(x);最后，在命令窗口输入以下代码，进行结果实现。test1(fun,1,4)10.1单变量函数的求解第10章【例10-1】使用二分法求解f(x)=x3cosx+2x2-2sinx在区间[1,4]内的根图10-1交点图得到以下结果。k=1c=-0.1074error=2.2148yc=0.2362ans=-0.107410.1单变量函数的求解10.1.2迭代法第10章在给定的实数域上的光滑实值函数ψ(x)及初值x0，定义数列： 其中数列xn(n=0,1,…)称为迭代函数ψ(x)的迭代序列。迭代过程的几何意义就是将求方程f(x)=0的根的问题转化为求两曲线的交点问题，交点的横坐标就是方程的根x*。10.1单变量函数的求解第10章迭代法的算法取初始点x0，最大迭代次数N和精度要求ε，令k=0计算xk+1=ψ(xk)若|xk+1-xk|ε，则停止计算若k=N，则停止计算；否则，令k=k+1，转第二步10.1单变量函数的求解10.1.3抛物线法第10章抛物线法使用的是3个点来进行准确根的逼近，该方法称为抛物线法。假定方程f(x)=0的根为x’，使用迭代计算的相邻的3个点的坐标分别为xk-2、xk-1和xk。接着将使用3点推导获得下个坐标xk+1。记f(xk-2)=fk-2，f(xk-1)=fk-1，f(xk)=fk，则过3点的抛物线可写成： 由此可知该式通过曲线上的3个点分别为(xk-2,fk-2)、(xk-1,fk-1)、(xk,fk)。10.1单变量函数的求解第10章在MATLAB中没有现成的函数可以实现抛物线法求解非线性方程根的解，可以自行编写进行实现。functionxr=test3(fun3,x0,x1,x2,D)ifnargin5D=1e-6;endak=inf;whileabs(ak)D;f2=feval(fun3,x2);f1=feval(fun3,x1);f0=feval(fun3,x0);ak=f2;ck=[(f0-f2)/(x0-x2)-(f0-f2)/(x0-x2)]/(x0-x2);bk=(f2-f1)/(x2-x1)+ck*(x2-x1);x0=x1;x1=x2;x2=x2-2*ak/[bk+sign(bk)*sqrt(bk^2-4*ak*ck)];endxr=x2;10.1单变量函数的求解第10章【例10-2】使用抛物线法，求解方程f(x)=x3cosx+2x2-2sinx在区间[1,3]内的解。创建M文件，文件名为fun3.m，输入以下代码，进行保存。functiony=fun3(x)y=x.^3.*cos(x)+2*x.^2-2*sin(x);在命令窗口输入以下代码，进行结果实现。a=1;b=2;x0=3;xr=test3(fun3,a,b,x0)xr=2.397810.1单变量函数的求解10.1.4牛顿法第10章设x0为方程f(x)=0的一个近似的根，将f(x)在x0点附近展开成泰勒级数： 取其线性部分作为非线性方程f(x)=0的近似方程，则有： 设，则其解为： 再把f(x)在x1附近展开成泰勒级数，取其线性部分作为f(x)=0的近似方程；得到牛顿法的迭代序列为： 10.1单变量函数的求解10.1.4牛顿法第10章牛顿法的算法取初始点x0，最大迭代次数为N，精度要求为ε，令k=0如果f’(xk)=0，则停止计算；否则计算：若，则停止计算若k=N，则停止计算；否则，令k=k+1，转到第二步10.1单变量函数的求解10.1.5正割法第10章用差商来代替牛