实验32方程（组）的数值解法.docVIP

实验3.2 方程（组）的数值解法 实验目的 （1）非线性方程（组）的数值解法 （2）线性方程组的数值解法 实验内容 1、掌握非线性方程与非线性方程组数值解的Matlab命令； 2、掌握线性方程组的数值解法。 1. 多项式的一般表达及非线性方程（组）的求根 （1）多项式在MATLAB中是由行向量来表达的，向量中的元素是多项式系数的降幂排列。其中最后一个元素代表多项式中的0幂项，即常数项，千万不可省略，否则MATLAB将无法识别这一项。 【例题3.8】MATLAB中输入四次多项式，只需在命令区中输入： p=[0.6955,0.436,0.668,1.35]; 这样就建立了这个多项式。 （2）多项式的求根命令： ① roots函数 【例题3.9】求例题3.2中多项式的根 解：输入：p=[0.6955,0.436,0.668,1.35]; x=roots(p) 回车后将得到如图3.3所示的结果。因为采用矩阵算法，所以MATLAB中的多项式和它的根都是向量。多项式为行向量，根为列向量。 图3.3 ②如果不用roots命令求根，还可以利用solve函数的命令。可以得到方程根的解析解或数值解，其命令格式为： solve(‘eqn1’,’eqn2’,…,’eqnN’); solve(‘eqn1’,’eqn2’,…,’eqnN’,’var1,var2,...,varN’); 【例题3.