138 ATLAB5手册 China-pub.com 下载 然而这样得到的特征值没有用MATLAB命令eig(A)得到的特征值的精度高，而且有效 性也差些： usedEig=eig(A) 给出： usedEig -6.0000 10.0000 4.0000 结果是一样的，但是顺序正好相反。 (g)对于所有的矩阵A都有：po1yva1m(po1y(A),A)。=0 这是Cayley-Hamilton法则。这个法则对于秩为5的方阵来说： Magical magic(5); AlmostZero polyvalm(poly(Magical),Magical) AlmostZero 1.0e-07* 0.2794 0.3551 0.1723 0.1770 0.2654 0.2765 0.2887 0.2049 0.2142 0.2561 0.1775 0.2468 0.2701 0.3073 0.2375 0.1942 0.2744 0.2759 0.2608 0.2282 0.2082 0.3120 0.2608 0.2515 0.2049 ■ 10.2 函数的零值 MATLAB的M文件可以表示数学函数：参见2.9节。函数： 5x-6.4 9x x-0.4 g(x)= x-1.3)2+0.002+x3+0.03-x-0.922+0.005 如果输入下面的M文件g.m,这个函数就可以在MATLAB中调用： function y =g(x) y=(5.*x-6.4)./((x-1.3).2+0.002)+.. (9.*x)./(x.3+0.03)-,., (x-0.4)./((x-0.92).2+0.005); 使用元素运算符*、人、A、+和一定义MATLAB函数g。结果是如果这个函数被一个向量调用， 那么得到的结果也是一个向量。本章中提到的所MATLAB函数需要以这种方式来定义数学函数。 用plot命令可以画出函数的图形： x=linspace(0,2); 号生成向量x plot(x,g(x)); 号画g(x)图形 grid; 冬画格栅 title('The g(x)function 告给出图标题 或者使用fplot命令： p1ot(g,【02]): 号画g(x)图形 grid; 号画格栅 title('The g(x)functidn 号给出图标题然而这样得到的特征值没有用 M AT L A B命令e i g ( A )得到的特征值的精度高，而且有效 性也差些： u s e d E i g = e i g ( A ) 给出： 结果是一样的，但是顺序正好相反。 (g) 对于所有的矩阵A都有：polyvalm(poly(A), A)=0 。 这是C a y l e y - H a m i l t o n法则。这个法则对于秩为5的方阵来说： ■ 10.2 函数的零值 M AT L A B的M文件可以表示数学函数；参见 2 . 9节。函数： 如果输入下面的M文件g . m，这个函数就可以在M AT L A B中调用： 使用元素运算符. *、. /、. ^、+和－定义M AT L A B函数g。结果是如果这个函数被一个向量调用， 那么得到的结果也是一个向量。本章中提到的所有M AT L A B函数需要以这种方式来定义数学函数。 用p l o t命令可以画出函数的图形： x=linspace(0, 2); % 生成向量x p l o t ( x , g ( x ) ) ; % 画g ( x )图形 g r i d ; % 画格栅 t i t l e ('The g(x) function') % 给出图标题 或者使用f p l o t命令： f p l o t ('g', [0 2]); % 画g ( x )图形 g r i d ; % 画格栅 t i t l e ('The g(x) function') % 给出图标题 1 3 8 M ATLAB 5 手册 下载