小木虫

当前位置：首页 > 程序语言 >MATLAB编程遇问题

MATLAB编程遇问题

作者 zhangyuque: 来源: 小木虫 400 8 举报帖子

如何将下面的程序用向量的方式表示，不用for循环。望大神指教！
for k1=1:M
for k2=1:N
      for k3=1:N1
      a2(k1,k2)=a2(k1,k2)+((x(k1)-r*cos(2*pi*k3/N1))^2+ (y(k2)-r*sin(2*pi*k3/N1))^2+Z^2)^(-11);
      end
end
end 返回小木虫查看更多

今日热帖

精华评论

zhangyuque

引用回帖:
2楼: Originally posted by super929 at 2018-01-05 16:02:24
k1=
? ? k2=
? ?? k3=
? ?? ???a2(k1,k2)=a2(k1,k2)+((x(k1)-r*cos(2*pi*k3/N1))^2+ (y(k2)-r*sin(2*pi*k3/N1))^2+Z^2)^(-11);?

矩阵维度不一致，在这里M=100;N=100;N1=4;
?�?懳实�

k1=1:100;
k2=1:100;
k3=1:4;
r=1;
[X,Y,K3] = meshgrid(k1,k2,k3);
a2=zeros(length(k1),length(k2),length(k3));
Z=ones(length(k1),length(k2),length(k3));
a1=a2+((X-r*cos(2*pi*K3/max(k3))).^2+(Y-r*cos(2*pi*K3/max(k3))).^2+(Z).^2).^-11;
a=sum(a1,3);

猜不出来你的确定的想法，所以就按我自己的想法给你写了一个作为参考吧！
zhangyuque

引用回帖:
5楼: Originally posted by ?�?懳实� at 2018-01-07 20:02:10
k1=1:100;
k2=1:100;
k3=1:4;
r=1;
= meshgrid(k1,k2,k3);
a2=zeros(length(k1),length(k2),length(k3));
Z=ones(length(k1),length(k2),length(k3));
a1=a2+((X-r*cos(2*pi*K3/max(k3))).^2+(Y-r*cos(2*pi ...

谢谢您啦
zhangyuque

引用回帖:
5楼: Originally posted by ?�?懳实� at 2018-01-07 20:02:10
k1=1:100;
k2=1:100;
k3=1:4;
r=1;
= meshgrid(k1,k2,k3);
a2=zeros(length(k1),length(k2),length(k3));
Z=ones(length(k1),length(k2),length(k3));
a1=a2+((X-r*cos(2*pi*K3/max(k3))).^2+(Y-r*cos(2*pi ...

倒数第二句的时候就说矩阵维度不一致
，