8ª Aula Prática de TFCOMP-II

Comecem por criar e guardar em disco um Notebook com o nome TFC8_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).


Problema 1

Usando os nomes de dois elementos do seu grupo, crie duas imagens com letras do tipo [Graphics:Images/Pratica8Sol_gr_1.gif] e transforme-as em duas matrizes quadradas de números reais (use as funcionalidades de [Graphics:Images/Pratica8Sol_gr_2.gif] e/ou [Graphics:Images/Pratica8Sol_gr_3.gif]). Usando [Graphics:Images/Pratica8Sol_gr_4.gif][Graphics:Images/Pratica8Sol_gr_5.gif],   represente com [Graphics:Images/Pratica8Sol_gr_6.gif] ou [Graphics:Images/Pratica8Sol_gr_7.gif] o efeito nas imagens de alterar a fase ou amplitude dos respectivos coeficientes de Fourier. Observe também a convolução e a correlação destas imagens uma com a outra. Que conclusões tira sobre a fase e a amplitude dos coeficientes de Fourier?



Problema 2

Determine a solução da equação de difusão [Graphics:Images/Pratica8Sol_gr_8.gif] com a condição inicial [Graphics:Images/Pratica8Sol_gr_9.gif] se [Graphics:Images/Pratica8Sol_gr_10.gif] e zero fora deste intervalo. Faça uma animação mostrando a evolução da solução para [Graphics:Images/Pratica8Sol_gr_11.gif].

Solução 2

Solução usando transformadas de Fourier

Para a solução deste problema veja o documento incluído na página da cadeira sob o título "Equações de Difusão e Fourier"

Condição inicial e sua transformada.

[Graphics:Images/Pratica8Sol_gr_12.gif]
[Graphics:Images/Pratica8Sol_gr_13.gif]
[Graphics:Images/Pratica8Sol_gr_14.gif]
Transformada da Equação de Fourier (atenção ao [Graphics:Images/Pratica8Sol_gr_15.gif]: sem ele o resultado é errado)
[Graphics:Images/Pratica8Sol_gr_16.gif]
[Graphics:Images/Pratica8Sol_gr_17.gif]
[Graphics:Images/Pratica8Sol_gr_18.gif]
[Graphics:Images/Pratica8Sol_gr_19.gif]
A solução obtém-se pela transformada inversa de [Graphics:Images/Pratica8Sol_gr_20.gif] quando se substitui [Graphics:Images/Pratica8Sol_gr_21.gif] pela transformada da condição inicial.
[Graphics:Images/Pratica8Sol_gr_22.gif]
[Graphics:Images/Pratica8Sol_gr_23.gif]
[Graphics:Images/Pratica8Sol_gr_28.gif]

Solução por convolução com o "Heat Kernel" [Graphics:Images/Pratica8Sol_gr_29.gif] da condição inicial [Graphics:Images/Pratica8Sol_gr_30.gif]

[Graphics:Images/Pratica8Sol_gr_31.gif]
[Graphics:Images/Pratica8Sol_gr_32.gif]
[Graphics:Images/Pratica8Sol_gr_33.gif]
[Graphics:Images/Pratica8Sol_gr_34.gif]
[Graphics:Images/Pratica8Sol_gr_39.gif]



Problema 3

Uma corda de guitarra está presa em [Graphics:Images/Pratica8Sol_gr_40.gif][Graphics:Images/Pratica8Sol_gr_41.gif]. Sugira condições iniciais [Graphics:Images/Pratica8Sol_gr_42.gif][Graphics:Images/Pratica8Sol_gr_43.gif] e, usando transformadas de Fourier e/ou de Laplace, resolva a equação  de ondas [Graphics:Images/Pratica8Sol_gr_44.gif], onde [Graphics:Images/Pratica8Sol_gr_45.gif][Graphics:Images/Pratica8Sol_gr_46.gif] a tensão na corda e [Graphics:Images/Pratica8Sol_gr_47.gif] a sua densidade linear de massa. Mostre a evolução da solução numa animação e use a função [Graphics:Images/Pratica8Sol_gr_48.gif] para ouvir o som gerado pela corda, e compare a intensidade do tom fundamental com os das harmónicas seguintes .

Solução 3

A solução geral de uma equação do tipo [Graphics:Images/Pratica8Sol_gr_49.gif] é [Graphics:Images/Pratica8Sol_gr_50.gif], onde [Graphics:Images/Pratica8Sol_gr_51.gif][Graphics:Images/Pratica8Sol_gr_52.gif] são funções arbitrárias, representando perfis que se deslocam para a direita e a esquerda com velocidade [Graphics:Images/Pratica8Sol_gr_53.gif]. De facto:
[Graphics:Images/Pratica8Sol_gr_54.gif]
[Graphics:Images/Pratica8Sol_gr_55.gif]
[Graphics:Images/Pratica8Sol_gr_56.gif]
[Graphics:Images/Pratica8Sol_gr_57.gif]
Quando se especificam condições iniciais, restringe-se a forma de [Graphics:Images/Pratica8Sol_gr_58.gif][Graphics:Images/Pratica8Sol_gr_59.gif].
[Graphics:Images/Pratica8Sol_gr_60.gif]
[Graphics:Images/Pratica8Sol_gr_61.gif]
[Graphics:Images/Pratica8Sol_gr_62.gif]
[Graphics:Images/Pratica8Sol_gr_63.gif]
[Graphics:Images/Pratica8Sol_gr_64.gif]
[Graphics:Images/Pratica8Sol_gr_65.gif]
[Graphics:Images/Pratica8Sol_gr_66.gif]
Vemos assim que a solução que verifica as condições iniciais é, àparte as constantes, [Graphics:Images/Pratica8Sol_gr_67.gif]
[Graphics:Images/Pratica8Sol_gr_68.gif]
[Graphics:Images/Pratica8Sol_gr_69.gif]
Frequentemente, as condições iniciais indicam que a corda deformada parte do repouso, i.e. [Graphics:Images/Pratica8Sol_gr_70.gif].
Neste caso [Graphics:Images/Pratica8Sol_gr_71.gif][Graphics:Images/Pratica8Sol_gr_72.gif].
No caso contrário a corda parte do equilíbrio [Graphics:Images/Pratica8Sol_gr_73.gif] com velocidade inicial dada por [Graphics:Images/Pratica8Sol_gr_74.gif], e então [Graphics:Images/Pratica8Sol_gr_75.gif], onde [Graphics:Images/Pratica8Sol_gr_76.gif].

Para além de condições iniciais pede-se ainda que se verifiquem condições fronteira, i.e. no caso [Graphics:Images/Pratica8Sol_gr_77.gif][Graphics:Images/Pratica8Sol_gr_78.gif].  Vamos considerar o caso de [Graphics:Images/Pratica8Sol_gr_79.gif].

[Graphics:Images/Pratica8Sol_gr_80.gif]
O método das imagens permite visualizar de forma simples a solução pretendida. Consideremos uma solução que é a sobreposição de três destas formas de onda: para [Graphics:Images/Pratica8Sol_gr_81.gif][Graphics:Images/Pratica8Sol_gr_82.gif], onde [Graphics:Images/Pratica8Sol_gr_83.gif] .  Esta forma propaga-se para a esquerda como [Graphics:Images/Pratica8Sol_gr_84.gif] e para a direita como [Graphics:Images/Pratica8Sol_gr_85.gif], e para garantir que [Graphics:Images/Pratica8Sol_gr_86.gif] é necessário adicionar à esquerda uma solução partindo de [Graphics:Images/Pratica8Sol_gr_87.gif] para [Graphics:Images/Pratica8Sol_gr_88.gif],  e à direita [Graphics:Images/Pratica8Sol_gr_89.gif]para [Graphics:Images/Pratica8Sol_gr_90.gif],  . Estas condições iniciais para segmentos virtuais da corda são imagens virtuais, com focos em [Graphics:Images/Pratica8Sol_gr_91.gif][Graphics:Images/Pratica8Sol_gr_92.gif],  da forma inicial da corda no intervalo [Graphics:Images/Pratica8Sol_gr_93.gif].
[Graphics:Images/Pratica8Sol_gr_94.gif]
Vamos  exemplificar com [Graphics:Images/Pratica8Sol_gr_98.gif][Graphics:Images/Pratica8Sol_gr_99.gif].
[Graphics:Images/Pratica8Sol_gr_100.gif]
[Graphics:Images/Pratica8Sol_gr_103.gif]

Infelizmente esta não é a solução final porque, embora se verifiquem as condições fronteira, apenas se observam duas reflexões.   Para se observar um movimento periódico tem que se  estender a definição de [Graphics:Images/Pratica8Sol_gr_104.gif] da seguinte forma

[Graphics:Images/Pratica8Sol_gr_105.gif]
Obtemos assim uma condição inicial periódica de período [Graphics:Images/Pratica8Sol_gr_106.gif] sobre toda a recta [Graphics:Images/Pratica8Sol_gr_107.gif], e s sua evolução temporal como solução da equação da corda vibrante dá-nos a solução pretendida no intervalo [Graphics:Images/Pratica8Sol_gr_108.gif]

[Graphics:Images/Pratica8Sol_gr_111.gif]

A seguinte animação mostra como [Graphics:Images/Pratica8Sol_gr_112.gif][Graphics:Images/Pratica8Sol_gr_113.gif].

[Graphics:Images/Pratica8Sol_gr_117.gif]

Olhemos agora para o movimento de um ponto fixo da corda. O som gerado pela corda tem mais a ver com este movimento temporal do que com a forma espacial da corda (o som corresponde a uma oscilação longitudinal do ar, enquanto a deformação da corda é transversal).  Para ouvir um som tem que existir uma oscilação de frequência [Graphics:Images/Pratica8Sol_gr_118.gif].  Quando [Graphics:Images/Pratica8Sol_gr_119.gif] é conhecido, podemos estimar que [Graphics:Images/Pratica8Sol_gr_120.gif]para formas iniciais muito simples.

[Graphics:Images/Pratica8Sol_gr_121.gif]
[Graphics:Images/Pratica8Sol_gr_122.gif]
[Graphics:Images/Pratica8Sol_gr_123.gif]
[Graphics:Images/Pratica8Sol_gr_124.gif]
[Graphics:Images/Pratica8Sol_gr_125.gif]
[Graphics:Images/Pratica8Sol_gr_126.gif]
[Graphics:Images/Pratica8Sol_gr_127.gif]

Compare este com um tom puro de [Graphics:Images/Pratica8Sol_gr_131.gif]:



Problema 4

A evolução de uma partícula num potencial [Graphics:Images/Pratica8Sol_gr_132.gif] é descrita em mecânica quântica por uma equação do tipo [Graphics:Images/Pratica8Sol_gr_133.gif]. Usando transformadas de Fourier mostre que, para a partícula livre ([Graphics:Images/Pratica8Sol_gr_134.gif]) esta equação é fácilmente resolúvel e construa o propagador [Graphics:Images/Pratica8Sol_gr_135.gif] usando uma discretização adequada para uma forma de onda inicial arbitrária [Graphics:Images/Pratica8Sol_gr_136.gif]. Represente a sua evolução temporal gráficamente.


Converted by Mathematica      June 22, 2001