作者maygood (我本闲人)
看板MATLAB
标题[其他] 请问一下plot的问题
时间Thu Jun 20 12:05:12 2013
我在写动作分析的角度问题时
我写了以下的程式
但是plot出来的线不连续
请问是哪里出问题
跪求高人指点
感激!!!!!!!
len=length(D_Lanklem(:,1));
for i=1:len
Shank_COM(i,:)=Ram(i,:)+0.567.*(Rkm(i,:)-Ram(i,:));
Foot_COM(i,:)=(D_Rtoe(i,:)+0.5.*(D_Rheel(i,:)-D_Rtoe(i,:)));
%% shank coordinate system
Z_shank(i,:)=(Rkm(i,:)-Ram(i,:))./(norm(Rkm(i,:)-Ram(i,:)));
X_shank(i,:)=cross(-VS(i,:),c(i,:))./norm(cross(-VS(i,:),c(i,:)));
Y_shank(i,:)=cross(Z_shank(i,:), X_shank(i,:));
R_shank(i*3-2:i*3,:)=[X_shank(i,:)' Y_shank(i,:)' Z_shank(i,:)'];
%% foot coordinate system
Z_foot(i,:)=(D_Rtoe(i,:)-Ram(i,:))./norm(D_Rtoe(i,:)-Ram(i,:)); % toe和ankle中点连线
X_foot(i,:)=cross(-Z_foot(i,:),e(i,:))./norm(cross(-Z_foot(i,:),e(i,:)));
Y_foot(i,:)=cross(Z_foot(i,:), X_foot(i,:));
R_foot(i*3-2:i*3,:)=[X_foot(i,:)' Y_foot(i,:)' Z_foot(i,:)'];
%% 计算角度
R_ankle(i*3-2:i*3,:)=R_shank(i*3-2:i*3,:)'*R_foot(i*3-2:i*3,:);
Beda(i,:)=asind(-R_ankle(i*3-1,3));
A=Beda(i,:);
alpha(i,:)=180/pi*atan2((R_ankle(i*3-2,3)/cos(A)),(R_ankle(i*3-1,3)/cos(A)));
gama(i,:)=180/pi*atan2((R_ankle(i*3-1,1)/cos(A)),(R_ankle(i*3-1,2)/cos(A)));
angle_ankle(i,:)=[alpha(i,:) Beda(i,:) gama(i,:)];
a=angle_ankle(i,:);
%% plot
subplot(2,2,1),plot(alpha,'d');
subplot(2,2,2),plot(Beda,'o');
subplot(2,2,3),plot(gama,'x');
end
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 163.15.154.17
1F:→ jatj:plot(alpha, 'd-'); 06/20 12:09
2F:嘘 jatj:也不回报收到了或公不工作... 06/22 23:54
3F:→ jatj:射後不理嘛 06/22 23:54
4F:→ maygood:楼上抱歉~~这几天在忙~没有上来看结果~感谢您的回覆 06/23 20:50