第二章 矩阵的数值运算4课时..pptVIP

4. conv多项式乘运算 例: a(x)=x2+2x+3; b(x)=4x2+5x+6; 求：c(x)=a(x)*b(x) a=[1 2 3];b=[4 5 6]; c=conv(a,b) %计算两多项式相乘后系数 c = 4 13 28 27 18 p=poly2str(c,x) p = 4 x^4 + 13 x^3 + 28 x^2 + 27 x + 18 5. deconv多项式除运算 例： p1=conv([1 0 2],conv([1 4],[1 1])); p2=[1 0 1 1]; [q r]=deconv(p1,p2) q = 1 5 r = 0 0 5 4 3 cq=‘商多项式’；cr=‘余多项式’； disp([cq,poly2str(q,’s’)]),disp([cr,poly2str(r,’s’)) 商多项式 s + 5 余多项式 5 s^2 + 4 s + 3 6. 多项式微分 matlab提供了polyder函数求多项式的微分。 命令格式： polyder(p): 求p的微分 例： a=[1 2 3 4 5]; poly2str(a,’x’) ans = x^4 + 2 x^3 + 3 x^2 + 4 x + 5 b=polyder(a) b = 4 6 6 4 poly2str(b,’x’) ans = 4 x^3 + 6 x^2 + 6 x + 4 2.5 函数的极值 1. 固定区间单变量函数的最小值 利用fminbnd函数求固定区间单变量函数的最小值。用法： fun：单变量函数的函数名，[x1 x2]为给定的计算区间，[x,fval]为最小值点的位置和最小值。 例：求 在区间[0，2]的最小值 f=inline(x.^3-2*x.^2-5,x);％采用内联函数表示测试函数 x=0:0.1:2;plot(x,f(x)) [x,y]=fminbnd(f,0,2) x = 1.3333 y = -6.1852 或：f=x.^3-2*x.^2-5 2. 固定区间多变量函数的最小值 利用fminsearch函数求固定区间多变量函数的最小值。用法： 从函数fun的给定点开始搜索附近的最小值点x，fval为最小值点对应的最小值。 例：求 在区间x=0.3,y=0.7附近的最小值点。 f=inline(‘p(1).^2+p(2).^2+7’,‘p’); ％采用内联函数表示测试函数 [minpoint,minval]=fminsearch(f,[0.3,0.7]) minpoint = 1.0e-004 * -0.1766 -0.1636 minval = 7.0000 2.6 代数方程求解 求非线性方程式的根的步骤： 定义方程式（ 形式），可用inline函数定义； 代入适当范围的x及其对应f(x)值，画出方程函数曲线，了解该方程式根的大致位置； 由图中决定f(x)与x轴相交的大致位置x0，通过fzero函数求出在 x0附近的根。 例：求方程 f=inline(2*exp(x).*sin(2*pi*x)-0.5,x)； x=0:0.1:2;y=f(x);plot(x,y) r=fzero(f,0.5) r = 0.4752 r=fzero(f,1) r = 1.0144 r=fzero(f,1.5) r = 1.4910 2.7 数值积分与微分 数值积分基本原理 求解定积分的数值方法多种多样，如简单的梯 形法、辛普生(Simpson)法、牛顿－柯特斯(Newton- Cotes)法等都是经常采用的方法。 它们的基本思想都是将整个积分区间[a,b]分成 n个子区间[xi,xi+1]，i=1,2,…,n，其中x1=a，xn+1=b。 这样求定积分问题就分解为求和问题。 数值积分的实现方法 变步长辛普生法 基于变步长辛普生法，MATLAB给出了quad函数来求定积分。该函数的调用格式为： [I n]=quad(fname,a,b) 其中fname是被积函数名。a和b分别是定积分的下限和上限。返回参数I即定积分值。n为被积函数的调用次数。 例：