实验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 先编制函数文件fu.m： function y=fu(x) y(1)=x(1)-0.7*sin(x(1))-0.2*cos(x(2)); y(2)=x(2)-0.7*cos(x(1))+0.2*sin(x(2)); y=[y(1), y(2)]; 在命令窗口调用fu运算： x1=fsolve(fu,[0.5,0.5]) x1= 0.5265 0.5079 （2）零点法 格式：fzero(‘fun’, x0) %求函数fun在x0附近的零点. 说明：估计值x0若为标量时，则在x0附近查找零点，x0=[x1, x2]向量时，则首先要求函数值fun(x1) fun(x2)0，然后将严格在[x1, x2]区间内寻找零点，若找不到，系统将给出提示. 【例1.74】 求函数f (x)=sinx2/x+xex－4的零点. fn=inline(sin(x^2)/x+x*exp(x)-4); x=fzero(fn,[1,2]) %这里的fn不要加单引号 x = 1.0748 注意：1.方程解的估计值x0可用fplot作图看出； 2.用function建立函数文件fn，求解调用时fn两边要加单引号，而用inline时fn两边不要加单引号； 3.这两种方法也可解线性方程组. 2. 代数方程的符号解 格式：g = solve(eq) %求解方程eq=0，输入参量eq可是符号表达式或字符串表达式. g = solve(eq,var) %对eq中指定的变量var求解方程eq(var)=0. g = solve(eq1,eq2,…,eqn) %求解方程组eq1=0,eq2=0,…,eqn=0. 【例】 solve(a*x^2 + b*x + c) solve(a*x^2 + b*x + c,b) [x,y]=solve(x + y = 1,x - 11*y = 5) [a,u,v] = solve(a*u^2 + v^2, u - v = 1, a^2 - 5*a +6) 计算结果为： ans = [ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))] ans = -(a*x^2+c)/x x = 4/3 y = -1/3 a = [ 2] [ 2] [ 3] [ 3] u = [ 1/3+1/3*i*2^(1/2)] [ 1/3-1/3*i*2^(1/2)] [ 1/4+1/4*i*3^(1/2)] [ 1/4-1/4*i*3^(1/2)] v = [ -2/3+1/3*i*2^(1/2)] [ -2/3-1/3*i*2^(1/2)] [ -3/4+1/4*i*3^(1/2)] [ -3/4-1/4*i*3^(1/2)] 注意：对于单个的方程或方程组，若不存在符号解，则返回方程（组）的数值解. 问题：用符号法求解问题1.28中的方程，结果不对，所以要验根，多用几种方法相互验证. 用符号法解方程3x2－ex=0，解的表达式不易懂，怎么办？ x=solve（3*x^2－exp（x）） x = [ -2*lambertw(-1/6*3^(1/2))] [ -2*lambertw(-1,-1/6*3^(1/2))] [ -2*lambertw