实验8 非线形程和常微分方程的解法.docVIP

实验8 非线性方程和常微分方程的解法 一、实验目的 求解非线性方程和常微分方程. 二、实验内容与要求 1. 非线性方程的数值解 （1）最小二乘法 格式：fsolve(‘fun’, x0) %求方程fun=0在估计值x0附近的近似解. 【例】 求方程x－e－x=0的解. fc=inline(x-exp(-x)); x1=fsolve(fc,0) x1= 0.5671 f(x)=5x2sinx－e-x图像 问题：求解方程5x2sinx－e-x=0，观察知有多解，如何求之？ 先用命令fplot(‘[5*x^2*sin(x)－exp(－x),0]’,[0,10])）作图1.13，注意，[5*x^2*sin(x)－exp(－x),0]中的“[…，0]”是作y=0直线，即x轴. 方程在[0,10]区间从图中可看出有4个解，分别在0，3，6，9附近，所以用命令： fun=inline(5*x.^2.*sin(x)-exp(-x))； fsolve(fun,[0,3,6,9],1e-6) 得出结果： ans = 0.5018 3.1407 6.2832 9.4248 【例】 求解方程组x－0.7sinx－0.2cosy=0 y－0.7cosx+0.2siny=0 先编制函数文件f