clc;clear all
A=rand(200,300);
[r1,r2]=size(A);
ku=rand(1,r2);
ko=rand(r1,1);
w=zeros(r1);
w1=zeros(r1);
tic
for i=1:r1
    for j=1:r1
        s=0;
        for k=1:r2
            s=s+(A(i,k)*A(j,k)/ku(k));
        end
        w(i,j)=s/ko(j);
    end
end
toc
tic
for i=1:r1
    for j=1:r1
        s=sum(A(i,:).*A(j,:)./ku);
        w1(i,j)=s/ko(j);
    end
end
toc
isequal(w,w1)
运行结果:
Elapsed time is 2.092196 seconds.
Elapsed time is 0.169971 seconds.
ans =
     1