Здравствуйте, столкнулся с такой проблемой, что после преобразования, лишь Y совпадает по графику с начальными данными. Были ли у кого-то ещё такие проблемы и как вы её решили?
Код MATLAB:
Код MATLAB:
Цитата |
---|
a = 6378245.0; % Большая (экваториальная) полуось b = 6356863.019; % Малая (полярная) полуось aj = 1/298.3; %сжатие эллипсоида e = (2*aj - aj^2)^2; for i=1:n siB(i) = sin(B(i)); siL(i) = sin(L(i)); coB(i) = cos(B(i)); coL(i) = cos(L(i)); siBd(i) = siB(i).^2; coz(i) = e*siBd(i); N(i) = a / sqrt(1 - coz(i)); N = N.'; coz = coz.'; x(i) = (N(i) + H(i))*coB(i)*coL(i); y(i) = (N(i) + H(i))*coB(i)*siL(i); z(i) = [(1 - e)*N(i) + H(i)]*(siB(i)); end |