资源描述:
《MATLAB--3数值计算》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、第三章数值计算1巨阵的变换;2线性方程(组)的求解;3非线性方程(组)的求解;4MATLAB函数文件;5函数的零点;6函数的极点;矩阵的运算矩阵的转置:’;矩阵的逆:A+=inv(A);矩阵的行列式:d=det(A);求矩阵的对角元素或对角化:x=diag(V,k),x=diag(V)V=diag(x,k),V=diag(x)稀疏矩阵优点:节约存储空间;生成S=sparse(A)S=sparse(i,j,s,m,n,nzmax)S=sparse(i,j,s,m,n)S=sparse(i,j,s)S=sparse(m,n)稀疏矩阵稀疏矩阵还原全矩阵:A=full(S),S为一稀
2、疏矩阵稀疏带状矩阵生成:[B,d]=spdiags(A)B=spdiags(A,d)A=spdiags(B,d,A)Aspdiags(B,d,m,n)矩阵的特征值和特征向量:d=eig(A);[V,D]=eig(A);求解线性方程组线性方程组:Ax=bX=Ab;X=inv(A)*b;求解非线性方程(组)非线性方程(组):f(x)=0函数:X=fsolve(fun,x0)注意:fun为MATLAB函数文件。函数文件创建函数:1函数定义行:function[out1,out2,…]=function_name(in)2函数帮助文件(%开头);3函数体:函数功能实现;function[s,
3、l]=circle(r)s=pi*r^2;l=2*pi*r;函数调用:可以在m文件和命令窗口中直接调用所定义的函数。调用格式与定义时相同:[out1,out2,…]=function_name(in)求解函数的零点一元函数的零点:1思路:先作图求出函数零点的大概区域,再使用函数fzero求解零点。2函数调用格式:x=fzero(fun,x0);[x,fval]=fzero(fun,x0)多元函数零点1思路:先做图确定零点范围,在反复计算逼近真实值;2函数调用格式:x=fsolve(fun,x0);[x,fval]=fsolve(fun,x0);x0是一向量。函数极点一元函数极点:1思路
4、:先做图确定极点大概位置,再用函数fminbnd确定具体位置;2调用格式:[x,fval]=fminbnd(fun,x1,x2);其中,(x1,x2)为区间范围;多元函数极点:无约束条件1单纯形法:[x,fval]=fminsearch(fun,x0);2拟牛顿法:[x,fval]=fminunc(fun,x0);x0为起始猜测点位置;有约束条件:[x,fval]=fmincon(fun,x0,A,b);x0为起始点,A,b为变换约束条件中的系数矩阵和常数矩阵: