for i=1:a;
for j=1:b;
m=i*Cos(i)+jSin(j);
A(i,j)=A(i,j)+m;
end
end
a=5;b=6;
A=zeros(5,6);
for i=1:a;
for j=1:b;
m(i,j)=i*cos(i)+j*sin(j);
A(i,j)=A(i,j)+m(i,j);
end
end
A
%向量化
[x,y]=meshgrid(1:a,1:b);
A=x.*cos(x)+y.*sin(y);
A=A'
for k=1;m
[i,j]=meshgrid(1:a,1:b);
L=(i-1)*cos(k)+(j-1)*sin(k);
n = fix(L);
cL=L-n;
A= A+(1-cL).*B(n)+cL.*B(n+1);
end