回復(fù)武漢某大學(xué)老師的一個(gè)問題:北太天元求根和范德蒙行列式

%北太天元求行列式給出的多項(xiàng)式的根
%?det(?[1 1 1; 2 4 x; 3 9 x^2 ]?) = 0 的根
% 實(shí)際上有范德蒙行列式的知識(shí),我們知道 x = 2 或者 3
% 但是我們假設(shè)沒有學(xué)過范德蒙矩陣的知識(shí),
A = [ 1 1 ;
?????????2?3 ;
?????????4?9 ; ]
???c = zeros(3,1);
???for j=1:3
??????if j==1
?????????ind = [2,3];
??????elseif j==2
???????????????ind = [1, 3];
??????else
????????????????ind = [1,2];
??????end
??????c(j) = (-1)^(j+3) * det( A( ind,:));
???end
???%p(x) = c(1) +c(2)*x+c(3)*x^2
??% 但是使用roots的時(shí)候要求輸入的系數(shù)是按降冪排列
???c = c(end:-1:1);
???[r,a]?= roots(c)
???% roots 是 a 的特征值
???[v,d] = eig(a);
???xx = diag(d)
標(biāo)簽: