EDA365电子论坛网
标题:
为什么循环的结果都一样
[打印本页]
作者:
flamesmnm
时间:
2023-1-12 10:52
标题:
为什么循环的结果都一样
G=6.67*10^(-11);
% l: ~7 Q8 Y; t# ], G' Z; D) V1 D
r=100;
4 D5 i- q' T9 f% Z6 V& u s) T
h=500;
7 ` C0 h3 l9 s$ l8 l1 |
m=0.1*(4/3)*pi*r^3;
) S) s7 ~$ l' V' x( A2 J3 a9 q: N
g=zeros(41,41);
8 ?6 O1 Q' G% c, G3 q1 S3 g
for x=-1000:50:1000
/ s# t" s' p8 T" `& ?) j( m
for y=-1000:50:1000
2 |. O% [7 {5 E/ P+ e
g(1:41,：)=G*m*h/((x^2+y^2+h^2)^(3/2));
! e( m# O1 K6 K& U
end
* d! A+ @: `7 a4 k& |
end
0 P5 B/ P4 l& w/ m' `' k5 u# t
作者:
modengxian111
时间:
2023-1-12 13:16
你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
8 Y7 J5 V3 d; e, X
x=-1000:50:1000;
/ m" ] a' d9 f7 |( D: l6 k
y=-1000:50:1000;
U% X" t- @. h" ^" v
for i=1:41
/ ]7 F" q+ l T) l
for j=1:41
5 e. u4 M+ W) P
g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
9 B+ l' Q7 `1 F4 D. G' r2 l
end
6 ^9 u5 r$ p7 V
end
t$ x1 c l5 s! b' o0 Y, v8 U6 K: m
复制代码
0 Y: {( m2 N X1 s; T3 a
( E& j. S# L" s* @
试试这个
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2