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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 简单的MATLAB几何建模问题
51/1返回列表
查看: 1136  |  回复: 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的回帖

碗碗儿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的回帖
查看全部 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 个回答
信息提示
关闭
请填处理意见
关闭