自己设计了一个求解，矩阵权向量的算法，希望能做到更优化

thinking_—zeng

clc
clear
9 J7 b: y5 G8 m4 `- ~a=input('请输入矩阵:');
! }/ a' B, k% A7 O' Ufor j=1:1:length(a)
9 s9 I/ L" |1 `4 N    b(j)=sum(a(:,j));
! @3 ^0 D' y5 f' h    w(:,j)=a(:,j)/b(j);
end
2 t2 n7 `0 p2 n0 |" R* J% h" Sfor i=1:numel(a)/length(a)
% E+ \' J1 R  o+ S7 E    c(i,1)=sum(w(i,:));
end
# E- S9 U9 y) [) E: x+ qd=sum(c);
0 a1 Q+ ]: C5 c+ i, N; ^for i=1:1:numel(a)/length(a)
6 d  M* r% U6 ]; Sw1(i)=c(i)/d;
end
m=w1'

madio

是采用的和法吗？可以使用幂法

thinking_—zeng

madio 发表于 2013-7-28 17:30
是采用的和法吗？可以使用幂法

7 Z& @; B* K+ L4 Y是的，不会冥法，