应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 简单的MATLAB几何建模问题
51/1返回列表
查看: 1244  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

碗碗儿2954

新虫 (小有名气)

[交流] 简单的MATLAB几何建模问题 已有2人参与

为了建立一系列随机分布的小球，编写程序如下，请懂的人帮忙看看，真是没辙了！   
    k=1;
    i=2:50;
    m=1:49;
    R(50)=ones;
    V(50)=ones;
    VSUM(50)=ones;
    [X,Y,Z]=sphere(50);
    R(1)=rand(1)/10;
    V(1)=(4/3)*pi*R(1)^3;
    originX=rand(1,50);
    originY=rand(1,50);
    originZ=rand(1,50);
    surf(X*R(1)+originX(1),Y*R(1)+originY(1),Z*R(1)+originZ(1));
    VSUM(1)=V(1);
for i=2:50
    R(i)=rand(1)/10;
    V(i)=(4/3)*pi*R(i)^3;
end
while (VSUM<v)
     for m=1:i-1
         A=R(i)+R(m);
         B=sqrt((originX(i)-originX(m))^2+(originY(i)-originY(m))^2+(originZ(i)-originZ(m))^2);
         if (A<B)
            k=k+1;
         end
         continue
     end
end
if (k==i)
      surf(X*R(i)+originX(i),Y*R(i)+originY(i),Z*R(i)+originZ(i));
      fprintf('The volume of the ith sphere is %12.5f\n',V(i));
      VSUM(i)=VSUM(i-1)+V(i);
      hold on
end
axis equal
回复此楼

» 猜你喜欢
努力！
1楼 2016-01-27 11:09:48
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

灿烂豆沙包

木虫 (小有名气)
厉害

[ 发自手机版 http://muchong.com/3g ]
回复此楼
keepcalmandcarryon..
8楼2016-03-05 00:16:19
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

碗碗儿2954

新虫 (小有名气)
每次输进去MATLAB就卡住了，难道是死循环了？
赞一下 回复此楼
努力！
2楼2016-01-27 11:10:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
碗碗儿2954: 金币+10 2016-01-28 09:35:50
简单写了下代码，楼主参考参考。
CODE:
function spheres

n = 100;
r = rand( 1, n ) / 10;

for i = 1 : n
    color_value = rand( 1, 3 ); % 颜色
    position = rand( 1, 3 );  % 球心坐标
     [ x, y, z ] = unitsphere;
     
     x =  r( i ) * x + position( 1 );
     y =  r( i ) * y + position( 2 );
     z =  r( i ) * z + position( 3 );
     
     h = surf( x, y, z);
     set( h, 'EdgeColor', 'none', ...
         'FaceColor', color_value  );
     hold on
end
axis equal

V = 4 * pi * r .^3 / 3;
VSUM = sum( V )

function [ x, y, z ] = unitsphere

r = 1; % 单位圆
n = 20;

theta = ( -n : 2 : n ) / n * pi;
phi = ( -n : 2 : n )' / n * pi / 2;

cosphi = cos( phi );
cosphi( 1 ) = 0;
cosphi( n + 1 ) = 0;
sintheta = sin( theta );
sintheta( 1 ) = 0;
sintheta( n + 1 ) = 0;

x = r * cosphi * cos( theta );
y = r * cosphi * sintheta;
z = r * sin( phi ) * ones( 1, n + 1 );
赞一下(1人) 回复此楼
明德厚学、求是创新
3楼2016-01-27 16:41:01
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

碗碗儿2954

新虫 (小有名气)
引用回帖:
3楼: Originally posted by xiegangmai at 2016-01-27 16:41:01
简单写了下代码，楼主参考参考。

function spheres

n = 100;
r = rand( 1, n ) / 10;

for i = 1 : n
    color_value = rand( 1, 3 ); % 颜色
    position = rand( 1, 3 );  % 球心坐标
      = unit ...

好心人您好，我看了您的程序是可行的，但是我想要通过比较第 i 个球和前面生成的 (i-1)的球半径和与球心距离的大小来避免球体之间的重合，在总体积小于比如 v 的情况下，由于position只是任意的1*3数组，并不能存储第 i 个球体的位置信息？我想出的以下语句好像也有问题
while (VSUM<v)
     for m=1:i-1
         A=R(i)+R(m);
         B=sqrt((positionX(i)-positionX(m))^2+(positionY(i)-positionY(m))^2+(positionZ(i)-positionZ(m))^2);
         if (A<B)
            k=k+1;
         end
         continue
     end
end
赞一下 回复此楼
努力！
4楼2016-01-28 09:14:49
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 山东高校教师考核超级无底线，员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 298求调剂 +9 钉叮咚冬瓜 2026-04-07 11/550 2026-04-11 09:35 by zhq0425
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 一志愿东北大学控制工程085406数二英二385，求调剂 +8 Ezra_Zhang 2026-04-09 8/400 2026-04-11 09:15 by 猪会飞
[考研] 342电子信息专硕求调剂 +9 你让我怎么荔枝 2026-04-10 10/500 2026-04-11 08:33 by zhq0425
[材料工程] 材料调剂推荐 +8 蛋糕x2 2026-04-07 8/400 2026-04-10 23:13 by Ftglcn90
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 材料专业344求调剂 +16 hualkop 2026-04-10 21/1050 2026-04-10 17:28 by laoshidan
[考研] 0703化学求调剂 +21 不知名的小卅 2026-04-08 21/1050 2026-04-09 18:55 by l_paradox
[考研] 材料调剂 +14 一样YWY 2026-04-06 14/700 2026-04-08 23:00 by 猪会飞
[考研] 软件工程求调剂22软工296分求调剂，接受跨调 +4 yangchen2017 2026-04-08 5/250 2026-04-08 21:56 by 土木硕士招生
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 专硕085403，291分，有两篇专利，一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
[考研] 一志愿北交大材料工程总分358求调剂 +10 cs0106 2026-04-05 12/600 2026-04-06 19:41 by 无际的草原
[考研] 材料调剂 +5 小刘同学吖吖 2026-04-06 5/250 2026-04-06 18:34 by sherry_1901
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
信息提示
关闭
请填处理意见
关闭