MATLAB程序代碼

1、MATLAB 程序代碼以及運(yùn)行結(jié)果function = xy_1( A % Detailed explanation goes herex0=653.779y0=604.47 %JD0的坐標(biāo)x1=757.119y1=569.527 %JD1的坐標(biāo)dx=x0-x1dy=y0-y1L=(dx2+dy20.5 %JD1到ID2的距離T=T1(12,28,37 %切線長xk0=T-Lyk0=0 %JD2的局部坐標(biāo)c=0.9473s=-0.3203 %預(yù)設(shè)cos和sin的值%求左端緩和曲線坐標(biāo)for l=0:10:40x=l-(l5/(40*(A2+l9/(3456*(A4%求左端緩和曲線X局部坐標(biāo)y=

2、l3/(6*A-(l7/(336*(A3%求左端緩和曲線Y局部坐標(biāo)dxk=x-xk0dyk=y-yk0B=x0;y0+c,-s;s,c*dxk;dyk%進(jìn)行坐標(biāo)換算endendfunction T1 = T1( a,b,c%求左端切線長% Detailed explanation goes hereA=a+b/60+c/3600r=750p1=p(40,750p2=p(30,750m1=m(40,750T1=(r+p2-(r+p1*cosd(A/sind(A+m1endfunction x = JZ1( %左端坐標(biāo)系坐標(biāo)轉(zhuǎn)換矩陣% Detailed explanation goes herex

3、0=653.779y0=604.47 %JD0的坐標(biāo)x1=757.119y1=569.527 %JD1的坐標(biāo)dx=x0-x1dy=y0-y1L=(dx2+dy20.5 %JD1到ID2的距離T=T1(12,28,37 %切線長xk0=T-Lyk0=0 %JD0的局部坐標(biāo)xk1=Tyk1=0 %JD1的局部坐標(biāo)dxk=xk0-xk1dyk=yk0-yk1A=dxk,-dyk;dyk,dxkb=dx,dyx=inv(A*b%依次輸出cos、sin 的值endxy_1(30000A =30000x0 =653.7790y0 =604.4700x1 =757.1190 y1 = 569.5270 dx

4、 =-103.3400 dy =34.9430 L =109.0879A =12.4769 r =750p =0.0889 p1 =0.08890.0500 p2 =0.0500 m =19.9995 m1 =19.9995 T1 = 101.8149 T =101.8149 xk0 =-7.2730 yk0 =c =0.9473s =-0.3203x =y =dxk =7.2730dyk =B =660.6687 602.1405x =10.0000y =0.0056dxk =17.2730 dyk = 0.0056 B= 670.1435 598.9427 x= 19.9999 y= 0.0444 dxk = 27.2729 dyk = 0.0444 B= 679.6289 595.7766 x= 29.9993 6 y= 0.1500 dxk = 37.2723 dyk = 0.1500 B= 689.1351 592.6738 x=