クロソイド曲線CASIO fx-4850P版 測量計算ソフト プログラムソース |
|
|
FX-603P版はこちらです→FX-603P クロソイド曲線
クロソイド曲線上の任意の点(P)の座標(Xp,Yp)及び、諸数値(R,τ,σ,So,H,Ho,Hp,Xm,Ym,Tk,TL)を計算します。
クロソイド曲線 【ファイル名:P15】Fixm:Defm 20:H=0Lbl 0:Cls:"<CLOTHOID>" {H} "LEFT_CURVE_⇒1" H"RIGHT_CURVE⇒2_" H=1⇒K=-1:Goto 1△ H=2⇒K=1:Goto 1△ Goto 0 Lbl 1:Cls:{A,B,C,D} A"X0_":B"Y0_" Cls C"XZ_":D"YZ_" T=A:U=B:V=C:W=D Prog "ZA.SUB" Z=999⇒Goto 1△ Z[14]=Y:Z[15]=X Lbl 2:Cls:{E} E"LZ_" E≦0⇒Prog "ER.SUB":Goto 2△ Lbl 3:Cls:{F} F"A_" F≦0⇒Prog "ER.SUB":Goto 3△ Z[10]=F^2÷E Z[11]=1-(E^2÷56÷Z[10]^2) Z[11]=Z[11]+(E^4÷7040÷Z[10]^4) Z[11]=Z[11]-(E^6÷161280÷Z[10]^6) Z[11]=Z[11]×(E^2÷6÷Z[10]) Z[12]=1-(E^2÷40÷Z[10]^2) Z[12]=Z[12]+(E^4÷3456÷Z[10]^4) Z[12]=Z[12]-(E^6÷599040÷Z[10]^6) Z[12]=Z[12]×E Z[13]=Atan (Z[11]÷Z[12]) L=Z[14]-Z[13]×K L<0⇒L=L+360△ L≧360⇒L=L-360△ Lbl 4:Cls:{G} G"L_" G≦0⇒Goto 4△ M=F^2÷G N=1-(G^2÷56÷M^2) N=N+(G^4÷7040÷M^4) N=N-(G^6÷1612800÷M^6) N=N×G^2÷6÷M O=1-(G^2÷40÷M^2) O=O+(G^4÷3456÷M^4) O=O-(G^6÷599040÷M^6) O=O×G Q=Atan (N÷O) R=√(N^2+O^2) P=G÷2÷M×180÷π Z[1]=L+P×K Z[1]<0⇒Z[1]=Z[1]+360△ Z[1]≧360⇒Z[1]=Z[1]-360△ Z[2]=L+Q×K Z[2]<0⇒Z[2]=Z[2]+360△ Z[2]≧360⇒Z[2]=Z[2]-360△ Z[3]=A+R×cos Z[2] Z[4]=B+R×sin Z[2] Z[6]=O−sin P×M Z[5]=N+cos P×M Z[7]=N÷sin P P=90⇒Z[8]=O:Goto 5△ Z[8]=O−N÷tan P Lbl 5:Cls "A=":F:Pause 0 "R=":M▲ "S0=":R:Pause 0 "H=":L→DMS▲ "H0=":Z[2]→DMS "HP=":Z[1]→DMS▲ "XP=":Z[3]:Pause 0 "YP=":Z[4]▲ Goto 8 Lbl 6:Cls "T=":P→DMS "S=":Q→DMS▲ "XM=":Z[5]:Pause 0 "YM=":Z[6]▲ "TK=":Z[7]:Pause 0 "TL=":Z[8]▲ Lbl 8:Cls:I=0 {I} "REPEAT⇒0_MORE⇒5" "CONTINUE_⇒1" I"TOP⇒2_END⇒9_" I=0⇒Goto 5△ I=1⇒Goto 4△ I=2⇒Goto 0△ I=5⇒Goto 6△ I≠9⇒Goto 8△ エラー表示サブルーチン 【ファイル名:ER.SUB】Z=999"DATA_ERROR" Pause 3 Cls 座標変換サブルーチン 【ファイル名:ZA.SUB】Fixm:Z=0(T-V)≠0⇒Goto 1△ (U-W)≠0⇒Goto 1△ Prog "ER.SUB":Goto 2 Lbl 1:X=W-U:Y=V-T Deg:Pol(Y,X):J≦0⇒J=J+360△ X=I:Y=J Lbl 2 ご注意下さい
|
|
■■■■
Copyright(C) 2005 AsakazeNet All rights reserved. 無断転載はご遠慮願います。 ■■■■
当サイトはあさかぜネットが運営しています。 | ||