18+ О жизни Артёма... дома и не только

Дабы не отставать от прогрессивной общественности, в копилку своих знаний этот пост забираю... Как раз у меня много чего можно мониторить между камеральными "кофе-брейками"...

ЗЫ... Осваиваюсь с облачным хранением медиа..

 

Бестолковый "инстаграм" какой-то: все квадратное и на рабочем ПК не открывается...

 

Ну, как бы намек на ограниченность ресурсов при фотоскринах уловил, а привычный импульс неправильно публикуется с точки зрения серийных снимков...

 

Учитывая колоритность участника, за развитием сюжета послежу.. Ну и помогу, ежели моих знаний хватит...

...

--- Сообщения объединены, 28 апр 2018, Оригинальное время сообщения: 28 апр 2018 ---

ЗЫ. До кучи притянул, хотя под ногами будет путаться "матлабовский" код метода четвертого порядка

 

Трофей дня...
Всё... Я сам себе всё доказал: при установке на пункт триангуляции второго класса "Гудзево" заранее откалиброванный на один пункт второразрядной полигонометрии в топоцентрической проекции МСК23 инструмент выдает такую вот поисковую точность для марки (включен режим выноса точки в РТК от ЭФТ-ПДБС г.Крымск)

 

Тут выход на цель в пять сантиметров плановые значения и сантиметров шестьдесят-семьдесят по высоте. Конкретно этот пункт имеет характеристику "тригонометрическое", т.е. высотная отметка сойтись и не должна была

 

По дороге домой помянул ГАИшников и новостей по существу начитался... Аминь...

 

Артем Скурихин а начитанное существо точно не утка? Вопрос конечно праздный - но всё таки первыми это увидят (и прочуствуют на себе) ежедневные пользователи ГНСС.

 

Он это по-русски написал?... Вроде, кирилицей...

 

Не ругайся) .. просто давно привык относиться к информации из медиа и интернета относиться критически и с известной долей недоверия.. если честно сильно цинично смотрю на подобные статьи

Оффтоп

меня просто много лет на полит занятиях просвещали - ну нету у нас союзников, тут те, тут эти а вот тут из-за угла вон те облизываются и ждут когда мыло уронят .. потому сильно сомневаюсь что МО РФ так уж сильно захочет делиться своими спутниками (пусть даже и с попутчиками)

вроде да ..

 

Никак не привыкну... Это где-то в шаговой доступности от моего дома (из маршрутного такси вылез - щёлкнул)...

 

Следующая цель:

 

Есть трофей!...

Выхожу по темноте... Домой добираться ещё ой-ой.. но как красив Новороссийск в ночных огнях!...

 

Из вчерашнего:

 

Общая атмосфера праздника...

 

Суточные "приключения" спутника ГЛОНАСС по версии частного кода алгоритма в среде Matlab

Код:

% ГЛОНАСС НКА координаты
clear
clc
% это все время
T = 1; % [ с ]
% Tk = 56640; % [ с ] % не используется
% Day = 1031 ; % [ день ] % не используется
Tb = 56700; % [ с ]
Gn = -4.5474735088646410e-012; % [ с/с ] % с / с - сокращаются и эта штука по идее безразмерная
Tn = -7.0374459028244020e-005; % [ с ]
% координаты
X = 1.7244073242187500e+004; % [ км ]
Y = -9.5787133789062500e+003; % [ км ]
Z = 1.6202982910156250e+004; % [ км ]
% скорость
Xt = -1.4764909744262700e+000; % [ км/с ]
Yt = 1.6627349853515620e+000; % [ км/с ]
Zt = 2.5435619354248050e+000; % [ км/с ]
% по размерности получается что это ускорение
Xtt= 0.0000000000000000e+000; % [ км/с**2 ]
Ytt= 0.0000000000000000e+000; % [ км/с**2 ]
Ztt= -9.3132257461547850e-010; % [ км/с**2 ]
% какое то тоже время
Tc = -3.8743019104003910e-007; % [ с ]
s(1:6) = 0;
Tmdv = T + Tn - Gn*(T - Tb) + Tc;
s(1) = X*1000;
s(2) = Y*1000;
s(3) = Z*1000;
s(4) = Xt*1000;
s(5) = Yt*1000;
s(6) = Zt*1000;
Xtt = Xtt*1000;
Ytt = Ytt*1000;
Ztt = Ztt*1000;
h = 1;
if Tmdv < Tb
h=-h;
end
hl = rem(Tmdv - Tb,h);
if hl==0
arrt = Tb+h:h:Tmdv;
ind = 0;
else
arrt = [Tb+h:h:Tmdv, Tmdv];
ind = 1;
end
s0 = s;
i=1;
back(1:3,1:length(arrt)+1)=0;
for t = arrt
if t==Tmdv && ind ~= 0
h = hl;
end
%display([t,h]); % эта штука выводит время и шаг. можно поставить % - и ничего не будет выводится
back(:,i) = s(1:3);
i=i+1;
arg = s;
k1 = h*fglon(arg,Xtt,Ytt,Ztt);
arg = s + k1/2;
k2 = h*fglon(arg,Xtt,Ytt,Ztt);
arg = s + k2/2;
k3 = h*fglon(arg,Xtt,Ytt,Ztt);
arg = s + k3;
k4 = h*fglon(arg,Xtt,Ytt,Ztt);
ds = (k1 + 2*k2 + 2*k3 + k4)/6;
s = s + ds;
end
back(1:3,length(arrt)+1)=s(1:3);
for i = 1:6
disp([num2str(s0(i)),' s(',num2str(i),') = ',num2str(s(i))]); % здесь выводится
% правая и левая колонки с координатами и скоростями
end
backa = reshape(back,[],1);
backb = reshape(backa,3,[]);
backc = rot90(backb);
backc = rot90(backb);
backc = rot90(backb);
% arrq = rot90(backc); % не забыть в экселе извлечь эту матрицу
 
plot = 'y'; % если стоит 'n' - график не выводится, если 'y' - график выводится
switch plot
case 'y'
plot3(back(1,:),back(2,:),back(3,:),'.r','LineWidth',3,'MarkerFaceColor','w')
hold on
plot3(0,0,0,'.r','LineWidth',3,'MarkerFaceColor','w')
[ex,ey,ez] = sphere(50);
rz=6378136;
load('topo.mat','topo','topomap1');
colormap(topomap1);
props.AmbientStrength = 0.8;
%props.DiffuseStrength = 1;
props.FaceColor= 'texture';
props.EdgeColor = 'none';
props.FaceLighting = 'phong';
props.Cdata = topo;
surface(-ex*rz,-ey*rz,ez*rz,props);
axis equal;
set(gca,'Color','none');
set(gcf,'Color','w');
%light('Position',[-1 0 0],'Style','infinite');
xlabel('X'); ylabel('Y'); zlabel('Z');
box on
view(s(1:3))
case 'n'
return
end
% это все время
T = 86400; % [ с ]
% Tk = 56640; % [ с ] % не используется
% Day = 1031 ; % [ день ] % не используется
Tb = 56700; % [ с ]
Gn = -4.5474735088646410e-012; % [ с/с ] % с / с - сокращаются и эта штука по идее безразмерная
Tn = -7.0374459028244020e-005; % [ с ]
% координаты
X = 1.7244073242187500e+004; % [ км ]
Y = -9.5787133789062500e+003; % [ км ]
Z = 1.6202982910156250e+004; % [ км ]
% скорость
Xt = -1.4764909744262700e+000; % [ км/с ]
Yt = 1.6627349853515620e+000; % [ км/с ]
Zt = 2.5435619354248050e+000; % [ км/с ]
% по размерности получается что это ускорение
Xtt= 0.0000000000000000e+000; % [ км/с**2 ]
Ytt= 0.0000000000000000e+000; % [ км/с**2 ]
Ztt= -9.3132257461547850e-010; % [ км/с**2 ]
% какое то тоже время
Tc = -3.8743019104003910e-007; % [ с ]
s(1:6) = 0;
Tmdv = T + Tn - Gn*(T - Tb) + Tc;
s(1) = X*1000;
s(2) = Y*1000;
s(3) = Z*1000;
s(4) = Xt*1000;
s(5) = Yt*1000;
s(6) = Zt*1000;
Xtt = Xtt*1000;
Ytt = Ytt*1000;
Ztt = Ztt*1000;
h = 1;
if Tmdv < Tb
h=-h;
end
hl = rem(Tmdv - Tb,h);
if hl==0
arrt = Tb+h:h:Tmdv;
ind = 0;
else
arrt = [Tb+h:h:Tmdv, Tmdv];
ind = 1;
end
s0 = s;
i=1;
forw(1:3,1:length(arrt)+1)=0;
for t = arrt
if t==Tmdv && ind ~= 0
h = hl;
end
%display([t,h]); % эта штука выводит время и шаг. можно поставить % - и ничего не будет выводится
forw(:,i) = s(1:3);
i=i+1;
arg = s;
k1 = h*fglon(arg,Xtt,Ytt,Ztt);
arg = s + k1/2;
k2 = h*fglon(arg,Xtt,Ytt,Ztt);
arg = s + k2/2;
k3 = h*fglon(arg,Xtt,Ytt,Ztt);
arg = s + k3;
k4 = h*fglon(arg,Xtt,Ytt,Ztt);
ds = (k1 + 2*k2 + 2*k3 + k4)/6;
s = s + ds;
end
forw(1:3,length(forw)+1)=s(1:3);
for i = 1:6
disp([num2str(s0(i)),' s(',num2str(i),') = ',num2str(s(i))]); % здесь выводится
% правая и левая колонки с координатами и скоростями
end
forwa = reshape(forw,[],1);
forwb = reshape(forwa,3,[]);
forwc = rot90(forwb);
forwc = rot90(forwc);
forwc = rot90(forwc);
forwc = fliplr(forwc)
% arrq = rot90(forwc); % не забыть в экселе извлечь эту матрицу
 
plot = 'y'; % если стоит 'n' - график не выводится, если 'y' - график выводится
switch plot
case 'y'
plot3(forw(1,:),forw(2,:),forw(3,:),'.r','LineWidth',3,'MarkerFaceColor','w')
hold on
plot3(0,0,0,'.r','LineWidth',3,'MarkerFaceColor','w')
[ex,ey,ez] = sphere(50);
rz=6378136;
load('topo.mat','topo','topomap1');
colormap(topomap1);
props.AmbientStrength = 0.8;
%props.DiffuseStrength = 1;
props.FaceColor= 'texture';
props.EdgeColor = 'none';
props.FaceLighting = 'phong';
props.Cdata = topo;
surface(-ex*rz,-ey*rz,ez*rz,props);
axis equal;
set(gca,'Color','none');
set(gcf,'Color','w');
%light('Position',[-1 0 0],'Style','infinite');
xlabel('X'); ylabel('Y'); zlabel('Z');
box on
view(s(1:3))
case 'n'
return
end
 

ЗЫ. А это фрагмент кода на фортране, куда внесена правка для возможной генерации ГЛОНАССовской части в условный проект "СИГОГ-ГЛОНАСС"

Код:

...
C fecha=creation date of RINEX files
C fileop=file with satellite-receiver distances for sat prn #
C idnval=array with available satellite at one observation epoch
C estac=receiver name
C trer=generate tropospheric error? S/N
C End definition variables
c=299792458.0d0
PI=4.0D0*DATAN(1.0D0)
C f1=1575.42D6
C f2=1227.60D6
f1=1602.01D6
f2=1246.01D6
 

C To generate different types of observables:
...
 
 

 

Старшенькую взял в горы посмотреть-проверить её выносливость... Запретил ей всякую хрень в рюкзак пихать... И что бы вы думали?...