102 lines
2.0 KiB
Scilab
102 lines
2.0 KiB
Scilab
res_log=read('../data/gpscntl_log_0.m',-1,11);
|
|
pvt_log=read('../data/gps_navfix_log.m',-1, 16);
|
|
[r_n,c_n]=size(res_log);
|
|
idx=[1:r_n];
|
|
|
|
|
|
|
|
//xset('color',3);
|
|
//xset('background', 1);
|
|
xbasc();
|
|
// plot I result
|
|
xset('window',0);
|
|
subplot(221);
|
|
plot2d(res_log(:,2));
|
|
xtitle("I Output",'ms','Mag');
|
|
xset('background', 4);
|
|
|
|
// plot Q result
|
|
subplot(223);
|
|
plot2d(res_log(:,3));
|
|
xtitle("Q Output",'ms','Mag');
|
|
xset('background', 4);
|
|
|
|
subplot(222);
|
|
plot2d(res_log(:,6)-(4.1304e6));
|
|
xtitle("Err of carrier freq", 'ms','Hz');
|
|
xset('background', 4);
|
|
|
|
subplot(224);
|
|
plot2d(res_log(:,8)/2-1.023e6);
|
|
xtitle("Err of code freq", 'ms', 'Hz');
|
|
|
|
xset('background', 4);
|
|
|
|
xset('window',1);
|
|
plot2d(idx, [res_log(:,9), res_log(:,10), res_log(:,11)], leg="Early@Prompt@Late");
|
|
xtitle("Correlation result for Early, Prompt, Late phase");
|
|
|
|
|
|
xset('window',2);
|
|
subplot(311);
|
|
plot2d(pvt_log(:,1));
|
|
xtitle("ECEF position x",'TIC','m');
|
|
xset('background',4);
|
|
|
|
subplot(312);
|
|
plot2d(pvt_log(:,2));
|
|
xtitle("ECEF position y",'TIC','m');
|
|
xset('background',4);
|
|
|
|
subplot(313);
|
|
plot2d(pvt_log(:,3));
|
|
xtitle("ECEF position z",'TIC','m');
|
|
xset('background',4);
|
|
|
|
xset('window',3);
|
|
subplot(311);
|
|
plot2d(pvt_log(:,8));
|
|
xtitle("ECEF velocity x",'TIC','m/s');
|
|
xset('background',4);
|
|
|
|
subplot(312);
|
|
plot2d(pvt_log(:,8));
|
|
xtitle("ECEF velocity y",'TIC','m/s');
|
|
xset('background',4);
|
|
|
|
subplot(313);
|
|
plot2d(pvt_log(:,10));
|
|
xtitle("ECEF velocity z",'TIC','m/s');
|
|
xset('background',4);
|
|
|
|
xset('window',4);
|
|
subplot(211);
|
|
plot2d(pvt_log(:,11));
|
|
xtitle("Clock drift rate", 'TIC', '1/s');
|
|
xset('background',4);
|
|
|
|
subplot(212);
|
|
plot2d(pvt_log(:,12));
|
|
xtitle("Clock bias", 'TIC', 's');
|
|
xset('background',4);
|
|
|
|
xset('window',5);
|
|
subplot(411);
|
|
plot2d(pvt_log(:,13));
|
|
xtitle("hdop", 'TIC', '1/s');
|
|
xset('background',4);
|
|
|
|
subplot(412);
|
|
plot2d(pvt_log(:,14));
|
|
xtitle("vdop", 'TIC', 's');
|
|
xset('background',4);
|
|
subplot(413);
|
|
plot2d(pvt_log(:,15));
|
|
xtitle("pdop", 'TIC', '1/s');
|
|
xset('background',4);
|
|
|
|
subplot(414);
|
|
plot2d(pvt_log(:,16));
|
|
xtitle("gdop", 'TIC', 's');
|
|
xset('background',4);
|