找回密码
 注册
关于网站域名变更的通知
查看: 458|回复: 1
打印 上一主题 下一主题

为什么循环的结果都一样

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2023-1-12 10:52 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
G=6.67*10^(-11);% {- @" z2 f1 H% r
r=100;/ r0 B! ^; S: }/ Z5 @' s8 X  G* b" X
h=500;
9 {2 b+ ^9 t0 X$ _m=0.1*(4/3)*pi*r^3;
9 _# x" y2 x* Z, Tg=zeros(41,41);) o' w+ a2 `7 l% L% G
for x=-1000:50:1000! i/ y0 o! I( }: |2 Z  U
    for y=-1000:50:1000* c, z+ L4 t4 Q! B! h0 }- h# k
        g(1:41,:)=G*m*h/((x^2+y^2+h^2)^(3/2));
$ t2 H+ H3 P6 }: m! d0 C    end3 P+ R' R9 i' n1 `; h, m- M" G
end
* X5 s% b( ^$ T# R
  • TA的每日心情
    开心
    2023-1-3 15:10
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
    发表于 2023-1-12 13:16 | 只看该作者
    你的赋值语句不对,循环中每次赋值都是给g矩阵的整体赋值,所以结果就是最后一次循环给g整体赋值的结果
    ; S" T- r9 F, `* O) Kx=-1000:50:1000;
    ! n; k- g6 R4 D4 Sy=-1000:50:1000;  s2 P0 D; b; b6 b
    for i=1:41; W, J8 w' e9 Y  e# Z! M
        for j=1:416 @4 {9 |- E0 a. q% t  N/ x, G3 N1 w
            g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));1 x) |* k7 H  R. D* u$ i; L0 L( ^
        end
    3 G5 q( S& M  G' _4 Hend
    : y% O8 S; r0 r. S复制代码
    9 P4 N  @  C+ I* f
    ( D* o7 h8 f* W' I' R9 N试试这个
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-11-23 15:04 , Processed in 0.156250 second(s), 23 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表