我用sym定義了變量x,然后將x和矩陣進(jìn)行數(shù)乘,但是報(bào)錯(cuò),coffes(det( A - x * eye(n) ), x)報(bào)錯(cuò)不支持該類型輸入
這個(gè)問題可能是由于符號表達(dá)式未正確展開或函數(shù)參數(shù)傳遞格式有問題導(dǎo)致的,但還需看具體代碼確定原因和給解決辦法,方便給個(gè)代碼嗎
引用:這個(gè)問題可能是由于符號表達(dá)式未正確展開或函數(shù)參數(shù)傳遞格式有問題導(dǎo)致的,但還需看具體代碼確定原因和給解決辦法,方便給個(gè)代碼嗎
function [V_Eig, D_Eig] = Eig(A)
n = size(A, 1);
lambda = sym('lambda');
coef = coeffs(det(A - lambda * eye(n)), lambda);
V_Eig = zeros(n, n);
D_Eig = roots(coef);
for i = 1:size(A)
V_Eig(:, i) = null(A - D_Eig(i) * eye(n));
end
end
不支持的輸入類型
位于文件 C:\baltamatica\腳本\Eig1.m (第 12 行)
coef = coeffs(det(A - lambda * eye(n)), lambda)
引用:function [V_Eig, D_Eig] = Eig(A)n = size(A, 1);lambda = sym('lambda');coef = coeffs(det(A - lambda * eye(n)), lambda);V_Eig = zeros(n, n);D_Eig = roots(coef);for i = 1:size(A)V_Eig(:, i) = nul...
下一個(gè)版本已經(jīng)重載了eig函數(shù),不需要自己寫Eig.