(3)多项式求值(polyval): 将多项式a中的自变量x 赋予值xv时，
该多项式的值可用
f=polyval (a，xv)
求得，其中xv可以是复数，而且可以是矩阵或 数组，此时polyval对输入变元作元素群运算， 这对于求线性系统的频率特性特别方便。 polyvalm 则对输入的变元阵方阵作矩阵多项式 运算。
【例4-3-2】 设a为系统分母系数向量，b为系统分子系数向量，求此系统的频率响应并画出频率特性。先令频率数组w取线性间隔: w=linspace(0,10);
%在w等于0~10之间按线性间隔取100点(默认值) a=polyval(a，j*w);
%分别求分母分子多项式的值(为复数数组) b=polyval(b，j*w); subplot(2,1,1); plot(w，abs(b./a))
%画两者元素群相除所得的幅频特性 subplot(2,1,2); plot(w，angle(b./a))
%画相频特性
频率特性通常在对数坐标中绘制。
因此输入频率数组取对数等间隔: w1=logspace(-1,1) %在w1从10-1~10之间，按对数分割为 50点(默认值) f=polyval(b,j*w1)./polyval(a，j*w1); %求出这些点上的频率响应(复数) subplot(2,1,1)，loglog(w1，abs(f)) %在双对数坐标中画出幅频特性 subplot(2,1,2); semilogx(w1，angle(f)) %在双对数坐标(x)中画出相频特性 所得曲线如图4-3所示。
