俩种方法tim1=0.007 tim2=2,用矩阵点乘比循环快。
w=rand(4,1000)t1=clock()w1=w.*ww1=sqrt(sum(w1))t2=clock()tim1=etime(t2,t1);shu=zeros([1 1000])t3=clock()for i=1:1000 shu(i)=sqrt(w(1:4,i)'*w(1:4,i))endt4=clock()tim2=etime(t4,t3);