polyval polyvalm

polyval 只是计算一个多项式，这个你应该明白吧，help里的例子很清楚：

 

就是把x值带进去。


polyvalm比较特别的。但是明白这个之前，你需要明白什么是 characteristic polynomial:

In  linear algebra , one associates a  polynomial  to every  square matrix , its  characteristic polynomial . This polynomial encodes several important properties of the  matrix , most notably its  eigenvalues , its  determinant  and its  trace .


Suppose we want to compute the characteristic polynomial of the matrix

  We have to compute the determinant of

  and this determinant is

  The latter is the characteristic polynomial of  A .

现在回到我们自己的问题， 对于Y = polyvalm(p,X)， 他的计算方法是：

Y = P(1)*X^N + P(2)*X^(N-1) + ... + P(N)*X + P(N+1)*I

如果写出for循环是：

for i = 1:np %多项式长度，m是X的维数
    Y = X * Y + diag(p(i) * ones(m,1));
end



这个循环比较简单，不要我解释吧？

举个例子给你看，我把循环拆开：

1. 
   
2. p=[1 2 3]
   
3. X=[1 2; 3 4]
   
4. np = length(p);
   
5. [m,n] = size(X);
   
6. 
   
7. Y = zeros(m,m)
   
8. i = 1
   
9.     Y = X * Y + diag(p(i) * ones(m,1))
   
10. 
    
11. i = 2
    
12.     Y = X * Y + diag(p(i) * ones(m,1))
    
13.     
    
14.     
    
15. i = 3
    
16.     Y = X * Y + diag(p(i) * ones(m,1))

复制代码

你再运行：

p=[1 2 3]
X=[1 2; 3 4]
polyvalm(p,X)

看看，结果是不是一样的。

总结一下：
polyval(p, t)是计算出来：
p(1)*t^n + p(2)*t^(n-1) + .... + p(n)
polyvalm(p, X)是计算出来：
p(1)*X^n + p(2)*X^(n-1) + .... + p(n)*I
这里要注意最后的I就是单位阵，不要忘了，否则会出错，比如：
>> p=[1 2 3];
>> X=[1 2; 3 4];
>> polyvalm(p, X)

ans =

    12    14
    21    33

>> p(1)*X*X+p(2)*X+p(3)

ans =

    12    17
    24    33

>> p(1)*X*X+p(2)*X+p(3)*eye(2, 2)

ans =

    12    14
    21    33

多谢各位，认真研究了下帖子，终于明白了，就是将矩阵X整体带入多项式中，求矩阵的多项式值:lol