Matlab圆弧插补程序+郑国安文档格式.docx
- 文档编号:17202742
- 上传时间:2022-11-28
- 格式:DOCX
- 页数:5
- 大小:14.96KB
Matlab圆弧插补程序+郑国安文档格式.docx
《Matlab圆弧插补程序+郑国安文档格式.docx》由会员分享,可在线阅读,更多相关《Matlab圆弧插补程序+郑国安文档格式.docx(5页珍藏版)》请在冰豆网上搜索。
NorF=input('
请选择圆心(1代表靠近原点2代表远离原点):
\n'
SorN=input('
请选择补走向(1代表顺时针/2代表逆时针):
h=input('
请输入步长\nh='
if((Xe==X0)&
(Ye==Y0))
x01=0;
y01=0;
x02=2*Xe;
y02=2*Ye;
R=sqrt(Xe^2+Ye^2);
mark=0;
else
mark=1;
k1=(Ye-Y0)/(Xe-X0);
k2=-1/k1;
Xz=(X0+Xe)/2;
Yz=(Y0+Ye)/2;
%两点中点坐标
L1=sqrt((X0-Xe)^2+(Y0-Ye)^2)/2;
%两点之间距离的一半
L2=sqrt(R^2-L1^2);
beta=atan(k2);
x01=Xz-L2*cos(beta);
y01=Yz-L2*sin(beta);
%靠近原点的圆心
x02=Xz+L2*cos(beta);
y02=Yz+L2*sin(beta);
%远离原点的圆心
end
if(NorF==1)%判断圆心位置
if((x01^2+y01^2-x02^2-y02^2)<
=0)
x0=x01;
y0=y01;
else
x0=x02;
y0=y02;
end
else
end
%画基准圆
alpha=0:
pi/20:
2*pi;
xx=R*cos(alpha)+x0;
yy=R*sin(alpha)+y0;
plot(xx,yy,'
g:
'
holdon;
axisequal;
Xm=X0;
Ym=Y0;
%NXY=(abs(Xe-X0)+abs(Ye-Y0))/h;
step=0;
Fm=0;
while((Xm-Xe)^2+(Ym-Ye)^2>
h*h/2|(step==0&
mark==0))
if((Xm-x0)>
0&
(Ym-y0)>
=0)XOY=1;
%判断动点所在象限
if((Xm-x0)<
=0&
0)XOY=2;
(Ym-y0)<
=0)XOY=3;
0)XOY=4;
switchXOY
case1
if(SorN==1)
if(Fm>
x1=[Xm,Xm];
y1=[Ym,Ym-h];
x1=[Xm,Xm+h];
y1=[Ym,Ym];
end
if(Fm<
y1=[Ym,Ym+h];
x1=[Xm,Xm-h];
case2
0)
case3
case4
step=step+1;
plot(x1,y1,'
r-'
%由此点和前一点坐标组成的2个向量画直线
Xm=x1
(2);
%保存此点坐标供下次作图和比较时使用
Ym=y1
(2);
Fm=(Xm-x0)^2+(Ym-y0)^2-R^2;
text((x1
(1)+x1
(2))/2,(y1
(1)+y1
(2))/2,[num2str(step)])
pause(0.5);
%延时程序形参为每走一步所用时间
xlabel('
X'
)
ylabel('
Y'
title(['
四象限圆弧插补'
])
holdoff;
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Matlab 圆弧 程序 郑国安
![提示](https://static.bdocx.com/images/bang_tan.gif)