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

为什么循环的结果都一样

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
G=6.67*10^(-11);/ A! m  C( }/ D
r=100;3 b- [% n1 j3 h3 y  }) H8 I
h=500;- ^5 c* T7 F* [( \
m=0.1*(4/3)*pi*r^3;& C# m% b# }) j, a2 W9 ^
g=zeros(41,41);
8 H; ]" T3 s- ]8 P  [4 m) [  }for x=-1000:50:10001 _. ^: {- o- ]8 Z  Z
    for y=-1000:50:1000
: Y  [4 P, p  @& d        g(1:41,:)=G*m*h/((x^2+y^2+h^2)^(3/2));0 h7 T( q9 v6 C$ r
    end
) T  _/ D* o% P5 I! l. Eend
) O  |9 v+ `( z1 N6 k2 Q5 l+ N$ V
  • TA的每日心情
    开心
    2023-1-3 15:10
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
    发表于 2023-1-12 13:16 | 只看该作者
    你的赋值语句不对,循环中每次赋值都是给g矩阵的整体赋值,所以结果就是最后一次循环给g整体赋值的结果
    - \. K3 w# d. i  J3 n- v, T" H; L* Ox=-1000:50:1000;
    4 g6 T# o6 [/ w, ~y=-1000:50:1000;  K  O2 e. Y8 C& ]) w" G( ^5 E
    for i=1:41/ \: S" x2 ?" X& s+ s/ c+ [
        for j=1:410 j9 l# l9 O- Y) W
            g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
    / f9 \! L; ]+ R5 O8 C/ z1 X    end
    9 D; m6 S' D$ hend
    * P5 B" X/ g. T3 |' F复制代码4 Z2 r- p; I6 B. B5 h

    * v4 g) T2 A9 M% H' P3 M试试这个
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-10-4 13:16 , Processed in 0.140625 second(s), 23 queries , Gzip On.

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

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

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