목차
1. Title
2. Introduction
3. Materials and Methods
4. Results
5. Discussion and Conclusion
2. Introduction
3. Materials and Methods
4. Results
5. Discussion and Conclusion
본문내용
= cp_data(1,3);
Zr= cp_data(1,4);
Tr=0; % initial value for time error
l=1;
while(l<=10)
%%%%%%%%%%%%%%%%%
% absolute posltlonlng
%%%%%%%%%%%%%%%%%
v1 = epoch_data(:,10);
v3 = sqrt((Xr - epoch_data(:,7)).^2+(Yr-epoch_data(:,8)).^2 + (Zr-epoch_data(:,9)).^2);
v2 = c* epoch_data(:,6);
K = v1+ v2-v3;
v21 = (Xr - epoch_data(:,7))./v3;
v22 = (Yr -epoch_data(:,8))./v3;
v23 = (Zr - epoch_data(:,9))./v3;
v24 = ones(k,1);
B = cat(2, v21,v22,v23,v24);
X = inv(B'*B)*B'*K;
Xr= Xr + X(1);
Yr= Yr + X(2);
Zr= Zr + X(3);
Tr= Tr + X(4)/c;
l = l+1;
end
dx(j) = cp_data(1,2)-Xr;
dy(j) = cp_data(1,3)-Yr;
dz(j) = cp_data(1,4)-Zr;
dt(j) = Tr;
BB = diag(inv(B'*B));
PDOP(j) = norm(BB(1:3));
l=l+1;
G(j) = PDOP(j)*10;
%%%%%%%%%%%%%%%%%
% PDOP
%%%%%%%%%%%%%%%%%
end
end
t=1:1:2880; % X-axls (tlme)
t=t.*30/60/60;
%%%%%%%%%%%%%%%%%
% plot the graphes
%%%%%%%%%%%%%%%%%
subplot(2,2,1);
plot(t,dx,'b');
hold on
plot(t,G-30,'r');
subplot(2,2,2);
plot(t,dy,'b');
hold on
plot(t,G-30,'r');
subplot(2,2,3);
plot(t,dz,'b');
hold on
plot(t,G-30,'r');
subplot(2,2,4);
plot(t,dt,'b');
hold on
plot(t,G-30,'r');
disp('lf you want to contlnue, press any number');
disp('Or you want to qult, press 0');
q=input('q= ');
end
Zr= cp_data(1,4);
Tr=0; % initial value for time error
l=1;
while(l<=10)
%%%%%%%%%%%%%%%%%
% absolute posltlonlng
%%%%%%%%%%%%%%%%%
v1 = epoch_data(:,10);
v3 = sqrt((Xr - epoch_data(:,7)).^2+(Yr-epoch_data(:,8)).^2 + (Zr-epoch_data(:,9)).^2);
v2 = c* epoch_data(:,6);
K = v1+ v2-v3;
v21 = (Xr - epoch_data(:,7))./v3;
v22 = (Yr -epoch_data(:,8))./v3;
v23 = (Zr - epoch_data(:,9))./v3;
v24 = ones(k,1);
B = cat(2, v21,v22,v23,v24);
X = inv(B'*B)*B'*K;
Xr= Xr + X(1);
Yr= Yr + X(2);
Zr= Zr + X(3);
Tr= Tr + X(4)/c;
l = l+1;
end
dx(j) = cp_data(1,2)-Xr;
dy(j) = cp_data(1,3)-Yr;
dz(j) = cp_data(1,4)-Zr;
dt(j) = Tr;
BB = diag(inv(B'*B));
PDOP(j) = norm(BB(1:3));
l=l+1;
G(j) = PDOP(j)*10;
%%%%%%%%%%%%%%%%%
% PDOP
%%%%%%%%%%%%%%%%%
end
end
t=1:1:2880; % X-axls (tlme)
t=t.*30/60/60;
%%%%%%%%%%%%%%%%%
% plot the graphes
%%%%%%%%%%%%%%%%%
subplot(2,2,1);
plot(t,dx,'b');
hold on
plot(t,G-30,'r');
subplot(2,2,2);
plot(t,dy,'b');
hold on
plot(t,G-30,'r');
subplot(2,2,3);
plot(t,dz,'b');
hold on
plot(t,G-30,'r');
subplot(2,2,4);
plot(t,dt,'b');
hold on
plot(t,G-30,'r');
disp('lf you want to contlnue, press any number');
disp('Or you want to qult, press 0');
q=input('q= ');
end
소개글