计算机辅助设计课程设计.docVIP

程序内容说明： 1.本程序是采用matlab7.01编译。 2.程序中的一些符号含义： A：原始数据矩阵； B：曲线要素矩阵； C：里程坐标矩阵； D：直线距离矩阵； E：外矢距； F，坐标转换后的矩阵； Alpha：转角； Azimuth：方位角； R：半径； T：切线长； L：曲线长； Beta：缓和曲线角； Theta:：切线角； Lc：里程； Lo：缓和曲线长； P：内移距； M：切垂距； X1,Y1：局部坐标系下的坐标。 原始程序和必要解释说明： clc; %清除屏幕之前打过的代码 clear;清除工作空间的变量 format; A=[1 2996 265 0 0 ;2 7514 14651 510 8000;3 805 28078 510 8000 ;4 5049 41093 0 0] %原始数据，依次为编号，X坐标，Y坐标，缓和曲线长，曲线半径。 n=size(A,1);%交点个数 i=1; pi=3.1415929;%求各个交点间直线的方位角 while in D(i,1)=sqrt((A(i+1,2)-A(i,2)).^2+(A(i+1,3)-A(i,3)).^2); if A(i+1,2)-A(i,2)0A(i+1,3)-A(i,3)0; azimuth(i,1)=atan((A(i+1,3)-A(i,3))./(A(i+1,2)-A(i,2)));%第一象限 elseif A(i+1,2)-A(i,2)=0A(i+1,3)-A(i,3)0; azimuth(i,1)=pi-atan(-(A(i+1,3)-A(i,3))./(A(i+1,2)-A(i,2)));%第二象限 elseif A(i+1,2)-A(i,2)=0A(i+1,3)-A(i,3)=0; azimuth(i,1)=atan((A(i+1,3)-A(i,3))./(A(i+1,2)-A(i,2)))+pi;%第三象限 elseif A(i+1,2)-A(i,2)=0A(i+1,3)-A(i,3)=0; azimuth(i,1)=2*pi-atan(-(A(i+1,3)-A(i,3))./(A(i+1,2)-A(i,2)));%第四象限 end i=i+1; end k=1; %求直线间的转角 while k=n-2; alpha(k,1)=azimuth(k+1)-azimuth(k); k=k+1; end alpha %计算各曲线要素： for j=1:n; if A(j,5)0; p(j,1)=A(j,4).^2./(24*A(j,5)); %内移距 m(j,1)=A(j,4)./2; %切垂距 e(j,1)=A(j,5).*(1./cos(abs(alpha(j-1,1))/2)-1); %外矢距 t(j,1)=(A(j,5)+p(j,1)).*tan(abs(alpha(j-1,1))./2)+m(j,1);%切线长 l(j,1)=A(j,5).*abs(alpha(j-1,1))+A(j,4); %曲线长 else p(j,1)=0; m(j,1)=0; e(j,1)=0; t(j,1)=0; l(j,1)=0; end end for s=2:n-1; bianhao(s-1,1)=s-1; R(s-1,1)=A(s,5); Lo(s-1,1)=A(s,4); T(s-1,1)=t(s); L(s-1,1)=l(s); E(s-1,1)=e(s); end %输出曲线要素点的坐标，一次为编号，转角，曲线半径，缓和曲线长，切线长，曲线长，外矢距 B=[ bianhao alpha R Lo T L E] %计算各曲线要素点的坐标 d1=D(1)-B(1,5); xZH1=A(1,2)+d1.*cos(azimuth(1)); yZH1=A(1,3)+d1.*sin(azimuth(1)); xZH1 yZH1 d2=D(2)-B(1,5)-B(2,5); xHZ1=A(2,2)+B(1,5).*cos(azimuth(2)); yHZ1=A(2,3)+B(1,5).*sin(azimuth(2)); xHZ1 yHZ1 ly=D(2)-