《MATLAB程序设计基础教程》课件第4章.ppt

6．乘法运算重载

多项式的乘法实际上可以表示为系数向量的卷积，可以由conv()函数直接获得。故可以如下重载定义多项式的乘法运算：

functionp=mtimes(a,b)

a=polynom(a);b=polynom(b);

p=polynom(conv(a.c,b.c));7．乘方运算重载

多项式的乘方运算只限于正整数乘方的运算，其n次方相当于将该多项式自乘n次。若n=0，则结果为1。这样我们就可以重载定义多项式的乘方运算：

functionp=mpower(a,n)

ifn=0,n=floor(n);a=polynom(a);p=1;

ifn=1,

fori=1:n,p=p*a;end

end

else,error(乘方运算只限于正整数！)

end8．多项式求值重载

可以对多项式求值函数polyval()进行重载定义：

functiony=polyval(a,x)

a=polynom(a);

y=polyval(a.c,x);9．根函数重载

可以对多项式求根函数roots()进行重载定义：

functionr=roots(a)

%roots(obj)returnsavectorcontainingtherootsofa

r=roots(a.c);

end10．微分函数重载

可以对多项式求根函数roots()进行重载定义：

functionr=roots(a)

functionq=diff(a)

%diff(a)isthederivativeoftheDocPolynoma

C=a.c;

d=length(C)-1;%degree

q=Polynom(a.c(1:d).*(d:-1:1));

end11．绘图函数重载

可以对多项式求绘图函数plot()进行重载定义：

functionplot(a)

%plot(obj)plotstheDocPolynomobj

r=max(abs(roots(a)));

x=(-1.1:0.01:1.1)*r;

y=polyval(a,x);

plot(x,y);

title([y=char(a)])

xlabel(X)

ylabel(Y,Rotation,0)

gridon

end将上述所有重载函数和引导函数polynom()都保存在@polynom子目录中。

使用methods()函数可以列出刚建的新类和已经定义的方法函数名。例如：

methods(polynom)

Methodsforclasspolynom:

chardisplayminusmtimespluspolyval

diffdoublempowerplotpolynomroots

定义了这些类之后，我们就可以方便地进行多项式处理了。例4-4-1构建重载函数的多项式对象。

建立两个多项式对象：P(s)=x^4-1.2x^3+3x^2-6和Q(s)=4x^4+3x^3-5x^2+6x+10。

其MATLAB语句如下：

P=polynom([1,-1.2,3,0,-6]),Q=polynom([4,3,-5,6,10])

P=x^4-1.2*x^3+3*x^2-6

Q=4*x^4+3*x^3-5*x^2+6*x+10

然后调用重载函数就可以得出相应的计算结果。在函数tic()中，变量TICTOC说明为全局的，因此它的值由调用函数clock()来设定。以后在函数toc()中，变量TICTOC也说明为全局的，让toc()访问存储在TICTOC中的值。利用这个值，toc()计算自执行函数tic()以来消失的时间。值得注意的是，变量TICTOC存在于tic()和toc()的工作空间，而不在MATLAB的工作空间。(2)实际编程中，无论什么时候应尽量避免使用全局变量。如果要使用全局变量，建议全局变量名要长，它应包含所有要使用该全局变量的函数名的大写字母，并选择以首次出现的M文件的名字开头。例如上面的M函数文件tic()首次出现，全局变量声明为：globalTICTOC。如果遵循这些建议