herryzl
2008-01-13, 21:23
比如用QR分解求解逆矩阵invmat
for i=1:3
init=zeros(3,1);
init(i)=1;
vec=r\(q'*init);
if i==1 invmat=vec;
else invmat=cat(2,invmat,vec);
end
end
中间须加一个判断i是否为1的命令，能不能直接创建一个空向量invmat，省掉这个判断语句？
thx!
for i=1:3
init=zeros(3,1);
init(i)=1;
vec=r\(q'*init);
if i==1 invmat=vec;
else invmat=cat(2,invmat,vec);
end
end
中间须加一个判断i是否为1的命令，能不能直接创建一个空向量invmat，省掉这个判断语句？
thx!