我一直得到这个错误消息在情节我的代码:
指数在位置2超出数组边界(不得超过1)。
错误CEE210_CP3(第283行)
p =情节(NodePoly (: 1), NodePoly (:, 3));
这是我的原始代码
% ....输入问题数据(点的半径,#,#的成员,θ)
R = 10;%半径
nNodesPoly = 6;
nMembersPoly = nNodesPoly - 1;
θ=π/ nMembersPoly;%一个半圆形的拱
% %……使用一个for循环计算x, y和z每个点的坐标
% ....创建一个矩阵来存储每一个点。
% ....您可以使用变量名“nodePoly”节点的坐标
% ....多边形。
为n = 1: nNodesPoly
x (n:) = [R * cosθ((n - 1) *)];% x R和θ的函数
y (n:) = (R * sinθ((n - 1) *)];% y R和θ的函数
:z (n) = [0];% z R和θ的函数
结束
NodePoly = [x, y, z];
% ....使用一个for循环计算成员或连接和单位
% ....向量的每个成员。
% ....您可以使用变量名“UnVecPoly”单位向量。
PosVecPoly = 0 (nMembersPoly, 3);%初始化内存存储
UnitVecPoly = 0 (nMembersPoly, 3);
OppPosVecPoly = 0 (nMembersPoly, 3);
n = 1: nMembersPoly;
:MemberPoly (n) = (n, n + 1);
SN = MemberPoly (n, 1);
EN = MemberPoly (n, 2);
PosVecPoly (n) = NodePoly (EN,:) -NodePoly (SN:);
长度(n:) =规范(PosVecPoly (n,:));
UnVecPoly (n:) = PosVecPoly长度(n:) / (n:);
结束
% %共创输出命令窗口
流(' % s \ n ', ' - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -”)
流(' % s \ n”、“多边形- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -”)
流(' % s \ n ', ' - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -”)
(i = 1: nNodesPoly)
流(' % s % 8我% 8.3 f % 8.3 f % 8.3 f \ n”、“节点:”,NodePoly(我:))
结束
流(“\ n”)
(i = 1: nMembersPoly)
流(' % s % 8我% 8.3 f % 8.3 f % 8.3 f \ n”,“单位向量:”,我,UnVecPoly(我:))
结束
流(“\ n \ n \ n”)
% %……创建一个阴谋
如果=图(3);clf;网格;轴相等;抓住;
包含(“X”);ylabel (Y);标题(“多边形桁架”);
p =情节(NodePoly (:, 2), NodePoly (:, 3));
集(p,“颜色”,“黑色”,“线宽”,2);
散射(NodePoly (:, 2), NodePoly(:, 3),“填满”,“绿色”);
%的程序
