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

为什么循环的结果都一样

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
G=6.67*10^(-11);
. ]1 m7 D4 T8 y$ ^' z7 ar=100;
1 B3 ~% _0 `+ J+ Dh=500;& R7 ?8 O3 e9 w; s8 K2 B: ]" S4 H% g
m=0.1*(4/3)*pi*r^3;
( E7 J$ h; G, n( t3 Ng=zeros(41,41);2 X' V9 Y) [7 C2 D9 ~1 `0 O( q
for x=-1000:50:1000
; O: r# u  k' `, O0 e) Q    for y=-1000:50:1000! n- N0 F  W9 Y' Y
        g(1:41,:)=G*m*h/((x^2+y^2+h^2)^(3/2));9 p1 M2 P* Y, J
    end
- x( w3 A  r( L4 d2 _6 Oend
4 M# S" p' A' ~+ T  X
  • TA的每日心情
    开心
    2023-1-3 15:10
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
    发表于 2023-1-12 13:16 | 只看该作者
    你的赋值语句不对,循环中每次赋值都是给g矩阵的整体赋值,所以结果就是最后一次循环给g整体赋值的结果
    # m* T1 }$ }' u# M% dx=-1000:50:1000;
    * J0 r, r4 G, o8 }1 u8 dy=-1000:50:1000;  \" c6 _  {* i& o8 b4 X
    for i=1:41
    / K# G7 a) N$ P) f/ A    for j=1:41) }2 j5 v4 }- u2 m( k
            g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));/ v& N" Z& ?7 _) y2 G. ^0 z
        end: _! e" k) ?1 R  H) D" \4 y
    end
    ; i! m" X8 O" q# r复制代码
    . p% y* ^) {/ b  B8 N1 \" ~9 M1 Y9 u
    试试这个
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-23 14:21 , Processed in 0.140625 second(s), 23 queries , Gzip On.

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

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

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