矩阵求极限算法是否正确

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

A矩阵求极限如下：
M=limn→∞(A)^n


算法这样弄可以吗？为什么我得不出来结果，一直在运行
clear
syms A n;
A=[0.1786,0.2205,0.2366,0.2205;
    0.2496,0.1985,0.2563,0.2567;
    0.3221,0.3243,0.2509,0.3243;
    0.2496,0.2567,0.2563,0.1985];
M=limit(A^n,n,Inf);

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：limit Lim CLE Lea IMI MATLAB 矩阵 极限

求特征值给出以下结果，说明这个极限不可能存在
ans =
   1.000008538586829
  -0.043025817257730
  -0.072282721329098
  -0.058200000000000
方法也不正确，矩阵幂的极限一般使用特征值来求，具体可参考常见线性代数书籍中关于矩阵特征值的部分，我编了以下函数实现

function M = matlim(A)
%compute limit of matrix exponential A^n by eig, make sure A be square
norma = norm(A);
if norma>1
   disp('norm larger than 1, no limit exist')
   return
end

[U,eiga] = eig(A);

eiglim = round(eiga);

Uinv = U\eye(size(A));

M = U*diag(eiglim)*Uinv;
end