主题:[讨论]哪位高手帮我看看这段程序,怎么修改,才不出错.谢谢!!
unnamed=[1,2,2,3;3,3,2,1]
err=0.0000001;
w=zeros(size(unnamed,1),1);
w(1)=1;
r=2;
k=1;
while(abs(abs(r)-1)>err)
k=k+1;
w(:,k)=unnamed*(tanh(w(:,k-1)'*unnamed))-tanh(w(:,k-1)'*unnamed)*w(:,k-1);
w(:,k)=w(:,k)/norm(w(:,k));
r=w(:,k)'*w(:,k-1);
end
q(:,1)=w(:,k);
k
[em18]
err=0.0000001;
w=zeros(size(unnamed,1),1);
w(1)=1;
r=2;
k=1;
while(abs(abs(r)-1)>err)
k=k+1;
w(:,k)=unnamed*(tanh(w(:,k-1)'*unnamed))-tanh(w(:,k-1)'*unnamed)*w(:,k-1);
w(:,k)=w(:,k)/norm(w(:,k));
r=w(:,k)'*w(:,k-1);
end
q(:,1)=w(:,k);
k
[em18]