gusucode.com > matlab通信工程仿真源码(张德丰等编著)程序书籍 > matlab_code/matlab通信工程仿真源码(张德丰等编著)/第3章/li3_38.m
A = rand(3); E = 0*A; F = eye(size(A)); N = 1; while norm(E+F-E)>0 E = E+F; F = A*F/N; N = N+1; end