mv-pereira / forensic-gpm-trajectories Goto Github PK

Ao invés de Gravar x,y,z é melhor gravar latitude/longitude e altura

altura = 60;
descendente = 0;
phi = 1M_PI/180;
gamma = 183M_PI/180;
v = 230;
m=10.240/1000.0;
a=M_PIpowl((8.82/1000.0),2)/4; // A divisão por quatro leva em conta o raio.
//a=log(M_PI/4)+2log(8.82/1000);
dextrogiro=1;
c = 0.235800;
vento = 10/3.6;
ang = (100+180.0)M_PI/180.0;
latitude = -8.127727M_PI/180.0;
longitude = -34.898383;
azimute = gamma;

a ver com #28

Deixar todos as constantes coerentes (inclusive entre as imagens).

Algum critério pra acelerar correção de θ, Azimute...

altura = 60;
descendente = 1;
phi = 0M_PI/180;
gamma = 183M_PI/180;
v = 230;
m=10.240/1000.0;
a=M_PIpowl((8.82/1000.0),2)/4; // A divisão por quatro leva em conta o raio.
//a=log(M_PI/4)+2log(8.82/1000);
dextrogiro=1;
c = 0.235800;
vento = 10/3.6;
ang = (100+180.0)M_PI/180.0;
latitude = -8.127727M_PI/180.0;
longitude = -34.898383;
azimute = gamma;

ao invés de ir de ±0.001 ou algo assim, avançar em termos de médias:
ex:
theta est = 5
phi simulado = 8
phi = 10

-> theta = theta+(phi-phisimu)/2 ou /x...

π/2 < ϕ < -π/2, mas se ϕ -> 0 pela esquerda, o programa não finaliza.

É possível que seja problema de aproximação/resolução.

transformar as 3 etapas em um grande while

• Testar ângulos diversos de impactação (phi).
• Trajetórias ascendentes e descendentes.
• Disparos descendentes.

Tiro completamente descendente, o programa transforma em ascendente pra dar o resultado.

Refatorar Cd como Trajetórias balísticas.

Inserir uma opção para marcar a localização de uma edificação entre a origem do disparo e a impactação. Com a localização (coordenadas) traçar se é possível o projétil ter sua origem nessa edificação, em algum pavimento.

Relacionado ao issue #26 .

Checar se há problema de descontinuidade em 0 e em, hipoteticamente, outros valores de theta.

para phi

== 0+: Origem (NMM) -8.121949, -34.898103

== 0-: Origem (NMM) -8.122646, -34.898132

considerando os parâmetros de debug.

Colocar prints do google earth, ensinado essa parte dele pra medir altura. Usar o Holliday como ponto de partida -8.1236011, -34.8981805

Função de Cd pelo número de Mach em G1, G2... G7. (McCoy, Pg 52)

Para as funções, laços, etc...

Apenas o y0 variaria, nada mais?

Fazer ou explicar o uso do programa Cd.

Azimute inicial pode ser simplesmente γ. γ é a melhor estimativa...

Se mudar isso, mudar o Readme.md.