资源描述：
《线非线性方程近似解》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、输入多项式f(x)的系数c(用数组给出，按降幕排列)，输出r为方程f(x)的全部根；输入f(x)=0的全部根r(用数组给出)，输出c为方程f(x)的系数(按降幕排列)；c=poly(r)实验项目3：线、非线性方程近似解一、实验目的用MATLAB数值求解线性代数方程组、数值求解非线性代数方程的迭代法和牛顿法并对作初步分析。实验要求掌握用MATLAB数值求解线性代数方程组、数值求解非线性代数方程的方法,并能解决简化的实际问题。三、实验内容1、主要命令和注意事项：fzero(@f,x0)用来求单变量方程f(x)=0的根
2、，最简单的一种调用方式;一般调用方式为：[x,fv,ef,out]=fzero(@f,xO,opt,pl,p2,...)；fsolve(@f,x0)用于非线性方程组的求解，最简单的一种调用方式；一般调用方式为:[x,fv,ef,out,jac]=fsolve(@f,xO,opt,pl,p2,・..)；2•—些矩阵运算的函数(用来求解线性方程组的数值解)[xD]=eig(A):同时给出A的特征根与特征向量，输出的x是特征根，输出的D是特征向量。[x,刃=lu(A):若A可逆且顺序主子式不为0,输出x为单位下三角阵L
3、,y为上三角阵U,使A=LU；若A可逆，x为一交换阵与单位下三角阵之积。[x,y,p]=lu(A):若A可逆，输岀x为单位下三角阵L,y为上三角阵U,p为一交换阵P,使PA=LU。对正定对称矩阵A的Cholesky分解，输出u为上三角阵U,u=chol(A)使A=U，U。n=norm(x)c=cond(x)c=rcond(x)v=diag(x)v=diag(diag(x))y=triu(x)y=tril(x)y=triu(x,l)y=tril(x,—1)输入x为向量或矩阵，输出为x的2—范数输入x为矩阵，输出为x
4、的2—条件数输入x为方阵，输出为x的条件数的倒数输入向量x,输出v是以x为对角元素的对角阵；输入矩阵X,输出V是X的对角元素构成的向量；输入矩阵X,输出V是X的对角元素构成的对角阵，可用于迭代法中从A中提取D；输入矩阵X,输出V是X的上三角阵；输入矩阵X,输出V是X的下三角阵；输入矩阵x,输出v是x的上三角阵，但对角元素为0,可用于迭代法中从A中提取U；输入矩阵x,输出v是x的下三角阵；但对角元素为0,可用于迭代法中从A中提取L；x,r=roots(c)rref(A)a=sparse(r,c,v,m,n)求极大无
5、关组在第r行、第c列输入数值v,矩阵供m行n列，输岀a为稀疏矩阵，只给出(r,c)及v；aa=full(a)输入稀疏矩阵a,输出aa为满矩阵(包含零元素)】3.实验练习：编程求解下列各题P63：2；3；5四、编写实验报告将实验练习屮的2、3题写在实验报告上，第5题用Word文档形式递交,程序用M■文件的形式..数计学院—学院实验报告专业名称数学与应用数学实验室2栋201实验课程数学建模实验名称线、非线性方程近似解姓名张宇学号100701188同组人员无实验日期2013年3月27日注：实验报告应包含(实验目的，实验
6、原理，主要仪器设备和材料，实验过程和步骤，实验原始数据记录和处理，实验结果和分析，成绩评定)等七项内容。具体内容可根据专业特点和实验性质略作调整，页面不够可附页。一、实验目的用MATLAB数值求解线性代数方程组、数值求解非线性代数方程的迭代法和牛顿法并对作初步分析。二、实验原理1、主要命令和注意事项：fzero(@f,xO)用来求单变量方程f(x)=O的根，最简单的一种调用方式；一般调用方式为:[x,fv,ef,out]=fzero(@f,x0,opt,pl,p2,...)；fsolve(@f,xO)用于非线性方
7、程组的求解，最简单的一种调用方式；—般调用方式为:[x,fv,ef,out,jac]二fsolve(@f,xO,opt,p1,p2,...)；r=roots(c)输入多项式f(x)的系数c(用数组给出,按降幕排列)，输出i•为方程f(x)的全部根；c=poly(r)输入f(x)=0的全部根r(用数组给出),输出c为方程f(x)的系数(按降幕排列);2•—些矩阵运算的函数(用来求解线性方程组的数值解)[xD]=eig(A):同时给出A的特征根与特征向量，输出的x是特征根，输出的D是特征向量。[x,y]=lu(A):
8、若A可逆且顺序主子式不为0,输出x为单位下三角阵L,y为上三角阵U,使A=LU；若A可逆，x为一交换阵与单位下三角阵之积。[x,y,pl=lu(A):若A可逆，输出x为单位下三角阵L,y为上三角阵U,p为一交换阵P,使PA=LUou=chol(A)对正定对称矩阵A的Cholesky分解，输出u为上三角阵U,使A=LFU。n=norm(x)输入x为向量或矩阵，输出为x的2—
