QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 1209|回复: 4
打印 上一主题 下一主题

程序出现“访问什么的时候指数必须是整数或逻辑 ”的错误

[复制链接]
字体大小: 正常 放大

18

主题

12

听众

44

积分

升级  41.05%

  • TA的每日心情
    慵懒
    2016-7-13 15:45
  • 签到天数: 6 天

    [LV.2]偶尔看看I

    自我介绍
    好还哦好海鸥
    跳转到指定楼层
    1#
    发表于 2016-7-5 10:29 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    clear all;clc
    load samp22.txt
    idx=find(samp22(:,4)==1);
    pointS=samp22(idx,1:3);
    data.x=pointS(:,1);
    data.y=pointS(:,2);
    data.z=pointS(:,3);
    nPoint=length(data.x);
    pointS=[data.x';data.y';data.z'];ch=convhull(data.x,data.y);
    area=polyarea(data.x(ch,1),data.y(ch,1));
    densityHorizontal=nPoint/area;
    xDis=1/sqrt(densityHorizontal);
    xDis=ceil(xDis*10)/10;
    yDis=xDis;   

    ch1=convhull(data.x,data.z);
    area1=polyarea(data.x(ch,1),data.z(ch,1));
    densityHorizontal1=nPoint/area1;
    xDis1=1/sqrt(densityHorizontal1);
    xDis1=ceil(xDis1*10)/10;

    ch2=convhull(data.y,data.z);
    area2=polyarea(data.y(ch,1),data.z(ch,1));
    densityHorizontal2=nPoint/area2;
    yDis1=1/sqrt(densityHorizontal2);
    yDis1=ceil(yDis1*10)/10;
    zDis=min(xDis1,yDis1);%%%取y,z和x,z之中比较小的
    min_x=min(pointS(1,);    max_x=max(pointS(1,);
    min_y=min(pointS(2,);    max_y=max(pointS(2,);
    min_z=min(pointS(3,);    max_z=max(pointS(3,);
    Column=ceil((max_x-floor(min_x))/xDis);%列(j)
    Row = ceil((max_y-floor(min_y))/yDis);%行 (i)
    Height= ceil((max_z-floor(min_z))/zDis);%高(k)
    elementMesh=zeros(Row,Column,Height);
    X=pointS(1,;  Y=pointS(2,;  Z=pointS(3,;
    >> for i=1:nPoint
        xColumn=ceil((X(i)-floor(min_x))/xDis);
        yRow=ceil((Y(i)-floor(min_y))/yDis);
        zHeight=ceil((Z(i)-floor(min_z))/zDis);
        elementMesh(yRow,xColumn,zHeight)=1;
    end
    clear i

    出现这样的错误:
    ??? Attempted to access elementMesh(2,0,21); index must be a positive integer or logical.


    这是由于什么引起的,对别的数据没有事,但是对这组数据就出现这样的问题,麻烦大神指导一下!!



    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

    2088

    主题

    31

    听众

    3238

    积分

  • TA的每日心情
    奋斗
    2016-11-10 17:20
  • 签到天数: 102 天

    [LV.6]常住居民II

    自我介绍
    很烦呀

    社区QQ达人

    群组2013年电工杯B题讨论群

    群组2013电工杯A题讨论群组

    群组数学中国第二期SAS培训

    群组2016美赛交流群组

    回复

    使用道具 举报

    2088

    主题

    31

    听众

    3238

    积分

  • TA的每日心情
    奋斗
    2016-11-10 17:20
  • 签到天数: 102 天

    [LV.6]常住居民II

    自我介绍
    很烦呀

    社区QQ达人

    群组2013年电工杯B题讨论群

    群组2013电工杯A题讨论群组

    群组数学中国第二期SAS培训

    群组2016美赛交流群组

    回复

    使用道具 举报

    16

    主题

    12

    听众

    494

    积分

    升级  64.67%

  • TA的每日心情
    奋斗
    2016-11-4 09:23
  • 签到天数: 36 天

    [LV.5]常住居民I

    自我介绍
    是我 就是我 你们的英雄 小哪吒

    群组2016国赛备战群组

    回复

    使用道具 举报

    3503

    主题

    538

    听众

    5990

    积分

  • TA的每日心情
    开心
    2017-2-7 15:12
  • 签到天数: 691 天

    [LV.9]以坛为家II

    社区QQ达人 元老勋章 发帖功臣 新人进步奖 优秀斑竹奖 金点子奖 原创写作奖 最具活力勋章 助人为乐奖 风雨历程奖

    群组2013年国赛赛前培训

    群组2014年地区赛数学建模

    群组数学中国第二期SAS培训

    群组物联网工程师考试

    群组2013年美赛优秀论文解

    由于xColumn = 0引起的
    你可以在你最后的for循环中,按照以下方式修改(添加红字部分):
    for i=1:nPoint
        xColumn=ceil((X(i)-floor(min_x))/xDis);
        yRow=ceil((Y(i)-floor(min_y))/yDis);
        zHeight=ceil((Z(i)-floor(min_z))/zDis);
        if (xColumn>0)
            elementMesh(yRow,xColumn,zHeight)=1;
        end
    end
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2025-9-12 07:28 , Processed in 0.531949 second(s), 73 queries .

    回顶部