5800計(jì)算器程序辛普森公式程序(通用)解析

1、V=19：FixmJ（與預(yù)留的自由變量數(shù)有關(guān)?。㎡：O“1.ZX 2.BX” JQ=0JLbI 1JOH仁QG :Q”L”:G”XJJD”（顯示“ L”輸入邊樁寬度， 均為正值，顯示“ XJJD”輸入邊樁測點(diǎn)與線路中線上點(diǎn)的連線 與線路中線上點(diǎn)的切線的夾角，如右法線為 90 度，左法線為 270 度或者-90 度，其余同理。）Prog”ZB ”:X二X+Rec（Q,T+G）: Y=Y+J:”X= ” X : Pause 0 Y= ” ： Y （X ”X= ” Pause 0 : YY二”：）”T= ” : T DMS （T”T二” ）Goto 1（注釋：括號內(nèi)寫法為 Casio4800P 計(jì)算

2、器顯示方式）2.ZB ”:坐標(biāo)計(jì)算子程序V=19|FixmJ（與預(yù)留的自由變量數(shù)有關(guān)?。￤=21J（與預(yù)留的自由變量數(shù)有關(guān)！）K:;LbI 2 J1.”XLZB ”坐標(biāo)計(jì)算主程序F=0:W=27=F=39 903 |: （注釋：方框內(nèi)數(shù)字 27第七段為斷鏈的段號；方框內(nèi)數(shù)字 39.903 為斷鏈值，短鏈輸“ + ”，長鏈輸“-”；如遇斷鏈，依此增加即可！)E= ZW+3a+ZW+F :K ZW+3a =K V=W 1 :工 W=W+1: Goto 2 ；K 24200=K=K F: (方框?yàn)閿噫滭c(diǎn)的后樁號，即去向樁 號，如遇斷鏈，依此增加即可！)S=K ZV+3a+1:P=(ZV+2a+1 Z

3、V+a+1) - ZV+1:匸P*S:C=ZV+4a+1:D=ZV+a+1:J二C+(l+2D)*S*90/n:L =1:X=ZV+5a+1+S*(CosC+CosJ+4刀(Cos(C+(L+0.5)l/B+2D)(L+0.5)*90*S/B/n(C+(L*l/B+2D)*L*90*S/B/n),L,1,(B-1)/6/B:),L,0,(B-1)+2E(CosY二ZV+6a+1+S*(SinC+SinJ+4刀(Sin(C+(L+0.5)I/B+2D)(L+0.5)*90*S/B/n(C+(L*I/B+2D)*L*90*S/B/n),L,1,(B-1)/6/B:),L,0,(B-1)+2E(Si

4、nT=J:3.SJK.DAT”:數(shù)據(jù)庫主程序Mcl:Defm(20+7a):Z21=?:Z21+3a=?:Z21+4a=?:Z21+5a=?:Z21+6a=?:Prog”SJK”:(輸入程序時(shí)，Defm 后須換成相應(yīng)數(shù)據(jù)，a 為曲線元段數(shù))4.SJK”： 數(shù)據(jù)庫計(jì)算子程序a二？:V=20:T=0:B二？（a:曲線元段數(shù)；B:等分?jǐn)?shù)，最少為 2, 般為 4 就可以滿足要求）Lbl 0JF=0:T=6=F=39 9031: （注釋：方框內(nèi)數(shù)字 6第七段減一為斷鏈的段號減去 1 所得；方框內(nèi)數(shù)字 39.903 為斷鏈值，短鏈輸“ + ”,長鏈輸“-”；如遇斷鏈，依此增加即可?。〦= ZV+3a+1+

5、ZV+1+F:ZV+3a+2=E :S=ZV+1:匸ZV+2a+1 ZV+a+1:C=ZV+4a+1: D=ZV+a+1:J= C+(I+2D)*S*90/n:ZV+4a+2=J:L=1 :ZV+5a+2=ZV+5a+1+S*(CosC+CosJ+4刀Cos(C+(L+0.5)I/B+2D)(L+0.5)*90*S/B/n), L,0,(B-1)Cos(C+(L*l/B+2D)*L*90*S/B/n),L,1,(B-1)/6/B:ZV+6a+2=ZV+6a+1+S*(SinC+SinJ+4刀Sin(C+(L+0.5)I/B+2D)(L+0.5)*90*S/B/n), L,0,(B-1)Sin(

6、C+(L*I/B+2D)*L*90*S/B/n),L,1,(B-1)/6/B:V=V+1:T=T+1:T a-2 =”END OF THE PROGRAM ” Goto 0為了使程序能夠計(jì)算連續(xù)多個(gè)曲線元的坐標(biāo)，減少計(jì)算過程中 的數(shù)據(jù)輸入，可利用 CASIOfx-4800/4850P 計(jì)算器提供的擴(kuò)展變量設(shè) 置曲線元數(shù)據(jù)庫， 假定某一段線路包含有 a 個(gè)曲線元， 則可設(shè)置曲線 元數(shù)據(jù)庫如下表所示：曲線元數(shù)據(jù)庫設(shè)置表名稱第一段第二段第 i 段第 a-1 段第 a 段曲線兀長度212220+i19+a20+a起點(diǎn)曲率21+a22+a20+a+i19+2a20+2a終點(diǎn)曲率21+2a22+2a20+

7、2a+i19+3a20+3a起點(diǎn)樁號21+3a起點(diǎn)方位角21+4a起點(diǎn)x坐標(biāo)21+5a起點(diǎn)y坐標(biāo)21+6a擴(kuò)展變量 Z1Z20用作自由變量，可以用來彌補(bǔ)字母變量的不足，此項(xiàng)設(shè)置可自由選擇！因此，數(shù)據(jù)庫的起點(diǎn)本程序從Z21開始，表中 Z21Z21+3a以及 Z21+4a、Z21+5a、Z21+6a需預(yù)先置入，表中空白部分由程序自動(dòng)計(jì)算生成。CASIOfx-4800/4850P計(jì)算器程序源程序該算法共有四個(gè)程序： 1.線路坐標(biāo)計(jì)算主程序（ XLZB ）；2.數(shù)據(jù)庫生成主程序（SJK.DAT）; 3坐標(biāo)計(jì)算子程序（ZB）; 4數(shù)據(jù)庫計(jì)算子程序（SJK）。程序使用說明不管是整體式還是分離式，主線還是匝

8、道，都能用同一個(gè) 計(jì)算器就可進(jìn)行各路段相關(guān)坐標(biāo)計(jì)算，只是分別采用了各自路 段的數(shù)據(jù)庫，如果再結(jié)合施工放樣正反算程序及高程程序、隧 道及路基開挖斷面程序，則可形成一個(gè)完善實(shí)用的施工測量程 序包，特別適合于公路、鐵路現(xiàn)場施工測量人員使用。在計(jì)算某一段線路時(shí)，首先必須運(yùn)行相應(yīng)路段數(shù)據(jù)庫生成主程序（SJK.DAT），然后運(yùn)行線路坐標(biāo)計(jì)算主程序（XLZB ）,按照程序提示輸入是計(jì)算中線還是邊線，若是計(jì)算邊線除需輸入 樁號外，還需輸入寬度及斜交角度。因此，根據(jù)線路實(shí)際情況 可編制相應(yīng)路段數(shù)據(jù)庫生成主程序即可完成相應(yīng)路段坐標(biāo)計(jì)算 工作。算例本文提供聶讓編著全站儀與高等級公路測量 306308 頁算例的 計(jì)算

9、結(jié)果如下：樁號原著坐標(biāo)本算法坐標(biāo)備注XYXY直O(jiān)O)K0+1161378.2142822.9501378.2142822.950起點(diǎn)樁號：kO+116,起點(diǎn)方位角200，線+1401355.6612814.7411355.6612814.742xq=1378.214,yq=2822.950A + 1501346.2642811.3211346.2642811.321回A +1501346.2642811.3211346.2642811.321pa=0pb=1/124旋+2001300.1422792.1211300.1422792.121線B +2241279.8452779.3631279.

10、8452779.364圓B +2241279.8452779.3631279.8452779.364pb=1/124pc=1/124曲+2801244.1792736.8071244.1792736.808線C +341.8401230.6822677.1131230.6822677.114回C +341.8401230.6822677.1131230.6822677.114pc=1/124pd=1/60旋+3801238.6942640.0771238.6762640.065線D+407.6501254.8242617.8411254.7852617.831圓D+407.6501254.82

11、42617.8411254.7852617.831pd=1/60pe=1/60曲+4601303.3722603.2091303.3172603.148線E+495.8261335.2752618.3081335.2362618.214回E+495.8261335.2752618.3081335.2362618.214pe=1/60pf=0旋+5401357.2722656.0001357.2732655.882線F+577.4931364.6202692.7311364.6592692.605直F+577.4931364.6202692.7311364.6592692.605pf=0pg=0

12、線+6001367.9752714.9871368.0372714.857G +6401373.9382754.5401374.0412754.404名稱第一段第二段第三段第四段第五段第六段第七段曲線兀長度3474117.8465.8188.17681.66762.507起點(diǎn)曲率001/1241/1241/601/600終點(diǎn)曲率01/1241/1241/601/6000起點(diǎn)樁號116起點(diǎn)方位角200起點(diǎn)x坐標(biāo)1378.214起點(diǎn)y坐標(biāo)2822.950段數(shù)12345678910長度202.5262.324202.5202.5216.244202.5744.7153.6526.26153.6位置Z

13、21222324252627282930起點(diǎn)曲率0-1/1000-1/100001/10001/100000-1/1500-1/1500位置Z31323334353637383940終點(diǎn)曲率-1/1000-1/100001/10001/100000-1/1500-1/15000位置Z41424344454647484950起點(diǎn)樁號22258.882位置Z51525354555657585960起點(diǎn)方位角239-11-15.18位置Z61626364656667686970起點(diǎn)X坐標(biāo)318841.194位置Z71727374757677787980起點(diǎn)Y坐標(biāo)493525.938位置Z818283

14、84858687888990說明：1、該線路在第七段線元有斷鏈，具體是 K24160.097二K24200，短鏈 39.903m;2、左轉(zhuǎn)曲率為“一” ；“ + ”；3、該數(shù)據(jù)庫對應(yīng)的程序見前。已知坐標(biāo)求該點(diǎn)里程并帶入線路高程及主洞開挖放樣的復(fù)合計(jì)算程序Lbi 7IJZ1=Z5+IcosJ （Z1= I 儀器直接顯示坐標(biāo)時(shí)用）Z2=Z6+ISinJ （Z2= J 儀器直接顯示坐標(biāo)時(shí)用）prog“ZB”Lbi 8Z10=（X Z1）Sin（U 90） （Y Z2）Cos（U 90）AbsZ10 K二K+Z10:Pol（Z1 X,Z2 Y）:Z4= : Z4= ICos（J U+90）“ K= ”

15、 : K Z8=1 = Prog “YXGC ”：半Prog “ZXGC”： Prog “ZDKW” :Goto 7:半K=K+Z10: prog “ZB” : Goto 8: 注釋：Z5 、Z6 ：測站坐標(biāo)（直接輸坐標(biāo)時(shí)不用）I、J:實(shí)測距離及方位角（或?qū)崪y坐標(biāo)）JSZCXZ1 、 Z2 ：實(shí)測點(diǎn)坐標(biāo)Z10 ：中間變量Z4 ：實(shí)測點(diǎn)到線路中線的距離，負(fù)數(shù)為線左，正數(shù)為線右“ GC” ：設(shè)計(jì)高程程序(子程序)“ ZDKW ”：開挖放樣子程序X、Y:實(shí)測點(diǎn)對應(yīng)里程的設(shè)計(jì)坐標(biāo)，由線路坐標(biāo)計(jì)算子程序“ ZB”帶入U(xiǎn) :實(shí)測點(diǎn)對應(yīng)里程的切線方位角三心圓開挖放樣程序“ZDKW”(也可用于同心圓) ：請自

16、己畫示意圖理解！CR :C “ H”：Z9+C A=V(Abs(Z4 士 Z8x?)+d)2+(Z9+C H h2)2) RB二Abs(Z4 士 Z8x?)+d V( (R+ r1) 2 (Z9+C H h2)2) V=Z9+C-H-h2-V(R+r1)2(Z4 士 Z8x?)2)工 Z9+C A=V(Abs(Z4 士 Z8x?) 2+(Z9+C H h1)2) RB=Abs(Z4 士 Z8x?) V(R2(Z9+CHh1)2)V=Z9+C-H-h1-V(R2(Z4 士 Z8x?)2)工 V=Z9+C H h1-V(R2 Abs(Z4 士 Z8X?) 2) 注釋：C（顯示 H）:實(shí)測點(diǎn)至儀器中心的高差（如直接輸入高程，則將C “H”改為 C “SCGC”，且程序中的 Z9+C 就直接用 C 代替）R :開挖半徑（小）Z9 :儀高（直接輸實(shí)測高程時(shí)不用）A :實(shí)測點(diǎn)半徑與設(shè)計(jì)半徑之差值（徑向） ， +超挖，欠挖B :實(shí)測點(diǎn)與設(shè)計(jì)點(diǎn)之差值（水平方向），+超挖，一欠挖V :實(shí)測點(diǎn)與設(shè)計(jì)點(diǎn)之豎向距離r1 :大半徑與小半徑