实验8非线性方程和微常分方程.docVIP

实验8 非线性方程和常微分方程的解法 一、实验目的 学会用MATLAB软件求解非线性方程和常微分方程。 二、实验内容与要求 1. 非线性方程的整值解 （1）最小二乘法 格式：fsolve（’fun’,）%求方程fun=0在估计值x附近的近似解。 [例1.72] 求方程的解。 fc=inline(‘x-exp(-x)’); xl=fsolve(fc,0) xl= 0.5671 问题1.28：求解方程,观察知有多解，如何求之？ 先用命令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],le-6) 得出结果： ans= 0.5918 3.1407 6.2832 9.4248 【例 1.73】求解方程组x-0.7sinx-0.2cosy=0 先编制函数文件 fu.m: function y=fu (x) y (1)=x (1) - 0.7*sin ( x (1) ) - 0.2*cos( x (2) ) ; y (2)=x (1