工程测量学 变形监测沉降曲线绘制 以及简单曲线元素的计算 MatLab 实现

一、沉降曲线

Average_Settlement = [[0,0.858750000000000,3.29250000000000,4.28250000000000,6.64375000000000,8.33000000000000,10.9912500000000,14.9975000000000]];
Average_Velocity = [0 0.053672 0.135208 0.070714 0.131181 0.099191 0.076036 0.072841];
Date = [68 84 102 116 134 151 186 241];

x = Date
y1 = Average_Settlement
y2 = Average_Velocity
[AX,H1,H2] = plotyy(x,y1,x,y2,'plot');
set(get(AX(1),'Ylabel'),'String','累计沉降 mm') 
set(get(AX(2),'Ylabel'),'String','沉降速度 mm/Day') 
xlabel('年积日/2013') 
title('沉降曲线') 
set(H1,'LineStyle','--','LineWidth',2)
set(H2,'LineStyle','-.','LineWidth',2)
legend('累计沉降','沉降速度')
grid on

二、简单圆曲线要素计算

%此小程序针对简单圆曲线的要素以及主点的计算
R = 490 %已知半径 m
alpha = (31/3600 + 58/60 + 26) * (pi / 180) %已知转角 弧度的形式
%Dis_Circlr = % 圆曲线细部点间距
Kjd = 523.6 % 已知里程点桩号（仅有数字字母自己加即可）
T = R * tan(alpha / 2);% 切线长度 tan(弧度 / 2)
L = R * alpha ;% 圆曲线长度 长度乘弧度的计算格式
E = R * (sec(alpha / 2) - 1);%外矢距
q = 2 * T - L; %切曲线

% 主点里程的计算
Kzy = Kjd - T; %直圆点
Kyz = Kzy + L; %圆直点 
Kqz = Kyz - L / 2;%曲中点

%校核
Kjd2 = Kqz + q / 2;%检查计算的与已知的误差