3matlab数值计算.ppt

【例】求f(x)?=?e^ (?2x?)??x?=?0在x0?=?0附近的根。 步骤如下： （1）建立函数文件funx.m。 function fx=funx(x) fx=exp(-2*x)-x; （2）调用fzero函数求根。 z=fzero(@funx,0.0) z = 0.4263 也可以直接使用表达式字符串求解： z=fzero(exp(-2*x)-x,0.0) 求解非线性方程组需要用到MATLAB的优化工具箱（Optimization Toolbox）。对于非线性方程组F(X)=0，用fsolve函数求其数值解。fsolve函数的调用格式为 X=fsolve(@fun,X0,option) 其中X为返回的解，fun是用于定义需求解的非线性方程组的函数文件名，X0是求根过程的初值，option为最优化工具箱的选项设定。最优化工具箱提供了20多个选项，用户可以使用optimset命令将它们显示出来。 如果想改变其中某个选项，则可以调用optimset()函数来完成。例如，Display选项决定函数调用时中间结果的显示方式，其中off为不显示，iter表示每步都显示，final则只显示最终结果。optimset(Display,off)将设定Display选项为off。 非线性方程组的求解 步骤如下： （1）建立函数文件myfun.m function q=myfun(p) x1=p(1); x2=p(2); q(1)=x1^2+x1-x2^2-1; q(2)=x2-sin(x1^2) （2）在给定的初值 [0.5,0.5]下,调用fsolve函数求根 X=fsolve(@myfun,[0.5,0.5],optimset(Display,off)) X = 0.7260 0.5029 本 章 小 结 特殊矩阵的生成 （1）通用的特殊矩阵（5个） （2)面向特定应用的特殊矩阵（6个） 矩阵分析 矩阵结构变换 对角阵(diag)、三角阵（triu,tril)、矩阵转置（‘）、 矩阵旋转（rot90(A,k))、矩阵的翻转（fliplr，flipud，flipdim） 矩阵求值 （1）方阵的行列式值（det) （2）矩阵的秩与迹（rank,trace) (3) 向量和矩阵的范数（norm(A,1),norm(A),norm(A,inf)) （4）矩阵的条件数（cond(A,1),cond(A),cond(A,inf)) 矩阵的特征值与特征向量（eig(A),[V,D]=eig(A),[V,D]=eig(A,’nobalance’)) * 本章小结 线性方程组求解 （1）矩阵求逆（inv）及利用矩阵求逆方法解线性方程组； （2）利用左除运算符求解线性方程组 矩阵分解 （1）LU分解； （2）QR分解； （3）Cholesky分解； （4）Schur分解 （5）矩阵的奇异值分解； （6）矩阵的QZ分解等 * 本 章 小 结 稀疏矩阵的处理 （1）矩阵的存储方式（完全存储，稀疏存储） （2）矩阵的稀疏存储方式（sparse(m,n);sparse(u,v,S);[u,v,S]=find(A);full(A)) （3）产生稀疏存储矩阵（spconvert(A)) （4）带状稀疏存储矩阵（spdiags(B,d,m,n)) （5）单位矩阵的稀疏存储（speye(m,n)) 数据统计处理 （1）最大值和最小值（max(A),[Y,U]=max(A),max(A,[],dim))，min函数与max函数调用格式相同； （2）两个向量或矩阵对应元素的比较（max(A,B);max(A,n)),min函数与max调用格式相同； * 本 章 小 结 （3）求和与求积（sum(A);sum(A,dim);prod(A),prod(A,dim)) （4）平均值与中值（mean(A),mean(A,dim);median同mean) （5）累加和与累乘积（cumsum(A),cumsum(A，dim),cumprod同cumsum) （6）标准方差与相关系数（std(A,flag,dim);corrcoef(X,Y),corrcoef(X)) （7）排序（sort(A,dim,mode)) 多项式计算 （1）多项式四则运算（加减，乘除（conv(p1,p2),[Q,r]=deconv(p1,p2)); （2）多项式的导数（p=polyder(p),p=polyder(P,Q),[p,q]=polyder(P,Q)) 本 章 小 结 （3