第3章 MATLAB计算概要.ppt

第3章 MATLAB计算 3.1方程组的求解 3.2 插值和拟合 3.3 函数的极值点 3.4 数值微积分 3.5 符号对象 3.6 符号微积分 3.7 符号方程的求解 3.1方程组的求解 3.1.1 多项式及其运算 多项式的四则运算 加减运算 对多项式加法，MATLAB不提供一个直接的函数。如果两个多项式向量大小相同，标准的数组加法有效。例 a=[1 2 3 4] ; b=[1 4 9 16]; d=a+b d = 2 6 12 20 结果是d(x)= 2x3＋6x2＋12x＋20。当两个多项式阶次不同，低阶的多项式必须用首零填补，使其与高阶多项式有同样的阶次。 p1=[2 0 -6 1]; p2=[1 -1]; p3=[0,0,p2] q=p1+p3 利用子函数mmpadd，适用于任意阶次间多项式相加 function p=mmpadd(a,b) % MMPADD Polynomial addition. % MMPADD(A,B) adds the polynomial A and B % Copyright (c) 1996 by Prentice Hall,Inc. if nargin2 error( Not enough input arguments ) end a=a(:). ; % make sure inputs are polynomial row vectors b=b(:). ; na=length(a) ; % find lengths of a and b nb=length(b) ; p=[zeros(1,nb-na) a]+[zeros(1,na-nb) b] ; % add zeros as necessary 乘法运算 利用函数conv对它们的系数作卷积来实现,p=conv(p1,p2) 除法运算 通过函数deconv进行系数解卷积来实现，[p,r]=deconv(p1,p2)，其中k是商，r是余数。 多项式求导 利用函数polyder来实现，语法格式： poly(p)：返回多项式p的导数。 poly(p1,p2)：返回多项式p1*p2的导数。 [q,d]=poly(p1,p2)：返回多项式p1/p2的导数，q是分子，d是分母。 多项式求值 y=polyval(p,x)：按数组运算规则计算多项式值。 y=polyvalm(p,X)：按矩阵运算规则计算多项式值,且X只能是方阵 Y = P(1)*X^N + P(2)*X^(N-1) + ... + P(N)*X + P(N+1)*I 例： 已知多项式， 分别计算当取1、取在[2,3]区间均匀分布的4个点、取2行4列的正态分布的随机阵、取3行3列均匀分布的随机阵时，多项式的值。 p=[2 0 -6 1]; x1=1;x2=linspace(2,3,4);x3=randn(2,4);x4=rand(3); 多项式求根---利用roots(p)函数来实现 例： 求多项式 的根 3.1.2 线性方程组的求解 可直接通过矩阵相除得到 例3-7 求解线性方程组 3.1.3 非线性方程组的求解 求解单变量非线性方程的根 利用fzero函数，基本的语法格式： x=fzero(fun,x0) 返回函数fun的是离x0最近的根。其中fun的输入方式可以是M函数文件的函数句柄、匿名函数或字符串 例3-8 求 距离最近的根。 求非线性方程组的解 利用fsolve函数，其基本语法格式：x=fsolve(fun,x0) 返回根初始猜测值向量为x0、非线性方程组fun的数值解。 fun可以以M函数文件的函数句柄、匿名函数或字符串方式输入。 求解时还可以在指令中加入options： x=fsolve(fun,x0,options) options是最优化工具箱的选项设定，通过调用optimset函数来设定最大优化代数、优化精度、中间结果的显示等等。 例如Display选项决定显示优化的中间结果的方式，其调用格式为： options=optimset(Display,off或iter或final) 其中off表示不显示，iter表示每步都显示，final表示仅显示最终结果。 例3-9 求解下列非线性方程组在初值为 的数值