为什么循环的结果都一样

aid4her

G=6.67*10^(-11);
r=100;
h=500;
$ g  B9 L( N4 k4 @6 f5 g1 g4 O; ym=0.1*(4/3)*pi*r^3;
6 M! E1 n& W/ L3 C9 b( Mg=zeros(41,41);
for x=-1000:50:1000
2 M2 q  V& {! r    for y=-1000:50:1000
, j' V& R; P! c! _% ^" C, I        g(1:41,: )=G*m*h/((x^2+y^2+h^2)^(3/2));
2 S6 I1 U& \, D8 c    end
end

R_myself

你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
# ?4 \7 Z- g9 }2 Q! S4 E8 F. Z2 Lx=-1000:50:1000;
* g- n: G' B' ^$ _y=-1000:50:1000;
' b3 }8 j- R; Tfor i=1:41
    for j=1:41
        g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
7 T7 g4 V( N( H" h2 ]1 q# {    end
: Y+ |+ N% ?( H) K2 ]  p: j; o7 B8 O! send
复制代码
看一下这是不是你要的结果，如果不是请自己理解并修改一下吧。

fuu65iwi

这严格来说就不是MATLAB问题，是你循环逻辑写错了，你在循环中没有控制变量，在不停地重复相同的赋值。打个比方，裁缝店来了高个子客人应该给长衣服，来了胖客人应该给宽衣服，裁缝店这个提供衣服的循环，应当根据控制变量也就是客人的具体身材，进行控制。可你现在这个循环，实质上是不管来了什么身材的客人，统统扔出来41*41套一模一样的衣服。
改成
xSet = -1000:50:1000;
ySet = -1000:500:1000;
% a$ @8 @+ y- k* H2 P0 U7 B& vfor index_row = 1 : 1 : numel( xSet )
    x = xSet( index_row );
    for index_column = 1 : 1 : numel( ySet )
1 @0 x$ Z, x; L% s  {2 X        y = ySet( index_column );
        g( index_row, index_column ) = G * m * h / ( ( x^2 + y^2 + h^2 ) ^ ( 3/2 ) );
    end
( [% F  g1 B0 |' Zend

I_believe

这么改也可以：
7 B% E& g0 H0 T4 _/ p+ Eclear;clc
4 w* T( P  z* s
G=6.67*10^(-11);
r=100;
h=500;
m=0.1*(4/3)*pi*r^3;
$ A$ R6 t+ _. K" ?' q" fg=zeros(41,41);
+ U, \/ Z9 n+ T' _) D( ^
i = 0;
for x=-1000:50:1000
6 B- ^; E( p& a9 t    i = i + 1;
    y = -1000:50:1000;
8 _0 m/ a. i  y0 q    g(i,: )=G*m*h./((x^2+y.^2+h^2).^(3/2));
end