第七章MATLAB解方程.pptVIP

第7章 MATLAB解方程与函数极值 7.1 线性方程组求解 7.2 非线性方程数值求解 7.3 常微分方程初值问题的数值解法 7.4 函数极值 7.1 线性方程组求解 7.1.1 直接解法 1、利用左除运算符的直接解法 对于线性方程组Ax=b，可以利用左除运算符“\”或pinv指令求解,调用格式如下: x=A\b 或x=pinv(A)*b 7.1.2 迭代解法 迭代解法非常适合求解大型系数矩阵的方程组。在数值分析中，迭代解法主要包括 Jacobi迭代法、Gauss-Serdel迭代法、超松弛迭代法和两步迭代法。 7.1.3 求线性方程组的通解 线性方程组的求解分为两类：一类是求方程组的惟一解即特解，另一类是求方程组的无穷解即通解。 7.2 非线性方程与最优化问题求解 7.2.1 非线性方程数值求解 7.2.2 无约束最优化问题求解 7.2.3 有约束最优化问题求解 7.2.2 非线性方程组的求解 对于非线性方程组F(X)=0，用fsolve函数求其数值解。fsolve函数的调用格式为： X=fsolve(fun,X0,option) 其中X为返回的解，fun是用于定义需求解的非线性方程组的函数文件名，X0是求根过程的初值，option为最优化工具箱的选项设定。 7.2.2 无约束最优化问题求解 无约束最优化问题的一般描述为: 所以,无约束最优化问题又称为最小化问题 7.2.3 有约束最优化问题求解 有约束最优化问题的一般描述为: 下标表示x要满足G(x)的约束条件 7.3 常微分方程的数值解法 7.3.1 龙格－库塔法简介 一次一阶常微分方程的初值问题的一般形式： y′=f(x,y) y(x0)=y0 在求解y时，y在x0点的值y(x0)=y0是已知的，由高数中值定理有： y(x0+h)=y1≈y0+hf(x0,y0),h0,为步长 y(x0+2h)=y2≈y1+hf(x1,y1),h0 7.4 函数极值 MATLAB提供了基于单纯形算法求解函数极值的函数fmin和fmins，它们分别用于单变量函数和多变量函数的最小值，其调用格式为： x=fmin(fname,x1,x2) x=fmins(fname,x0) 这两个函数的调用格式相似。其中fmin函数用于求单变量函数在区间(a,b)上最小值点。fname是被最小化的目标函数名，x1和x2限定自变量的取值范围。 fmins函数用于求多变量函数在估计值x0附近的最小值点。 例7-9 设有初值问题， 试求其数值解，并与精确解相比较 (精确解为 ）。 (1) 建立函数文件funt.m。 function yp=funt(t,y) yp=(y^2-t-2)/4/(t+1); (2) 求解微分方程。 t0=0;tf=10; y0=2; [t,y]=ode23(funt,[t0,tf],y0); %求数值解 y1=sqrt(t+1)+1; %求精确解 t y y1 y为数值解，y1为精确值，显然两者近似。 例7-10 求微分方程 y’=-2(y/x)+4x y(1)=2 初值问题在[1,3]区间内的数值解,并将结果与解析解(y=x2+1/x2)进行比较. (1)建立函数文件 fxy1.m: function f=f(x,y) f=-2.*y./x+4*x; %注意使用点运算符号 return (2)命令窗口输入: [X,Y]=ode45(‘fxy1’,[1,3],2); X’ ; %显示自变量的一组采样点 Y’ ; %显示求解函数与采样点对应的一组数值解 (x.^2+1./x.^2)’ %显示求解函数与采样点对应的一 组解析解 例7-11 求微分方程 y〞+xy′-x2+5=0 y(0)=5 y′(0)=6 初值问题在[0,2]区间内的数值解 分析:这是一个二阶微分方程,ode45只能求解一阶微分方程和方程组,故应该先将它转换为一阶微分方程 令u=y,v= y′得到: v′=-xv+x2-5 u′=v u(0)=5 v(0)=6 (1)建立函数文件 fxy