第4讲多项式与符号运算资料.ppt

第4讲 多项式与符号运算;提 纲;数据统计;3 求和与求积 数据序列求和与求积函数是sum和prod，其使用方法类似。设x是一个向量，A是一个矩阵，函数的调用格式为： sum(x)：返回向量x各元素之和。 prod(x)：返回向量x各元素的乘积。 sum(A)：返回一个行向量，其第i个元素是A的第i列的元素之和。 prod(A)：返回一个行向量，其第i个元素是A的第i列的元素乘积。 sum(A，dim)：当dim为1时，该函数等同于sum(A)；当dim为2时，返回一个列向量，其第i个元素是A的第i行的元素之和。 prod(A，dim)：当dim为1时，该函数等同于prod(A)；当dim为2时，返回一个列向量，其第i个元素是A的第i行的元素乘积。;4 平均值、标准方差;MATLAB提供了mean，std函数来计算平均值、标准方差或方差。这些函数的调用方法如下： mean(x)：返回向量x的算术平均值。 std(x)：返回向量x的标准方差。 对于矩阵A，mean函数的一般调用格式为： y=mean(A，dim) 这里，dim取1或2。当dim=1时，返回一个行向量y，y的第i个元素是A的第i列元素的平均值；当dim=2时，返回一个列向量y，y的第i个元素是A的第i行元素的平均值。 对于矩阵A，std函数的一般调用格式为： y=std(A，flag，dim) 这里，dim取1或2。当dim=1时，求各列元素的标准方差；当dim=2时，求各行元素的标准方差。flag取0或1，当flag=0时，按计算标准方差；当flag=1时，按计算方差。缺省flag=0，dim=1。;6 排序 对向量元素的进行排序是一种经常性的操作，MATLAB提供了sort函数对向量x进行排序。 y=sort(x)：返回一个对x中元素按升序排列后的向量y。 [y，i]=sort(x)：返回一个对x中的元素按升序排列的向量y，而i记录y中元素在x中的位置。;多项式;多项式的四则运算;2 多项式的乘法运算;3 多项式除法;多项式的求导;多项式的求值和求根;roots函数用来求代数多项式的根， 其调用格式为： x=roots(p) 如果x为向量，则p=poly(x)可以建立一个以x为其根的多项式。;MATLAB还提供了一个fzero函数，可以用来求单变量非线性方程的求根。该函数的调用格式为： z=fzero(‘fname’,x0) 其中fname是待求根的函数文件名，x0为搜索的起点。一个函数可能有多个根，但fzero函数只能给出离x0最近的那个根。;函数极值;建立函数文件f.m。 function w=f(p) x=p(1); y=p(2); z=p(3); w=x+y^2/(4*x)+z^2/y+2/z; 调用fminsearch函数求多元函数在[1/2,1/2,1/2]附近的最小值点。 w=fminsearch(f ,[1/2,1/2,1/2]) w = 0.5000 1.0000 1.0000 计算多元函数的最小值。 f(w) ans = 4.0000;MATLAB符号计算;应用sym函数还可以定义符号变量，使用符号变量进行代数运算时和数值常量进行的运算不同。下面的命令用于比较符号常量与数值常量在代数运算时的差别 ;（2）syms函数 函数sym一次只能定义一个符号变量，使用不方便。MATLAB提供了另一个函数syms，一次可以定义多个符号变量。syms函数的一般调用格式为： syms 符号变量名1 符号变量2 … 符号变量n 用这种格式定义符号变量时，变量间用空格而不要用逗号分隔。例如，用syms函数定义4个符号变量a,b，命令如下： syms a b;2 建立符号表达式;（3）利用已经定义的符号变量组成符号表达式。 例如 syms x y; z=3*x^2-5*y+2*x*y+6 z = 3*x^2-5*y+2*x*y+6;3 符号表达式中变量的确定 利用函数findsym(s)可以确定符号表达式s中的全部符号变量。例如： syms a b x y; %定义4个符号变量 c=sym(3); %定义1个符号常量 s=3*x+y; findsym(s) ans = x, y findsym(5*x+2) ans = x findsym(a*x+b*y+c) %符号变量c不会出现在结果中 ans = a, b, x, y;符号表达式运算;2 符号表达式的因式分解与展开 符号表达式的因式分解和展开运算，可用函数factor和expand来实现，其调用格式为： factor(s)：对符号表达式s分解因式。 expand(s)：对符号表达式s进行展开。 例如： syms x y; s1=x^3-6*x^