matlab根据给定3点画圆弧
matlab根据给定3点坐标画圆弧
实现原理看代码及注释
主程序arcPlot.m
function arcPlot(pos,colour)
p1=pos(1,:);
p2=pos(2,:);
p3=pos(3,:);
p1p2=p2-p1;
u1=p1p2/sum(p1p2.^2);
p2p3=p3-p2;
u2=unitVec(p2p3);
normal=cross(p1p2,p2p3);
normal=unitVec(normal);
per2=cross(normal,u2);
per2=unitVec(per2);
mid1=(p1+p2)/2;
mid2=(p2+p3)/2;
dis2=sqrt(sum(p2p3.^2));
L=(vecDot(mid1,u1)-vecDot(mid2,u1))/vecDot(per2,u1);
arc_R=sqrt(L^2+1.0/4*dis2^2);
C=mid2+L*per2;
vec1=p1-C;
vec1=unitVec(vec1);
vec2=p3-C;
vec2=unitVec(vec2);
theta=vecDot(vec1,vec2);
normal2=cross(vec1,vec2);
theta=acos(theta);
if(sqrt(sum((normal-normal2).^2))&1.0e-5)
theta=2*pi-
vx=[1,0,0];
vy=[0,1,0];
vz=[0,0,1];
v2=cross(v3,v1);
v2=unitVec(v2);
r11=vecDot(v1,vx);r12=vecDot(v2,vx);r13=vecDot(v3,vx);
r21=vecDot(v1,vy);r22=vecDot(v2,vy);r23=vecDot(v3,vy);
r31=vecDot(v1,vz);r32=vecDot(v2,vz);r33=vecDot(v3,vz);
Tr=[r11,r12,r13,0;
r21,r22,r23,0;
r31,r32,r33,0;
Tt=[1,0,0,C(1);
0,1,0,C(2);
0,0,1,C(3);
t=0:(theta)/100:
x1=arc_R*cos(t);
y1=arc_R*sin(t);
pt=[x1;y1;z1;ones(size(x1))];
plot3(pt(1,:),pt(2,:),pt(3,:),colour);
xlabel('x')
ylabel('y')
zlabel('z')
view(-37.5,30)
axis equal
自定义函数vecDot()
function z=vecDot(x,y)
z=sum(x.*y);
自定义函数unitVec()
function vector=unitVec(vec)
vector=vec/sqrt(sum(vec.^2));
pos=[0 0 0;2 2 0;4 0 0];
arcPlot(pos,'r');
结果如下：
祝您生活愉快！
本文已收录于以下专栏：
相关文章推荐
已经知道圆心和两点，画出两点间的圆弧？思路，先分别求出原点到两点向量和x轴间的夹角余弦值，然后求反余弦得出角度，这里要注意的是，求出的向量的夹角是按几何坐标系中定义的逆时针，然后用DrawArc()函...
原文链接：http://blog.csdn.net/lgm252008/article/details/、
文中的图很重要！！！
UIBezierPath通过
- (void...
之前在微博上看到一个用弧形展示的百分比图，感觉很有创意，就想使用R来实现一样的功能，效果如下：
//wbphotos/lar...
package com.example.
import java.util.ArrayL
import android.os.B
import android.app...
emgucv 画圆弧
他的最新文章
讲师：汪剑
讲师：刘道宽
您举报文章：
举报原因：
原文地址：
原因补充：
(最多只允许输入30个字)}