Upload
buitram
View
217
Download
1
Embed Size (px)
Citation preview
Caracterização da Qualidade Geométrica das Vias-Férreas
CARACTERIZAÇÃO DA QUALIDADE GEOMÉTRICA DE VIAS FERROVIÁRIAS
MARIA CLARA RIBEIRO DE MOURA
Dissertação submetida para satisfação parcial dos requisitos do grau de
MESTRE EM ENGENHARIA CIVIL — ESPECIALIZAÇÃO EM VIAS DE COMUNICAÇÃO
Orientador: Professor Doutor António José Fidalgo do Couto
Co-Orientador: Professora Doutora Cecília Maria Nogueira Alvarenga Santos do Vale
JUNHO DE 2010
Caracterização da Qualidade Geométrica das Vias-Férreas
iii
MESTRADO INTEGRADO EM ENGENHARIA CIVIL 2009/2010
DEPARTAMENTO DE ENGENHARIA CIVIL
Tel. +351-22-508 1901
Fax +351-22-508 1446
Editado por
FACULDADE DE ENGENHARIA DA UNIVERSIDADE DO PORTO
Rua Dr. Roberto Frias
4200-465 PORTO
Portugal
Tel. +351-22-508 1400
Fax +351-22-508 1440
� http://www.fe.up.pt
Reproduções parciais deste documento serão autorizadas na condição que seja mencionado o Autor e feita referência a Mestrado Integrado em Engenharia Civil - 2009/2010 - Departamento de Engenharia Civil, Faculdade de Engenharia da Universidade do Porto, Porto, Portugal, 2009.
As opiniões e informações incluídas neste documento representam unicamente o ponto de vista do respectivo Autor, não podendo o Editor aceitar qualquer responsabilidade legal ou outra em relação a erros ou omissões que possam existir.
Este documento foi produzido a partir de versão electrónica fornecida pelo respectivo Autor.
Caracterização da Qualidade Geométrica das Vias-Férreas
i
AGRADECIMENTOS
Aos meus pais, António e Clara, pelo esforço e dedicação, por terem sido os principais responsáveis pelo que hoje sou e por tudo o que consegui alcançar. Obrigada!
Á Professora Doutora Cecília Maria Nogueira Alvarenga dos Santos Vale, pela orientação neste trabalho, pelo rigor, detalhe que depositou neste projecto, pelo material facultado, pelo acompanhamento, disponibilidade e pela atenção dispensada.
Ao Professor Doutor Rui Artur Bártolo Calçada pela disponibilidade, simpatia e sugestão do tema.
Aos professores da secção de Vias de Comunicação, pela amizade que criaram entre alunos e professores. A todos os amigos de curso que de uma forma ou de outra contribuíram para a realização desta tese.
Às amigas, Lara, Ana e Marta pela amizade e companheirismo.
Ao Manuel pelo interesse, dedicação e apoio.
Ao Zé pela inestimável amizade e apoio.
Caracterização da Qualidade Geométrica das Vias-Férreas
iii
RESUMO
O tema do presente trabalho, “Caracterização da qualidade geométrica de vias ferroviárias”, foi proposto no âmbito da realização da dissertação, do curso de Engenharia Civil, da Faculdade de Engenharia da Universidade do Porto.
Este trabalho tem como objectivo principal identificar metodologias para a caracterização da qualidade geométrica de vias-férreas.
Neste trabalho, caracteriza-se um pequeno troço de via da Linha do Norte à luz dos parâmetros geométricos para avaliação da qualidade geométrica das vias-férreas segundo a norma europeia EN 13848-5 2008 “Railway application – Track - Track geometry quality – Part 5: Geometric quality levels”[1].
Por fim são apresentadas conclusões e desenvolvimentos para trabalhos futuros.
Palavras-chave
Via-férrea, Qualidade geométrica, Alinhamento.
Caracterização da Qualidade Geométrica das Vias-Férreas
v
ABSTRACT
The theme of the present work, “Caracterização da qualidade geométrica de vias ferroviárias”, was proposed in order to realize this dissertation, of the Civil Engineering course of the Faculdade de Engenharia da Universidade do Porto.
This has as primary objective to identify methodologies in order to characterize the geometric quality of railways.
On the present work is characterized a small section of railway, belonging to the Linha do Norte , under the geometric parameters of railways quality evaluation according to the European norm EN 13848-5 2008 “Railway application – Track - Track geometry quality – Part 5: Geometric quality levels”[1].
For last, conclusions and developments are presented in order to be used for further studies.
Keywords
Railway, Geometric quality, Alignment.
Caracterização da Qualidade Geométrica das Vias-Férreas
vii
ÍNDICE GERAL
1 INTRODUÇÃO ..................................................................................................................... 1
1.1 Enquadramento Histórico .......................................................................................................... 1
1.2 Objectivos ...................................................................................................................................... 3
1.3 Estrutura do Trabalho ................................................................................................................ 3
2 Qualidade Geométrica da Via-férrea .......................................................... 5
2.1 Sumário .......................................................................................................................................... 5
2.2 Inspecção da Qualidade Geométrica da Via ......................................................................... 5
2.3 Parâmetros para a caracterização da qualidade d a via-férrea ........................................ 6
2.3.2. NIVELAMENTO LONGITUDINAL..................................................................................................... 7
2.3.3 NIVELAMENTO TRANSVERSAL ...................................................................................................... 7
2.3.4 ALINHAMENTO ............................................................................................................................... 8
2.3.5 EMPENO ........................................................................................................................................ 9
2.4 Caracterização da qualidade geométrica da Via-f érrea segundo a norma EN 13848-5 2008 “Railway application – Track - Track geometr y quality – Part 5: Geometric quality levels” ...................................................................................................................................... 9
2.4.1 LIMITES DE SEGURANÇA ............................................................................................................. 11
2.4.2 LIMITES DE INTERVENÇÃO .......................................................................................................... 13
2.4.3 LIMITES DE ALERTA .................................................................................................................... 15
3 Caso de Estudo: Avaliação da qualidade geométrica de um troço de via da Linha do Norte ................................................................. 19
3.1 Sumário ........................................................................................................................................ 19
3.2 Caracterização do troço ferroviário em estudo ................................................................. 19
3.3 Avaliação da qualidade Geométrica do troço em e studo ............................................... 19
3.3.1 INFLUÊNCIA DO TRAÇADO EM PLANTA ....................................................................................... 20
3.3.1.1. Segmento A .............................................................................................................................. 20
3.3.1.2. Segmento b ............................................................................................................................... 24
3.3.1.3 Comparação entre os segmentos A e B ..................................................................................... 28
Caracterização da Qualidade Geométrica das Vias-Férreas
viii
3.4 EVOLUÇÃO DO ALINHAMENTO AO LONGO DO TEMPO .................................................................. 28
3.4.1 METODOLOGIA PARA A ANÁLISE DA EVOLUÇÃO DO DESVIO PADRÃO DO ALINHAMENTO NO
TEMPO ................................................................................................................................................... 30
3.4.2 ANÁLISE E DISCUSSÃO DE RESULTADOS ................................................................................... 33
4 Conclusões e Desenvolvimentos Futuros ........................................ 45
Caracterização da Qualidade Geométrica das Vias-Férreas
ix
ÍNDICE DE FIGURAS
Figura 1-veículo de inspecção Plasser&Theurer EM120 ........................................................................ 6
Fonte: REFER ......................................................................................................................................... 6
Figura 2-Bitola ......................................................................................................................................... 7
Figura 3-NivelamentoTransversal ........................................................................................................... 8
Figura 4-Alinhamento .............................................................................................................................. 8
Figura 5-Empeno ..................................................................................................................................... 9
Figura 6- Função sinusoidal .................................................................................................................. 10
Figura 7- Bitola ...................................................................................................................................... 20
Figura 8-Nivelamento Longitudinal D1 Esquerdo e Direito ................................................................... 21
Figura 9-Detalhe do Nivelamento Longitudinal D1 Esquerdo e Direito ................................................ 21
Figura 10-Nivelamento Longitudinal D2 Esquerdo e Direito ................................................................. 22
Figura 11-Nivelamento Transversal ...................................................................................................... 22
Figura 12-Alinhamento D1 Esquerdo e Direito ..................................................................................... 23
Figura 13-Detalhe do Alinhamento D1 Esquerdo e Direito ................................................................... 23
Figura 14-Bitola ..................................................................................................................................... 24
Figura 15-Nivelamento Longitudinal D1 Esquerdo e Direito ................................................................. 25
Figura 16-Nivelamento Longitudinal D2 Esquerdo e Direito ................................................................. 25
Figura 17-Nivelamento Transversal ...................................................................................................... 26
Figura 18-Nivelamento transversal e irregularidades de nivelamento transversal ............................... 26
Figura 19-Nivelamento transversal, comparação com a norma ........................................................... 26
Figura 20-Alinhamento D1 Esquerdo e Direito ................................................................................... 287
Figura 21-Variação do Alinhamento do carril direito no tempo ............................................................. 29
Figura 22-Variação do Alinhamento do carril esquerdo no tempo...................................................... 309
Figura 23- 2007 Dez Vs 2009 Jan, sem ajuste temporal ...................................................................... 30
Figura 24 - 2007 Dez Vs 2009 Jan, com ajuste temporal ..................................................................... 31
Figura 25 - variação do desvio padrão no tempo para o trecho 57 ...................................................... 34
Figura 26-variação do desvio padrão no tempo para o trecho 60 ........................................................ 34
Figura 27 - variação do desvio padrão no tempo para o trecho 121 .................................................... 35
Caracterização da Qualidade Geométrica das Vias-Férreas
x
Figura 28 - variação do desvio padrão no tempo para o trecho 126 .................................................... 35
Figura 29-variação do desvio padrão no tempo para o trecho 5 .......................................................... 36
Figura 30-variação do desvio padrão no tempo para o trecho 6 .......................................................... 36
Figura 31-variação do desvio padrão no tempo para o trecho 7 .......................................................... 36
Figura 32-variação do desvio padrão no tempo para o trecho 12 ........................................................ 37
Figura 33-variação do desvio padrão no tempo para o trecho 23 ........................................................ 37
Figura 34-variação do desvio padrão no tempo para o trecho 32 ........................................................ 38
Figura 35-variação do desvio padrão no tempo para o trecho 24 ........................................................ 38
Figura 36-variação do desvio padrão no tempo para o trecho 39 ........................................................ 39
Figura 37-variação do desvio padrão no tempo para o trecho 41 ........................................................ 39
Figura 38variação do desvio padrão no tempo para o trecho 68 .......................................................... 40
Figura 39-variação do desvio padrão no tempo para o trecho 69 ........................................................ 40
Figura 40-variação do desvio padrão no tempo para o trecho 92 ........................................................ 41
Figura 41-variação do desvio padrão no tempo para o trecho 107 ...................................................... 41
Figura 42-variação do desvio padrão no tempo para o trecho 121 ...................................................... 42
Figura 43-variação do desvio padrão no tempo para o trecho 129 ...................................................... 42
Figura 44-variação do desvio padrão no tempo para o trecho 140 ...................................................... 43
Caracterização da Qualidade Geométrica das Vias-Férreas
xi
INDICE DE TABELAS
Tabela 1-Bitola - defeitos isolados ........................................................................................................ 11
Tabela 2-Bitola ...................................................................................................................................... 12
Tabela 3-Nivelamento Longitudinal – defeitos isolados – Diferença entre a média e o valor de pico . 12
Tabela 4-Nivelamento Transversal – defeitos isolados - Diferença entre a média e o valor de pico ... 12
Tabela 5-Alinhamento – defeitos isolados - Diferença entre a média e o valor de pico ....................... 13
Tabela 6-Bitola - defeitos isolados ........................................................................................................ 13
Tabela 7-Bitola ...................................................................................................................................... 14
Tabela 8-Nivelamento Longitudinal - defeitos isolados – Diferença entre a média e o valor de pico .. 14
Tabela 9-Nivelamento Transversal - defeitos isolados – Diferença entre a média e o valor de pico ... 14
Tabela 10-Alinhamento – defeitos isolados - Diferença entre a média e o valor de pico ..................... 15
Tabela 11-Bitola - defeitos isolados ...................................................................................................... 15
Tabela 12 - Bitola .................................................................................................................................. 16
Tabela 13-Nivelamento Longitudinal - defeitos isolados – Diferença entre a média e o valor de pico 16
Tabela 14-Nivelamento Transversal – Desvio Padrão - Diferença entre a média e o valor de pico .... 16
Tabela 15-Alinhamento - Defeitos Isolados - Diferença entre a média e o valor de pico ..................... 17
Tabela 16 – Alinhamento - Desvio Padrão ........................................................................................... 17
Tabela 17 – Classificação dos trechos ...…………………………………………………………………....43
Tabela 18 – Taxa de Degradação ……………………………………………………………………………44
Caracterização da Qualidade Geométrica das Vias-Férreas
1
1 INTRODUÇÃO
1.1 ENQUADRAMENTO HISTÓRICO
Na segunda metade do século XIX, era discutida a melhor forma de modernizar o país por elementos da elite intelectual, política e económica. Muitos defendiam a necessidade de construção de vias de comunicação, mas Portugal ainda não se tinha recuperado das convulsões políticas e das guerras civis que enfrentou, não conseguindo, por isso, obter o capital necessário para tal investimento.
Em Dezembro de 1844 foi fundada a Companhia das Obras Públicas em Portugal, cujo principal objectivo era a construção do caminho-de-ferro de Lisboa até à fronteira espanhola, com a intenção de ligar o país à Europa. Não se obtendo nenhum resultado prático, só depois de 1851 foram reunidas as condições necessárias para iniciar a construção da rede ferroviária nacional.[2]
A linha de caminho-de-ferro entre Lisboa (Santa Apolónia) e Carregado foi inaugurada a 28 de Outubro de 1856 [2], mas só em Março de 1857 se deu a dissolução da Companhia Central Peninsular dos Caminhos-de-Ferro de Portugal. Sendo que a primeira Companhia Real dos Caminhos-de-ferro Portugueses foi formalmente constituída e declarada em 1860, após a realização de 40 por cento do seu capital, podendo ser considerada como o primeiro antepassado da CP – Comboios de Portugal. [3]
Entretanto, nasce a Companhia dos Caminhos-de-Ferro ao Sul, que se propõe explorar uma linha entre Aldeia Galega e Vendas Novas. Posteriormente a origem passa para o Barreiro e são previstas extensões para Setúbal, Évora e Beja. Em 1854 foi aberto o concurso para a construção das linhas, e em Maio de 1858, os comboios da Companhia já circulavam até Bombel com bitola de 1,44 metros [3]. A 16 de Julho de 1859 foi apresentado pelo governo o programa de concurso para o prolongamento da construção até Évora e Beja. Não havendo interessados, a 3 de Janeiro de 1860, o governo contratou a Companhia dos Caminhos-de-Ferro do Sueste para a construção e exploração dos prolongamentos. [4]
A Linha do Norte, que teve a sua origem ao quilómetro 106,3 da linha do Leste, ficou completa até às Devesas, na margem esquerda do rio Douro, a 7 de Julho de 1864, estabelecendo assim a ligação directa de Lisboa com a “capital” do Norte – o Porto. [4]
Caracterização da Qualidade Geométrica das Vias-Férreas
2
A Linha do Douro, de via larga, ligava Ermesinde, na Linha do Minho, a Barca D´Alva, numa extensão total de 212 km. Na segunda metade do século XIX a região do Douro surgia como uma importante fonte de riqueza para o país. O aumento da produção de Vinho do Porto e a necessidade de melhorar o seu transporte, tanto em termos económicos como em segurança e rapidez, nomeadamente entre a Régua e o Porto, levaram a que se começasse a pensar no caminho-de-ferro como alternativa aos meios de transporte existentes. Em 1862, o Governo incumbiu Sousa Brandão de fazer o reconhecimento da região e proceder aos primeiros estudos de terreno entre Porto e a Régua [3]. Em 1864 [4] os trabalhos de Sousa Brandão incluíam já o anteprojecto do prolongamento da linha da Régua até Barca D´Alva, com ligação a Salamanca e Fuentes d´Oñoro, em Espanha. O estudo ficou concluído em Julho de 1867 [4] e o governo decidiu-se pela construção de duas linhas, a do Minho e a do Douro. As linhas a norte do Douro faziam assim a ligação do Porto a Braga e à fronteira da Galiza e uma outra ao longo do Rio Douro, que ligava o Porto a Pinhão.
É de salientar que uma das maiores obras ferroviárias portuguesas do século XIX, foi a construção da linha de ferro que ligava Barca D´Alva a La Frageneda e a Salamanca. Esta linha internacional era de via larga e ligava a linha do Douro à rede ferroviária espanhola. Os trabalhos foram iniciados no ano de 1882 [3], depois de muitos avanços e recuos, nomeadamente em relação à localização da ponte internacional sobre o rio Águeda. Foi inaugurada a 9 de Dezembro de 1887 [5] com a abertura à exploração pública do serviço directo Porto – Salamanca. Devido às difíceis características do terreno próximo da fronteira, na margem direita do Rio Águeda, foi necessário executar inúmeras obras de arte: 20 pontes (13 das quais nos 17 km de linha mais próximos da fronteira) e 19 túneis entre Barca D´Alva e La Frageneda, e tudo financiado e construído pelos portugueses. Até ficar concluída a rede ferroviária nacional, foi preciso mais de meio século.
Ao longo da segunda metade do século XIX, várias linhas iam sendo estudadas, construídas e exploradas, por iniciativa do Estado ou de empresas particulares. Como já foi referido, o caminho-de-ferro português surgiu como um empreendimento privado e depois de um período em que coexistiram empresas públicas e privadas, a integração ocorreu em 1951. Em 1975 foi criada a CP – Caminhos-de-ferro Portugueses, E.P., como empresa pública responsável pela gestão do caminho-de-ferro em Portugal. [4]
A CP – Comboios de Portugal é actualmente o operador público de transporte ferroviário em Portugal. Devido às profundas alterações introduzidas no sector ferroviário português e à luz das orientações comunitárias, nomeadamente com o surgimento da empresa pública gestora da infra-estrutura ferroviária (REFER - Rede Ferroviária Nacional, EP) a CP concentra-se hoje na exploração do transporte de passageiros e mercadorias nas linhas férreas e ramais que integram a rede ferroviária nacional, com excepção do transporte suburbano no eixo ferroviário Norte/Sul da região de Lisboa.
Ainda em crescimento, a rede ferroviária nacional carece de acções de manutenção para que continuem a ser circuláveis as redes existentes. Assim tona-se necessário conhecer os parâmetros e o modo de os estudar e analisar por forma a serem identificados os meios de actuação em caso de necessidade de intervenção. Posto isto, são apresentados neste trabalho os parâmetros de avaliação da qualidade geométrica da via-férrea bem como um método de avaliação determinístico para a análise desses parâmetros aplicado a um troço de aproximadamente 33Km situado na Linha do Norte.
A avaliação destes parâmetros foi feita com recurso a um programa, também ele desenvolvido para esta tese, em linguagem Matlab cujo código integral surge em anexo.
Caracterização da Qualidade Geométrica das Vias-Férreas
3
1.2 OBJECTIVOS
Tal como foi dito anteriormente o tema do presente trabalho, “Caracterização da qualidade geométrica de vias ferroviárias”, foi proposto no âmbito da realização da dissertação, do curso de Engenharia Civil, da Faculdade de Engenharia da Universidade do Porto.
A análise da qualidade geométrica de vias-férreas visa, além do cumprimento dos demais objectivos do trabalho apresentados seguidamente, aprofundar os conhecimentos em engenharia ferroviária aliando o presente estudo aos conhecimentos adquiridos na disciplina de “caminho de ferro” leccionada no referido curso.
O presente trabalho tem os seguintes objectivos:
• estudar o campo de aplicação da norma europeia EN 13848-5 2008 “Railway application – Track - Track geometry quality – Part 5: Geometric quality levels”[1];
• aplicar a norma europeia EN 13848-5 2008 “Railway application – Track - Track geometry quality – Part 5: Geometric quality levels”[1] a um troço de aproximadamente 31 Km, da via ascendente, da Linha do Norte;
• compreender a influência do traçado em planta da via na qualidade geométrica da via;
• avaliar a taxa de degradação do desvio padrão do alinhamento no troço de aproximadamente 31 km;
• apresentar as conclusões e desenvolvimentos futuros do trabalho.
1.3 ESTRUTURA DO TRABALHO
Este trabalho está organizado em 4 Capítulos.
No presente capítulo, Introdução, faz-se uma breve abordagem ao estado da arte onde se apresenta de forma sintética a origem e história dos caminhos-de-ferro, no contexto nacional. Apresenta-se também os objectivos e a estrutura do trabalho.
No capítulo 2, Qualidade Geométrica da Via-Férrea, começa-se por se apresentar de modo sumário o método de recolha de dados do veículo de inspecção da via. Seguidamente apresenta-se a definição dos parâmetros de avaliação da qualidade geométrica de linhas ferroviárias e faz-se uma descrição da norma europeia EN 13848-5 2008 “Railway application – Track - Track geometry quality – Part 5: Geometric quality levels”[1], mostrando-se o âmbito e campo de aplicação.
No capítulo 3, Caso de Estudo: avaliação da qualidade geométrica de um troço de via da Linha do Norte, fazem-se duas análises. A primeira tem em vista a avaliação da importância da geometria da via na manutenção dos parâmetros de qualidade geométrica. A segunda tem como objectivo avaliar a evolução temporal do alinhamento num espaço de tempo de 4 anos.
Por fim, no capítulo 4, Conclusões e Desenvolvimentos Futuros, apresentam-se as conclusões alcançadas e perspectivam-se possíveis melhorias desta tese.
Caracterização da Qualidade Geométrica das Vias-Férreas
5
2 Qualidade Geométrica da Via-férrea
2.1 SUMÁRIO
Neste capítulo aborda-se a sumariamente a temática da recolha dos parâmetros geométricos das vias-férreas.
Apresentam-se também os parâmetros geométricos de avaliação da qualidade das vias-férreas, demonstrando-se, quando oportuno, o respectivo método de cálculo.
Na parte final deste capítulo apresentam-se os limites alarme preconizados norma europeia EN 13848-5 2008 “Railway application – Track - Track geometry quality – Part 5: Geometric quality levels”[1].
2.2 INSPECÇÃO DA QUALIDADE GEOMÉTRICA DA VIA
Muitas das linhas férreas integradas na rede ferroviária nacional estão em operação há mais de cem anos, apesar de a sua construção ter sido feita com baixos recursos tecnológicos e a sua concepção baseada em critérios empíricos, especialmente no que respeita aos aspectos geotécnicos.
No âmbito da modernização dos caminhos-de-ferro, particularmente quando é necessário proceder-se à renovação de plataformas ferroviárias, reaproveitando, portanto, as vias já existentes é importante ter indicadores de desempenho que permitam avaliar a qualidade da via-férrea para que se possa obter uma comparação antes a após a remodelação.
Existem cinco parâmetros ou indicadores para a análise da qualidade geométrica das vias-férreas, nomeadamente:
• Bitola;
• Nivelamento longitudinal;
• Nivelamento transversal;
• Alinhamento;
• Empeno (ou escala).
Os dados relativos aos parâmetros geométricos de qualidade da via-férrea analisados neste trabalho foram recolhidos pelo veículo de inspecção Plasser & Theurer EM 120 (figura 1).
Caracterização da Qualidade Geométrica das Vias-Férreas
6
Figura 1-veículo de inspecção Plasser&Theurer EM120
Fonte: REFER
Com o auxílio deste veículo os parâmetros geométricos calculam-se com base nos dados sobre as sucessivas posições dos carris.
Neste caso, não são realizadas medições directas em pontos fixos nem existe qualquer contacto físico com os carris durante o processo de medição.
Este sistema torna-se muito versátil para a recolha de dados relativos a parâmetros geométricos definidos por comprimentos de onda, uma vez que permite escolher e, portanto, filtrar os comprimentos de onda desejados.
2.3 PARÂMETROS PARA A CARACTERIZAÇÃO DA QUALIDADE DA VIA -FÉRREA
2.3.1 Bitola
A bitola é “o parâmetro de maior importância na definição das características geométricas, exactamente por definir a base de uma ferrovia.”(Lima 1998)
A bitola é a distância entre dois carris medida 15 mm abaixo da cabeça do carril, medidos em esquadria com os mesmos, como se mostra na figura 2.
Caracterização da Qualidade Geométrica das Vias-Férreas
7
Figura 2-Bitola
Define-se desvio de bitola como se mostra em (1):
������������ ���������� � ������������(1)
2.3.2. NIVELAMENTO LONGITUDINAL
Na via-férrea a variação do alinhamento longitudinal dos carris é um dos principais indicadores de desempenho e define-se como o desvio na direcção z, perpendicular ao plano de rolamento, em consecutivas posições, do eixo de cada carril, em relação a uma linha de referência paralela ao plano de rolamento.
2.3.3 NIVELAMENTO TRANSVERSAL
O nivelamento transversal é definido como a diferença de altura entre dois carris, como mostra na figura 3 e na expressão 2.
Bitola Actual
Bitola Nominal
Caracterização da Qualidade Geométrica das Vias-Férreas
8
Figura 3-NivelamentoTransversal
Nivelamento transversal = Escala Actual – Escala Desejada (2)
2.3.4 ALINHAMENTO
O alinhamento é definido como a média da posição lateral de cada carril. O alinhamento é medido separadamente para os carris da esquerda e da direita. A figura 4 mostra a definição de alinhamento.
Figura 4-Alinhamento
Escala
Actual Escala
Desejada
Alinhamento do carril direito
Alinhamento do carril esquerdo
Eixo nominal da via (projecto) Actual eixo da via
Alinhamento médio
Caracterização da Qualidade Geométrica das Vias-Férreas
9
2.3.5 EMPENO
Considerando-se quatro pontos sobre a mesa de rolamento dos carris, dois sobre cada carril, formando um rectângulo, define-se como empeno, a distância vertical de um dos pontos ao plano formado pelos outros três. Dito de outra forma, o valor do empeno corresponde à diferença de dois nivelamentos transversais numa determinada base de medição que habitualmente é 3 m, como se mostra na figura 5.
Figura 5-Empeno
2.4 CARACTERIZAÇÃO DA QUALIDADE GEOMÉTRICA DA VIA-FÉRREA SEGUNDO A NORMA EN
13848-5 2008 “R AILWAY APPLICATION – TRACK - TRACK GEOMETRY QUALITY – PART 5: GEOMETRIC QUALITY LEVELS ”[1]
Tal como foi dito anteriormente neste trabalho, no que respeita à avaliação da qualidade geométrica da via-férrea, dá-se grande ênfase, à norma europeia EN 13848-5 2008 “Railway application – Track - Track geometry quality – Part 5: Geometric quality levels”[1]. Nesse sentido, apresenta-se de seguida uma breve abordagem ao conteúdo da mesma.
A norma europeia EN 13848-5 2008 “Railway application – Track - Track geometry quality – Part 5: Geometric quality levels”[1] estabelece valores máximos e mínimos admissíveis para o que considera serem os parâmetros de caracterização da qualidade geométrica da via-férrea já referidos anteriormente neste trabalho. Assim, preconiza a existência de três indicadores descritivos da qualidade da via-férrea, nomeadamente:
• Valores extremos de defeitos isolados;
• Desvio padrão num comprimento de 200m;
• Média.
Considera ainda três níveis ou limites de alarme, estudados com mais detalhe em 2.4.1, 2.4.2 e 2.4.3:
Caracterização da Qualidade Geométrica das Vias-Férreas
10
• Limite de Segurança: refere-se ao valor que, quando excedido, requer tomada de medidas imediata o que pode levar à diminuição da velocidade máxima dos comboios ou mesmo ao encerramento temporário da linha até que o defeito seja corrigido;
• Limite de Intervenção: refere-se ao valor que, quando excedido, requer medidas correctivas de manutenção para que não seja atingido o limite de segurança até à próxima inspecção;
• Limite de Alerta: refere-se ao valor que, quando excedido, requer que as condições geometria da via sejam analisadas e consideradas nas operações da manutenção previstas.
Estes valores são dados em função da velocidade de circulação da via que é um factor fundamental para a avaliação da qualidade da via.
Os limites de Segurança e Intervenção são dados apenas para defeitos isolados enquanto os limites de alerta são dados para os demais defeitos isolados e irregularidades de desvio padrão.
Ao contrário dos limites de segurança, que têm em conta parâmetros relacionados com interacção via/veículo, bem como o risco de situações inesperadas, os limites de Intervenção e Alerta estão essencialmente relacionados com a política de manutenção.
Nesta norma, os limites de nivelamento e de alinhamento são definidos para duas gamas de comprimentos de onda denominados de D1 e D2 a que correspondem, respectivamente, comprimentos de onda entre 3 mm e 25 mm e entre 25 mm e 70 mm.
Define-se comprimento de onda como a distância entre valores repetidos num padrão de onda, usualmente representado pela letra grega lambda (λ).
Numa onda sinusoidal, o comprimento de onda é a distância entre picos (ou máximos), tal como se mostra na figura 6.
No eixo x representa-se a distância e no eixo y a amplitude da onda sinusóide. Na figura, pode também ver-se o modo de medição do comprimento de onda. Tal como explicado anteriormente, pode ver-se no gráfico o comprimento de onda marcada entre dois máximos da função.
Figura 6- Função sinusoidal
Caracterização da Qualidade Geométrica das Vias-Férreas
11
A política de manutenção pode ser direccionada apenas para aspectos de segurança por si só ou pode ter em consideração parâmetros relacionados com a qualidade da viagem ou com a prestação de serviços mais atractivos como viagens a velocidades superiores aliados, inevitavelmente, à segurança. Os limites e Alerta e Intervenção estabelecidos pela entidade responsável pela infra-estrutura serão definidos de modo a assegurar, no mínimo, a segurança embora possam ser alargados para que se atinja um determinado nível de qualidade na circulação.
A média e o desvio padrão dos parâmetros que caracterizam a qualidade geométrica da via (para o limite a alerta são avaliados em trechos de 200m ou, no caso da bitola, em trechos de 100 m.
2.4.1 LIMITES DE SEGURANÇA
Os valores para o limite de Segurança apresentados na referida norma derivam da experiência e de considerações teóricas sobre a interacção roda – carril bem como de testes realizados com diferentes veículos até ao ponto de descarrilamento. Exceder estes valores limite implica a aplicação de medidas específicas para a redução do risco de descarrilamento, ou outros riscos, até um nível aceitável.
Nas tabelas de 1 a 5, que se seguem, apresentam-se os valores indicados pela norma europeia EN 13848-5 2008 “Railway application – Track - Track geometry quality – Part 5: Geometric quality levels”[1] para o Limite de Segurança, relativamente aos parâmetros geométricos bitola, nivelamento longitudinal, nivelamento transversal e alinhamento.
Tabela 1-Bitola - defeitos isolados
Velocidade (Km/h) Valor de Pico para a Bitola Nominal (mm)
Mínimo Máximo
V≤80 -13 +35
80<V≤120 -13 +35
120<V≤160 -13 +35
160<V≤220 -7 +27
220<V≤300 -7 +27
Caracterização da Qualidade Geométrica das Vias-Férreas
12
Tabela 2-Bitola
Velocidade (Km/h) Diferença entre a Bitola nominal e a bitola média em 100 m (mm)
Mínimo Máximo
V≤80 -9 +32
80<V≤120 -7 +27
120<V≤160 -5 +20
160<V≤220 -5 +20
220<V≤300 -5 +20
Nas tabelas 1 e 2 apresentam-se valores mínimos e máximos para a bitola. Como é possível verificar pelas tabelas seguintes o mesmo não acontece para os restantes parâmetros. Esta distinção deriva dos possíveis defeitos de bitola tais como o estreitamento e o alargamento. Também se verifica que os limites mínimos são, em valor absoluto, menores que os máximos, facto este que se deve à gravidade acrescida de um estreitamento em relação a um alargamento de bitola.
Tabela 3-Nivelamento Longitudinal – defeitos isolados – Diferença entre a média e o valor de pico
Velocidade (Km/h) Comprimento de Onda (mm)
D1 D2
V≤80 29 N/A
80<V≤120 26 N/A
120<V≤160 24 N/A
160<V≤220 20 33
220<V≤300 17 28
Tabela 4-Nivelamento Transversal – defeitos isolados - Diferença entre a média e o valor de pico
Velocidade (Km/h) Diferença entre a média e o valor de pico (mm)
V≤80 16
80<V≤120 14
120<V≤160 12
160<V≤220 11
220<V≤300 10
Caracterização da Qualidade Geométrica das Vias-Férreas
13
Tabela 5-Alinhamento – defeitos isolados - Diferença entre a média e o valor de pico
Velocidade (Km/h) Comprimento de Onda (mm)
D1 D2
V≤80 22 N/A
80<V≤120 17 N/A
120<V≤160 14 N/A
160<V≤220 12 24
220<V≤300 10 20
A norma contempla ainda os limites para o empeno que não são apresentados nesta tese dado este aspecto não é considerado no âmbito deste trabalho.
2.4.2 LIMITES DE INTERVENÇÃO
Nas tabelas de 6 a 10, apresentadas de seguida, constam os Limites de Intervenção preconizados pela norma para os parâmetros geométricos bitola, nivelamento longitudinal, nivelamento transversal e alinhamento.
Tabela 6-Bitola - defeitos isolados
Velocidade (Km/h) Valor de Pico para a Bitola Nominal (mm)
Mínimo Máximo
V≤80 -11 +30
80<V≤120 -11 +30
120<V≤160 -11 +30
160<V≤220 -5 +22
220<V≤300 -5 +22
Caracterização da Qualidade Geométrica das Vias-Férreas
14
Tabela 7-Bitola
Velocidade (Km/h) Diferença entre a Bitola nominal e a bitola média em 100m (mm)
Mínimo Máximo
V≤80 -8 +28
80<V≤120 -6 +25
120<V≤160 -4 +18
160<V≤220 -4 +18
220<V≤300 -4 +18
Tabela 8-Nivelamento Longitudinal - defeitos isolados – Diferença entre a média e o valor de pico
Velocidade (Km/h) Comprimento de Onda (mm)
D1 D2
V≤80 16-20 N/A
80<V≤120 12-18 N/A
120<V≤160 10-17 N/A
160<V≤220 9-14 18-23
220<V≤300 8-12 16-20
Tabela 9-Nivelamento Transversal - defeitos isolados – Diferença entre a média e o valor de pico
Velocidade (Km/h) Diferença entre a média e o valor de pico (mm)
V≤80 12-14
80<V≤120 10-12
120<V≤160 8-10
160<V≤220 7-9
220<V≤300 6-8
Caracterização da Qualidade Geométrica das Vias-Férreas
15
Tabela 10-Alinhamento – defeitos isolados - Diferença entre a média e o valor de pico
Velocidade (Km/h) Comprimento de Onda (mm)
D1 D2
V≤80 14-16 N/A
80<V≤120 10-12 N/A
120<V≤160 8-10 N/A
160<V≤220 7-9 14-17
220<V≤300 6-8 12-14
2.4.3 LIMITES DE ALERTA
Nas tabelas de 11 a 16, apresentadas de seguida, constam os Limites de Alerta preconizados pela norma para os parâmetros geométricos bitola, nivelamento longitudinal, nivelamento transversal e alinhamento.
Tabela 11-Bitola - defeitos isolados
Velocidade (Km/h) Valor de Pico para a Bitola Nominal (mm)
Mínimo Máximo
V≤80 -9 +25
80<V≤120 -9 +25
120<V≤160 -9 +25
160<V≤220 -4 +20
220<V≤300 -4 +20
Caracterização da Qualidade Geométrica das Vias-Férreas
16
Tabela 12 - Bitola
Velocidade (Km/h) Diferença entre a Bitola nominal e a bitola média em 100m (mm)
Mínimo Máximo
V≤80 -7 +25
80<V≤120 -5 +22
120<V≤160 -3 +16
160<V≤220 -3 +16
220<V≤300 -3 +16
Tabela 13-Nivelamento Longitudinal - defeitos isolados – Diferença entre a média e o valor de pico
Velocidade (Km/h) Comprimento de Onda (mm)
D1 D2
V≤80 12-18 N/A
80<V≤120 10-16 N/A
120<V≤160 8-15 N/A
160<V≤220 7-12 14-20
220<V≤300 6-10 12-18
Tabela 14-Nivelamento Transversal – Desvio Padrão - Diferença entre a média e o valor de pico
Velocidade (Km/h) Diferença entre a média e o valor de pico (mm)
V≤80 10-13
80<V≤120 8-11
120<V≤160 6-9
160<V≤220 6-8
220<V≤300 5-7
Caracterização da Qualidade Geométrica das Vias-Férreas
17
Tabela 15-Alinhamento - Defeitos Isolados - Diferença entre a média e o valor de pico
Velocidade (Km/h) Comprimento de Onda (mm)
D1 D2
V≤80 12-15 N/A
80<V≤120 8-11 N/A
120<V≤160 6-9 N/A
160<V≤220 5-8 10-15
220<V≤300 4-7 8-13
Tabela 16 – Alinhamento - Desvio Padrão
Velocidade (Km/h) Comprimento de Onda (mm)
D1
V≤80 1,5-1,8
80<V≤120 1,2-1,5
120<V≤160 1,0-1,3
160<V≤220 0,8-1,1
220<V≤300 0,7-1,0
Em síntese, pode dizer-se que a norma europeia EN 13848-5 2008 “Railway application – Track - Track geometry quality – Part 5: Geometric quality levels”[1] contempla os aspectos relevantes relacionados com a qualidade geométrica de vias ferroviárias, tendo em conta a velocidade de circulação na via, os comprimentos de onda das irregularidades, os defeitos pontuais e distribuídos e ainda os diversos parâmetros caracterizadores da qualidade da via: bitola, nivelamento longitudinal, nivelamento transversal, alinhamento e empeno.
Caracterização da Qualidade Geométrica das Vias-Férreas
19
3 Caso de Estudo: Avaliação da qualidade geométrica d e um
troço de via da Linha do Norte
3.1 SUMÁRIO
Neste capítulo apresenta-se a avaliação de qualidade geométrica de um troço de via da Linha do Norte de aproximadamente 31 km de extensão.
Nesse sentido, faz-se uma análise primordial onde foram analisados os parâmetros geométricos da qualidade da via-férrea, para dois segmentos de via com características geométricas diferentes (recta e curva), com o objectivo de avaliar a influência do traçado na qualidade geométrica da via-férrea.
Posteriormente, numa segunda análise, a partir dos registos medidos nas campanhas de inspecção entre 2005 e 2009, faz-se uma tentativa de avaliação da evolução do desvio padrão do alinhamento ao longo do tempo.
3.2 CARACTERIZAÇÃO DO TROÇO FERROVIÁRIO EM ESTUDO
Tal como foi dito no capítulo 1 este trabalho contempla a avaliação e caracterização de acordo com os parâmetros da qualidade, de um troço da via ascendente de aproximadamente 31 km, situado na Linha do Norte, cujas características geométricas do traçado se apresentam no anexo3.
O troço em estudo é composto por rectas e curvas compostas (curvas de transição, de entrada e saída, e curva circular central).
A velocidade de circulação varia entre 50 km/h e 140 km/h.
No troço em análise o carril usado é do tipo 54 E1.
3.3 AVALIAÇÃO DA QUALIDADE GEOMÉTRICA DO TROÇO EM ESTUDO
Dado que os dados recolhidos correspondem a medições de 0,25 m em 0,25 m consegue perceber-se, facilmente, que para cada campanha e para 31 km os dados apresentam 124000 registos por carril. Deste modo a estratégia usada para a primeira análise dos dados foi definir dois segmentos, um em recta (segmento A) e outro em curva (segmento B), onde a velocidade fosse a mesma (V=120 km/h) a fim de verificar a existência, ou não, de defeitos e ainda de avaliar a influência do traçado na qualidade
Caracterização da Qualidade Geométrica das Vias-Férreas
20
geométrica da via-férrea. Esta análise contempla apenas os registos respeitantes à campanha de Janeiro de 2009.
3.3.1 INFLUÊNCIA DO TRAÇADO EM PLANTA
3.3.1.1. SEGMENTO A
Na figura 7, apresenta-se a variação da bitola em função da posição no traçado para o segmento A que apresenta aproximadamente 700 m de extensão.
Após a aplicação da norma, é perceptível que ao nível da bitola e para este segmento não são ultrapassados os limites de alerta, intervenção e segurança, representados no gráfico com as cores verde, amarelo e vermelho, respectivamente.
Verificam-se ainda pequenas e expectáveis variações do parâmetro ao longo do troço.
Figura 7- Bitola
Na figura 8 mostra-se a variação do nivelamento longitudinal na gama de comprimentos de onda D1 em função da posição no traçado. Após a análise da mesma figura verifica-se que, tal como para a bitola, no parâmetro Nivelamento Longitudinal D1 os limites da norma são claramente cumpridos.
Veja-se ainda o comportamento do parâmetro para o carril esquerdo, apresentado a azul, e do carril direito, apresentado a ciano, cujas variações de nivelamento são pouco significativas tal como se vê na figura 9.
200 200.1 200.2 200.3 200.4 200.5 200.6 200.7 200.8 200.9 201-15
-10
-5
0
5
10
15
20
25
30
35
40
Posição (Km)
Bito
la (
mm
)
BitolaLim. Segurança Min.
Lim. Segurança Max.
Lim. Intervenção Min.
Lim. Intervenção Max.
Lim. Alerta Min.Lim. Alerta Min.
Caracterização da Qualidade Geométrica das Vias-Férreas
21
Figura 8-Nivelamento Longitudinal D1 Esquerdo e Direito
Figura 9-Detalhe do Nivelamento Longitudinal D1 Esquerdo e Direito
Relativamente ao Nivelamento longitudinal na gama de comprimentos de onda D2 representado na figura 10, em função da posição, verifica-se uma evolução semelhante nos carris esquerdo e direito.
Note-se que o nivelamento Longitudinal na gama D2 só apresenta limites de segurança, intervenção e alerta para velocidades superiores a 160 km/h não se aplicando ao presente estudo porque a norma não estabelece limites, para este parâmetro, para a velocidade em estudo.
200 200.2 200.4 200.6 200.8 201 201.2-30
-20
-10
0
10
20
30
Posição (Km)
Niv
elam
ento
Lon
gitu
dina
l D1
(mm
)
Nivelamento Longitudinal Esquerdo D1
Nivelamento Longitudinal Direito D1
Limite de SegurançaLimite de Intervenção
Limite de Alerta
200.02 200.04 200.06 200.08 200.1 200.12 200.14
-6
-4
-2
0
2
4
Posiçao (Km)
Niv
elam
ento
Lon
gitu
dina
l D1
(mm
)
Nivelamento Longitudinal Esquerdo D1
Nivelamento Longitudinal Direito D1
Caracterização da Qualidade Geométrica das Vias-Férreas
22
Figura 10-Nivelamento Longitudinal D2 Esquerdo e Direito
Na figura 11 mostra-se a variação do nivelamento transversal em função da posição. Tal como anteriormente os limites impostos pela norma são cumpridos ainda que em determinados pontos do segmento estejam muito próximos do limite de Alerta.
Figura 11-Nivelamento Transversal
Na figura 12 apresenta-se a variação do alinhamento na gama D1 em função da posição. Analisando a respectiva figura observa-se que ao nível do alinhamento os limites da norma são também francamente cumpridos.
200 200.1 200.2 200.3 200.4 200.5 200.6 200.7 200.8-10
-8
-6
-4
-2
0
2
4
6
8
10
Posiçao (Km)
Niv
elam
ento
Lon
gitu
dina
l D2
(mm
)
Niv. Long. Esq. D2
Niv. Long. Dir. D2
200 200.1 200.2 200.3 200.4 200.5 200.6 200.7 200.8 200.9 201-15
-10
-5
0
5
10
15
Posição (Km)
Niv
elam
ento
Tra
nsve
rsal
(m
m)
Nivelamento Transversal
Limite de SegurançaLimite de Intervenção
Limite de Alerta
Caracterização da Qualidade Geométrica das Vias-Férreas
23
Figura 12-Alinhamento D1 Esquerdo e Direito
Na figura 13 apresenta-se um detalhe da figura 12 do alinhamento para a gama D1 de ambos os carris onde se pode verificar que não a evolução não é semelhante. Esta variação pode dever-se ao movimento de lacete que também pode ocorrer em alinhamento recto.
O movimento de lacete traduz-se numa série de oscilações do eixo da via que se tornam possíveis porque a distância entre os bordos interiores das duas rodas de um mesmo eixo é inferior à bitola da via.
Figura 13-Detalhe do Alinhamento D1 Esquerdo e Direito
200 200.1 200.2 200.3 200.4 200.5 200.6 200.7 200.8 200.9 201-20
-15
-10
-5
0
5
10
15
20
Posição (Km)
Alin
ham
ento
D1
(mm
)
Alinhamento Esquerdo D1
Alinhamento Direito D1
Limite de SegurançaLimite de Intervenção
Limite de Alerta
200.44 200.45 200.46 200.47 200.48 200.49 200.5 200.51 200.52 200.53 200.54-3
-2
-1
0
1
2
3
Posiçao (Km)
Alin
ham
ento
D1
(mm
)
Alinhamento Esquerdo D1Alinhamento Direito D1
Caracterização da Qualidade Geométrica das Vias-Férreas
24
Relativamente ao Alinhamento na gama D2 os ficheiros de dados estudados não apresentam valores para este segmento em recta pelo que não será apresentado o estudo análogo aos anteriores para este parâmetro. No entanto importa referir que tal como para o Nivelamento Longitudinal da gama D2 e para a velocidade de 120 km/h, não são impostos, pela norma, quaisquer valores limite de segurança, alerta ou intervenção.
3.3.1.2. SEGMENTO B
Faz-se agora uma análise idêntica à apresentada para o segmento A, mas agora para um segmento de via em curva, denominado segmento B.
Na figura 14 apresenta-se a variação da bitola em função da posição no traçado. Observando a figura é perceptível que ao nível da bitola e para este segmento não existem defeitos ao nível da segurança, alerta ou intervenção.
Verificam-se ainda pequenas e expectáveis variações do parâmetro ao longo do troço.
Figura 14-Bitola
Na figura 15 apresenta-se a variação do nivelamento longitudinal na gama D1 (comprimento de onda de 3 a 25 mm) em função da posição na via.
Pela análise da figura verifica-se que, tal como para a bitola, o parâmetro nivelamento longitudinal cumpre os limites indicados na norma.
Constata-se também que existem para o carril esquerdo, apresentado a azul, e para o direito, apresentado a ciano, variações de nivelamento que se justificam com questões inerentes à geometria do traçado em curva. Importa referir que as maiores diferenças entre o carril esquerdo e direito verificam-se nas curvas de transição de entrada e curva circular central.
201 201.1 201.2 201.3 201.4 201.5 201.6-15
-10
-5
0
5
10
15
20
25
30
35
40
Posição (Km)
Bito
la (
mm
)
BitolaLimite de Segurança Min
Limite de Segurança Max
Limite de Intervenção Min
Limite de Intervenção Max
Limite de Alerta MinLimite de Alerta Min
Caracterização da Qualidade Geométrica das Vias-Férreas
25
Figura 15-Nivelamento Longitudinal D1 Esquerdo e Direito
Relativamente ao Nivelamento longitudinal na gama D2 representado na figura 16, em função da posição, verifica-se um comportamento semelhante entre os carris esquerdo e direito.
Figura 16-Nivelamento Longitudinal D2 Esquerdo e Direito
Na figura 17 mostra-se a variação do nivelamento transversal em função da posição. Repare-se, para este parâmetro, na forma do gráfico que se explica-se com a própria geometria da curva estudada que sendo uma curva composta apresenta clotóide de entrada, representada na parte inicial da figura, uma curva circular central, evidenciada na zona central e uma clotóide de saída, representada na parte final da figura.
201 201.1 201.2 201.3 201.4 201.5 201.6 201.7-30
-20
-10
0
10
20
30
Posição (Km)
Niv
elam
ento
Lon
gitu
dina
l D1
(mm
)
Nivelamento Longitudinal Esquerdo D1
Nivelamento Longitudinal Direito D1
Limite de SegurançaLimite de Intervenção
Limite de Alerta
200.95 201 201.05 201.1 201.15 201.2 201.25 201.3 201.35 201.4-6
-4
-2
0
2
4
6
Posiçao (Km)
Niv
elam
ento
Lon
gitu
dina
l D2
(mm
)
Niv. Long. Esq. D2Niv. Long. Dir. D2
Caracterização da Qualidade Geométrica das Vias-Férreas
26
Figura 17-Nivelamento Transversal
Figura 18 - Nivelamento transversal e irregularidades de nivelamento transversal
Na figura 18 mostra-se, para além das irregularidades no nivelamento transversal, o valor deste nivelamento ao longo do troço, devido à geometria do troço. Para poder aplicar-se a norma, deve subtrair-se o nivelamento transversal do troço ao perfil de irregularidades (indicado na figura 17), obtendo-se a figura 19.
201 201.1 201.2 201.3 201.4 201.5-100
-80
-60
-40
-20
0
20
Posição (Km)
Niv
elam
ento
Tra
nsve
rsal
(m
m)
Nivelamento Transversal
Limite de SegurançaLimite de Intervenção
Limite de Alerta
y = -766,2x + 153972
R² = 0,9581
y = 777,82x - 156670
R² = 0,9492-100
-80
-60
-40
-20
0
20
200,95 201,15 201,35 201,55
Niv
ela
me
nto
tra
nsv
ers
al
(mm
)
Posição (Km)
Nivelamento Transversal
Linha média
clotoideE
clotoideS
Linear (clotoideE)
Linear (clotoideS)
Caracterização da Qualidade Geométrica das Vias-Férreas
27
Figura 19- Nivelamento transversal, comparação com os limites da norma
Na figura 20 observa-se a variação do alinhamento na gama D1 em função da posição. Analisando a respectiva figura observa-se que os limites da norma são francamente cumpridos. Notam-se ainda algumas diferenças na variação deste parâmetro geométrico nos dois carris da via. Como pode ver-se pela Figura, o carril esquerdo apresenta maiores irregularidades o que é um facto expectável uma vez que para uma linha de tráfego misto, como é o caso, os carris sofrem desgaste relativo ao excesso e insuficiência de escala o que provoca, dependendo do tipo de tráfego mais frequente, maiores irregularidades num dos carris.
Caracterização da Qualidade Geométrica das Vias-Férreas
28
Figura 20-Alinhamento D1 Esquerdo e Direito
Relativamente ao Alinhamento na gama D2 os ficheiros estudados não apresentam valores para este segmento em curva pelo que não será apresentado o estudo análogo aos anteriores para este parâmetro.
3.3.1.3 COMPARAÇÃO ENTRE OS SEGMENTOS A E B
Pela análise das figuras 7 e 14 verifica-se que a variação da bitola, em curva (segmento B) e para a mesma velocidade, é ligeiramente maior do que no segmento A, dado existirem picos mais acentuados naquele segmento.
Pela análise da figura 8, relativa ao segmento em recta, e da figura 15, segmento em curva, é clara a diferença no parâmetro alinhamento longitudinal D1 sendo que em curva e apara a mesma velocidade a via apresenta variações significativamente menores ao nível deste parâmetro.
O nivelamento transversal é o parâmetro onde a diferença entre o troço em recta e o troço em curva é mais significativa. Com é evidenciado pelas figuras 11 e 17 o troço em recta apresenta uma forma aproximadamente linear enquanto o troço em curva apresenta uma forma parabólica que traduz a variação da curvatura ao longo do segmento: curva de transição – curva circular - curva de transição.
3.4 EVOLUÇÃO DO ALINHAMENTO AO LONGO DO TEMPO
Nesta fase do trabalho, tal como foi dito anteriormente, foi estudado um troço de aproximadamente 31Km na qual circulam comboios desde os 50 km/h aos 140 km/h.
Por uma questão de tempo limitado, nesta análise, optou-se por analisar apenas a evolução do alinhamento do carril esquerdo.
Para isso foram investigadas 9 campanhas fim de se conseguir detectar uma tendência no comportamento da via:
• Março de 2005;
201 201.1 201.2 201.3 201.4 201.5 201.6-20
-15
-10
-5
0
5
10
15
20
Posição (Km)
Alin
ham
ento
D1
(mm
)
Alinhamento Esquerdo D1
Alinhamento Direito D1
Limite de SegurançaLimite de Intervenção
Limite de Alerta
Caracterização da Qualidade Geométrica das Vias-Férreas
29
• Dezembro de 2006;
• Março de 2007;
• Junho de 2007;
• Setembro de 2007;
• Dezembro de 2007;
• Março de 2008;
• Junho de 2008;
• Janeiro de 2009.
Nas figuras 19 e 20 apresenta-se a variação do alinhamento ao longo do troço para as nove campanhas, para o carril direito e esquerdo, respectivamente.
Figura 21-Variação do Alinhamento do carril direito no tempo
2005 2005.5 2006 2006.5 2007 2007.5 2008 2008.5 2009 2009.5
200
210
220
230
-50
0
50
Campanha
Direita
Posiçao (Km)
Alin
ham
ento
(m
m)
Caracterização da Qualidade Geométrica das Vias-Férreas
30
Figura 22-Variação do Alinhamento do carril esquerdo no tempo
Da análise dos gráficos observa-se que aposição dos picos do alinhamento se mantém ao longo dos anos nas campanhas consideradas.
3.4.1 METODOLOGIA PARA A ANÁLISE DA EVOLUÇÃO DO DESVIO PADRÃO DO ALINHAMENTO NO TEMPO
Como já foi referido anteriormente, cada campanha apresenta aproximadamente 124000 registos. Dividindo cada campanha em trechos de 200 m obtém-se 159 trechos. Deste modo, tornou-se importante o uso de uma ferramenta versátil e com potencialidades de programação – O Matlab. Toda a análise gráfica e numérica foi feita com recurso a esta ferramenta aliada, em alguns casos, ao Excel.
A metodologia para a análise da evolução do desvio padrão do alinhamento no tempo consistiu em várias etapas:
1. ajuste dos dados das várias campanhas;
2. divisão do techo de 31 Km, para cada campanha, em trechos de 200 m;
3. cálculo do desvio padrão do alinhamento em cada trecho de 200 m e em cada campanha de inspecção.
No que diz respeito à primeira tarefa, a campanha de referência usada para analisar a evolução do alinhamento no tempo foi a de Janeiro de 2009, ou seja, cada uma das campanhas acima indicadas foi comparada graficamente com a campanha de 2009 e posteriormente ajustada, visualmente, a fim de se determinar para cada uma o desvio padrão em trechos de 200 m.
O ajuste temporal fez-se visualmente, isto é, após a elaboração do gráfico encontraram-se os picos de cada uma das curvas que, visualmente, se conseguiam identificar como sendo o mesmo. Identificados os picos, mediu-se a distância entre eles e graficamente procedeu-se ao ajuste retirando ou somando o valor da diferença à curva da campanha de referência (Janeiro de 2009).
20052005.5
20062006.5
20072007.5
20082008.5
20092009.5
200
205
210
215
220
225
230
235
-50
0
50
Campanha
Esquerda
Posiçao (Km)
Alin
ham
ento
(m
m)
Caracterização da Qualidade Geométrica das Vias-Férreas
31
A figura 21 apresenta a variação do alinhamento do carril esquerdo, num troço de meio quilómetro, em função da posição no traçado.
Figura 23- 2007 Dez Vs 2009 Jan, sem ajuste temporal
Na figura 22 pode ver-se o mesmo trecho acima descrito após o ajuste temporal.
Figura 24 - 2007 Dez Vs 2009 Jan, com ajuste temporal
Observando a figura 20 constata-se o claro ajuste dos picos e por conseguinte, o ajuste das campanhas.
200 200.05 200.1 200.15 200.2 200.25 200.3 200.35 200.4 200.45 200.5-4
-3
-2
-1
0
1
2
3
4
Posição (Km)
Alin
ham
ento
(m
m)
2007 Dez VS 2009, 0,5km, Esquerda
Alin. D1 2007 Dez 0,5km
Alin. D1 2009 0,5km
200 200.05 200.1 200.15 200.2 200.25 200.3 200.35 200.4 200.45 200.5-4
-3
-2
-1
0
1
2
3
4
Posição (Km)
Alin
ham
ento
(m
m)
2007 Dez VS 2009, 0,5 km, com ajuste, carril esquerdo
Alin. D1 E1 2007 Dez 0,5 km ajuste
Alin. D1 2009 0,5 km ajuste
Caracterização da Qualidade Geométrica das Vias-Férreas
32
O próximo passo na análise da evolução do desvio padrão do alinhamento ao longo do tempo é a divisão do techo de 31 Km, para cada campanha, em trechos de 200m recomendado pela norma europeia EN 13848-5 2008 “Railway application – Track - Track geometry quality – Part 5: Geometric quality levels”[1].
De seguida apresenta-se o código de programação usado para a divisão da campanha de Janeiro de 2009, ajustadas em trechos de 200 m e respectivo cálculo do desvio padrão.
No Anexo A3 consta o restante código para as demais campanhas.
• Definição dos trechos de via em 200 m: carril esquerdo e carril direito
• Cálculo do desvio padrão do alinhamento para o carril esquerdo
- Definição dos trechos de via em 200 m: carril esquerdo e carril direito
data_09=dlmread('campanha2009pronto.txt'); % ler ficheiro de dados
km_09=data_09(:,1); % lê a coluna dos kilómetros no ficheiro de 2009
loc_09=data_09(:,2)*10^-3; % lê a coluna da localização no ficheiro de 2009
x_09=[loc_09(1: 132329)+km_09(1: 132329)] %define a posição
ali_09_E=data_09(:,4); % alinhamento carril esquerdo D1
ali_09_D=data_09(:,5); %alinhamento carril direito D1
n_09=fix(size(x_09)/800); % define trechos de 200 m (800 registos)
Caracterização da Qualidade Geométrica das Vias-Férreas
33
3.4.2 ANÁLISE E DISCUSSÃO DE RESULTADOS
Depois de feita a análise descrita anteriormente em Matlab, obteve-se o desvio padrão do alinhamento em 159 trechos, tanto para o carril E como para o D.
% cálculo do desvio padrão do alinhamento no carril esquerdo
for i=1:n_09
n1_09(i)=1+800*(i-1);
n2_09(i)=800*i;
dp9(i)=[std([ali_09_E([n1_09(i):n2_09(i)],:)])];
dp_09=dp9';
end
% cálculo do desvio padrão do alinhamento no carril direito
for i=1:n_09
n1_09(i)=1+800*(i-1);
n2_09(i)=800*i;
dp9r(i)=[std([ali_09_D([n1_09(i):n2_09(i)],:)])];
dp_09r=dp9r';
end
Caracterização da Qualidade Geométrica das Vias-Férreas
34
No Anexo A1 apresentam-se os valores obtidos para o desvio padrão do alinhamento em cada um dos 159 trechos e para as 9 campanhas de inspecção analisadas. As campanhas encontram-se por ordem cronológica de 2005 até 2009.
No anexo A2 mostra-se, graficamente, a evolução desvio padrão do alinhamento ao longo do tempo, em dias, para os 159 trechos.
Nas Figuras 23, 24,25 e 26 mostra-se para quatro trechos de via a evolução do desvio padrão do alinhamento no carril esquerdo. Nos trechos 25 e 26, verifica-se uma tendência de evolução ao contrário do que acontece nos trechos 23 e 24.
Figura 25 - variação do desvio padrão no tempo para o trecho 57
Figura 26-variação do desvio padrão no tempo para o trecho 60
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 57
trecho 57
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 60
trecho 60
Caracterização da Qualidade Geométrica das Vias-Férreas
35
Figura 27 - variação do desvio padrão no tempo para o trecho 121
Figura 28 - variação do desvio padrão no tempo para o trecho 126
No âmbito deste trabalho, foi feita ainda uma tentativa de avaliar a taxa de degradação do desvio padrão do alinhamento, a partir de regressões lineares. No entanto, dado o número limitado de campanhas investigadas, só foi possível definir essa taxa (declive da recta de regressão linear) para 16 trechos como se mostra nas figuras seguintes.
y = 0,0028x - 0,5929
R² = 0,9727
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo(dias )
trecho 121
série 121
série 121'
Linear (série
121')
y = 0,0017x - 0,1856
R² = 0,9946
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 129
série 129
série 129'
Linear (série
129')
Caracterização da Qualidade Geométrica das Vias-Férreas
36
Figura 29-variação do desvio padrão no tempo para o trecho 5
Figura 30-variação do desvio padrão no tempo para o trecho 6
Figura 31-variação do desvio padrão no tempo para o trecho 7
y = 0,0004x + 0,6616
R² = 0,681
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 5
série5
série 5'
Linear (série 5')
y = 0,001x + 0,5856
R² = 0,7434
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 6
série6
série 6'
Linear (série 6')
y = 0,0003x + 1,0731
R² = 0,7577
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
De
svio
pa
drã
o
tempo (dias)
trecho 7
série 7
série 7'
Linear (série 7')
Caracterização da Qualidade Geométrica das Vias-Férreas
37
Figura 32-variação do desvio padrão no tempo para o trecho 12
Figura 33-variação do desvio padrão no tempo para o trecho 23
y = 0,001x + 0,05
R² = 0,6987
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 12
série 12
série 12'
Linear (série 12')
y = 0,0005x + 0,6094
R² = 0,647
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 23
série 23
série 23'
Linear (série 23')
Caracterização da Qualidade Geométrica das Vias-Férreas
38
Figura 34-variação do desvio padrão no tempo para o trecho 32
Figura 35-variação do desvio padrão no tempo para o trecho 24
y = 0,0006x + 0,1918
R² = 0,7886
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 32
série 32
série 32'
Linear (série 32')
y = 0,0006x + 0,209
R² = 0,7651
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 34
série 34
série 34'
Linear (série 34')
Caracterização da Qualidade Geométrica das Vias-Férreas
39
Figura 36-variação do desvio padrão no tempo para o trecho 39
Figura 37-variação do desvio padrão no tempo para o trecho 41
y = 0,0004x + 0,2902
R² = 0,6174
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 39
série 39
série 39'
y = 0,0003x + 0,5001
R² = 0,8374
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 41
série 41
série 41'
Linear (série 41')
Caracterização da Qualidade Geométrica das Vias-Férreas
40
Figura 38variação do desvio padrão no tempo para o trecho 68
Figura 39-variação do desvio padrão no tempo para o trecho 69
y = 0,0003x + 0,5922
R² = 0,8934
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 68
série 68
série 68'
Linear (série 68')
y = 0,0003x + 0,6947
R² = 0,9107
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 69
trecho 69
série 69'
Linear (série 69')
Caracterização da Qualidade Geométrica das Vias-Férreas
41
Figura 40-variação do desvio padrão no tempo para o trecho 92
Figura 41-variação do desvio padrão no tempo para o trecho 107
y = 0,0009x + 0,7288
R² = 0,918
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
Trecho 92
série 92
série 92'
Linear (série 92')
y = 0,0009x + 0,1982
R² = 0,956
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
de
pa
drã
o
campanhas (dias)
Trecho 107
série 107
série 107'
Linear (série 107')
Caracterização da Qualidade Geométrica das Vias-Férreas
42
Figura 42-variação do desvio padrão no tempo para o trecho 121
Figura 43-variação do desvio padrão no tempo para o trecho 129
y = 0,0028x - 0,5929
R² = 0,9727
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo(dias )
trecho 121
série 121
série 121'
Linear (série 121')
y = 0,0017x - 0,1856
R² = 0,9946
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 129
série 129
série 129'
Linear (série 129')
Caracterização da Qualidade Geométrica das Vias-Férreas
43
Figura 44-variação do desvio padrão no tempo para o trecho 140
No cálculo da taxa de degradação do desvio padrão do alinhamento, foram considerados os valores (indicados a vermelho), que não incluem acções de manutenção.
Pela análise das Figuras, nota-se que R2 varia entre 0,617e 0,994, o que parece ser insatisfatório. No entanto, prevê-se que a análise de mais campanhas de inspecção possa fazer melhorar a qualidade do ajuste.
Para efeitos Para efeito de cálculo da taxa de degradação só foram tidas em conta regressões com classificação “significativo” de acordo com a seguinte convenção:
• 0 ≤ R2 ≤ 0.50 , Insignificante
• 0.50 < R2 ≤ 0.85, Pouco significativo
• 0.85 < R2 ≤ 1.0 , Significativo
Feita esta avaliação para todos os gráficos obtidos, que são apresentados em anexo, observa-se e determina-se:
Tabela 17- Classificação de resultados
y = 0,0016x - 0,245
R² = 0,9336
0,0000
0,5000
1,0000
1,5000
2,0000
2,5000
3,0000
3,5000
0,0 500,0 1000,0 1500,0
de
svio
pa
drã
o
tempo (dias)
trecho 140
série 140
série 140'
Linear (série 140')
Caracterização da Qualidade Geométrica das Vias-Férreas
44
Para os únicos três gráficos classificados com “significativo” apresenta-se seguidamente a respectiva taxa de degradação.
Tabela 18 – Taxa de degradação
Caracterização da Qualidade Geométrica das Vias-Férreas
45
4 Conclusões e Desenvolvimentos Futuros
De um modo geral com este trabalho aprofundaram-se os conhecimentos sobre a qualidade geométrica da via.
A análise de evolução temporal dos parâmetros geométricos permite identificar tanto as secções de fraco desempenho - defeitos pontuais ou de pico - como a variação das condições de uma determinada secção.
A perda de regularidade geométrica depende do comportamento da super-estrutura e infra-estrutura e geralmente ocorre como resultado da degradação da via ferroviária imputada à passagem sucessiva de comboios e às variações sazonais das condições meteorológicas.
Como desenvolvimento futuro proponho:
• O estudo de mais campanhas;
• Análise de todos os parâmetros de qualidade geométrica;
• Uso de métodos de análise probabilísticos em vez dos comuns métodos determinísticos.
Com as sugestões enunciadas anteriormente pode fazer-se uma melhor análise da qualidade geométrica das vias-férreas e possivelmente poderia determinar-se um valor para a taxa de degradação das vias ferroviárias com determinadas características.
Caracterização da Qualidade Geométrica das Vias-Férreas
47
BIBLIOGRAFIA
[1]CEN 2008 “EN 13848-5 2008 “Railway application – Track - Track geometry quality – Part 5: Geometric quality levels”[1] European Committee for Standardization.
[2]http://www.cp.pt/cp/displayPage.do?vgnextoid=1d94079ffa057010VgnVCM1000007b01a8c0RCRD. Abril de 2010.
[3]http://www.cp.pt/cp/displayPage.do?vgnextoid=2fdfcef780e5c010VgnVCM1000007b01a8c0RCRD#2. Abril de 2010.
[4]http://www.cp.pt/cp/displayPage.do?vgnextoid=d5433cddefcb7010VgnVCM1000007b01a8c0RCRD#1. Abril de 2010.
[5] http://pt.wikipedia.org/wiki/Cronologia_do_Caminho-de-ferro_em_Portugal. Abril de 2010.
4
1 2 3 4 5 6 7 8 9
0,9658 1,1396 0,9434 1,1324 1,0106 1,1145 1,7194 0,7615 0,7702
0,9658 1,1396 0,9434 1,1324 1,0106 1,1145 1,7194 0,7615 0,7702
0,1104 0,9106 0,7499 1,0107 0,9598 1,3829 1,2995 0,8921 1,1934
0,1080 0,7865 0,6894 0,9998 0,9750 1,3795 1,2816 0,8256 1,0619
0,1144 1,0254 0,7762 1,0257 1,0731 1,5111 1,3692 0,8866 1,2110
0,0985 0,8359 0,7041 0,8870 1,0629 1,6324 1,3510 0,9306 1,2114
0,1000 0,9227 0,8173 1,0059 1,0146 1,3460 1,3980 0,9036 1,1500
0,1248 1,4990 1,0557 1,3269 1,2574 1,8489 1,9734 0,7797 1,2733
0,1064 1,1817 0,7736 1,0501 1,2013 1,9186 1,5870 0,9551 1,2328
10 11 12 13 14 15 16 17 18
0,7677 0,7075 1,8161 1,2224 1,3268 1,2457 0,8566 0,8863 1,3677
0,7677 0,7075 1,8161 1,2224 1,3268 1,2457 0,8566 0,8863 1,3677
0,8308 0,9702 0,9390 1,2550 1,6547 1,3202 0,9701 0,9861 1,1245
0,6435 0,6897 0,7964 1,1055 1,2571 1,1736 0,7314 0,9254 1,0375
0,9774 1,4802 1,1376 1,3691 1,5559 1,5358 1,1747 0,9659 1,1153
0,6864 1,5652 0,9975 1,4322 1,6068 1,1164 0,8855 0,8650 1,0449
0,9730 0,7484 0,9159 1,1086 1,2909 1,1980 1,2052 0,9985 1,1080
1,9972 1,5761 1,3288 1,9935 2,1681 1,2646 1,1426 1,2304 1,4779
1,3732 2,4019 1,6100 1,6299 1,6654 1,8468 1,1655 1,0325 1,1267
19 20 21 22 23 24 25 26 27
1,2105 1,2550 1,2108 1,2871 0,8694 1,3649 1,4201 1,3818 1,0122
1,2105 1,2550 1,2108 1,2871 0,8694 1,3649 1,4201 1,3818 1,0122
1,1050 1,0413 1,4131 1,5541 1,0566 1,2346 1,4005 0,6954 0,6264
1,0050 0,9527 1,2596 1,3800 0,9140 0,9359 1,1790 0,6627 0,6259
1,3286 1,1653 1,4563 1,5679 1,1536 1,1313 1,4003 0,8181 0,6334
0,9882 0,9850 1,3863 1,4438 0,9672 1,2175 1,2238 0,6962 0,6130
1,2542 1,2992 1,2419 1,6703 1,1566 1,1298 1,4580 0,8298 0,7327
1,2763 1,4310 1,3843 1,3070 1,1242 1,1056 1,9112 1,1141 1,0113
1,4848 1,2967 1,2751 1,7119 1,2925 1,0116 1,4383 0,9045 0,7417
5
28 29 30 31 32 33 34 35 36
1,4075 1,4252 1,5783 0,9277 0,8889 0,9183 0,9073 0,9313 1,0948
1,4075 1,4252 1,5783 0,9277 0,8889 0,9183 0,9073 0,9313 1,0948
0,9106 1,0063 1,4506 0,7574 0,6681 0,7397 0,7062 0,9007 0,8981
0,9448 0,9976 1,2805 0,5777 0,6032 0,7121 0,5890 0,8387 0,8334
0,8866 0,9416 1,5320 0,7277 0,6403 0,8051 0,8356 0,9069 0,8965
0,9424 1,0561 1,3561 0,7524 0,7341 0,7068 0,7113 0,9110 0,7447
0,9709 0,9929 1,4733 0,8311 0,8292 0,9071 0,8150 1,0345 0,9539
1,1386 1,2976 1,3820 1,0793 0,9705 1,0878 0,9196 1,1699 1,0334
0,9316 1,0023 1,5418 0,8637 0,9322 0,9136 1,0482 1,0687 0,9189
37 38 39 40 41 42 43 44 45
1,0843 1,1154 0,9960 1,1099 1,1661 1,1473 1,2598 1,4992 1,5123
1,0843 1,1154 0,9960 1,1099 1,1661 1,1473 1,2598 1,4992 1,5123
0,8896 0,6827 0,6322 0,6652 0,7299 0,8506 1,3281 1,7564 1,8587
0,6449 0,6164 0,5938 0,6104 0,7820 0,8866 1,4077 1,7824 1,8417
0,9642 0,7444 0,6582 0,7050 0,7568 0,8700 1,4057 1,6884 1,9592
0,7562 0,5981 0,6243 0,5680 0,7834 0,7147 1,3391 1,7143 1,7917
0,9624 0,7626 0,7638 0,7667 0,8071 0,9285 1,4448 1,6414 1,8527
1,1864 0,8993 0,9227 0,9601 0,9171 1,0947 1,3585 1,7307 1,7595
1,1074 0,8445 0,8105 0,8303 0,9254 0,9853 1,4673 1,6783 1,8503
46 47 48 49 50 51 52 53 54
1,7294 1,1989 1,0518 1,5608 1,1583 0,8578 0,8845 0,7649 0,5905
1,7294 1,1989 1,0518 1,5608 1,1583 0,8578 0,8845 0,7649 0,5905
2,0825 1,1772 1,0055 1,4392 1,0765 0,9617 0,9275 0,6635 0,6880
1,8397 1,0758 0,8621 1,4763 1,0209 0,9334 0,8216 0,5530 0,6909
1,9878 1,1748 0,9852 1,5008 1,1175 0,9874 0,8998 0,6120 0,5390
2,1004 1,0607 0,8493 1,2822 0,9074 0,8322 0,8077 0,5714 0,7886
1,9380 1,1224 0,9877 1,4432 1,1423 1,0316 1,0100 0,7515 0,6391
1,8276 1,1396 1,0224 1,5246 1,1078 0,9291 0,9736 0,6394 0,7955
2,0906 1,2070 1,0207 1,4837 1,1828 1,0875 1,0879 0,6900 0,7463
6
55 56 57 58 59 60 61 62 63
0,6481 0,8934 1,4298 0,5982 0,8466 1,2398 0,8609 0,8789 1,3764
0,6481 0,8934 1,4298 0,5982 0,8466 1,2398 0,8609 0,8789 1,3764
0,8197 1,0847 0,8495 0,7451 0,6807 0,6548 0,7065 0,9749 2,0049
0,7730 1,0105 1,3559 0,6665 0,5616 1,0338 0,6949 1,1568 1,6900
0,6697 1,0407 0,8930 0,7302 1,0162 0,8904 0,7208 0,9504 1,8494
0,6754 1,0676 0,7857 0,6646 1,0225 0,9069 0,6637 1,0238 1,6946
0,6701 1,0233 0,8569 0,6230 0,9122 0,8574 0,7515 0,8014 1,3900
0,8744 0,9843 1,5735 0,6663 0,9793 1,1446 1,1007 1,2904 1,8302
0,8127 1,0685 0,8583 0,7020 0,9887 0,9613 0,7851 0,8278 1,4440
64 65 66 67 68 69 70 71 72
1,2695 1,5715 2,0800 0,8184 0,8233 0,7360 0,5248 0,4940 0,6046
1,2695 1,5715 2,0800 0,8184 0,8233 0,7360 0,5248 0,4940 0,6046
1,3604 1,7089 1,8198 1,6632 1,0206 0,9977 0,6791 0,6134 0,9538
1,4164 1,7570 1,6712 1,1347 0,8636 0,9460 0,6130 0,5967 0,8993
1,3018 1,7081 1,7545 1,1461 0,9056 0,9664 1,0617 1,1164 1,1270
1,4199 2,3220 1,6934 0,9612 0,9117 0,9843 1,0625 1,1493 0,9666
1,3249 1,5691 2,0619 1,2501 1,0017 1,0580 1,0312 1,1163 1,0508
1,6003 1,9978 1,6201 1,2147 0,9631 1,0400 0,6552 0,5796 0,9631
1,4369 1,5824 2,0589 1,2384 1,0653 1,1013 1,1123 1,2223 1,1318
73 74 75 76 77 78 79 80 81
0,6975 0,6330 0,4680 0,5631 0,7652 0,8465 0,7605 0,9816 1,0563
0,6975 0,6330 0,4680 0,5631 0,7652 0,8465 0,7605 0,9816 1,0563
0,8063 0,7467 0,5889 0,7174 0,6279 0,8698 0,7023 1,0275 1,0146
0,7523 0,7048 0,3896 0,6021 0,5176 0,7274 0,6926 1,0047 1,1117
0,9244 0,8351 0,5417 0,6701 0,6547 0,8411 0,6842 1,0345 1,1902
0,7794 0,7107 0,4612 0,6882 0,6224 0,7591 0,6006 0,9868 1,1310
0,9607 0,8094 0,5541 0,7647 0,7004 0,8534 0,7229 0,9302 1,1532
0,8227 0,8135 0,5473 0,6838 0,8016 0,9694 0,8452 1,0542 1,0765
0,9751 0,8263 0,5549 0,7832 0,7011 0,8774 0,7085 0,9723 1,1677
7
82 83 84 85 86 87 88 89 90
0,8765 0,7032 1,9970 2,0674 2,2154 1,0271 1,6258 3,3036 3,9019
0,8765 0,7032 1,9970 2,0674 2,2154 1,0271 1,6258 3,3036 3,9019
0,9790 0,7493 1,2667 2,9114 3,3391 1,4893 1,9027 3,3428 3,4492
0,7860 0,6814 1,9178 2,5821 2,7944 1,5965 1,7982 3,1985 3,5744
1,0318 0,7385 1,9231 2,7540 2,6882 1,3542 1,9493 3,0557 4,0398
0,8614 0,7511 2,0094 1,7119 1,1571 1,0370 1,7898 4,4970 3,5220
1,0074 0,7017 1,3976 2,7006 3,3985 1,4713 1,7555 3,3110 3,5580
0,8335 0,8602 2,1216 2,5528 2,7345 1,2871 2,0975 3,1452 3,3289
1,0051 0,7008 1,3607 3,2191 2,7750 1,6054 2,1734 3,2400 3,5579
91 92 93 94 95 96 97 98 99
0,8951 1,1900 1,3776 1,3368 0,7933 0,8142 0,9733 0,9162 1,1216
0,8951 1,1900 1,3776 1,3368 0,7933 0,8142 0,9733 0,9162 1,1216
0,8907 1,5257 1,8652 1,1819 0,9118 0,7993 0,9050 1,4910 0,9188
0,8526 1,3671 1,7031 0,8075 0,6741 1,3731 1,3788 1,6289 0,9144
0,8717 1,6266 1,8564 1,2907 0,9003 0,7452 0,8577 0,8807 1,3303
0,9410 1,6314 1,6517 1,0103 0,8480 0,8364 1,0439 1,3714 0,8842
0,9474 1,6882 1,7287 1,0677 0,8310 0,9322 0,9546 0,9850 1,4608
0,9387 1,7659 1,7684 0,9617 0,9864 1,3972 1,2440 1,4304 0,8217
3,4569 2,0566 1,8241 1,3289 0,8889 1,1086 1,0191 1,1875 1,3802
100 101 102 103 104 105 106 107 108
0,8681 0,8850 1,0294 1,4486 0,9915 0,6074 0,7932 0,7688 1,2034
0,8681 0,8850 1,0294 1,4486 0,9915 0,6074 0,7932 0,7688 1,2034
0,9606 1,2437 0,8780 1,3585 0,8517 0,7195 0,6988 0,8943 1,4299
0,9999 1,0470 1,1586 1,3886 0,5766 0,6432 0,6558 0,9386 1,2257
1,0325 1,0911 0,8482 1,3895 0,8674 0,5936 0,6814 1,0136 1,2069
0,9017 0,9726 0,8329 1,4669 0,9037 0,7605 0,7864 1,0966 1,3382
0,9677 1,2612 0,9626 1,5254 1,0408 0,6639 0,7893 0,8571 1,2778
1,0470 0,9918 1,2098 2,3799 0,9561 0,9562 0,9923 1,1567 1,2519
1,0812 1,2552 1,0005 1,4529 1,0409 0,6985 0,8478 0,9254 1,3282
8
109 110 111 112 113 114 115 116 117
1,1721 1,1531 0,9999 0,7409 1,8263 0,7229 0,6228 0,9118 0,5549
1,1721 1,1531 0,9999 0,7409 1,8263 0,7229 0,6228 0,9118 0,5549
1,0155 1,3541 1,3308 1,2685 2,0787 0,8786 0,7255 1,0321 0,6937
0,9092 1,1396 1,1098 1,0625 1,9140 0,4065 0,6035
1,2143 1,2319 1,3832 1,2348 2,1029 0,9172 0,5946 1,2217 0,5780
1,0211 1,3268 1,1485 1,0786 2,2807 0,6595 0,6900 0,9641 0,6242
1,2595 1,2361 1,1974 1,1174 1,6268 0,9808 0,6399 0,7210 0,5278
1,1244 1,2715 1,0617 1,0822 1,8909 0,7896 0,7285
1,2293 1,2286 1,2839 1,1339 1,5918 0,9796 0,6716 0,8881 0,6880
118 119 120 121 122 123 124 125 126
0,7712 0,6905 0,6228 1,2910 1,3170 0,8150 1,6221 2,8231 1,6568
0,7712 0,6905 0,6228 1,2910 1,3170 0,8150 1,6221 2,8231 1,6568
0,9101 0,8614 0,8169 1,8484 1,0935 1,1873 2,0491 3,2148 1,1227
1,3141 1,7513 1,0908 1,2285 1,8724 2,6598 0,6602
0,8599 0,6926 0,8640 1,9837 0,8355 1,3575 2,0405 3,2606 1,4039
0,7396 0,7568 0,9103 2,1630 1,0250 1,2709 2,0192 3,1814 0,9759
0,8039 0,6864 0,8414 2,5490 1,0396 1,1798 1,8698 2,6375 1,4633
1,2539 1,6998 1,1149 1,4278 2,1281 3,0025 1,6762
0,8500 0,6867 0,9068 2,3594 0,9526 1,2313 1,9984 2,6012 1,5346
127 128 129 130 131 132 133 134 135
1,7120 1,1270 1,1834 1,2986 1,4299 2,1638 1,6757 0,7074 0,5637
1,7120 1,1270 1,1834 1,2986 1,4299 2,1638 1,6757 0,7074 0,5637
2,6014 1,2760 1,6910 1,2314 1,3110 1,9648 1,0396 0,7477 0,8830
1,9042 1,5781 1,3930 1,0917 1,7509 2,2428 0,7460 0,5258 0,6446
1,8187 1,2443 1,3603 1,2440 1,2320 1,9908 1,2004 0,6176 0,8673
1,9223 1,1081 1,5440 1,1326 1,5110 2,0744 0,9794 0,8755 0,9303
1,6698 1,2908 1,6931 1,4497 1,2884 1,9458 1,3263 0,6740 0,9810
2,1693 1,6384 1,8283 1,6196 2,0139 2,1965 1,6621 1,0824 0,8890
1,7096 1,2788 1,6670 1,3723 1,4658 2,1938 1,2577 0,6368 1,0281
9
136 137 138 139 140 141 142 143 144
0,6951 0,7346 1,5059 0,8645 0,8593 1,3097 1,3951 1,7821 1,4264
0,6951 0,7346 1,5059 0,8645 0,8593 1,3097 1,3951 1,7821 1,4264
0,8882 1,2523 1,0402 0,8530 1,0166 1,2620 1,2400 1,8588 1,2117
0,8100 0,8449 0,8406 0,6606 0,9794 0,9480 1,2343 1,6803 0,9681
0,9368 1,1042 1,0948 0,8890 1,3178 1,2676 1,2656 1,9257 1,4739
1,0141 1,1232 1,0153 1,0773 1,3657 1,2556 1,4349 1,9263 1,3149
0,9203 1,1258 1,2004 0,9509 1,5767 1,2539 1,3938 2,1173 1,4749
0,9908 0,9333 0,8980 1,7480 1,3450 1,3026 2,0144 2,0151 1,2504
0,9957 1,2005 1,1479 0,9673 1,8219 1,3095 1,4806 2,2491 1,5025
145 146 147 148 149 150 151 152 153
1,0878 0,8069 0,7622 1,0054 1,1459 1,4901 0,7716 0,7602 0,7471
1,0878 0,8069 0,7622 1,0054 1,1459 1,4901 0,7716 0,7602 0,7471
1,0416 1,0016 0,9303 1,0337 1,6301 1,6222 1,1487 1,0553 1,1678
0,7356 1,1936 1,4401 1,0459 2,2607 1,2544 0,8563 0,7588 1,9081
1,0112 1,0316 1,2681 0,9181 1,8826 1,6282 1,2959 1,3360 1,1035
0,8877 1,0069 1,1530 1,0962 1,5746 1,6336 1,1729 0,9395 1,1284
1,0991 1,0993 1,3466 1,1707 1,8189 1,4477 1,3415 1,3989 1,2619
0,9289 1,2994 1,4745 1,1420 2,0171 1,1666 0,8538 0,7535 1,6791
1,0263 1,3403 1,8232 1,1904 1,8256 1,4769 1,3969 1,5489 1,1159
154 155 156 157 158 159
2,4194 1,3462 1,9727 0,9950 1,2923 1,5808
2,4194 1,3462 1,9727 0,9950 1,2923 1,5808
2,7880 3,0556 1,2122 0,9569 0,6622 3,2846
2,2447 2,4643 1,1624 1,3930 1,5783 3,3522
2,8993 2,0337 1,3745 1,1206 1,1910 2,2313
2,5373 2,9041 1,0521 0,7604 0,9153 2,4186
1,6285 1,9574 1,6319 1,1683 1,3028 2,1418
2,8749 2,1213 1,1869 1,2080 1,5938 2,5977
1,1949 2,1226 1,4226 0,9076 1,3174 2,2575
94
• Análise dos segmentos A e B
%Análise de um trecho em Curva (V=120Km/h) data_curva=dlmread( 'curva.txt' ); data_seguranca=dlmread( 'limites.txt' ); %posição curva km_curva=data_curva(:,1); loc_curva=data_curva(:,2); x_curva=(km_curva*1000+loc_curva)/1000; %posição em Km %LIMITES %Bitola bitola_curva=data_curva(:,4); %limites de segurança lsmin=data_seguranca(1:1881,1); lsmax=data_seguranca(1:1881,2); limin=data_seguranca(1:1881,3); limax=data_seguranca(1:1881,4); lamin=data_seguranca(1:1881,5); lamax=data_seguranca(1:1881,6); %Gráfico Bitola figure(1) plot(x_curva,bitola_curva, 'b' ,x_curva,lsmin, 'g' ,x_curva,lsmax, 'g' ,x_curva,limin, 'y' ,x_curva,limax, 'y' ,x_curva,lamin, 'r' ,x_curva,lamax, 'r' ) legend ( 'Bitola' , 'Limite de Segurança Min' , 'Limite de Segurança Max' , 'Limite de Intervenção Min' , 'Limite de Intervenção Max' , 'Limite de Alerta Min' , 'Limite de Alerta Min' ) xlabel( 'Posiçao (Km)' ) ylabel( 'Bitola (mm)' ) %Nivelamento Transversal nivt_curva=data_curva(1:1881,5); %limites de segurança ls_nivt=data_seguranca(1:1881,13); li_nivt=data_seguranca(1:1881,14); la_nivt=data_seguranca(1:1881,15); %Gráfico Nivelamento Transversal figure(2)
95
plot(x_curva,nivt_curva, 'b' ,x_curva,ls_nivt, 'g' ,x_curva,li_nivt, 'y' ,x_curva,la_nivt, 'r' ) legend ( 'Nivelamento Transversal' , 'Limite de Segurança' , 'Limite de Intervenção' , 'Limite de Alerta' ) xlabel( 'Posiçao (Km)' ) ylabel( 'Nivelamento Transversal (mm)' ) %Nivelamento Longitudinal D1 nivled1_curva=data_curva(1:1881,6); nivldd1_curva=data_curva(1:1881,7); %limites de segurança ls_nivld1=data_seguranca(1:1881,10); li_nivld1=data_seguranca(1:1881,11); la_nivld1=data_seguranca(1:1881,12); %Gráfico Nivelamento Longitudinal D1 figure(3) plot(x_curva,nivled1_curva, 'b' ,x_curva,nivldd1_curva, 'c' ,x_curva,ls_nivld1,'g' ,x_curva,li_nivld1, 'y' ,x_curva,la_nivld1, 'r' ) legend ( 'Nivelamento Longitudinal Esquerdo D1' , 'Nivelamento Longitudinal Direito D1' , 'Limite de Segurança' , 'Limite de Intervenção' , 'Limite de Alerta' ) xlabel( 'Posiçao (Km)' ) ylabel( 'Nivelamento Longitudinal D1 (mm)' ) %Nivelamento Longitudinal D2 nivled2_curva=data_curva(1:1881,8); nivldd2_curva=data_curva(1:1881,9); %limites de segurança (Não é Aplicável) %Gráfico Nivelamento Longitudinal D2 figure(4) plot(x_curva,nivled2_curva, 'b' ,x_curva,nivldd2_curva, 'c' ) legend ( 'Nivelamento Longitudinal Esquerdo D2' , 'Nivelamento Longitudinal Direito D2' ) xlabel( 'Posiçao (Km)' ) ylabel( 'Nivelamento Longitudinal D2 (mm)' ) %Alinhamento D1
96
alied1_curva=data_curva(1:1881,10); alidd1_curva=data_curva(1:1881,11); %limites de segurança ls_alid1=data_seguranca(1:1881,7); li_alid1=data_seguranca(1:1881,8); la_alid1=data_seguranca(1:1881,9); %Gráfico Alinhamento D1 figure(5) plot(x_curva,alied1_curva, 'b' ,x_curva,alidd1_curva, 'c' ,x_curva,ls_alid1, 'g',x_curva,li_alid1, 'y' ,x_curva,la_alid1, 'r' ) legend ( 'Alinhamento Esquerdo D1' , 'Alinhamento Direito D1' , 'Limite de Segurança' , 'Limite de Intervenção' , 'Limite de Alerta' ) xlabel( 'Posiçao (Km)' ) ylabel( 'Alinhamento D1 (mm)' ) %Alinhamento D2 alied2_curva=data_curva(1:1881,12); alidd2_curva=data_curva(1:1881,13); %limites de segurança (Não é Aplicável) %Gráfico Alinhamento D2 figure(6) plot(x_curva,alied2_curva, 'b' ,x_curva,alidd2_curva, 'c' ) legend ( 'Alinhamento Esquerdo D2' , 'alinhamento Direito D2' ) xlabel( 'Posiçao (Km)' ) ylabel( 'Alinhamento D2 (mm)' )
%Análise de um trecho em recta (V=120Km/h) data_recta=dlmread( 'recta.txt' ); data_seguranca=dlmread( 'limites.txt' );
97
%posição recta km_recta=data_recta(:,1); loc_recta=data_recta(:,2); x_recta=(km_recta*1000+loc_recta)/1000; %posição em Km %LIMITES %Bitola bitola_recta=data_recta(:,4); %limites de segurança lsmin=data_seguranca(:,1); lsmax=data_seguranca(:,2); limin=data_seguranca(:,3); limax=data_seguranca(:,4); lamin=data_seguranca(:,5); lamax=data_seguranca(:,6); %Gráfico Bitola figure(1) plot(x_recta,bitola_recta, 'b' ,x_recta,lsmin, 'g' ,x_recta,lsmax, 'g' ,x_recta,limin, 'y' ,x_recta,limax, 'y' ,x_recta,lamin, 'r' ,x_recta,lamax, 'r' ) legend ( 'Bitola' , 'Lim. Segurança Min.' , 'Lim. Segurança Max.' , 'Lim. Intervenção Min.' , 'Lim. Intervenção Max.' , 'Lim. Alerta Min.' , 'Lim. Alerta Min.' ) xlabel( 'Posiçao (Km)' ) ylabel( 'Bitola (mm)' ) AXIS([200 201.4 -15 40]) %Nivelamento Transversal nivt_recta=data_recta(:,5); %limites de segurança ls_nivt=data_seguranca(:,13); li_nivt=data_seguranca(:,14); la_nivt=data_seguranca(:,15); %Gráfico Nivelamento Transversal figure(8) plot(x_recta,nivt_recta, 'b' ,x_recta,ls_nivt, 'g' ,x_recta,li_nivt, 'y' ,x_recta,la_nivt, 'r' ) legend ( 'Nivelamento Transversal' , 'Limite de Segurança' , 'Limite de Intervenção' , 'Limite de Alerta' ) xlabel( 'Posiçao (Km)' ) ylabel( 'Nivelamento Transversal (mm)' )
98
%Nivelamento Longitudinal D1 nivled1_recta=data_recta(:,6); nivldd1_recta=data_recta(:,7); %limites de segurança ls_nivld1=data_seguranca(:,10); li_nivld1=data_seguranca(:,11); la_nivld1=data_seguranca(:,12); %Gráfico Nivelamento Longitudinal D1 figure(9) plot(x_recta,nivled1_recta, 'b' ,x_recta,nivldd1_recta, 'c' ,x_recta,ls_nivld1,'g' ,x_recta,li_nivld1, 'y' ,x_recta,la_nivld1, 'r' ) legend ( 'Nivelamento Longitudinal Esquerdo D1' , 'Nivelamento Longitudinal Direito D1' , 'Limite de Segurança' , 'Limite de Intervenção' , 'Limite de Alerta' ) xlabel( 'Posiçao (Km)' ) ylabel( 'Nivelamento Longitudinal D1 (mm)' ) %Nivelamento Longitudinal D2 nivled2_recta=data_recta(:,8); nivldd2_recta=data_recta(:,9); %limites de segurança (Não é Aplicável) %Gráfico Nivelamento Longitudinal D2 figure(10) plot(x_recta,nivled2_recta, 'b' ,x_recta,nivldd2_recta, 'c' ) legend ( 'Nivelamento Longitudinal Esquerdo D2' , 'Nivelamento Longitudinal Direito D2' ) xlabel( 'Posiçao (Km)' ) ylabel( 'Nivelamento Longitudinal D2 (mm)' ) %Alinhamento D1 alied1_recta=data_recta(:,10); alidd1_recta=data_recta(:,11); %limites de segurança ls_alid1=data_seguranca(:,7); li_alid1=data_seguranca(:,8); la_alid1=data_seguranca(:,9);
99
%Gráfico Alinhamento D1 figure(11) plot(x_recta,alied1_recta, 'b' ,x_recta,alidd1_recta, 'c' ,x_recta,ls_alid1, 'g',x_recta,li_alid1, 'y' ,x_recta,la_alid1, 'r' ) legend ( 'Alinhamento Esquerdo D1' , 'Alinhamento Direito D1' , 'Limite de Segurança' , 'Limite de Intervenção' , 'Limite de Alerta' ) xlabel( 'Posiçao (Km)' ) ylabel( 'Alinhamento D1 (mm)' ) %Alinhamento D2 alied2_recta=data_recta(:,12); alidd2_recta=data_recta(:,13); %limites de segurança (Não é Aplicável) %Gráfico Alinhamento D2 figure(12) plot(x_recta,alied2_recta, 'b' ,x_recta,alidd2_recta, 'c' ) legend ( 'Alinhamento Esquerdo D2' , 'alinhamento Direito D2' ) xlabel( 'Posiçao (Km)' ) ylabel( 'Alinhamento D2 (mm)' )
• Sincronização temporal das campanhas de inspecção de 2005 a 2009
%Universidade do Porto %Faculdade de Engenharia %DEC 2010 %Clara Moura %2009 (REF) and 2008 Junho % troco total %loading data % 2008 data_08=dlmread( 'campanha2008pronto.txt' ); % 2009 REF data_09=dlmread( 'campanha2009pronto.txt' ); %posicao 2008 km_08j=data_08(:,1); loc_08j=data_08(:,2)*10^-3;
100
x_08j=[loc_08j(1:4000)+km_08j(1:4000) loc_08j(4001:8000)+km_08j(4001:8000) loc_08j(8001:12000)+km_08j(8001:12000) loc_08j(12001:16000)+km_08j(12001:16000) loc_08j(16001:20000)+km_08j(16001:20000) loc_08j(20001:24062)+km_08j(20001:24062) loc_08j(24063:28062)+km_08j(24063:28062) loc_08j(28063:32087)+km_08j(28063:32087) loc_08j(32088:36087)+km_08j(32088:36087) loc_08j(36088:40087)+km_08j(36088:40087) loc_08j(40088:44087)+km_08j(40088:44087) loc_08j(44088:48087)+km_08j(44088:48087) loc_08j(48088:52087)+km_08j(48088:52087) loc_08j(52088:56087)+km_08j(52088:56087) loc_08j(56088:60087)+km_08j(56088:60087) loc_08j(60088:64087)+km_08j(60088:64087) loc_08j(64088:68087)+km_08j(64088:68087) loc_08j(68088:72149)+km_08j(68088:72149) loc_08j(72150:76149)+km_08j(72150:76149) loc_08j(76150:80149)+km_08j(76150:80149) loc_08j(80150:84149)+km_08j(80150:84149) loc_08j(84150:88149)+km_08j(84150:88149) loc_08j(88150:92149)+km_08j(88150:92149) loc_08j(92150:96149)+km_08j(92150:96149) loc_08j(96150:100149)+km_08j(96150:100149) loc_08j(100150:104149)+km_08j(100150:104149) loc_08j(104150:108149)+km_08j(104150:108149) loc_08j(108150:112267)+km_08j(108150:112267) loc_08j(112268:116267)+km_08j(112268:116267) loc_08j(116268:120267)+km_08j(116268:120267) loc_08j(120268:124267)+km_08j(120268:124267) loc_08j(124268:128284)+km_08j(124268:128284)]; %posicao 2009 km_09=data_09(:,1); loc_09=data_09(:,2)*10^-3; x_09=[loc_09(1:4000)+km_09(1:4000) loc_09(4001:8000)+km_09(4001:8000) loc_09(8001:12000)+km_09(8001:12000) loc_09(12001:16000)+km_09(12001:16000) loc_09(16001:20057)+km_09(16001:20057) loc_09(20058:24057)+km_09(20058:24057) loc_09(24058:28057)+km_09(24058:28057) loc_09(28058:32057)+km_09(28058:32057) loc_09(32058:36057)+km_09(32058:36057) loc_09(36058:40095)+km_09(36058:40095) loc_09(40096:44095)+km_09(40096:44095) loc_09(44096:48095)+km_09(44096:48095) loc_09(48096:52095)+km_09(48096:52095) loc_09(52096:56095)+km_09(52096:56095) loc_09(56096:60155)+km_09(56096:60155)
101
loc_09(60156:64155)+km_09(60156:64155) loc_09(64156:68155)+km_09(64156:68155) loc_09(68156:72155)+km_09(68156:72155) loc_09(72156:76155)+km_09(72156:76155) loc_09(76156:80155)+km_09(76156:80155) loc_09(80156:84155)+km_09(80156:84155) loc_09(84156:88227)+km_09(84156:88227) loc_09(88228:92227)+km_09(88228:92227) loc_09(92228:96227)+km_09(92228:96227) loc_09(96228:100227)+km_09(96228:100227) loc_09(100228:104227)+km_09(100228:104227) loc_09(104228:108227)+km_09(104228:108227) loc_09(108228:112300)+km_09(108228:112300) loc_09(112301:116300)+km_09(112301:116300) loc_09(116301:120300)+km_09(116301:120300) loc_09(120301:124300)+km_09(120301:124300) loc_09(124301:128329)+km_09(124301:128329) loc_09(128330:132329)+km_09(128330:132329) loc_09(132330:136329)+km_09(132330:136329)]; %ajuste em 2008 total, esquerdo x_08_ajuste_E=x_08j-1.333*10^-3; %ajuste em 2008 total, direito x_08_ajuste_D=x_08j-1.7167*10^-3; %Ali carril esquerdo %ali E d1 2008 total ali_08_E=data_08(:,4); %ali E d1 2009 total ali_09_E=data_09(:,4); %Ali carril direito %ali D d1 2008 total ali_08_D=data_08(:,5); %ali D d1 2009 total ali_09_D=data_09(:,5); % 2008, total, Esquerdo, teste figure(1) plot(x_08j,ali_08_E, 'b' ) legend ( 'Ali E d1 2008 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '% 2008, total, Esquerdo, teste' ) % 2009, total, Esquerdo, teste figure(2)
102
plot(x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2009 Jun total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '% 2009, total, Esquerdo, teste' ) % 2008 VS 2009, total, Esquerdo figure(3) plot(x_08j,ali_08_E, 'b' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2008 total' , 'Ali E d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 VS 2009, total, Esquerdo' ) % 2008 VS 2009, total, Direito figure(4) plot(x_08j,ali_08_D, 'b' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2008 total' , 'Ali D d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 VS 2009, total, Direito' ) % 2008 VS 2009, total, Esquerdo, ajuste figure(5) plot(x_08_ajuste_E,ali_08_E, 'b' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2008 total ajuste' , 'Ali E d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 VS 2009, total, Esquerdo, ajuste' ) % 2008 VS 2009, total, Direito, ajuste figure(6) plot(x_08_ajuste_D,ali_08_D, 'b' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2008 total ajuste' , 'Ali D d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 VS 2009, total, Direito, ajuste' ) % calculo do desvio padrao em segmentos de 200 m, E squerda, 05, ajustado n=fix(size(x_08_ajuste_E)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_08_E([n1(i):n2(i)],:)])];
103
dp_08pronto_E=dp'; end % calculo do desvio padrao em segmentos de 200 m, D ireita, 05, ajustado n=fix(size(x_08_ajuste_D)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_08_D([n1(i):n2(i)],:)])]; dp_08pronto_D=dp'; end %Sincronizacao 1km %Parametro D1 % %1 km (1,4000) %posicao 2008 1km km_um_08=data_08(1:4000,1); loc_um_08=loc_08j(1:4000); x_um_08=km_um_08+loc_um_08; %posicao 2009 1km km_um_09=data_09(1:4000,1); loc_um_09=loc_09(1:4000); x_um_09=km_um_09+loc_um_09; %0,5 km (1,2000) %posicao 2008 0,5km km_meio_08=data_08(1:2000,1); loc_meio_08=loc_08j(1:2000); x_meio_08=km_meio_08+loc_meio_08; %ajuste em 2008 para 0,5 km esquerdo x_meio_08_ajuste_E=x_meio_08-1.333*10^-3
104
%ajuste em 2008 para 0,5 km direito x_meio_08_ajuste_D=x_meio_08-1.7167*10^-3 %posicao 2009 0,5km km_meio_09=data_09(1:2000,1); loc_meio_09=loc_09(1:2000); x_meio_09=km_meio_09+loc_meio_09; %Ali carril esquerdo %ali E D1 2008 1km ali_um_08_E=data_08(1:4000,4); ali_meio_08_E=data_08(1:2000,4); %ali E D1 2009 1km REF ali_um_09_E=data_09(1:4000,4); ali_meio_09_E=data_09(1:2000,4); % 2008 VS 2009, 1km, Esquerda figure(7) plot(x_um_08,ali_um_08_E, 'b' ,x_um_09,ali_um_09_E, 'g' ) legend ( 'Ali E1 2008 1km' , 'Ali E1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 VS 2009, 1km, Esquerda' ) % 2008 VS 2009, 0,5km, Esquerda figure(8) plot(x_meio_08,ali_meio_08_E, 'b' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2008 0,5km' , 'Ali E1 2009 0,5km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 VS 2009, 0,5km, Esquerda' ) % 2008 VS 2009, 0,5km, c/ ajuste figure(9) plot(x_meio_08_ajuste_E,ali_meio_08_E, 'b' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2008 0,5 km ajuste' , 'Ali E1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 VS 2009, 0,5km, c/ ajuste, Esquerda' ) %Ali Carril Direito %ali D D1 2008 1km ali_um_08_D=data_08(1:4000,5); ali_meio_08_D=data_08(1:2000,5);
105
%ali D D1 2009 1km REF ali_um_09_D=data_09(1:4000,5); ali_meio_09_D=data_09(1:2000,5); % 2008 VS 2009, 1km figure(10) plot(x_um_08,ali_um_08_D, 'b' ,x_um_09,ali_um_09_D, 'g' ) legend ( 'Ali D d1 2008 1km' , 'Ali D d1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 VS 2009, 1km, Direito' ) % 2008 VS 2009, 0,5km figure(11) plot(x_meio_08,ali_meio_08_D, 'b' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2008 0,5 km' , 'Ali D d1 2009 0,5 km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 VS 2009, 0,5km, Direito' ) % 2008 VS 2009, 0,5km, c/ ajuste, direito figure(12) plot(x_meio_08_ajuste_D,ali_meio_08_D, 'b' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2008 0,5 km ajuste' , 'Ali D d1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 VS 2009, 0,5km, c/ ajuste, direito' ) %Universidade do Porto %Faculdade de Engenharia %DEC 2010 %Clara Moura %2009 (REF) and 2008 Março % troco total %loading data % 2008 data_08=dlmread( 'campanha2008mar.txt' ); % 2009 REF data_09=dlmread( 'campanha2009pronto.txt' ); %posicao 2008 km_08=data_08(:,1); loc_08=data_08(:,2)*10^-3; x_08=[loc_08(1:4000)+km_08(1:4000)
106
loc_08(4001:8000)+km_08(4001:8000) loc_08(8001:12000)+km_08(8001:12000) loc_08(12001:16000)+km_08(12001:16000) loc_08(16001:20000)+km_08(16001:20000) loc_08(20001:24000)+km_08(20001:24000) loc_08(24001:28000)+km_08(24001:28000) loc_08(28001:32288)+km_08(28001:32288) loc_08(32289:36288)+km_08(32289:36288) loc_08(36289:40288)+km_08(36289:40288) loc_08(40289:44288)+km_08(40289:44288) loc_08(44289:48288)+km_08(44289:48288) loc_08(48289:52288)+km_08(48289:52288) loc_08(52289:56288)+km_08(52289:56288) loc_08(56289:60288)+km_08(56289:60288) loc_08(60289:64288)+km_08(60289:64288) loc_08(64289:68288)+km_08(64289:68288) loc_08(68289:72288)+km_08(68289:72288) loc_08(72289:76288)+km_08(72289:76288) loc_08(76289:80288)+km_08(76289:80288) loc_08(80289:84288)+km_08(80289:84288) loc_08(84289:88288)+km_08(84289:88288) loc_08(88289:92426)+km_08(88289:92426) loc_08(92427:96426)+km_08(92427:96426) loc_08(96427:100426)+km_08(96427:100426) loc_08(100427:104426)+km_08(100427:104426) loc_08(104427:108426)+km_08(104427:108426) loc_08(108427:112426)+km_08(108427:112426) loc_08(112427:116426)+km_08(112427:116426) loc_08(116427:120426)+km_08(116427:120426) loc_08(120427:124426)+km_08(120427:124426) loc_08(124427:128114)+km_08(124427:128114)]; %posicao 2009 km_09=data_09(:,1); loc_09=data_09(:,2)*10^-3; x_09=[loc_09(1:4000)+km_09(1:4000) loc_09(4001:8000)+km_09(4001:8000) loc_09(8001:12000)+km_09(8001:12000) loc_09(12001:16000)+km_09(12001:16000) loc_09(16001:20057)+km_09(16001:20057) loc_09(20058:24057)+km_09(20058:24057) loc_09(24058:28057)+km_09(24058:28057) loc_09(28058:32057)+km_09(28058:32057) loc_09(32058:36057)+km_09(32058:36057) loc_09(36058:40095)+km_09(36058:40095) loc_09(40096:44095)+km_09(40096:44095) loc_09(44096:48095)+km_09(44096:48095) loc_09(48096:52095)+km_09(48096:52095) loc_09(52096:56095)+km_09(52096:56095) loc_09(56096:60155)+km_09(56096:60155) loc_09(60156:64155)+km_09(60156:64155) loc_09(64156:68155)+km_09(64156:68155)
107
loc_09(68156:72155)+km_09(68156:72155) loc_09(72156:76155)+km_09(72156:76155) loc_09(76156:80155)+km_09(76156:80155) loc_09(80156:84155)+km_09(80156:84155) loc_09(84156:88227)+km_09(84156:88227) loc_09(88228:92227)+km_09(88228:92227) loc_09(92228:96227)+km_09(92228:96227) loc_09(96228:100227)+km_09(96228:100227) loc_09(100228:104227)+km_09(100228:104227) loc_09(104228:108227)+km_09(104228:108227) loc_09(108228:112300)+km_09(108228:112300) loc_09(112301:116300)+km_09(112301:116300) loc_09(116301:120300)+km_09(116301:120300) loc_09(120301:124300)+km_09(120301:124300) loc_09(124301:128329)+km_09(124301:128329) loc_09(128330:132329)+km_09(128330:132329) loc_09(132330:136329)+km_09(132330:136329)]; %ajuste em 2008 total, esquerdo x_08_ajuste_E=x_08-4.5*10^-3; %ajuste em 2008 total, direito x_08_ajuste_D=x_08-2.4167*10^-3; %Ali carril esquerdo %ali E d1 2008 total ali_08_E=data_08(:,4); %ali E d1 2009 total ali_09_E=data_09(:,4); %Ali carril direito %ali D d1 2008 total ali_08_D=data_08(:,5); %ali D d1 2009 total ali_09_D=data_09(:,5); % 2008 Mar VS 2009, total, Esquerdo, teste figure(1) plot(x_08,ali_08_E, 'b' ) legend ( 'Ali E d1 2008 Mar total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) % 2008 Mar VS 2009, total, Esquerdo figure(3) plot(x_08,ali_08_E, 'b' ,x_09,ali_09_E, 'c' ) legend ( 'Ali E d1 2008 Mar total' , 'Ali E d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' )
108
title( '2008 Mar VS 2009, total, Esquerdo' ) % 2008 Mar VS 2009, total, Direito figure(4) plot(x_08,ali_08_D, 'b' ,x_09,ali_09_D, 'c' ) legend ( 'Ali D d1 2008 Mar total' , 'Ali D d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 Mar VS 2009, total, Direito' ) % 2008 Mar VS 2009, total, Esquerdo, ajuste figure(5) plot(x_08_ajuste_E,ali_08_E, 'b' ,x_09,ali_09_E, 'c' ) legend ( 'Ali E d1 2008 Jun total ajuste' , 'Ali E d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 Mar VS 2009, total, Esquerdo, ajuste' ) % 2008 Mar VS 2009, total, Direito, ajuste figure(6) plot(x_08_ajuste_D,ali_08_D, 'b' ,x_09,ali_09_D, 'c' ) legend ( 'Ali D d1 2008 Mar total ajuste' , 'Ali D d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 Mar VS 2009, total, Direito, ajuste' ) % calculo do desvio padrao em segmentos de 200 m, E squerda, 05, ajustado n=fix(size(x_08_ajuste_E)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_08_E([n1(i):n2(i)],:)])]; dp_08mar_E=dp'; end % calculo do desvio padrao em segmentos de 200 m, D ireita, 05, ajustado n=fix(size(x_08_ajuste_D)/800); for i=1:n
109
n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_08_D([n1(i):n2(i)],:)])]; dp_08mar_D=dp'; end %Sincronizacao 1km %Parametro D1 %loading data % 1km 2008 data_um_08=dlmread( 'campanha2008jun1kmi.txt' ); % 1km 2009 REF data_um_09=dlmread( 'campanha2009pronto1kmi.txt' ); %1 km (1,4000) %posicao 2008 1km km_um_08=data_um_08(:,1); loc_um_08=data_um_08(:,2); x_um_08=km_um_08+loc_um_08; %posicao 2009 1km km_um_09=data_um_09(:,1); loc_um_09=data_um_09(:,2); x_um_09=km_um_09+loc_um_09; %0,5 km (1,2000) %posicao 2008 0,5km km_meio_08=data_um_08(1:2000,1); loc_meio_08=data_um_08(1:2000,2); x_meio_08=km_meio_08+loc_meio_08*10^-3; %ajuste em 2008 para 0,5 km esquerdo x_meio_08_ajuste_E=x_meio_08-4.333*10^-3; %ajuste em 2008 para 0,5 km direito x_meio_08_ajuste_D=x_meio_08-51*10^-3 %posicao 2009 0,5km
110
km_meio_09=data_um_09(1:2000,1); loc_meio_09=data_um_09(1:2000,2); x_meio_09=km_meio_09+loc_meio_09*10^-3; %Ali carril esquerdo %ali E D1 2008 1km ali_um_08_E=data_um_08(:,4); ali_meio_08_E=data_um_08(1:2000,4); %ali E D1 2009 1km REF ali_um_09_E=data_um_09(:,4); ali_meio_09_E=data_um_09(1:2000,4); % 2008 Jun VS 2009, 1km Esquerda figure(7) plot(x_um_08,ali_um_08_E, 'b' ,x_um_09,ali_um_09_E, 'm' ) legend ( 'Ali E1 2008 Mar 1km' , 'Ali E1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 Mar VS 2009, 1km, Esquerda' ) % 2008 Jun VS 2009, 0,5km , Esquerda figure(8) plot(x_meio_08,ali_meio_08_E, 'b' ,x_meio_09,ali_meio_09_E, 'm' ) legend ( 'Ali E1 2008 Mar 0,5km' , 'Ali E1 2009 0,5km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 Mar VS 2009, 0,5km, Esquerda' ) % 2008 Mar VS 2009, 0,5km, c/ ajuste figure(9) plot(x_meio_08_ajuste_E,ali_meio_08_E, 'b' ,x_meio_09,ali_meio_09_E, 'm' ) legend ( 'Ali E1 2008 Mar 0,5 km ajuste' , 'Ali E1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 Mar VS 2009, 0,5km, c/ ajuste, Esquerda' ) %Ali Carril Direito %ali D D1 2008 1km ali_um_08_D=data_um_08(:,5); ali_meio_08_D=data_um_08(1:2000,5); %ali D D1 2009 1km REF ali_um_09_D=data_um_09(:,5); ali_meio_09_D=data_um_09(1:2000,5); % 2008 Mar VS 2009, 1km figure(10) plot(x_um_08,ali_um_08_D, 'b' ,x_um_09,ali_um_09_D, 'g' ) legend ( 'Ali D d1 2008 Mar 1km' , 'Ali D d1 2009 1km' )
111
xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 Mar VS 2009, 1km, Direito' ) % 2008 Mar VS 2009, 0,5km figure(11) plot(x_meio_08,ali_meio_08_D, 'b' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2008 Mar 0,5 km' , 'Ali D d1 2009 0,5 km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 Mar VS 2009, 0,5km, Direito' ) % 2008 Mar VS 2009, 0,5km, c/ ajuste, direito figure(12) plot(x_meio_08_ajuste_D,ali_meio_08_D, 'b' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2008 Mar 0,5 km ajuste' , 'Ali D d1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2008 Mar VS 2009, 0,5km, c/ ajuste, direito' ) %Universidade do Porto %Faculdade de Engenharia %DEC 2010 %Clara Moura %2009 (REF) and 2007 Setembro % troco total %loading data % 2007 data_07=dlmread( 'campanha2007set.txt' ); % 2009 REF data_09=dlmread( 'campanha2009pronto.txt' ); %posicao 2007 km_07=data_07(:,1); loc_07=data_07(:,2)*10^-3; x_07=[loc_07(1:4000)+km_07(1:4000) loc_07(4001:8000)+km_07(4001:8000) loc_07(8001:12057)+km_07(8001:12057) loc_07(12058:16057)+km_07(12058:16057) loc_07(16058:20057)+km_07(16058:20057) loc_07(20058:24057)+km_07(20058:24057) loc_07(24058:28057)+km_07(24058:28057) loc_07(28058:32092)+km_07(28058:32092)
112
loc_07(32093:36092)+km_07(32093:36092) loc_07(36093:40092)+km_07(36093:40092) loc_07(40093:44092)+km_07(40093:44092) loc_07(44093:48124)+km_07(44093:48124) loc_07(48125:52124)+km_07(48125:52124) loc_07(52125:56124)+km_07(52125:56124) loc_07(56125:60124)+km_07(56125:60124) loc_07(60125:64124)+km_07(60125:64124) loc_07(64125:68124)+km_07(64125:68124) loc_07(68125:72124)+km_07(68125:72124) loc_07(72125:76124)+km_07(72125:76124) loc_07(76125:80124)+km_07(76125:80124) loc_07(80125:84124)+km_07(80125:84124) loc_07(84125:88124)+km_07(84125:88124) loc_07(88125:92194)+km_07(88125:92194) loc_07(92195:96194)+km_07(92195:96194) loc_07(96195:100194)+km_07(96195:100194) loc_07(100195:104194)+km_07(100195:104194) loc_07(104195:108117)+km_07(104195:108117) loc_07(108118:112262)+km_07(108118:112262) loc_07(112263:116243)+km_07(112263:116243) loc_07(116244:120243)+km_07(116244:120243) loc_07(120244:124243)+km_07(120244:124243) loc_07(124244:128152)+km_07(124244:128152)]; %posicao 2009 km_09=data_09(:,1); loc_09=data_09(:,2)*10^-3; x_09=[loc_09(1:4000)+km_09(1:4000) loc_09(4001:8000)+km_09(4001:8000) loc_09(8001:12000)+km_09(8001:12000) loc_09(12001:16000)+km_09(12001:16000) loc_09(16001:20057)+km_09(16001:20057) loc_09(20058:24057)+km_09(20058:24057) loc_09(24058:28057)+km_09(24058:28057) loc_09(28058:32057)+km_09(28058:32057) loc_09(32058:36057)+km_09(32058:36057) loc_09(36058:40095)+km_09(36058:40095) loc_09(40096:44095)+km_09(40096:44095) loc_09(44096:48095)+km_09(44096:48095) loc_09(48096:52095)+km_09(48096:52095) loc_09(52096:56095)+km_09(52096:56095) loc_09(56096:60155)+km_09(56096:60155) loc_09(60156:64155)+km_09(60156:64155) loc_09(64156:68155)+km_09(64156:68155) loc_09(68156:72155)+km_09(68156:72155) loc_09(72156:76155)+km_09(72156:76155) loc_09(76156:80155)+km_09(76156:80155) loc_09(80156:84155)+km_09(80156:84155) loc_09(84156:88227)+km_09(84156:88227) loc_09(88228:92227)+km_09(88228:92227) loc_09(92228:96227)+km_09(92228:96227)
113
loc_09(96228:100227)+km_09(96228:100227) loc_09(100228:104227)+km_09(100228:104227) loc_09(104228:108227)+km_09(104228:108227) loc_09(108228:112300)+km_09(108228:112300) loc_09(112301:116300)+km_09(112301:116300) loc_09(116301:120300)+km_09(116301:120300) loc_09(120301:124300)+km_09(120301:124300) loc_09(124301:128329)+km_09(124301:128329) loc_09(128330:132329)+km_09(128330:132329) loc_09(132330:136329)+km_09(132330:136329)]; % %ajuste em 2007 total, esquerdo x_07_ajuste_E=x_07-1.0*10^-3 %ajuste em 2007 total, direito x_07_ajuste_D=x_07-2.0*10^-3 %Ali carril esquerdo %ali E d1 2007 total ali_07_E=data_07(:,4); %ali E d1 2009 total ali_09_E=data_09(:,4); %Ali carril direito %ali D d1 2007 total ali_07_D=data_07(:,5); %ali D d1 2009 total ali_09_D=data_09(:,5); % 2007, total, Esquerdo, teste figure(1) plot(x_07,ali_07_E, 'b' ) legend ( 'Ali E d1 2007 Set total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '% 2007, total, Esquerdo, teste' ) % 2009, total, Esquerdo, teste figure(2) plot(x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2009 Jun total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '% 2009, total, Esquerdo, teste' ) % 2007 Set VS 2009, total, Esquerdo figure(3) plot(x_07,ali_07_E, 'b' ,x_09,ali_09_E, 'g' )
114
legend ( 'Ali E d1 2007 Set total' , 'Ali E d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Set VS 2009, total, Esquerdo' ) % 2007 Set VS 2009, total, Direito figure(4) plot(x_07,ali_07_D, 'b' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2007 Set total' , 'Ali D d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Set VS 2009, total, Direito' ) % 2007 Set VS 2009, total, Esquerdo, ajuste figure(5) plot(x_07_ajuste_E,ali_07_E, 'b' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2007 Set total ajuste' , 'Ali E d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Set VS 2009, total, Esquerdo, ajuste' ) % 2007 Set VS 2009, total, Direito, ajuste figure(6) plot(x_07_ajuste_D,ali_07_D, 'b' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2007 Set total ajuste' , 'Ali D d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Set VS 2009, total, Direito, ajuste' ) % calculo do desvio padrao em segmentos de 200 m, E squerda, 05, ajustado n=fix(size(x_07_ajuste_E)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_07_E([n1(i):n2(i)],:)])]; dp_07set_E=dp'; end % calculo do desvio padrao em segmentos de 200 m, D ireita, 05, ajustado n=fix(size(x_07_ajuste_D)/800);
115
for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_07_D([n1(i):n2(i)],:)])]; dp_07set_D=dp'; end %Sincronizacao 1km %Parametro D1 % %1 km (1,4000) %posicao 2007 1km km_um_07=data_07(1:4000,1); loc_um_07=loc_07(1:4000); x_um_07=km_um_07+loc_um_07; %posicao 2009 1km km_um_09=data_09(1:4000,1); loc_um_09=loc_09(1:4000); x_um_09=km_um_09+loc_um_09; %0,5 km (1,2000) %posicao 2007 0,5km km_meio_07=data_07(1:2000,1); loc_meio_07=loc_07(1:2000); x_meio_07=km_meio_07+loc_meio_07; %ajuste em 2007 para 0,5 km esquerdo x_meio_07_ajuste_E=x_meio_07-1.0*10^-3 %ajuste em 2007 para 0,5 km direito x_meio_07_ajuste_D=x_meio_07-2.0*10^-3 %posicao 2009 0,5km km_meio_09=data_09(1:2000,1); loc_meio_09=loc_09(1:2000); x_meio_09=km_meio_09+loc_meio_09;
116
%Ali carril esquerdo %ali E D1 2007 1km ali_um_07_E=data_07(1:4000,4); ali_meio_07_E=data_07(1:2000,4); %ali E D1 2009 1km REF ali_um_09_E=data_09(1:4000,4); ali_meio_09_E=data_09(1:2000,4); % 2007 Set VS 2009, 1km, Esquerda figure(7) plot(x_um_07,ali_um_07_E, 'b' ,x_um_09,ali_um_09_E, 'g' ) legend ( 'Ali E1 2007 Set 1km' , 'Ali E1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Set VS 2009, 1km, Esquerda' ) % 2007 Set VS 2009, 0,5km, Esquerda figure(8) plot(x_meio_07,ali_meio_07_E, 'b' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2007 Set 0,5km' , 'Ali E1 2009 0,5km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Set VS 2009, 0,5km, Esquerda' ) % 2007 Set VS 2009, 0,5km, c/ ajuste figure(9) plot(x_meio_07_ajuste_E,ali_meio_07_E, 'b' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2007 Set 0,5 km ajuste' , 'Ali E1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Set VS 2009, 0,5km, c/ ajuste, Esquerda' ) %Ali Carril Direito %ali D D1 2007 1km ali_um_07_D=data_07(1:4000,5); ali_meio_07_D=data_07(1:2000,5); %ali D D1 2009 1km REF ali_um_09_D=data_09(1:4000,5); ali_meio_09_D=data_09(1:2000,5); % 2007 Set VS 2009, 1km figure(10) plot(x_um_07,ali_um_07_D, 'b' ,x_um_09,ali_um_09_D, 'g' ) legend ( 'Ali D d1 2007 Set 1km' , 'Ali D d1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' )
117
title( '2007 Set VS 2009, 1km, Direito' ) % 2007 Set VS 2009, 0,5km figure(11) plot(x_meio_07,ali_meio_07_D, 'b' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2007 Set 0,5 km' , 'Ali D d1 2009 0,5 km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Set VS 2009, 0,5km, Direito' ) % 2007 Set VS 2009, 0,5km, c/ ajuste, direito figure(12) plot(x_meio_07_ajuste_D,ali_meio_07_D, 'b' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2007 Set 0,5 km ajuste' , 'Ali D d1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Set VS 2009, 0,5km, c/ ajuste, direito' ) %Universidade do Porto %Faculdade de Engenharia %DEC 2010 %Clara Moura %2009 (REF) and 2007 Março % troco total %loading data % 2007 data_07=dlmread( 'campanha07mar.txt' ); % 2009 REF data_09=dlmread( 'campanha2009pronto.txt' ); %posicao 2007 km_07m=data_07(:,1); loc_07m=data_07(:,2)*10^-3; x_07m=[loc_07m(1:4000)+km_07m(1:4000) loc_07m(4001:8000)+km_07m(4001:8000) loc_07m(8001:11991)+km_07m(8001:11991) loc_07m(11992:15991)+km_07m(11992:15991) loc_07m(15992:19991)+km_07m(15992:19991) loc_07m(19992:23980)+km_07m(19992:23980) loc_07m(23981:27980)+km_07m(23981:27980) loc_07m(27981:31971)+km_07m(27981:31971) loc_07m(31972:35971)+km_07m(31972:35971) loc_07m(35972:39971)+km_07m(35972:39971) loc_07m(39972:43971)+km_07m(39972:43971)
118
loc_07m(43972:47971)+km_07m(43972:47971) loc_07m(47972:51971)+km_07m(47972:51971) loc_07m(51972:55956)+km_07m(51972:55956) loc_07m(55957:59956)+km_07m(55957:59956) loc_07m(59957:63946)+km_07m(59957:63946) loc_07m(63947:67946)+km_07m(63947:67946) loc_07m(67947:71930)+km_07m(67947:71930) loc_07m(71931:75930)+km_07m(71931:75930) loc_07m(75931:79930)+km_07m(75931:79930) loc_07m(79931:83930)+km_07m(79931:83930) loc_07m(83931:87926)+km_07m(83931:87926) loc_07m(87927:91926)+km_07m(87927:91926) loc_07m(91927:95926)+km_07m(91927:95926) loc_07m(95927:99926)+km_07m(95927:99926) loc_07m(99927:103926)+km_07m(99927:103926) loc_07m(103927:107926)+km_07m(103927:107926) loc_07m(107927:111931)+km_07m(107927:111931) loc_07m(111932:115887)+km_07m(111932:115887) loc_07m(115888:119887)+km_07m(115888:119887) loc_07m(119888:123887)+km_07m(119888:123887) loc_07m(123888:127896)+km_07m(123888:127896)]; %posicao 2009 km_09=data_09(:,1); loc_09=data_09(:,2)*10^-3; x_09=[loc_09(1:4000)+km_09(1:4000) loc_09(4001:8000)+km_09(4001:8000) loc_09(8001:12000)+km_09(8001:12000) loc_09(12001:16000)+km_09(12001:16000) loc_09(16001:20057)+km_09(16001:20057) loc_09(20058:24057)+km_09(20058:24057) loc_09(24058:28057)+km_09(24058:28057) loc_09(28058:32057)+km_09(28058:32057) loc_09(32058:36057)+km_09(32058:36057) loc_09(36058:40095)+km_09(36058:40095) loc_09(40096:44095)+km_09(40096:44095) loc_09(44096:48095)+km_09(44096:48095) loc_09(48096:52095)+km_09(48096:52095) loc_09(52096:56095)+km_09(52096:56095) loc_09(56096:60155)+km_09(56096:60155) loc_09(60156:64155)+km_09(60156:64155) loc_09(64156:68155)+km_09(64156:68155) loc_09(68156:72155)+km_09(68156:72155) loc_09(72156:76155)+km_09(72156:76155) loc_09(76156:80155)+km_09(76156:80155) loc_09(80156:84155)+km_09(80156:84155) loc_09(84156:88227)+km_09(84156:88227) loc_09(88228:92227)+km_09(88228:92227) loc_09(92228:96227)+km_09(92228:96227) loc_09(96228:100227)+km_09(96228:100227) loc_09(100228:104227)+km_09(100228:104227) loc_09(104228:108227)+km_09(104228:108227)
119
loc_09(108228:112300)+km_09(108228:112300) loc_09(112301:116300)+km_09(112301:116300) loc_09(116301:120300)+km_09(116301:120300) loc_09(120301:124300)+km_09(120301:124300) loc_09(124301:128329)+km_09(124301:128329) loc_09(128330:132329)+km_09(128330:132329) loc_09(132330:136329)+km_09(132330:136329)]; %ajuste em 2007 total, esquerdo x_07_ajuste_E=x_07m+1*10^-3 %ajuste em 2007 total, direito x_07_ajuste_D=x_07m+2*10^-3 %Ali carril esquerdo %ali E d1 2007 total ali_07_E=data_07(:,4); %ali E d1 2009 total ali_09_E=data_09(:,4); %Ali carril direito %ali D d1 2007 total ali_07_D=data_07(:,5); %ali D d1 2009 total ali_09_D=data_09(:,5); % 2007, total, Esquerdo, teste figure(1) plot(x_07m,ali_07_E, 'b' ) legend ( 'Ali E d1 2007 Mar total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007, total, Esquerdo, teste' ) % 2009, total, Esquerdo, teste figure(2) plot(x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2009 Jun total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2009, total, Esquerdo, teste' ) % 2007 Mar VS 2009, total, Esquerdo figure(3) plot(x_07m,ali_07_E, 'b' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2007 Mar total' , 'Ali E d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' )
120
title( '2007 Mar VS 2009, total, Esquerdo' ) % 2007 Mar VS 2009, total, Direito figure(4) plot(x_07m,ali_07_D, 'b' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2007 Mar total' , 'Ali D d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Mar VS 2009, total, Direito' ) % 2007 Mar VS 2009, total, Esquerdo, ajuste figure(5) plot(x_07_ajuste_E,ali_07_E, 'b' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2007 Mar total ajuste' , 'Ali E d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Mar VS 2009, total, Esquerdo, ajuste' ) % 2007 Mar VS 2009, total, Direito, ajuste figure(6) plot(x_07_ajuste_D,ali_07_D, 'b' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2007 Mar total ajuste' , 'Ali D d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Mar VS 2009, total, Direito, ajuste' ) % calculo do desvio padrao em segmentos de 200 m, E squerda, 05, ajustado n=fix(size(x_07_ajuste_E)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_07_E([n1(i):n2(i)],:)])]; dp_07mar_E=dp'; end % calculo do desvio padrao em segmentos de 200 m, D ireita, 05, ajustado n=fix(size(x_07_ajuste_D)/800); for i=1:n
121
n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_07_D([n1(i):n2(i)],:)])]; dp_07mar_D=dp'; end %Sincronizacao 1km %Parametro D1 % %1 km (1,4000) %posicao 2007 1km km_um_07=data_07(1:4000,1); loc_um_07=loc_07m(1:4000); x_um_07=km_um_07+loc_um_07; %posicao 2009 1km km_um_09=data_09(1:4000,1); loc_um_09=loc_09(1:4000); x_um_09=km_um_09+loc_um_09; %0,5 km (1,2000) %posicao 2007 0,5km km_meio_07=data_07(1:2000,1); loc_meio_07=loc_07m(1:2000); x_meio_07=km_meio_07+loc_meio_07; %ajuste em 2007 para 0,5 km esquerdo x_meio_07_ajuste_E=x_meio_07+1*10^-3 %ajuste em 2007 para 0,5 km direito x_meio_07_ajuste_D=x_meio_07+2*10^-3 %posicao 2009 0,5km km_meio_09=data_09(1:2000,1); loc_meio_09=loc_09(1:2000); x_meio_09=km_meio_09+loc_meio_09; %Ali carril esquerdo %ali E D1 2007 1km ali_um_07_E=data_07(1:4000,4);
122
ali_meio_07_E=data_07(1:2000,4); %ali E D1 2009 1km REF ali_um_09_E=data_09(1:4000,4); ali_meio_09_E=data_09(1:2000,4); % 2007 Mar VS 2009, 1km, Esquerda figure(7) plot(x_um_07,ali_um_07_E, 'b' ,x_um_09,ali_um_09_E, 'g' ) legend ( 'Ali E1 2007 Mar 1km' , 'Ali E1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Mar VS 2009, 1km, Esquerda' ) % 2007 Mar VS 2009, 0,5km, Esquerda figure(8) plot(x_meio_07,ali_meio_07_E, 'b' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2007 Mar 0,5km' , 'Ali E1 2009 0,5km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Mar VS 2009, 0,5km, Esquerda' ) % 2007 Mar VS 2009, 0,5km, c/ ajuste figure(9) plot(x_meio_07_ajuste_E,ali_meio_07_E, 'b' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2007 Mar 0,5 km ajuste' , 'Ali E1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Mar VS 2009, 0,5km, c/ ajuste, Esquerda' ) %Ali Carril Direito %ali D D1 2007 1km ali_um_07_D=data_07(1:4000,5); ali_meio_07_D=data_07(1:2000,5); %ali D D1 2009 1km REF ali_um_09_D=data_09(1:4000,5); ali_meio_09_D=data_09(1:2000,5); % 2007 Mar VS 2009, 1km figure(10) plot(x_um_07,ali_um_07_D, 'b' ,x_um_09,ali_um_09_D, 'g' ) legend ( 'Ali D d1 2007 Mar 1km' , 'Ali D d1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Mar VS 2009, 1km, Direito' ) % 2007 Mar VS 2009, 0,5km figure(11) plot(x_meio_07,ali_meio_07_D, 'b' ,x_meio_09,ali_meio_09_D, 'g' )
123
legend ( 'Ali D d1 2007 Mar 0,5 km' , 'Ali D d1 2009 0,5 km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Mar VS 2009, 0,5km, Direito' ) % 2007 Mar VS 2009, 0,5km, c/ ajuste, direito figure(12) plot(x_meio_07_ajuste_D,ali_meio_07_D, 'b' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2007 Mar 0,5 km ajuste' , 'Ali D d1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Mar VS 2009, 0,5km, c/ ajuste, direito' ) %Universidade do Porto %Faculdade de Engenharia %DEC 2010 %Clara Moura %2009 (REF) and 2007 Junho % troco total %loading data % 2007 data_07=dlmread( 'campanha2007jun.txt' ); % 2009 REF data_09=dlmread( 'campanha2009pronto.txt' ); %posicao 2007 km_07j=data_07(:,1); loc_07j=data_07(:,2)*10^-3; x_07j=[loc_07j(1:4000)+km_07j(1:4000) loc_07j(4001:8000)+km_07j(4001:8000) loc_07j(8001:12000)+km_07j(8001:12000) loc_07j(12001:16000)+km_07j(12001:16000) loc_07j(16001:20000)+km_07j(16001:20000) loc_07j(20001:24000)+km_07j(20001:24000) loc_07j(24001:28000)+km_07j(24001:28000) loc_07j(28001:32000)+km_07j(28001:32000) loc_07j(32001:36000)+km_07j(32001:36000) loc_07j(36001:40000)+km_07j(36001:40000) loc_07j(40001:44000)+km_07j(40001:44000) loc_07j(44001:48000)+km_07j(44001:48000) loc_07j(48001:52000)+km_07j(48001:52000) loc_07j(52001:56000)+km_07j(52001:56000) loc_07j(56001:60000)+km_07j(56001:60000) loc_07j(60001:64000)+km_07j(60001:64000) loc_07j(64001:68000)+km_07j(64001:68000) loc_07j(68001:72000)+km_07j(68001:72000) loc_07j(72001:76000)+km_07j(72001:76000)
124
loc_07j(76001:80000)+km_07j(76001:80000) loc_07j(80001:84000)+km_07j(80001:84000) loc_07j(84001:87941)+km_07j(84001:87941) loc_07j(87942:91941)+km_07j(87942:91941) loc_07j(91942:95941)+km_07j(91942:95941) loc_07j(95942:99941)+km_07j(95942:99941) loc_07j(99942:103941)+km_07j(99942:103941) loc_07j(103942:107941)+km_07j(103942:107941) loc_07j(107942:111941)+km_07j(107942:111941) loc_07j(111942:115941)+km_07j(111942:115941) loc_07j(115942:119941)+km_07j(115942:119941) loc_07j(119942:123941)+km_07j(119942:123941) loc_07j(123942:127941)+km_07j(123942:127941)]; %posicao 2009 km_09=data_09(:,1); loc_09=data_09(:,2)*10^-3; x_09=[loc_09(1:4000)+km_09(1:4000) loc_09(4001:8000)+km_09(4001:8000) loc_09(8001:12000)+km_09(8001:12000) loc_09(12001:16000)+km_09(12001:16000) loc_09(16001:20057)+km_09(16001:20057) loc_09(20058:24057)+km_09(20058:24057) loc_09(24058:28057)+km_09(24058:28057) loc_09(28058:32057)+km_09(28058:32057) loc_09(32058:36057)+km_09(32058:36057) loc_09(36058:40095)+km_09(36058:40095) loc_09(40096:44095)+km_09(40096:44095) loc_09(44096:48095)+km_09(44096:48095) loc_09(48096:52095)+km_09(48096:52095) loc_09(52096:56095)+km_09(52096:56095) loc_09(56096:60155)+km_09(56096:60155) loc_09(60156:64155)+km_09(60156:64155) loc_09(64156:68155)+km_09(64156:68155) loc_09(68156:72155)+km_09(68156:72155) loc_09(72156:76155)+km_09(72156:76155) loc_09(76156:80155)+km_09(76156:80155) loc_09(80156:84155)+km_09(80156:84155) loc_09(84156:88227)+km_09(84156:88227) loc_09(88228:92227)+km_09(88228:92227) loc_09(92228:96227)+km_09(92228:96227) loc_09(96228:100227)+km_09(96228:100227) loc_09(100228:104227)+km_09(100228:104227) loc_09(104228:108227)+km_09(104228:108227) loc_09(108228:112300)+km_09(108228:112300) loc_09(112301:116300)+km_09(112301:116300) loc_09(116301:120300)+km_09(116301:120300) loc_09(120301:124300)+km_09(120301:124300) loc_09(124301:128329)+km_09(124301:128329) loc_09(128330:132329)+km_09(128330:132329) loc_09(132330:136329)+km_09(132330:136329)];
125
%ajuste em 2007 total, esquerdo x_07_ajuste_E=x_07j+1*10^-3; %ajuste em 2007 total, direito x_07_ajuste_D=x_07j+2*10^-3; %Ali carril esquerdo %ali E d1 2007 total ali_07_E=data_07(:,4); %ali E d1 2009 total ali_09_E=data_09(:,4); %Ali carril direito %ali D d1 2007 total ali_07_D=data_07(:,5); %ali D d1 2009 total ali_09_D=data_09(:,5); % 2007, total, Esquerdo, teste figure(1) plot(x_07j,ali_07_E, 'b' ) legend ( 'Ali E d1 2007 Jun total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '% 2007, total, Esquerdo, teste' ) % 2009, total, Esquerdo, teste figure(2) plot(x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2009 Jun total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '% 2009, total, Esquerdo, teste' ) % 2007 Jun VS 2009, total, Esquerdo figure(3) plot(x_07j,ali_07_E, 'b' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2007 Jun total' , 'Ali E d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Jun VS 2009, total, Esquerdo' ) % 2007 Jun VS 2009, total, Direito figure(4) plot(x_07j,ali_07_D, 'b' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2007 Jun total' , 'Ali D d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' )
126
title( '2007 Jun VS 2009, total, Direito' ) % 2007 Jun VS 2009, total, Esquerdo, ajuste figure(5) plot(x_07_ajuste_E,ali_07_E, 'b' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2007 Jun total ajuste' , 'Ali E d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Jun VS 2009, total, Esquerdo, ajuste' ) % 2007 Jun VS 2009, total, Direito, ajuste figure(6) plot(x_07_ajuste_D,ali_07_D, 'b' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2007 Jun total ajuste' , 'Ali D d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Jun VS 2009, total, Direito, ajuste' ) % calculo do desvio padrao em segmentos de 200 m, E squerda, 05, ajustado n=fix(size(x_07_ajuste_E)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_07_E([n1(i):n2(i)],:)])]; dp_07jun_E=dp'; end % calculo do desvio padrao em segmentos de 200 m, D ireita, 05, ajustado n=fix(size(x_07_ajuste_D)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_07_D([n1(i):n2(i)],:)])]; dp_07jun_D=dp'; end
127
%Sincronizacao 1km %Parametro D1 % %1 km (1,4000) %posicao 2007 1km km_um_07=data_07(1:4000,1); loc_um_07=loc_07j(1:4000); x_um_07=km_um_07+loc_um_07; %posicao 2009 1km km_um_09=data_09(1:4000,1); loc_um_09=loc_09(1:4000); x_um_09=km_um_09+loc_um_09; %0,5 km (1,2000) %posicao 2007 0,5km km_meio_07=data_07(1:2000,1); loc_meio_07=loc_07j(1:2000); x_meio_07=km_meio_07+loc_meio_07; %ajuste em 2007 para 0,5 km esquerdo x_meio_07_ajuste_E=x_meio_07+1*10^-3 %ajuste em 2007 para 0,5 km direito x_meio_07_ajuste_D=x_meio_07+2*10^-3 %posicao 2009 0,5km km_meio_09=data_09(1:2000,1); loc_meio_09=loc_09(1:2000); x_meio_09=km_meio_09+loc_meio_09; %Ali carril esquerdo %ali E D1 2007 1km ali_um_07_E=data_07(1:4000,4); ali_meio_07_E=data_07(1:2000,4); %ali E D1 2009 1km REF ali_um_09_E=data_09(1:4000,4); ali_meio_09_E=data_09(1:2000,4); % 2007 Jun VS 2009, 1km, Esquerda figure(7) plot(x_um_07,ali_um_07_E, 'b' ,x_um_09,ali_um_09_E, 'g' ) legend ( 'Ali E1 2007 Jun 1km' , 'Ali E1 2009 1km' )
128
xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Jun VS 2009, 1km, Esquerda' ) % 2007 Jun VS 2009, 0,5km, Esquerda figure(8) plot(x_meio_07,ali_meio_07_E, 'b' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2007 Jun 0,5km' , 'Ali E1 2009 0,5km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Jun VS 2009, 0,5km, Esquerda' ) % 2007 Jun VS 2009, 0,5km, c/ ajuste figure(9) plot(x_meio_07_ajuste_E,ali_meio_07_E, 'b' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2007 Jun 0,5 km ajuste' , 'Ali E1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Jun VS 2009, 0,5km, c/ ajuste, Esquerda' ) %Ali Carril Direito %ali D D1 2007 1km ali_um_07_D=data_07(1:4000,5); ali_meio_07_D=data_07(1:2000,5); %ali D D1 2009 1km REF ali_um_09_D=data_09(1:4000,5); ali_meio_09_D=data_09(1:2000,5); % 2007 Jun VS 2009, 1km figure(10) plot(x_um_07,ali_um_07_D, 'b' ,x_um_09,ali_um_09_D, 'g' ) legend ( 'Ali D d1 2007 Jun 1km' , 'Ali D d1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Jun VS 2009, 1km, Direito' ) % 2007 Jun VS 2009, 0,5km figure(11) plot(x_meio_07,ali_meio_07_D, 'b' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2007 Jun 0,5 km' , 'Ali D d1 2009 0,5 km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Jun VS 2009, 0,5km, Direito' ) % 2007 Jun VS 2009, 0,5km, c/ ajuste, direito figure(12) plot(x_meio_07_ajuste_D,ali_meio_07_D, 'b' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2007 Jun 0,5 km ajuste' , 'Ali D d1 2009 0,5 km ajuste' ) xlabel( 'Posicao' )
129
ylabel( 'Alinhamento' ) title( '2007 Jun VS 2009, 0,5km, c/ ajuste, direito' ) %Universidade do Porto %Faculdade de Engenharia %DEC 2010 %Clara Moura %2009 (REF) and 2007 Dezembro % troco total %loading data % 2007 data_07=dlmread( 'campanha2007dez.txt' ); % 2009 REF data_09=dlmread( 'campanha2009pronto.txt' ); %posicao 2007 km_07=data_07(:,1); loc_07=data_07(:,2)*10^-3; x_07=[loc_07(1:4000)+km_07(1:4000) loc_07(4001:8000)+km_07(4001:8000) loc_07(8001:12057)+km_07(8001:12057) loc_07(12058:16057)+km_07(12058:16057) loc_07(16058:20057)+km_07(16058:20057) loc_07(20058:24057)+km_07(20058:24057) loc_07(24058:28057)+km_07(24058:28057) loc_07(28058:32092)+km_07(28058:32092) loc_07(32093:36092)+km_07(32093:36092) loc_07(36093:40092)+km_07(36093:40092) loc_07(40093:44092)+km_07(40093:44092) loc_07(44093:48124)+km_07(44093:48124) loc_07(48125:52124)+km_07(48125:52124) loc_07(52125:56124)+km_07(52125:56124) loc_07(56125:60124)+km_07(56125:60124) loc_07(60125:64124)+km_07(60125:64124) loc_07(64125:68124)+km_07(64125:68124) loc_07(68125:72124)+km_07(68125:72124) loc_07(72125:76124)+km_07(72125:76124) loc_07(76125:80124)+km_07(76125:80124) loc_07(80125:84124)+km_07(80125:84124) loc_07(84125:88124)+km_07(84125:88124) loc_07(88125:92194)+km_07(88125:92194) loc_07(92195:96194)+km_07(92195:96194) loc_07(96195:100194)+km_07(96195:100194) loc_07(100195:104194)+km_07(100195:104194) loc_07(104195:108117)+km_07(104195:108117) loc_07(108118:112262)+km_07(108118:112262) loc_07(112263:116243)+km_07(112263:116243)
130
loc_07(116244:120243)+km_07(116244:120243) loc_07(120244:124243)+km_07(120244:124243) loc_07(124244:128281)+km_07(124244:128281)]; %posicao 2009 km_09=data_09(:,1); loc_09=data_09(:,2)*10^-3; x_09=[loc_09(1:4000)+km_09(1:4000) loc_09(4001:8000)+km_09(4001:8000) loc_09(8001:12000)+km_09(8001:12000) loc_09(12001:16000)+km_09(12001:16000) loc_09(16001:20057)+km_09(16001:20057) loc_09(20058:24057)+km_09(20058:24057) loc_09(24058:28057)+km_09(24058:28057) loc_09(28058:32057)+km_09(28058:32057) loc_09(32058:36057)+km_09(32058:36057) loc_09(36058:40095)+km_09(36058:40095) loc_09(40096:44095)+km_09(40096:44095) loc_09(44096:48095)+km_09(44096:48095) loc_09(48096:52095)+km_09(48096:52095) loc_09(52096:56095)+km_09(52096:56095) loc_09(56096:60155)+km_09(56096:60155) loc_09(60156:64155)+km_09(60156:64155) loc_09(64156:68155)+km_09(64156:68155) loc_09(68156:72155)+km_09(68156:72155) loc_09(72156:76155)+km_09(72156:76155) loc_09(76156:80155)+km_09(76156:80155) loc_09(80156:84155)+km_09(80156:84155) loc_09(84156:88227)+km_09(84156:88227) loc_09(88228:92227)+km_09(88228:92227) loc_09(92228:96227)+km_09(92228:96227) loc_09(96228:100227)+km_09(96228:100227) loc_09(100228:104227)+km_09(100228:104227) loc_09(104228:108227)+km_09(104228:108227) loc_09(108228:112300)+km_09(108228:112300) loc_09(112301:116300)+km_09(112301:116300) loc_09(116301:120300)+km_09(116301:120300) loc_09(120301:124300)+km_09(120301:124300) loc_09(124301:128329)+km_09(124301:128329) loc_09(128330:132329)+km_09(128330:132329) loc_09(132330:136329)+km_09(132330:136329)]; %ajuste em 2007 total, esquerdo x_07_ajuste_E=x_07-5.1667*10^-3; %ajuste em 2007 total, direito x_07_ajuste_D=x_07-5.3750*10^-3; %Ali carril esquerdo
131
%ali E d1 2007 total ali_07_E=data_07(:,4); %ali E d1 2009 total ali_09_E=data_09(:,4); %Ali carril direito %ali D d1 2007 total ali_07_D=data_07(:,5); %ali D d1 2009 total ali_09_D=data_09(:,5); % 2007, total, Esquerdo, teste figure(1) plot(x_07,ali_07_E, 'b' ) legend ( 'Ali E d1 2007 Dez total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '% 2007, total, Esquerdo, teste' ) % 2009, total, Esquerdo, teste figure(2) plot(x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2009 Jun total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '% 2009, total, Esquerdo, teste' ) % 2007 Dez VS 2009, total, Esquerdo figure(3) plot(x_07,ali_07_E, 'b' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2007 Dez total' , 'Ali E d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Dez VS 2009, total, Esquerdo' ) % 2007 Dez VS 2009, total, Direito figure(4) plot(x_07,ali_07_D, 'b' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2007 Dez total' , 'Ali D d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Dez VS 2009, total, Direito' ) % 2007 Dez VS 2009, total, Esquerdo, ajuste figure(5) plot(x_07_ajuste_E,ali_07_E, 'b' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2007 Dez total ajuste' , 'Ali E d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Dez VS 2009, total, Esquerdo, ajuste' )
132
% 2007 Dez VS 2009, total, Direito, ajuste figure(6) plot(x_07_ajuste_D,ali_07_D, 'b' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2007 Dez total ajuste' , 'Ali D d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Dez VS 2009, total, Direito, ajuste' ) % calculo do desvio padrao em segmentos de 200 m, E squerda, 05, ajustado n=fix(size(x_07_ajuste_E)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_07_E([n1(i):n2(i)],:)])]; dp_07dez_E=dp'; end % calculo do desvio padrao em segmentos de 200 m, D ireita, 05, ajustado n=fix(size(x_07_ajuste_D)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_07_D([n1(i):n2(i)],:)])]; dp_07dez_D=dp'; end %Sincronizacao 1km %Parametro D1 % %1 km (1,4000) %posicao 2007 1km
133
km_um_07=data_07(1:4000,1); loc_um_07=loc_07(1:4000); x_um_07=km_um_07+loc_um_07; %posicao 2009 1km km_um_09=data_09(1:4000,1); loc_um_09=loc_09(1:4000); x_um_09=km_um_09+loc_um_09; %0,5 km (1,2000) %posicao 2007 0,5km km_meio_07=data_07(1:2000,1); loc_meio_07=loc_07(1:2000); x_meio_07=km_meio_07+loc_meio_07; %ajuste em 2007 para 0,5 km esquerdo x_meio_07_ajuste_E=x_meio_07-5.1667*10^-3 %ajuste em 2007 para 0,5 km direito x_meio_07_ajuste_D=x_meio_07-5.3750*10^-3 %posicao 2009 0,5km km_meio_09=data_09(1:2000,1); loc_meio_09=loc_09(1:2000); x_meio_09=km_meio_09+loc_meio_09; %Ali carril esquerdo %ali E D1 2007 1km ali_um_07_E=data_07(1:4000,4); ali_meio_07_E=data_07(1:2000,4); %ali E D1 2009 1km REF ali_um_09_E=data_09(1:4000,4); ali_meio_09_E=data_09(1:2000,4); % 2007 Dez VS 2009, 1km, Esquerda figure(7) plot(x_um_07,ali_um_07_E, 'b' ,x_um_09,ali_um_09_E, 'g' ) legend ( 'Ali E1 2007 Dez 1km' , 'Ali E1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Dez VS 2009, 1km, Esquerda' ) % 2007 Dez VS 2009, 0,5km, Esquerda figure(8) plot(x_meio_07,ali_meio_07_E, 'b' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2007 Dez 0,5km' , 'Ali E1 2009 0,5km' )
134
xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Dez VS 2009, 0,5km, Esquerda' ) % 2007 Dez VS 2009, 0,5km, c/ ajuste figure(9) plot(x_meio_07_ajuste_E,ali_meio_07_E, 'b' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2007 Dez 0,5 km ajuste' , 'Ali E1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Dez VS 2009, 0,5km, c/ ajuste, Esquerda' ) %Ali Carril Direito %ali D D1 2007 1km ali_um_07_D=data_07(1:4000,5); ali_meio_07_D=data_07(1:2000,5); %ali D D1 2009 1km REF ali_um_09_D=data_09(1:4000,5); ali_meio_09_D=data_09(1:2000,5); % % 2007 Dez VS 2009, 1km figure(10) plot(x_um_07,ali_um_07_D, 'b' ,x_um_09,ali_um_09_D, 'g' ) legend ( 'Ali D d1 2007 Dez 1km' , 'Ali D d1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Dez VS 2009, 1km, Direito' ) % 2007 Dez VS 2009, 0,5km figure(11) plot(x_meio_07,ali_meio_07_D, 'b' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2007 Dez 0,5 km' , 'Ali D d1 2009 0,5 km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Dez VS 2009, 0,5km, Direito' ) % 2007 Dez VS 2009, 0,5km, c/ ajuste, direito figure(12) plot(x_meio_07_ajuste_D,ali_meio_07_D, 'b' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2007 Dez 0,5 km ajuste' , 'Ali D d1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2007 Dez VS 2009, 0,5km, c/ ajuste, direito' )
135
%Universidade do Porto %Faculdade de Engenharia %DEC 2010 %Clara Moura
%2009 (REF) and 2006 Dez % troco total %loading data % 2006 data_06=dlmread( 'campanha2006dez.txt' ); % 2009 REF data_09=dlmread( 'campanha2009pronto.txt' ); %posicao 2006 km_06=data_06(:,1); loc_06=data_06(:,2)*10^-3; x_06=[loc_05(1:4000)+km_06(1:4000) loc_06(4001:8000)+km_06(4001:8000) loc_06(8001:12000)+km_06(8001:12000) loc_06(12001:16000)+km_06(12001:16000) loc_06(16001:20000)+km_06(16001:20000) loc_06(20001:23987)+km_06(20001:23987) loc_06(23988:27987)+km_06(23988:27987) loc_06(27988:31987)+km_06(27988:31987) loc_06(31988:35987)+km_06(31988:35987) loc_06(35988:39987)+km_06(35988:39987) loc_06(39988:43987)+km_06(39988:43987) loc_06(43988:47987)+km_06(43988:47987) loc_06(47988:51987)+km_06(47988:51987) loc_06(51988:55987)+km_06(51988:55987) loc_06(55988:59987)+km_06(55988:59987) loc_06(59988:63965)+km_06(59988:63965) loc_06(63966:67965)+km_06(63966:67965) loc_06(67966:71965)+km_06(67966:71965) loc_06(71966:75965)+km_06(71966:75965) loc_06(75966:79965)+km_06(75966:79965) loc_06(79966:83965)+km_06(79966:83965) loc_06(83966:87936)+km_06(83966:87936) loc_06(87937:91936)+km_06(87937:91936) loc_06(91937:95936)+km_06(91937:95936) loc_06(95937:99936)+km_06(95937:99936) loc_06(99937:103936)+km_06(99937:103936) loc_06(103937:107936)+km_06(103937:107936) loc_06(107937:111939)+km_06(107937:111939) loc_06(111940:115939)+km_06(111940:115939) loc_06(115940:119939)+km_06(115940:119939) loc_06(119940:123939)+km_06(119940:123939) loc_06(123940:127915)+km_06(123940:127915)];
136
%posicao 2009 km_09=data_09(:,1); loc_09=data_09(:,2)*10^-3; x_09=[loc_09(1:4000)+km_09(1:4000) loc_09(4001:8000)+km_09(4001:8000) loc_09(8001:12000)+km_09(8001:12000) loc_09(12001:16000)+km_09(12001:16000) loc_09(16001:20057)+km_09(16001:20057) loc_09(20058:24057)+km_09(20058:24057) loc_09(24058:28057)+km_09(24058:28057) loc_09(28058:32057)+km_09(28058:32057) loc_09(32058:36057)+km_09(32058:36057) loc_09(36058:40095)+km_09(36058:40095) loc_09(40096:44095)+km_09(40096:44095) loc_09(44096:48095)+km_09(44096:48095) loc_09(48096:52095)+km_09(48096:52095) loc_09(52096:56095)+km_09(52096:56095) loc_09(56096:60155)+km_09(56096:60155) loc_09(60156:64155)+km_09(60156:64155) loc_09(64156:68155)+km_09(64156:68155) loc_09(68156:72155)+km_09(68156:72155) loc_09(72156:76155)+km_09(72156:76155) loc_09(76156:80155)+km_09(76156:80155) loc_09(80156:84155)+km_09(80156:84155) loc_09(84156:88227)+km_09(84156:88227) loc_09(88228:92227)+km_09(88228:92227) loc_09(92228:96227)+km_09(92228:96227) loc_09(96228:100227)+km_09(96228:100227) loc_09(100228:104227)+km_09(100228:104227) loc_09(104228:108227)+km_09(104228:108227) loc_09(108228:112300)+km_09(108228:112300) loc_09(112301:116300)+km_09(112301:116300) loc_09(116301:120300)+km_09(116301:120300) loc_09(120301:124300)+km_09(120301:124300) loc_09(124301:128329)+km_09(124301:128329) loc_09(128330:132329)+km_09(128330:132329) loc_09(132330:136329)+km_09(132330:136329)]; %ajuste em 2006 dezembro total, esquerdo x_06_ajuste_E=x_06-6.75*10^-3 %ajuste em 2006 dezembro total, direito x_06_ajuste_D=x_06+23*10^-3 %Ali carril esquerdo %ali E d1 2006 dezembro total ali_06_E=data_06(:,4); %ali E d1 2009 total ali_09_E=data_09(:,4);
137
%Ali carril direito %ali D d1 2006 dezembro total ali_06_D=data_06(:,5); %ali D d1 2009 total ali_09_D=data_09(:,5); % 2006 dezembro, total, Esquerdo, teste figure(1) plot(x_06,ali_06_E, 'y' ) legend ( 'Ali E d1 2006 dezembro total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro, total, Esquerdo, teste' ) % 2009, total, Esquerdo, teste figure(2) plot(x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2009 Jun total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2009, total, Esquerdo, teste' ) % 2006 dezembro VS 2009, total, Esquerdo figure(3) plot(x_06,ali_06_E, 'y' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2006 dezembro total' , 'Ali E d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, total, Esquerdo' ) % 2006 dezembro VS 2009, total, Direito figure(4) plot(x_06,ali_06_D, 'y' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2006 dezembro total' , 'Ali D d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, total, Direito' ) % 2006 dezembro VS 2009, total, Esquerdo, ajuste figure(5) plot(x_06_ajuste_E,ali_06_E, 'y' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2008 total ajuste' , 'Ali E d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, total, Esquerdo, ajuste' ) % 2006 dezembro VS 2009, total, Direito, ajuste figure(6) plot(x_06_ajuste_D,ali_06_D, 'y' ,x_09,ali_09_D, 'g' )
138
legend ( 'Ali D d1 2006 dezembro total ajuste' , 'Ali D d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, total, Direito, ajuste' ) % calculo do desvio padrao em segmentos de 200 m, E squerda, 05, ajustado n=fix(size(x_06_ajuste_E)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_06_E([n1(i):n2(i)],:)])]; dp_06_E=dp'; end % calculo do desvio padrao em segmentos de 200 m, D ireita, 05, ajustado n=fix(size(x_06_ajuste_D)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_06_D([n1(i):n2(i)],:)])]; dp_06_D=dp'; end %Sincronizacao 1km %Parametro D1 % %1 km (1,4000) %posicao 2006 dezembro 1km km_um_06=data_06(1:4000,1); loc_um_06=loc_06(1:4000); x_um_06=km_um_06+loc_um_06; %posicao 2009 1km km_um_09=data_09(1:4000,1);
139
loc_um_09=loc_09(1:4000); x_um_09=km_um_09+loc_um_09; %0,5 km (1,2000) %posicao 2006 dezembro 0,5km km_meio_06=data_06(1:2000,1); loc_meio_06=loc_06(1:2000); x_meio_06=km_meio_06+loc_meio_06; %ajuste em 2006 dezembro para 0,5 km esquerdo x_meio_06_ajuste_E=x_meio_06-6.75*10^-3 %ajuste em 2006 dezembro para 0,5 km direito x_meio_06_ajuste_D=x_meio_06+38*10^-3 %posicao 2009 0,5km km_meio_09=data_09(1:2000,1); loc_meio_09=loc_09(1:2000); x_meio_09=km_meio_09+loc_meio_09; %Ali carril esquerdo %ali E D1 2006 dezembro 1km ali_um_06_E=data_06(1:4000,4); ali_meio_06_E=data_06(1:2000,4); %ali E D1 2009 1km REF ali_um_09_E=data_09(1:4000,4); ali_meio_09_E=data_09(1:2000,4); % 2006 dezembro VS 2009, 1km, Esquerda figure(7) plot(x_um_06,ali_um_06_E, 'y' ,x_um_09,ali_um_09_E, 'g' ) legend ( 'Ali E1 2006 dezembro 1km' , 'Ali E1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, 1km, Esquerda' ) % 2006 dezembro VS 2009, 0,5km, Esquerda figure(8) plot(x_meio_06,ali_meio_06_E, 'y' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2006 dezembro 0,5km' , 'Ali E1 2009 0,5km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, 0,5km, Esquerda' ) % 2006 dezembro VS 2009, 0,5km, c/ ajuste
140
figure(9) plot(x_meio_06_ajuste_E,ali_meio_06_E, 'y' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2006 dezembro 0,5 km ajuste' , 'Ali E1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, 0,5km, c/ ajuste, Esquerda' ) %Ali Carril Direito %ali D D1 2006 dezembro 1km ali_um_06_D=data_06(1:4000,5); ali_meio_06_D=data_06(1:2000,5); %ali D D1 2009 1km REF ali_um_09_D=data_09(1:4000,5); ali_meio_09_D=data_09(1:2000,5); % 2006 dezembro VS 2009, 1km figure(10) plot(x_um_06,ali_um_06_D, 'y' ,x_um_09,ali_um_09_D, 'g' ) legend ( 'Ali D d1 2006 dezembro 1km' , 'Ali D d1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, 1km, Direito' ) % 2006 dezembro VS 2009, 0,5km figure(11) plot(x_meio_06,ali_meio_06_D, 'y' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2006 dezembro 0,5 km' , 'Ali D d1 2009 0,5 km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, 0,5km, Direito' ) % 2006 dezembro VS 2009, 0,5km, c/ ajuste, direito figure(12) plot(x_meio_06_ajuste_D,ali_meio_06_D, 'y' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2006 dezembro 0,5 km ajuste' , 'Ali D d1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, 0,5km, c/ ajuste, direito' )
141
%2009 (REF) and 2006 Dez % troco total %loading data % 2006 data_06=dlmread( 'campanha2006dez.txt' ); % 2009 REF data_09=dlmread( 'campanha2009pronto.txt' ); %posicao 2006 km_06=data_06(:,1); loc_06=data_06(:,2)*10^-3; x_06=[loc_05(1:4000)+km_06(1:4000) loc_06(4001:8000)+km_06(4001:8000) loc_06(8001:12000)+km_06(8001:12000) loc_06(12001:16000)+km_06(12001:16000) loc_06(16001:20000)+km_06(16001:20000) loc_06(20001:23987)+km_06(20001:23987) loc_06(23988:27987)+km_06(23988:27987) loc_06(27988:31987)+km_06(27988:31987) loc_06(31988:35987)+km_06(31988:35987) loc_06(35988:39987)+km_06(35988:39987) loc_06(39988:43987)+km_06(39988:43987) loc_06(43988:47987)+km_06(43988:47987) loc_06(47988:51987)+km_06(47988:51987) loc_06(51988:55987)+km_06(51988:55987) loc_06(55988:59987)+km_06(55988:59987) loc_06(59988:63965)+km_06(59988:63965) loc_06(63966:67965)+km_06(63966:67965) loc_06(67966:71965)+km_06(67966:71965) loc_06(71966:75965)+km_06(71966:75965) loc_06(75966:79965)+km_06(75966:79965) loc_06(79966:83965)+km_06(79966:83965) loc_06(83966:87936)+km_06(83966:87936) loc_06(87937:91936)+km_06(87937:91936) loc_06(91937:95936)+km_06(91937:95936) loc_06(95937:99936)+km_06(95937:99936) loc_06(99937:103936)+km_06(99937:103936) loc_06(103937:107936)+km_06(103937:107936) loc_06(107937:111939)+km_06(107937:111939) loc_06(111940:115939)+km_06(111940:115939) loc_06(115940:119939)+km_06(115940:119939) loc_06(119940:123939)+km_06(119940:123939) loc_06(123940:127915)+km_06(123940:127915)]; %posicao 2009 km_09=data_09(:,1); loc_09=data_09(:,2)*10^-3;
142
x_09=[loc_09(1:4000)+km_09(1:4000) loc_09(4001:8000)+km_09(4001:8000) loc_09(8001:12000)+km_09(8001:12000) loc_09(12001:16000)+km_09(12001:16000) loc_09(16001:20057)+km_09(16001:20057) loc_09(20058:24057)+km_09(20058:24057) loc_09(24058:28057)+km_09(24058:28057) loc_09(28058:32057)+km_09(28058:32057) loc_09(32058:36057)+km_09(32058:36057) loc_09(36058:40095)+km_09(36058:40095) loc_09(40096:44095)+km_09(40096:44095) loc_09(44096:48095)+km_09(44096:48095) loc_09(48096:52095)+km_09(48096:52095) loc_09(52096:56095)+km_09(52096:56095) loc_09(56096:60155)+km_09(56096:60155) loc_09(60156:64155)+km_09(60156:64155) loc_09(64156:68155)+km_09(64156:68155) loc_09(68156:72155)+km_09(68156:72155) loc_09(72156:76155)+km_09(72156:76155) loc_09(76156:80155)+km_09(76156:80155) loc_09(80156:84155)+km_09(80156:84155) loc_09(84156:88227)+km_09(84156:88227) loc_09(88228:92227)+km_09(88228:92227) loc_09(92228:96227)+km_09(92228:96227) loc_09(96228:100227)+km_09(96228:100227) loc_09(100228:104227)+km_09(100228:104227) loc_09(104228:108227)+km_09(104228:108227) loc_09(108228:112300)+km_09(108228:112300) loc_09(112301:116300)+km_09(112301:116300) loc_09(116301:120300)+km_09(116301:120300) loc_09(120301:124300)+km_09(120301:124300) loc_09(124301:128329)+km_09(124301:128329) loc_09(128330:132329)+km_09(128330:132329) loc_09(132330:136329)+km_09(132330:136329)]; %ajuste em 2006 dezembro total, esquerdo x_06_ajuste_E=x_06-6.75*10^-3 %ajuste em 2006 dezembro total, direito x_06_ajuste_D=x_06+23*10^-3 %Ali carril esquerdo %ali E d1 2006 dezembro total ali_06_E=data_06(:,4); %ali E d1 2009 total ali_09_E=data_09(:,4); %Ali carril direito %ali D d1 2006 dezembro total
143
ali_06_D=data_06(:,5); %ali D d1 2009 total ali_09_D=data_09(:,5); % 2006 dezembro, total, Esquerdo, teste figure(1) plot(x_06,ali_06_E, 'y' ) legend ( 'Ali E d1 2006 dezembro total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro, total, Esquerdo, teste' ) % 2009, total, Esquerdo, teste figure(2) plot(x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2009 Jun total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2009, total, Esquerdo, teste' ) % 2006 dezembro VS 2009, total, Esquerdo figure(3) plot(x_06,ali_06_E, 'y' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2006 dezembro total' , 'Ali E d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, total, Esquerdo' ) % 2006 dezembro VS 2009, total, Direito figure(4) plot(x_06,ali_06_D, 'y' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2006 dezembro total' , 'Ali D d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, total, Direito' ) % 2006 dezembro VS 2009, total, Esquerdo, ajuste figure(5) plot(x_06_ajuste_E,ali_06_E, 'y' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2008 total ajuste' , 'Ali E d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, total, Esquerdo, ajuste' ) % 2006 dezembro VS 2009, total, Direito, ajuste figure(6) plot(x_06_ajuste_D,ali_06_D, 'y' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2006 dezembro total ajuste' , 'Ali D d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, total, Direito, ajuste' )
144
% calculo do desvio padrao em segmentos de 200 m, E squerda, 05, ajustado n=fix(size(x_06_ajuste_E)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_06_E([n1(i):n2(i)],:)])]; dp_06_E=dp'; end % calculo do desvio padrao em segmentos de 200 m, D ireita, 05, ajustado n=fix(size(x_06_ajuste_D)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_06_D([n1(i):n2(i)],:)])]; dp_06_D=dp'; end %Sincronizacao 1km %Parametro D1 % %1 km (1,4000) %posicao 2006 dezembro 1km km_um_06=data_06(1:4000,1); loc_um_06=loc_06(1:4000); x_um_06=km_um_06+loc_um_06; %posicao 2009 1km km_um_09=data_09(1:4000,1); loc_um_09=loc_09(1:4000); x_um_09=km_um_09+loc_um_09;
145
%0,5 km (1,2000) %posicao 2006 dezembro 0,5km km_meio_06=data_06(1:2000,1); loc_meio_06=loc_06(1:2000); x_meio_06=km_meio_06+loc_meio_06; %ajuste em 2006 dezembro para 0,5 km esquerdo x_meio_06_ajuste_E=x_meio_06-6.75*10^-3 %ajuste em 2006 dezembro para 0,5 km direito x_meio_06_ajuste_D=x_meio_06+38*10^-3 %posicao 2009 0,5km km_meio_09=data_09(1:2000,1); loc_meio_09=loc_09(1:2000); x_meio_09=km_meio_09+loc_meio_09; %Ali carril esquerdo %ali E D1 2006 dezembro 1km ali_um_06_E=data_06(1:4000,4); ali_meio_06_E=data_06(1:2000,4); %ali E D1 2009 1km REF ali_um_09_E=data_09(1:4000,4); ali_meio_09_E=data_09(1:2000,4); % 2006 dezembro VS 2009, 1km, Esquerda figure(7) plot(x_um_06,ali_um_06_E, 'y' ,x_um_09,ali_um_09_E, 'g' ) legend ( 'Ali E1 2006 dezembro 1km' , 'Ali E1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, 1km, Esquerda' ) % 2006 dezembro VS 2009, 0,5km, Esquerda figure(8) plot(x_meio_06,ali_meio_06_E, 'y' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2006 dezembro 0,5km' , 'Ali E1 2009 0,5km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, 0,5km, Esquerda' ) % 2006 dezembro VS 2009, 0,5km, c/ ajuste figure(9) plot(x_meio_06_ajuste_E,ali_meio_06_E, 'y' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2006 dezembro 0,5 km ajuste' , 'Ali E1 2009 0,5 km ajuste' )
146
xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, 0,5km, c/ ajuste, Esquerda' ) %Ali Carril Direito %ali D D1 2006 dezembro 1km ali_um_06_D=data_06(1:4000,5); ali_meio_06_D=data_06(1:2000,5); %ali D D1 2009 1km REF ali_um_09_D=data_09(1:4000,5); ali_meio_09_D=data_09(1:2000,5); % 2006 dezembro VS 2009, 1km figure(10) plot(x_um_06,ali_um_06_D, 'y' ,x_um_09,ali_um_09_D, 'g' ) legend ( 'Ali D d1 2006 dezembro 1km' , 'Ali D d1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, 1km, Direito' ) % 2006 dezembro VS 2009, 0,5km figure(11) plot(x_meio_06,ali_meio_06_D, 'y' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2006 dezembro 0,5 km' , 'Ali D d1 2009 0,5 km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, 0,5km, Direito' ) % 2006 dezembro VS 2009, 0,5km, c/ ajuste, direito figure(12) plot(x_meio_06_ajuste_D,ali_meio_06_D, 'y' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2006 dezembro 0,5 km ajuste' , 'Ali D d1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2006 dezembro VS 2009, 0,5km, c/ ajuste, direito' )
%Universidade do Porto %Faculdade de Engenharia %DEC 2010 %Clara Moura
%2009 (REF) and 2005 Marco % troco total
147
%loading data % 2005 data_05=dlmread( 'campanha2005mar.txt' ); % 2009 REF data_09=dlmread( 'campanha2009pronto.txt' ); %posicao 2005 km_05=data_05(:,1); loc_05=data_05(:,2)*10^-3; x_05=[loc_05(1:4000)+km_05(1:4000) loc_05(4001:8000)+km_05(4001:8000) loc_05(8001:12081)+km_05(8001:12081) loc_05(12082:16081)+km_05(12082:16081) loc_05(16082:20081)+km_05(16082:20081) loc_05(20082:24081)+km_05(20082:24081) loc_05(24082:28081)+km_05(24082:28081) loc_05(28082:32081)+km_05(28082:32081) loc_05(32082:36081)+km_05(32082:36081) loc_05(36082:40144)+km_05(36082:40144) loc_05(40145:44144)+km_05(40145:44144) loc_05(44145:48144)+km_05(44145:48144) loc_05(48145:52144)+km_05(48145:52144) loc_05(52145:56144)+km_05(52145:56144) loc_05(56145:60207)+km_05(56145:60207) loc_05(60208:64207)+km_05(60208:64207) loc_05(64208:68207)+km_05(64208:68207) loc_05(68208:72207)+km_05(68208:72207) loc_05(72208:76207)+km_05(72208:76207) loc_05(76208:80207)+km_05(76208:80207) loc_05(80208:84207)+km_05(80208:84207) loc_05(84208:88294)+km_05(84208:88294) loc_05(88295:92294)+km_05(88295:92294) loc_05(92295:96294)+km_05(92295:96294) loc_05(96295:100294)+km_05(96295:100294) loc_05(100295:104294)+km_05(100295:104294) loc_05(104295:108294)+km_05(104295:108294) loc_05(108295:112384)+km_05(108295:112384) loc_05(112385:116384)+km_05(112385:116384) loc_05(116385:120393)+km_05(116385:120393) loc_05(120394:124393)+km_05(120394:124393) loc_05(124394:128393)+km_05(124394:128393)]; %posicao 2009 km_09=data_09(:,1); loc_09=data_09(:,2)*10^-3; x_09=[loc_09(1:4000)+km_09(1:4000) loc_09(4001:8000)+km_09(4001:8000) loc_09(8001:12000)+km_09(8001:12000) loc_09(12001:16000)+km_09(12001:16000) loc_09(16001:20057)+km_09(16001:20057)
148
loc_09(20058:24057)+km_09(20058:24057) loc_09(24058:28057)+km_09(24058:28057) loc_09(28058:32057)+km_09(28058:32057) loc_09(32058:36057)+km_09(32058:36057) loc_09(36058:40095)+km_09(36058:40095) loc_09(40096:44095)+km_09(40096:44095) loc_09(44096:48095)+km_09(44096:48095) loc_09(48096:52095)+km_09(48096:52095) loc_09(52096:56095)+km_09(52096:56095) loc_09(56096:60155)+km_09(56096:60155) loc_09(60156:64155)+km_09(60156:64155) loc_09(64156:68155)+km_09(64156:68155) loc_09(68156:72155)+km_09(68156:72155) loc_09(72156:76155)+km_09(72156:76155) loc_09(76156:80155)+km_09(76156:80155) loc_09(80156:84155)+km_09(80156:84155) loc_09(84156:88227)+km_09(84156:88227) loc_09(88228:92227)+km_09(88228:92227) loc_09(92228:96227)+km_09(92228:96227) loc_09(96228:100227)+km_09(96228:100227) loc_09(100228:104227)+km_09(100228:104227) loc_09(104228:108227)+km_09(104228:108227) loc_09(108228:112300)+km_09(108228:112300) loc_09(112301:116300)+km_09(112301:116300) loc_09(116301:120300)+km_09(116301:120300) loc_09(120301:124300)+km_09(120301:124300) loc_09(124301:128329)+km_09(124301:128329) loc_09(128330:132329)+km_09(128330:132329) loc_09(132330:136329)+km_09(132330:136329)]; %ajuste em 2005 total, esquerdo x_05_ajuste_E=x_05+25*10^-3; %ajuste em 2005 total, direito x_05_ajuste_D=x_05-13.5*10^-3; %Ali carril esquerdo %ali E d1 2005 total ali_05_E=data_05(:,4); %ali E d1 2009 total ali_09_E=data_09(:,4); %Ali carril direito %ali D d1 2005 total ali_05_D=data_05(:,5); %ali D d1 2009 total ali_09_D=data_09(:,5); % 2005, total, Esquerdo, teste
149
figure(1) plot(x_05,ali_05_E, 'm' ) legend ( 'Ali E d1 2005 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2005, total, Esquerdo, teste' ) % 2009, total, Esquerdo, teste figure(2) plot(x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2009 Jun total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2009, total, Esquerdo, teste' ) % 2005 VS 2009, total, Esquerdo figure(3) plot(x_05,ali_05_E, 'm' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2005 total' , 'Ali E d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2005 VS 2009, total, Esquerdo' ) % 2005 VS 2009, total, Direito figure(4) plot(x_05,ali_05_D, 'm' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2005 total' , 'Ali D d1 2009 total' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2005 VS 2009, total, Direito' ) % 2005 VS 2009, total, Esquerdo, ajuste figure(5) plot(x_05_ajuste_E,ali_05_E, 'm' ,x_09,ali_09_E, 'g' ) legend ( 'Ali E d1 2005 total ajuste' , 'Ali E d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2005 VS 2009, total, Esquerdo, ajuste' ) % 2005 VS 2009, total, Direito, ajuste figure(6) plot(x_05_ajuste_D,ali_05_D, 'm' ,x_09,ali_09_D, 'g' ) legend ( 'Ali D d1 2005 total ajuste' , 'Ali D d1 2009 total ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2005 VS 2009, total, Direito, ajuste' ) % calculo do desvio padrao em segmentos de 200 m, E squerda, 05, ajustado n=fix(size(x_05_ajuste_E)/800);
150
for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_05_E([n1(i):n2(i)],:)])]; dp_05_E=dp'; end % calculo do desvio padrao em segmentos de 200 m, D ireita, 05, ajustado n=fix(size(x_05_ajuste_D)/800); for i=1:n n1(i)=1+800*(i-1); n2(i)=800*i; dp(i)=[std([ali_05_D([n1(i):n2(i)],:)])]; dp_05_D=dp'; end %Sincronizacao 1km %Parametro D1 % %1 km (1,4000) %posicao 2005 1km km_um_05=data_05(1:4000,1); loc_um_05=loc_05(1:4000); x_um_05=km_um_05+loc_um_05; %posicao 2009 1km km_um_09=data_09(1:4000,1); loc_um_09=loc_09(1:4000); x_um_09=km_um_09+loc_um_09; %0,5 km (1,2000) %posicao 2005 0,5km
151
km_meio_05=data_05(1:2000,1); loc_meio_05=loc_05(1:2000); x_meio_05=km_meio_05+loc_meio_05; %ajuste em 2005 para 0,5 km esquerdo x_meio_05_ajuste_E=x_meio_05+25*10^-3; %ajuste em 2005 para 0,5 km direito x_meio_05_ajuste_D=x_meio_05-13.5*10^-3; %posicao 2009 0,5km km_meio_09=data_09(1:2000,1); loc_meio_09=loc_09(1:2000); x_meio_09=km_meio_09+loc_meio_09; %Ali carril esquerdo %ali E D1 2005 1km ali_um_05_E=data_05(1:4000,4); ali_meio_05_E=data_05(1:2000,4); %ali E D1 2009 1km REF ali_um_09_E=data_09(1:4000,4); ali_meio_09_E=data_09(1:2000,4); % 2005 VS 2009, 1km, Esquerda figure(7) plot(x_um_05,ali_um_05_E, 'm' ,x_um_09,ali_um_09_E, 'g' ) legend ( 'Ali E1 2005 1km' , 'Ali E1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2005 VS 2009, 1km, Esquerda' ) % 2005 VS 2009, 0,5km, Esquerda figure(8) plot(x_meio_05,ali_meio_05_E, 'm' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2005 0,5km' , 'Ali E1 2009 0,5km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2005 VS 2009, 0,5km, Esquerda' ) % 2005 VS 2009, 0,5km, c/ ajuste figure(9) plot(x_meio_05_ajuste_E,ali_meio_05_E, 'm' ,x_meio_09,ali_meio_09_E, 'g' ) legend ( 'Ali E1 2005 0,5 km ajuste' , 'Ali E1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2005 VS 2009, 0,5km, c/ ajuste, Esquerda' )
152
%Ali Carril Direito %ali D D1 2005 1km ali_um_05_D=data_05(1:4000,5); ali_meio_05_D=data_05(1:2000,5); %ali D D1 2009 1km REF ali_um_09_D=data_09(1:4000,5); ali_meio_09_D=data_09(1:2000,5); % 2005 VS 2009, 1km figure(10) plot(x_um_05,ali_um_05_D, 'm' ,x_um_09,ali_um_09_D, 'g' ) legend ( 'Ali D d1 2005 1km' , 'Ali D d1 2009 1km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2005 VS 2009, 1km, Direito' ) % 2005 VS 2009, 0,5km figure(11) plot(x_meio_05,ali_meio_05_D, 'm' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2005 0,5 km' , 'Ali D d1 2009 0,5 km' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2005 VS 2009, 0,5km, Direito' ) % 2005 VS 2009, 0,5km, c/ ajuste, direito figure(12) plot(x_meio_05_ajuste_D,ali_meio_05_D, 'm' ,x_meio_09,ali_meio_09_D, 'g' ) legend ( 'Ali D d1 2005 0,5 km ajuste' , 'Ali D d1 2009 0,5 km ajuste' ) xlabel( 'Posicao' ) ylabel( 'Alinhamento' ) title( '2005 VS 2009, 0,5km, c/ ajuste, direito' )
• cálculo do desvio padrão em segmentos de 200 m
% calculo do desvio padrao em segmentos de 200 m data_09=dlmread( 'campanha2009pronto.txt' ); %2009 %posicao 2009 km_09=data_09(:,1); loc_09=data_09(:,2)*10^-3; x_09=[loc_09(1:4000)+km_09(1:4000) loc_09(4001:8000)+km_09(4001:8000) loc_09(8001:12000)+km_09(8001:12000) loc_09(12001:16000)+km_09(12001:16000)
153
loc_09(16001:20057)+km_09(16001:20057) loc_09(20058:24057)+km_09(20058:24057) loc_09(24058:28057)+km_09(24058:28057) loc_09(28058:32057)+km_09(28058:32057) loc_09(32058:36057)+km_09(32058:36057) loc_09(36058:40095)+km_09(36058:40095) loc_09(40096:44095)+km_09(40096:44095) loc_09(44096:48095)+km_09(44096:48095) loc_09(48096:52095)+km_09(48096:52095) loc_09(52096:56095)+km_09(52096:56095) loc_09(56096:60155)+km_09(56096:60155) loc_09(60156:64155)+km_09(60156:64155) loc_09(64156:68155)+km_09(64156:68155) loc_09(68156:72155)+km_09(68156:72155) loc_09(72156:76155)+km_09(72156:76155) loc_09(76156:80155)+km_09(76156:80155) loc_09(80156:84155)+km_09(80156:84155) loc_09(84156:88227)+km_09(84156:88227) loc_09(88228:92227)+km_09(88228:92227) loc_09(92228:96227)+km_09(92228:96227) loc_09(96228:100227)+km_09(96228:100227) loc_09(100228:104227)+km_09(100228:104227) loc_09(104228:108227)+km_09(104228:108227) loc_09(108228:112300)+km_09(108228:112300) loc_09(112301:116300)+km_09(112301:116300) loc_09(116301:120300)+km_09(116301:120300) loc_09(120301:124300)+km_09(120301:124300) loc_09(124301:128329)+km_09(124301:128329) loc_09(128330:132329)+km_09(128330:132329) loc_09(132330:136329)+km_09(132330:136329)]; %Carril Esquerdo %ali E d1 2009 total ali_09_E=data_09(:,4); %ali D d1 2009 total ali_09_D=data_09(:,5); n_09=fix(size(x_09)/800); %170 grafs; for i=1:n_09 n1_09(i)=1+800*(i-1); n2_09(i)=800*i; dp9(i)=[std([ali_09_E([n1_09(i):n2_09(i)],:)])] ; dp_09=dp9'; end
154
for i=1:n_09 n1_09(i)=1+800*(i-1); n2_09(i)=800*i; dp9r(i)=[std([ali_09_D([n1_09(i):n2_09(i)],:)]) ]; dp_09r=dp9r'; end % 2007 Setembro data_07=dlmread( 'campanha2007set.txt' ); %posicao 2007 km_07=data_07(:,1); loc_07=data_07(:,2)*10^-3; x_07=[loc_07(1:4000)+km_07(1:4000) loc_07(4001:8000)+km_07(4001:8000) loc_07(8001:12057)+km_07(8001:12057) loc_07(12058:16057)+km_07(12058:16057) loc_07(16058:20057)+km_07(16058:20057) loc_07(20058:24057)+km_07(20058:24057) loc_07(24058:28057)+km_07(24058:28057) loc_07(28058:32092)+km_07(28058:32092) loc_07(32093:36092)+km_07(32093:36092) loc_07(36093:40092)+km_07(36093:40092) loc_07(40093:44092)+km_07(40093:44092) loc_07(44093:48124)+km_07(44093:48124) loc_07(48125:52124)+km_07(48125:52124) loc_07(52125:56124)+km_07(52125:56124) loc_07(56125:60124)+km_07(56125:60124) loc_07(60125:64124)+km_07(60125:64124) loc_07(64125:68124)+km_07(64125:68124) loc_07(68125:72124)+km_07(68125:72124) loc_07(72125:76124)+km_07(72125:76124) loc_07(76125:80124)+km_07(76125:80124) loc_07(80125:84124)+km_07(80125:84124) loc_07(84125:88124)+km_07(84125:88124) loc_07(88125:92194)+km_07(88125:92194) loc_07(92195:96194)+km_07(92195:96194) loc_07(96195:100194)+km_07(96195:100194) loc_07(100195:104194)+km_07(100195:104194) loc_07(104195:108117)+km_07(104195:108117) loc_07(108118:112262)+km_07(108118:112262) loc_07(112263:116243)+km_07(112263:116243) loc_07(116244:120243)+km_07(116244:120243) loc_07(120244:124243)+km_07(120244:124243) loc_07(124244:128152)+km_07(124244:128152)];
155
%ajuste em 2007 total, esquerdo x_07_ajuste_E=x_07-4.5*10^-3; %ajuste em 2007 total, direito x_07_ajuste_D=x_07-2.4167*10^-3; %ali E d1 2007 total Setembro ali_07_E=data_07(:,4); %ali D d1 2007 total Setembro ali_07_D=data_07(:,5); n_07=fix(size(x_07_ajuste_E)/800); %160 grafs; for j=1:n_07 n1_07(j)=1+800*(j-1); n2_07(j)=800*j; dp7(j)=[std([ali_07_E([n1_07(j):n2_07(j)],:)])] ; dp_07=dp7'; end for j=1:n_07 n1_07(j)=1+800*(j-1); n2_07(j)=800*j; dp7r(j)=[std([ali_07_D([n1_07(j):n2_07(j)],:)]) ]; dp_07r=dp7r'; end % 2008 marco data_08=dlmread( 'campanha2008mar.txt' ); %posicao 2008 km_08=data_08(:,1); loc_08=data_08(:,2)*10^-3; x_08=[loc_08(1:4000)+km_08(1:4000) loc_08(4001:8000)+km_08(4001:8000) loc_08(8001:12000)+km_08(8001:12000) loc_08(12001:16000)+km_08(12001:16000)
156
loc_08(16001:20000)+km_08(16001:20000) loc_08(20001:24000)+km_08(20001:24000) loc_08(24001:28000)+km_08(24001:28000) loc_08(28001:32288)+km_08(28001:32288) loc_08(32289:36288)+km_08(32289:36288) loc_08(36289:40288)+km_08(36289:40288) loc_08(40289:44288)+km_08(40289:44288) loc_08(44289:48288)+km_08(44289:48288) loc_08(48289:52288)+km_08(48289:52288) loc_08(52289:56288)+km_08(52289:56288) loc_08(56289:60288)+km_08(56289:60288) loc_08(60289:64288)+km_08(60289:64288) loc_08(64289:68288)+km_08(64289:68288) loc_08(68289:72288)+km_08(68289:72288) loc_08(72289:76288)+km_08(72289:76288) loc_08(76289:80288)+km_08(76289:80288) loc_08(80289:84288)+km_08(80289:84288) loc_08(84289:88288)+km_08(84289:88288) loc_08(88289:92426)+km_08(88289:92426) loc_08(92427:96426)+km_08(92427:96426) loc_08(96427:100426)+km_08(96427:100426) loc_08(100427:104426)+km_08(100427:104426) loc_08(104427:108426)+km_08(104427:108426) loc_08(108427:112426)+km_08(108427:112426) loc_08(112427:116426)+km_08(112427:116426) loc_08(116427:120426)+km_08(116427:120426) loc_08(120427:124426)+km_08(120427:124426) loc_08(124427:128114)+km_08(124427:128114)]; %ajuste em 2008 total, esquerdo x_08_ajuste_E=x_08-4.5*10^-3; %ajuste em 2008 total, direito x_08_ajuste_D=x_08-2.4167*10^-3; %ali E d1 2008 total marco ali_08_E=data_08(:,4); %ali D d1 2008 total marco ali_08_D=data_08(:,5); n_08=fix(size(x_08_ajuste_E)/800); %160 grafs; for h=1:n_08 n1_08(h)=1+800*(h-1); n2_08(h)=800*h; dp8(h)=[std([ali_08_E([n1_07(h):n2_07(h)],:)])] ;
157
dp_08=dp8'; end for h=1:n_08 n1_08(h)=1+800*(h-1); n2_08(h)=800*h; dp8r(h)=[std([ali_08_D([n1_07(h):n2_07(h)],:)]) ]; dp_08r=dp8r'; end % 2005 data_05=dlmread( 'campanha2005mar.txt' ); %posicao 2005 km_05=data_05(:,1); loc_05=data_05(:,2)*10^-3; x_05=[loc_05(1:4000)+km_05(1:4000) loc_05(4001:8000)+km_05(4001:8000) loc_05(8001:12081)+km_05(8001:12081) loc_05(12082:16081)+km_05(12082:16081) loc_05(16082:20081)+km_05(16082:20081) loc_05(20082:24081)+km_05(20082:24081) loc_05(24082:28081)+km_05(24082:28081) loc_05(28082:32081)+km_05(28082:32081) loc_05(32082:36081)+km_05(32082:36081) loc_05(36082:40144)+km_05(36082:40144) loc_05(40145:44144)+km_05(40145:44144) loc_05(44145:48144)+km_05(44145:48144) loc_05(48145:52144)+km_05(48145:52144) loc_05(52145:56144)+km_05(52145:56144) loc_05(56145:60207)+km_05(56145:60207) loc_05(60208:64207)+km_05(60208:64207) loc_05(64208:68207)+km_05(64208:68207) loc_05(68208:72207)+km_05(68208:72207) loc_05(72208:76207)+km_05(72208:76207) loc_05(76208:80207)+km_05(76208:80207) loc_05(80208:84207)+km_05(80208:84207) loc_05(84208:88294)+km_05(84208:88294) loc_05(88295:92294)+km_05(88295:92294) loc_05(92295:96294)+km_05(92295:96294) loc_05(96295:100294)+km_05(96295:100294) loc_05(100295:104294)+km_05(100295:104294) loc_05(104295:108294)+km_05(104295:108294) loc_05(108295:112384)+km_05(108295:112384) loc_05(112385:116384)+km_05(112385:116384)
158
loc_05(116385:120393)+km_05(116385:120393) loc_05(120394:124393)+km_05(120394:124393) loc_05(124394:128393)+km_05(124394:128393)]; %ajuste em 2005 total, esquerdo x_05_ajuste_E=x_05+25*10^-3; %ajuste em 2005 total, direito x_05_ajuste_D=x_05-13.5*10^-3; %ali E d1 2008 total marco ali_05_E=data_05(:,4); %ali D d1 2005 total marco ali_05_D=data_05(:,5); n_05=fix(size(x_05_ajuste_E)/800); %160 grafs; for i=1:n_08 n1_05(i)=1+800*(i-1); n2_05(i)=800*i; dp5(i)=[std([ali_05_E([n1_05(i):n2_05(i)],:)])] ; dp_05=dp5'; end for i=1:n_08 n1_05(i)=1+800*(i-1); n2_05(i)=800*i; dp5r(i)=[std([ali_05_D([n1_05(i):n2_05(i)],:)]) ]; dp_05r=dp5r'; end % 2006 data_06=dlmread( 'campanha2006dez.txt' ); %posicao 2006 km_06=data_06(:,1); loc_06=data_06(:,2)*10^-3;
159
x_06=[loc_05(1:4000)+km_06(1:4000) loc_06(4001:8000)+km_06(4001:8000) loc_06(8001:12000)+km_06(8001:12000) loc_06(12001:16000)+km_06(12001:16000) loc_06(16001:20000)+km_06(16001:20000) loc_06(20001:23987)+km_06(20001:23987) loc_06(23988:27987)+km_06(23988:27987) loc_06(27988:31987)+km_06(27988:31987) loc_06(31988:35987)+km_06(31988:35987) loc_06(35988:39987)+km_06(35988:39987) loc_06(39988:43987)+km_06(39988:43987) loc_06(43988:47987)+km_06(43988:47987) loc_06(47988:51987)+km_06(47988:51987) loc_06(51988:55987)+km_06(51988:55987) loc_06(55988:59987)+km_06(55988:59987) loc_06(59988:63965)+km_06(59988:63965) loc_06(63966:67965)+km_06(63966:67965) loc_06(67966:71965)+km_06(67966:71965) loc_06(71966:75965)+km_06(71966:75965) loc_06(75966:79965)+km_06(75966:79965) loc_06(79966:83965)+km_06(79966:83965) loc_06(83966:87936)+km_06(83966:87936) loc_06(87937:91936)+km_06(87937:91936) loc_06(91937:95936)+km_06(91937:95936) loc_06(95937:99936)+km_06(95937:99936) loc_06(99937:103936)+km_06(99937:103936) loc_06(103937:107936)+km_06(103937:107936) loc_06(107937:111939)+km_06(107937:111939) loc_06(111940:115939)+km_06(111940:115939) loc_06(115940:119939)+km_06(115940:119939) loc_06(119940:123939)+km_06(119940:123939) loc_06(123940:127915)+km_06(123940:127915)]; %ajuste em 2006 dezembro total, esquerdo x_06_ajuste_E=x_06-6.75*10^-3; %ajuste em 2006 dezembro total, direito x_06_ajuste_D=x_06+23*10^-3; %ali E d1 2006 dezembro total dez ali_06_E=data_06(:,4); %ali D d1 2006 dezembro total dez ali_06_D=data_06(:,5); n_06=fix(size(x_06_ajuste_E)/800); %159 grafs; for i=1:n_06 n1_06(i)=1+800*(i-1);
160
n2_06(i)=800*i; dp6(i)=[std([ali_05_E([n1_06(i):n2_06(i)],:)])] ; dp_06=dp6'; end for i=1:n_06 n1_06(i)=1+800*(i-1); n2_06(i)=800*i; dp6r(i)=[std([ali_05_D([n1_06(i):n2_06(i)],:)]) ]; dp_06r=dp6r'; end % 2007 Dezembro data_07d=dlmread( 'campanha2007dez.txt' ); %posicao 2007 km_07d=data_07d(:,1); loc_07d=data_07d(:,2)*10^-3; x_07d=[loc_07d(1:4000)+km_07d(1:4000) loc_07d(4001:8000)+km_07d(4001:8000) loc_07d(8001:12057)+km_07d(8001:12057) loc_07d(12058:16057)+km_07d(12058:16057) loc_07d(16058:20057)+km_07d(16058:20057) loc_07d(20058:24057)+km_07d(20058:24057) loc_07d(24058:28057)+km_07d(24058:28057) loc_07d(28058:32092)+km_07d(28058:32092) loc_07d(32093:36092)+km_07d(32093:36092) loc_07d(36093:40092)+km_07d(36093:40092) loc_07d(40093:44092)+km_07d(40093:44092) loc_07d(44093:48124)+km_07d(44093:48124) loc_07d(48125:52124)+km_07d(48125:52124) loc_07d(52125:56124)+km_07d(52125:56124) loc_07d(56125:60124)+km_07d(56125:60124) loc_07d(60125:64124)+km_07d(60125:64124) loc_07d(64125:68124)+km_07d(64125:68124) loc_07d(68125:72124)+km_07d(68125:72124) loc_07d(72125:76124)+km_07d(72125:76124) loc_07d(76125:80124)+km_07d(76125:80124) loc_07d(80125:84124)+km_07d(80125:84124) loc_07d(84125:88124)+km_07d(84125:88124) loc_07d(88125:92194)+km_07d(88125:92194) loc_07d(92195:96194)+km_07d(92195:96194)
161
loc_07d(96195:100194)+km_07d(96195:100194) loc_07d(100195:104194)+km_07d(100195:104194) loc_07d(104195:108117)+km_07d(104195:108117) loc_07d(108118:112262)+km_07d(108118:112262) loc_07d(112263:116243)+km_07d(112263:116243) loc_07d(116244:120243)+km_07d(116244:120243) loc_07d(120244:124243)+km_07d(120244:124243) loc_07d(124244:128281)+km_07d(124244:128281)]; %ajuste em 2007 total, esquerdo x_07d_ajuste_E=x_07d-5.1667*10^-3; %ajuste em 2007 total, direito x_07d_ajuste_D=x_07d-5.3750*10^-3; %ali E d1 2006 dezembro total ali_07d_E=data_07d(:,4); %ali D d1 2007 dezembro total ali_07d_D=data_07d(:,5); n_07d=fix(size(x_07d_ajuste_E)/800); %159 grafs; for i=1:n_07d n1_07d(i)=1+800*(i-1); n2_07d(i)=800*i; dp7d(i)=[std([ali_07d_E([n1_07d(i):n2_07d(i)],: )])]; dp_07d=dp7d'; end for i=1:n_07d n1_07d(i)=1+800*(i-1); n2_07d(i)=800*i; dp7dr(i)=[std([ali_07d_D([n1_07d(i):n2_07d(i)], :)])]; dp_07dr=dp7dr'; end % 2007 Junho
162
data_07j=dlmread( 'campanha2007jun.txt' ); %posicao 2007 Junho km_07j=data_07j(:,1); loc_07j=data_07j(:,2)*10^-3; x_07j=[loc_07j(1:4000)+km_07j(1:4000) loc_07j(4001:8000)+km_07j(4001:8000) loc_07j(8001:12000)+km_07j(8001:12000) loc_07j(12001:16000)+km_07j(12001:16000) loc_07j(16001:20000)+km_07j(16001:20000) loc_07j(20001:24000)+km_07j(20001:24000) loc_07j(24001:28000)+km_07j(24001:28000) loc_07j(28001:32000)+km_07j(28001:32000) loc_07j(32001:36000)+km_07j(32001:36000) loc_07j(36001:40000)+km_07j(36001:40000) loc_07j(40001:44000)+km_07j(40001:44000) loc_07j(44001:48000)+km_07j(44001:48000) loc_07j(48001:52000)+km_07j(48001:52000) loc_07j(52001:56000)+km_07j(52001:56000) loc_07j(56001:60000)+km_07j(56001:60000) loc_07j(60001:64000)+km_07j(60001:64000) loc_07j(64001:68000)+km_07j(64001:68000) loc_07j(68001:72000)+km_07j(68001:72000) loc_07j(72001:76000)+km_07j(72001:76000) loc_07j(76001:80000)+km_07j(76001:80000) loc_07j(80001:84000)+km_07j(80001:84000) loc_07j(84001:87941)+km_07j(84001:87941) loc_07j(87942:91941)+km_07j(87942:91941) loc_07j(91942:95941)+km_07j(91942:95941) loc_07j(95942:99941)+km_07j(95942:99941) loc_07j(99942:103941)+km_07j(99942:103941) loc_07j(103942:107941)+km_07j(103942:107941) loc_07j(107942:111941)+km_07j(107942:111941) loc_07j(111942:115941)+km_07j(111942:115941) loc_07j(115942:119941)+km_07j(115942:119941) loc_07j(119942:123941)+km_07j(119942:123941) loc_07j(123942:127941)+km_07j(123942:127941)]; %ajuste em 2007 total, esquerdo x_07j_ajuste_E=x_07j-1*10^-3; %ajuste em 2007 total, direito x_07j_ajuste_D=x_07j+2*10^-3; %ali E d1 2006 junho total ali_07j_E=data_07j(:,4); %ali D d1 2007 junho total ali_07j_D=data_07j(:,5);
163
n_07j=fix(size(x_07j_ajuste_E)/800); %159 grafs; for i=1:n_07j n1_07j(i)=1+800*(i-1); n2_07j(i)=800*i; dp7j(i)=[std([ali_07j_E([n1_07j(i):n2_07j(i)],: )])]; dp_07j=dp7j'; end for i=1:n_07j n1_07j(i)=1+800*(i-1); n2_07j(i)=800*i; dp7jr(i)=[std([ali_07j_D([n1_07j(i):n2_07j(i)], :)])]; dp_07jr=dp7jr'; end % 2007 Marco data_07m=dlmread( 'campanha07mar.txt' ); %posicao 2007 Marco km_07m=data_07m(:,1); loc_07m=data_07m(:,2)*10^-3; x_07m=[loc_07m(1:4000)+km_07m(1:4000) loc_07m(4001:8000)+km_07m(4001:8000) loc_07m(8001:11991)+km_07m(8001:11991) loc_07m(11992:15991)+km_07m(11992:15991) loc_07m(15992:19991)+km_07m(15992:19991) loc_07m(19992:23980)+km_07m(19992:23980) loc_07m(23981:27980)+km_07m(23981:27980) loc_07m(27981:31971)+km_07m(27981:31971) loc_07m(31972:35971)+km_07m(31972:35971) loc_07m(35972:39971)+km_07m(35972:39971) loc_07m(39972:43971)+km_07m(39972:43971) loc_07m(43972:47971)+km_07m(43972:47971) loc_07m(47972:51971)+km_07m(47972:51971) loc_07m(51972:55956)+km_07m(51972:55956) loc_07m(55957:59956)+km_07m(55957:59956)
164
loc_07m(59957:63946)+km_07m(59957:63946) loc_07m(63947:67946)+km_07m(63947:67946) loc_07m(67947:71930)+km_07m(67947:71930) loc_07m(71931:75930)+km_07m(71931:75930) loc_07m(75931:79930)+km_07m(75931:79930) loc_07m(79931:83930)+km_07m(79931:83930) loc_07m(83931:87926)+km_07m(83931:87926) loc_07m(87927:91926)+km_07m(87927:91926) loc_07m(91927:95926)+km_07m(91927:95926) loc_07m(95927:99926)+km_07m(95927:99926) loc_07m(99927:103926)+km_07m(99927:103926) loc_07m(103927:107926)+km_07m(103927:107926) loc_07m(107927:111931)+km_07m(107927:111931) loc_07m(111932:115887)+km_07m(111932:115887) loc_07m(115888:119887)+km_07m(115888:119887) loc_07m(119888:123887)+km_07m(119888:123887) loc_07m(123888:127896)+km_07m(123888:127896)]; %ajuste em 2007 total, esquerdo x_07m_ajuste_E=x_07m+1*10^-3; %ajuste em 2007 total, direito x_07m_ajuste_D=x_07m+2*10^-3; %ali E d1 2007 total marco ali_07m_E=data_07m(:,4); %ali D d1 2007 total marco ali_07m_D=data_07m(:,5); n_07m=fix(size(x_07m_ajuste_E)/800); %159 grafs; for i=1:n_07m n1_07m(i)=1+800*(i-1); n2_07m(i)=800*i; dp7m(i)=[std([ali_07m_E([n1_07m(i):n2_07m(i)],: )])]; dp_07m=dp7m'; end for i=1:n_07m n1_07m(i)=1+800*(i-1); n2_07m(i)=800*i;
165
dp7mr(i)=[std([ali_07m_D([n1_07m(i):n2_07m(i)], :)])]; dp_07mr=dp7mr'; end % 2008 Junho data_08j=dlmread( 'campanha2008pronto.txt' ); %posicao 2008 Junho km_08j=data_08j(:,1); loc_08j=data_08j(:,2)*10^-3; x_08j=[loc_08j(1:4000)+km_08j(1:4000) loc_08j(4001:8000)+km_08j(4001:8000) loc_08j(8001:12000)+km_08j(8001:12000) loc_08j(12001:16000)+km_08j(12001:16000) loc_08j(16001:20000)+km_08j(16001:20000) loc_08j(20001:24062)+km_08j(20001:24062) loc_08j(24063:28062)+km_08j(24063:28062) loc_08j(28063:32087)+km_08j(28063:32087) loc_08j(32088:36087)+km_08j(32088:36087) loc_08j(36088:40087)+km_08j(36088:40087) loc_08j(40088:44087)+km_08j(40088:44087) loc_08j(44088:48087)+km_08j(44088:48087) loc_08j(48088:52087)+km_08j(48088:52087) loc_08j(52088:56087)+km_08j(52088:56087) loc_08j(56088:60087)+km_08j(56088:60087) loc_08j(60088:64087)+km_08j(60088:64087) loc_08j(64088:68087)+km_08j(64088:68087) loc_08j(68088:72149)+km_08j(68088:72149) loc_08j(72150:76149)+km_08j(72150:76149) loc_08j(76150:80149)+km_08j(76150:80149) loc_08j(80150:84149)+km_08j(80150:84149) loc_08j(84150:88149)+km_08j(84150:88149) loc_08j(88150:92149)+km_08j(88150:92149) loc_08j(92150:96149)+km_08j(92150:96149) loc_08j(96150:100149)+km_08j(96150:100149) loc_08j(100150:104149)+km_08j(100150:104149) loc_08j(104150:108149)+km_08j(104150:108149) loc_08j(108150:112267)+km_08j(108150:112267) loc_08j(112268:116267)+km_08j(112268:116267) loc_08j(116268:120267)+km_08j(116268:120267) loc_08j(120268:124267)+km_08j(120268:124267) loc_08j(124268:128284)+km_08j(124268:128284)]; %ajuste em 2008 total, esquerdo x_08j_ajuste_E=x_08j-1.333*10^-3; %ajuste em 2008 total, direito x_08j_ajuste_D=x_08j-1.7167*10^-3;
166
%ali E d1 2008 total ali_08j_E=data_08j(:,4); %ali D d1 2008 total ali_08j_D=data_08j(:,5); n_08j=fix(size(x_08j_ajuste_E)/800); %159 grafs; for i=1:n_08j n1_08j(i)=1+800*(i-1); n2_08j(i)=800*i; dp8j(i)=[std([ali_08j_E([n1_08j(i):n2_08j(i)],: )])]; dp_08j=dp8j'; end for i=1:n_08j n1_08j(i)=1+800*(i-1); n2_08j(i)=800*i; dp8jr(i)=[std([ali_08j_D([n1_08j(i):n2_08j(i)], :)])]; dp_08jr=dp8jr'; end z=[x_06(1:800:127115), dp_05(1:159), dp_06, dp_07(1 :159), dp_07j(1:159), dp_07d(1:159), dp_08(1:159), dp_09(1:159), dp_07m(1 :159), dp_08j(1:159)]; zr=[x_06(1:800:127115), dp_05r(1:159), dp_06r, dp_0 7r(1:159), dp_07jr(1:159), dp_07dr(1:159), dp_08r(1:159), dp_0 9r(1:159), dp_07mr(1:159), dp_08jr(1:159)]; k=[2005.25, 2006.92, 2007.75, 2007.92, 2007.5, 2008 .25, 2009.25, 2007.25, 2008.5]; % Ciclo que faz os graficos do desvio padrao em fun cao dos anos %Esquerda for e=1:2 %159
167
figure(e) plot(k, z(e,2:10), '*' ) xlabel( 'Campanha' ) ylabel( 'Desvio Padrao' ) title( 'Desvio padrão carril Esquerdo' ) end %Direita for y=1:2 %159 figure(e+y) plot(k, zr(y,2:10), '+' ) xlabel( 'Campanha' ) ylabel( 'Desvio Padrao' ) title( 'Desvio padrão carril Direito' ) end %3D Esquerda x_3_d=[x_05_ajuste_E(1:127896), x_06_ajuste_E(1:127 896), x_07_ajuste_E(1:127896), x_07j_ajuste_E(1:127896), x_07d_ajuste_E(1:127896), x_08_ajuste_E(1:127896), x_09(1:127896), x_07m_ajuste_E, x_08j_ajuste_E(1:127896)]; ali_3_d=[ali_05_E(1:127896), ali_06_E(1:127896), al i_07_E(1:127896), ali_07j_E(1:127896), ali_07d_E(1:127896), ali_08_E( 1:127896), ali_09_E(1:127896), ali_07m_E, ali_08j_E(1:127896)] ; %3D Direita x_3_dr=[x_05_ajuste_D(1:127896), x_06_ajuste_D(1:12 7896), x_07_ajuste_D(1:127896), x_07j_ajuste_D(1:127896), x_07d_ajuste_D(1:127896), x_08_ajuste_D(1:127896), x_09(1:127896), x_07m_ajuste_D, x_08j_ajuste_D(1:127896)]; ali_3_dr=[ali_05_D(1:127896), ali_06_D(1:127896), a li_07_D(1:127896), ali_07j_D(1:127896), ali_07d_D(1:127896), ali_08_D( 1:127896), ali_09_D(1:127896), ali_07m_D, ali_08j_D(1:127896)] ; figure(e+y+1) for r=1:9 A=[x_3_d(:,r) ali_3_d(:,r)]; s=ones(size(A(:,1))); line(s*k(1,r),A(:,1),A(:,2)) xlabel( 'Campanha' ) ylabel( 'Posicao' ) zlabel( 'Alinhamento' ) title( 'Esquerda' ) end figure(e+y+2)