为什么循环的结果都一样

aid4her

G=6.67*10^(-11);
# b% D, y) `; o  b! w" l3 nr=100;
h=500;
m=0.1*(4/3)*pi*r^3;
! k. Z/ J! L5 Vg=zeros(41,41);
for x=-1000:50:1000
7 {1 V& W4 h' f# \9 c: n    for y=-1000:50:1000
- ~* c. S* c/ f% y2 D3 V  N, Q        g(1:41,: )=G*m*h/((x^2+y^2+h^2)^(3/2));
1 @5 y, ^; b2 H& j' K0 Q& ~, v# W    end
end
. b3 v7 r, u4 |+ b

R_myself

你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
x=-1000:50:1000;
y=-1000:50:1000;
for i=1:41
    for j=1:41
- n/ E) Q2 y  B0 e        g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
5 p0 L6 S. G+ I7 `, i    end
end
3 }$ g& c* K5 ]$ _" e% x5 z复制代码
看一下这是不是你要的结果，如果不是请自己理解并修改一下吧。

fuu65iwi

这严格来说就不是MATLAB问题，是你循环逻辑写错了，你在循环中没有控制变量，在不停地重复相同的赋值。打个比方，裁缝店来了高个子客人应该给长衣服，来了胖客人应该给宽衣服，裁缝店这个提供衣服的循环，应当根据控制变量也就是客人的具体身材，进行控制。可你现在这个循环，实质上是不管来了什么身材的客人，统统扔出来41*41套一模一样的衣服。
+ j7 g! f8 l1 o# |改成
8 `) r1 L# I4 T5 i5 }* H. a! oxSet = -1000:50:1000;
ySet = -1000:500:1000;
) ]; r3 `9 g, Z0 X* q) v; jfor index_row = 1 : 1 : numel( xSet )
    x = xSet( index_row );
    for index_column = 1 : 1 : numel( ySet )
        y = ySet( index_column );
- ]0 ]% P' T1 q% `4 w        g( index_row, index_column ) = G * m * h / ( ( x^2 + y^2 + h^2 ) ^ ( 3/2 ) );
. l- B, J2 {+ V# A    end
end

I_believe

这么改也可以：
clear;clc
4 O" I. Z4 v6 Z( [( @
5 T  P3 E& a1 ]5 c) M, z1 K/ QG=6.67*10^(-11);
r=100;
h=500;
m=0.1*(4/3)*pi*r^3;
3 q( Y; L/ D+ I  _g=zeros(41,41);
9 ^+ t7 N! f5 N! E( n4 G6 w3 |
8 D; A7 y# \7 T1 P5 Mi = 0;
' i1 w/ N5 h+ Z# J$ O! Nfor x=-1000:50:1000
    i = i + 1;
4 O- n% ~8 x  y    y = -1000:50:1000;
    g(i,: )=G*m*h./((x^2+y.^2+h^2).^(3/2));
end