24小时热门版块排行榜    

查看: 1138  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

碗碗儿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
回复此楼
努力!
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

碗碗儿2954

新虫 (小有名气)

每次输进去MATLAB就卡住了,难道是死循环了?
努力!
2楼2016-01-27 11:10:51
已阅   关注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的回帖

碗碗儿2954

新虫 (小有名气)

送红花一朵
引用回帖:
5楼: Originally posted by xiegangmai at 2016-01-28 12:06:46
以下代码应该能实现楼主的要求。
算法不一定简练,如
1)判断任意球心距离与相应球半径和大小,可采用数组矢量计算避免循环;
2)比较总体积与设定体积大小,这里需要讨论,有可能这次while循环中随机生成的半径 ...

您的程序写的真好,首先通过判断逻辑变量来决定是否生成小球,然后m=m+1避免了我的两个 for 语句嵌套操作,再有break和continue是我们非专业人士所困惑的,最后像 : 这类的符号我也是很佩服,说再多还是感谢您的宝贵时间!
努力!
6楼2016-01-29 09:46:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

碗碗儿2954

新虫 (小有名气)

引用回帖:
5楼: Originally posted by xiegangmai at 2016-01-28 12:06:46
以下代码应该能实现楼主的要求。
算法不一定简练,如
1)判断任意球心距离与相应球半径和大小,可采用数组矢量计算避免循环;
2)比较总体积与设定体积大小,这里需要讨论,有可能这次while循环中随机生成的半径 ...

我是学材料的,不知道有什么能帮得上您我有正版的MATLAB不知道您有没有需要?
努力!
7楼2016-01-29 09:49:08
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 碗碗儿2954 的主题更新
信息提示
请填处理意见