怎么编制卡西欧4850计算器计算曲线和高程的程序??

终于编出来了纯锋，累死我了，代码如下：

程序名:ZOU BIAO

M"JD":G"J(X)":D"J(Y)":F:R:A"J(＜)":I"LS":E"E(1,-1)"

T=(R+I2÷24÷R)tan(A÷2)+(I÷2-I3÷240÷R2):L=R(A-I÷R×180÷∏)×∏÷180

Z"ZH"=M-T:Z[3]"ZH(X)"=G-TcosF:Z[4]"ZH(Y)"=D-TsinF:Q"HY"=Z+I:S"做慎晌YH"=Q+L:O"HZ"=S+I:H=F+AE:Z[5]"HZ(X)"=G+TcosH:Z[6]"HZ(Y)"=D+TsinH:

LbI O

{K}

K≤Z L=Z-K:U=F:X"X"孝羡=Z[3]-LcosF⊿ Y"Y"=Z[4]-LsinF⊿ Goto 6

K≤Q Goto 2

K≤S Goto 3

K≤O Goto 5

K＞O L=K-O:N=H:U=N:Goto 1

LbI 1

X"X"=Z[5]+LcosN⊿ Y"Y"=Z[6]+LsinN⊿ Goto 6

LbI 2

L=K-Z:B=L-L5÷(40R2I2):C=L3÷(6RI):U=F+L2÷(2RI)×180÷∏×E: Goto 4

LbI 3

L=K-Q:V=(L÷R+I÷2÷R)×180÷∏:B=RsinV+I÷2-I3÷(240R2):C=R(1-cosV)+I2÷(24R):U=F+(I÷(2R)+L÷R)×180÷∏×E: Goto 4

LbI 4

X"X"=Z[3]+BcosF-ECsinF⊿ Y"Y"=Z[4]+BsinF+ECcosF⊿ Goto 6

LbI 5

L=O-K: B=L-L5÷(40R2I2):C=L3÷(6RI):U=H-L2÷(2RI)×180÷∏×E: X"X"=Z[5]-BcosH-ECsinH⊿ Y"Y"=Z[6]-BsinH+ECcosH⊿ Goto 6

LbI 6

{J}:J"J(0,1)"

J=0 Goto 0:≠ J=1 Goto 7:

LbI 7

{W,V}

Z[7]"Xz"=X+Wcos(U+V-180)⊿ Z[8]"Yz"=Y+Wsin(U+V-180)⊿

Z[9]"Xy"=X+Wcos(U+V)⊿ Z[10]"Yy"=Y+Wsin(U+V)⊿

Goto 0

计算带有缓和曲线的路线的中、边桩坐标

M…转点桩号 G、D…转点坐标

F…方位角 R…半径

A…转角 I…缓和曲线长

E…计算条件，左转为-1右转为1

T…切线长

K…计算桩号

前直线段

第一缓和曲线段

圆曲线段

第二缓和曲线段

后直线段

J…计算条件0转入中桩1转入边桩

W…边桩与中桩之间的距离

如何用matlab计算两条曲线的相关性？这个问题我们可碧腊以按下列方法来处理：

第一步，将第一条曲线的逗行数据赋值给A，将第二条曲线的数据赋值给B

第二步，计算A 的均值 μA和标准差σA，计算B 的均值 μB和标准差σB

第三步，利用Pearson 相关系数公式，计算A和B的相关性，计算方法如下

第四步，判断其相关性是正相关，还是负相关

计算实例：悔指滑已知：

t=[15 20 25 30 35 40 45 50 55 60 65 70]

A=[0.3 0.4 0.7 1.1 1.7 2.3 3 3.6 4.1 4.4 4.7 4.8]

B=[0.2854 0.4540 0.7080 1.0725 1.5604 2.1561 2.8068 3.4369 3.9795 4.4015 4.7045 4.9098]

求其相关系数

计算结果为

ρ(A,B)=0.99844，正相关

(CASIO-4850)竖曲线标高计算程序--全线贯通

(CASIO-4850)竖曲线标高计算程序

1,加了水准点数据.

2,加了竖曲线数据库(实现了全线贯通)

3,加了测量仪器的选择功能

4,能进行现场高程放样工作

(能计算除了超高闷纤段以外的任意桩号的中桩及任意距离边桩高程)

[SJG]------------------主程序

E“LGP”:N“JGC”:Lbi 0:A“BM”:O“HS”: “1,SZY”: “-1,QZY”:K“SZY,QZY”:{LM}:L“K”:M“ZBJ”:Prog“SQX”⊿T=RAbs(J-I)/2:C=B-L:F=1:I>J=>F=-1←┘

L ≤0=>{AO}:Goto 0 : ≠>L<B-T=>Z=0:P=I: ≠>L<B=>Z=1:P=I :≠>L<B+T=>Z=1:P=J: ≠>Z=0:P=J⊿⊿⊿⊿ ←┘

Fix 3:G=H-N-CP+ZF(T-AbsC) 2/(2R)+ME:Goto 1←┘

Lbi 1:{D}:D“QS”:Fix 3:X=A+KO-KD:G:Pause 0:X: Pause 0:Q=G-X: “+T,-W=”:Q◢

Goto 0←┘

[SQX]------------------子程序

L≤932.25(第二竖曲线起点桩号)=>H=179.17(第一竖曲线交点高程):B=160(第一竖曲线交点桩号):R=7000(第一竖陆罩宴曲线半径):I=-0.003(第一竖曲线前坡度):J=0.005(第一竖曲线后坡度):Goto 1⊿

L≤1173.96(第三竖曲线起点桩号)=>H=183.17(第二竖曲线交点高程):B=960(第二竖曲线交点桩号):R=18000(第二竖曲线半径):I=0.005(第二竖曲线前坡度):J=0.00808(第二竖曲线后坡度):Goto 1⊿

Lbi 1

…………………………继续添加竖曲线要素(坡度上+下早银-),括号内的数据供测试程序用

说明:

选择主程序进入

1, LGP?----------输入路拱坡度,负值 (按EXE键)

2, JGC?----------输入结构层厚度(按EXE键)

3, BM?-----------输入水准点高程(按EXE键)

4, HS?------------输入后视读数(全站仪测量时带+,-号)(按EXE键)

5, 1,SZY

-1,QZY

SZY,QZY?----选择测量仪器:输入1用水准仪,输入-1用全站仪 (按EXE键)

6, K≤0 BREAK

K?--------------输入放样桩号, K值≤0时返回程序起始位置 (按EXE键)

7, ZBJ?-----------输入边桩至中桩的距离,中桩放样时输入0 (按EXE键)

8, QS?-------------输入前视读数(全站仪测量时带+,-号) (按EXE键)

显示:

9,12.3-------------设计高程

11.2-------------实测高程

+T,-W

1.1(填挖量,+填-挖)

(按EXE键 重复6,7,8,9步,进行下一个桩点的放样工作.

---