为什么循环的结果都一样

aid4her

G=6.67*10^(-11);
r=100;
h=500;
0 ?6 D6 V9 c8 \/ K! Ym=0.1*(4/3)*pi*r^3;
7 _; Z! K# C9 _/ ag=zeros(41,41);
6 A  F+ l# H8 B1 \: j! Rfor x=-1000:50:1000
' I+ w! _; Y8 G4 Y2 k. |    for y=-1000:50:1000
& J) J% S! X+ N6 `2 q7 k        g(1:41,: )=G*m*h/((x^2+y^2+h^2)^(3/2));
    end
end

R_myself

你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
. u+ y* D+ S. c/ W9 B1 R" T* T" Qx=-1000:50:1000;
7 O/ C- _6 D- D& W( Ay=-1000:50:1000;
% ]* Q0 H* H# R8 P4 K8 G* Efor i=1:41
" {; V! d6 Y5 Z( F. l. G    for j=1:41
" {2 N% V" a% O( M9 Z" T6 x        g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
/ T4 g0 o1 n5 n% ~3 C. Y8 ^2 _) z3 ]    end
end
8 W) v% S( _, A; C# }复制代码
: L' \/ K1 e2 m: F看一下这是不是你要的结果，如果不是请自己理解并修改一下吧。

fuu65iwi

这严格来说就不是MATLAB问题，是你循环逻辑写错了，你在循环中没有控制变量，在不停地重复相同的赋值。打个比方，裁缝店来了高个子客人应该给长衣服，来了胖客人应该给宽衣服，裁缝店这个提供衣服的循环，应当根据控制变量也就是客人的具体身材，进行控制。可你现在这个循环，实质上是不管来了什么身材的客人，统统扔出来41*41套一模一样的衣服。
, C6 M, H, v  R" w改成
xSet = -1000:50:1000;
ySet = -1000:500:1000;
  f% ]! [7 ]) w4 y- f6 M7 q( ]for index_row = 1 : 1 : numel( xSet )
    x = xSet( index_row );
    for index_column = 1 : 1 : numel( ySet )
8 }0 l# [- E0 _$ ]* s+ C; r8 {( S- d% {/ N        y = ySet( index_column );
/ o7 b& a) e* I1 A" [7 s        g( index_row, index_column ) = G * m * h / ( ( x^2 + y^2 + h^2 ) ^ ( 3/2 ) );
) k7 K  Y7 b1 e/ V! P    end
  o2 u8 K2 t% G; c( o3 @1 _end

I_believe

这么改也可以：
clear;clc
3 H, Q- O, E4 j4 X: K
G=6.67*10^(-11);
- h4 `8 Q; p. l$ m$ s/ w) or=100;
# Y8 b% N, M, [& t- E: A1 d3 `h=500;
m=0.1*(4/3)*pi*r^3;
. k% |$ @+ m0 @  U$ Z, @g=zeros(41,41);

! v, l2 N- y+ G% L* Gi = 0;
for x=-1000:50:1000
    i = i + 1;
: Q" v- K" K" v" `% |" V# E6 a8 b1 u    y = -1000:50:1000;
8 Z! G3 g$ [+ a% h& e& g    g(i,: )=G*m*h./((x^2+y.^2+h^2).^(3/2));
end