AULA PRÁTICA DE TFCOMP-II

Solução


                 ©Amaro R. Silva, Prof                     http://centra.ist.utl.pt/~amaro 


Comecem por criar e guardar em disco um Notebook com o nome TFC2_grupo.nb, onde 'grupo' é o nome do vosso grupo de trabalho. Neste Notebook criem uma célula de título identificando os elementos constituintes do grupo (nome, número e contactos).
Use a função nativa [Graphics:Images/Pratica2_gr_1.gif] para gerar uma função [Graphics:Images/Pratica2_gr_2.gif] que imprima no Notebook a data de cada dia em que for avaliada, no formato compatível com, por exemplo: [Graphics:Images/Pratica2_gr_3.gif]. Guarde essa definição num documento "data.m" donde possa carregar a definição em qualquer altura que o deseje.
[Graphics:Images/Pratica2_gr_4.gif]
[Graphics:Images/Pratica2_gr_5.gif]
[Graphics:Images/Pratica2_gr_6.gif]



Redefina a função nativa [Graphics:Images/Pratica2_gr_7.gif] de forma que, quando avaliada, execute as regras usuais, por exemplo [Graphics:Images/Pratica2_gr_8.gif]. Como faria para implementar a transformação inversa agora?
[Graphics:Images/Pratica2_gr_9.gif]
[Graphics:Images/Pratica2_gr_10.gif]
[Graphics:Images/Pratica2_gr_11.gif]
[Graphics:Images/Pratica2_gr_12.gif]
[Graphics:Images/Pratica2_gr_13.gif]
[Graphics:Images/Pratica2_gr_14.gif]



Crie versores [Graphics:Images/Pratica2_gr_15.gif] formando uma base ortogonal para o sistema de coordenadas polares esféricas em [Graphics:Images/Pratica2_gr_16.gif]. A partir destes defina a matriz [Graphics:Images/Pratica2_gr_17.gif] que transforma os versores [Graphics:Images/Pratica2_gr_18.gif][Graphics:Images/Pratica2_gr_19.gif][Graphics:Images/Pratica2_gr_20.gif]do sistema cartesiano de coordenadas nos versores das coordenadas polares esféricas. Mostre que esta matriz [Graphics:Images/Pratica2_gr_21.gif] é ortogonal, i.e. [Graphics:Images/Pratica2_gr_22.gif](matriz transposta). Verifique que o determinante de [Graphics:Images/Pratica2_gr_23.gif] é unidade e que as linhas e colunas formam vectores ortogonais entre si.
[Graphics:Images/Pratica2_gr_24.gif]
[Graphics:Images/Pratica2_gr_25.gif]
[Graphics:Images/Pratica2_gr_26.gif]
[Graphics:Images/Pratica2_gr_27.gif]
[Graphics:Images/Pratica2_gr_28.gif]
[Graphics:Images/Pratica2_gr_29.gif]
[Graphics:Images/Pratica2_gr_30.gif]
[Graphics:Images/Pratica2_gr_31.gif]
[Graphics:Images/Pratica2_gr_32.gif]
[Graphics:Images/Pratica2_gr_33.gif]
[Graphics:Images/Pratica2_gr_34.gif]
[Graphics:Images/Pratica2_gr_35.gif]
[Graphics:Images/Pratica2_gr_36.gif]
[Graphics:Images/Pratica2_gr_37.gif]
[Graphics:Images/Pratica2_gr_38.gif]
[Graphics:Images/Pratica2_gr_39.gif]
[Graphics:Images/Pratica2_gr_40.gif]
[Graphics:Images/Pratica2_gr_41.gif]
[Graphics:Images/Pratica2_gr_42.gif]
[Graphics:Images/Pratica2_gr_43.gif]
[Graphics:Images/Pratica2_gr_44.gif] [Graphics:Images/Pratica2_gr_45.gif] [Graphics:Images/Pratica2_gr_46.gif]
[Graphics:Images/Pratica2_gr_47.gif] [Graphics:Images/Pratica2_gr_48.gif] [Graphics:Images/Pratica2_gr_49.gif]
[Graphics:Images/Pratica2_gr_50.gif] 0 [Graphics:Images/Pratica2_gr_51.gif]
[Graphics:Images/Pratica2_gr_52.gif]
[Graphics:Images/Pratica2_gr_53.gif]



Mostre explícitamente que o ângulo formado por dois vectores [Graphics:Images/Pratica2_gr_54.gif][Graphics:Images/Pratica2_gr_55.gif] em [Graphics:Images/Pratica2_gr_56.gif], bem como a magnitude de cada, é invariante para a transformação [Graphics:Images/Pratica2_gr_57.gif] definida no problema anterior (de facto isso é uma propriedade de todas as transformações ortogonais!)
[Graphics:Images/Pratica2_gr_58.gif]
Portanto:
[Graphics:Images/Pratica2_gr_59.gif]
[Graphics:Images/Pratica2_gr_60.gif]



Use as funções matriciais nativas do Mathematica para determinar, dada qualquer matriz [Graphics:Images/Pratica2_gr_61.gif]não singular, a matriz de transformação [Graphics:Images/Pratica2_gr_62.gif] que diagonaliza [Graphics:Images/Pratica2_gr_63.gif], i.e. [Graphics:Images/Pratica2_gr_64.gif] é diagonal. Crie uma função [Graphics:Images/Pratica2_gr_65.gif] para calcular a exponencial duma matriz a partir desta transformação (note que no referencial emque [Graphics:Images/Pratica2_gr_66.gif] é diagonal a exponenciação é trivial). Confirme o resultado obtido com a função nativa [Graphics:Images/Pratica2_gr_67.gif].
[Graphics:Images/Pratica2_gr_68.gif]
[Graphics:Images/Pratica2_gr_69.gif]
[Graphics:Images/Pratica2_gr_70.gif]
[Graphics:Images/Pratica2_gr_71.gif]
[Graphics:Images/Pratica2_gr_72.gif]
[Graphics:Images/Pratica2_gr_73.gif]
A matriz [Graphics:Images/Pratica2_gr_74.gif]é não-singular
[Graphics:Images/Pratica2_gr_75.gif]
[Graphics:Images/Pratica2_gr_76.gif]
[Graphics:Images/Pratica2_gr_77.gif] representa o conjunto [Graphics:Images/Pratica2_gr_78.gif] para [Graphics:Images/Pratica2_gr_79.gif]
[Graphics:Images/Pratica2_gr_80.gif]
[Graphics:Images/Pratica2_gr_81.gif]
A matriz [Graphics:Images/Pratica2_gr_82.gif] cujas colunas são os vectores próprios de [Graphics:Images/Pratica2_gr_83.gif] transforma uma base [Graphics:Images/Pratica2_gr_84.gif]numa em que os vectores [Graphics:Images/Pratica2_gr_85.gif]são esses vectores próprios de [Graphics:Images/Pratica2_gr_86.gif].
[Graphics:Images/Pratica2_gr_87.gif]
[Graphics:Images/Pratica2_gr_88.gif]
Os  valores próprios de [Graphics:Images/Pratica2_gr_89.gif] também são dados por [Graphics:Images/Pratica2_gr_90.gif]
[Graphics:Images/Pratica2_gr_91.gif]
[Graphics:Images/Pratica2_gr_92.gif]
A matriz diagonal [Graphics:Images/Pratica2_gr_93.gif] cujas entradas são os valores próprios [Graphics:Images/Pratica2_gr_94.gif] obtém-se simplesmente de
[Graphics:Images/Pratica2_gr_95.gif]
[Graphics:Images/Pratica2_gr_96.gif]
Invertendo esta relação obtém-se que
[Graphics:Images/Pratica2_gr_97.gif]
[Graphics:Images/Pratica2_gr_98.gif]
Note agora que [Graphics:Images/Pratica2_gr_99.gif], por exemplo
[Graphics:Images/Pratica2_gr_100.gif]
A função [Graphics:Images/Pratica2_gr_101.gif] tem um desenvolvimento em série muito familiar
[Graphics:Images/Pratica2_gr_102.gif]
[Graphics:Images/Pratica2_gr_103.gif]
Se [Graphics:Images/Pratica2_gr_104.gif] for substituído por [Graphics:Images/Pratica2_gr_105.gif], então é fácil de ver que
[Graphics:Images/Pratica2_gr_106.gif]
[Graphics:Images/Pratica2_gr_107.gif]
ou seja, factorizando [Graphics:Images/Pratica2_gr_108.gif][Graphics:Images/Pratica2_gr_109.gif]
[Graphics:Images/Pratica2_gr_110.gif]
Assim [Graphics:Images/Pratica2_gr_111.gif] reduz-se a calcular [Graphics:Images/Pratica2_gr_112.gif] para a matriz diagonal [Graphics:Images/Pratica2_gr_113.gif], o que é ainda uma matriz diagonal com entradas iguais à exponencial das de [Graphics:Images/Pratica2_gr_114.gif], e multplicar tudo à direita e à esquerda por [Graphics:Images/Pratica2_gr_115.gif][Graphics:Images/Pratica2_gr_116.gif].
[Graphics:Images/Pratica2_gr_117.gif]
[Graphics:Images/Pratica2_gr_118.gif]
[Graphics:Images/Pratica2_gr_119.gif]
[Graphics:Images/Pratica2_gr_120.gif]
[Graphics:Images/Pratica2_gr_121.gif]
[Graphics:Images/Pratica2_gr_122.gif]

Converted by Mathematica      April 18, 2001