trace_acc.m

北京到广州铁路路线仿真

function [long,lat,v1]=trace_acc(long0,lat0,theta,t0)
%加速过程
Re=6378137;
f=1/298.157;
trans=2*pi/360;
v=50;
v0=0;
dt=10;
Tacc=600;
a=v/Tacc;
i=1;
while v0<v
    Rm=Re*(1+2*f-3*f*sin(trans*lat0)^2);%该点的子午圈半径
    Rn=Re/((1-f)^2*sin(trans*lat0)^2+cos(trans*lat0)^2);%该点的卯酉圈半径   
    deltalong=(v0*dt+0.5*a*dt)*cos(theta)/(Rm*cos(trans*lat0)*trans);%十秒间隔内经度的增加量
    deltalat=(v0*dt+0.5*a*dt)*sin(theta)/(Rn*trans);%10秒间隔内纬度的增加量
    long(i)=long0+deltalong;%实时计算每个点的经纬度
    lat(i)=lat0+deltalat;
    t0=t0+dt; 
    v0=v0+a*dt;%
    v1(i)=v0;
    long0=long(i);
    lat0=lat(i);
    i=i+1;
end