GPS数据解析

clear; close all;



x1 = [3427819.8077

3427819.5156

3427819.5964

3427819.5864

3427819.3460

3427818.9751

3427819.1378

3427819.0409

3427819.3261

3427819.1508

3427819.1148 

3427819.1545 

3427819.5146

3427819.1058

3427819.1268

3427819.3561

3427819.2898

3427819.1150 

3427819.4789

3427818.9605

3427819.6698

3427819.4966

3427819.0015];





y1 = [ 603663.8520

603663.9057

 603663.8764

 603663.7779

 603663.6453

 603663.8508

 603663.7367

 603663.8545

 603664.0401

 603663.9542

603663.8564 

603663.7138 

603663.7253

603664.0669

603664.2620

603664.1093 

 603664.1499

603664.1756 

 603664.2031

 603664.0892

 603664.0161

 603663.8148

 603664.0675];

 

z1 = [  5326881.9105

 5326881.5277

  5326881.2313

  5326880.6163

  5326881.7420

  5326881.5032

  5326880.9846

  5326880.1988

 5326880.6466

  5326880.8853

 5326880.9158 

5326880.9822 

 5326881.4205

 5326881.0524

 5326881.3719 

 5326882.1826 

 5326881.9534

5326880.9812 

 5326881.4709

  5326880.9796

 5326881.5869

 5326881.4339

  5326880.4310];



 [lat1,lon1,h1] = ecef2geodetic(referenceEllipsoid('wgs84'),mean(x1),mean(y1),mean(z1));

  

  

 [xEast,yNorth,zUp] = ecef2enu(x1,y1,z1,lat1,lon1,h1,referenceEllipsoid('wgs84'));

  

figure;

subplot(3, 1, 1);

plot(xEast, 'b.-'); axis tight; grid on;

subplot(3, 1, 2);

plot(yNorth, 'b.-'); axis tight; grid on;

subplot(3, 1, 3);

plot(zUp, 'b.-'); axis tight; grid on;



 range(x1),  range(y1),  range(z1)

 range(xEast),  range(yNorth), range(zUp)

GPS数据解析，再配合函数lla2ecef使用。

我转载的这个链接

http://lordleid.lofter.com/post/7546f_a54ad1c

这个人写了很多关于matlab的东西，可以看看。