Corda Vibrante: mudanças entre as edições

De Física Computacional
Ir para navegação Ir para pesquisar
corrigida a constante de normalização na análise espectral
análise espectral: bases enumeráveis
Linha 26: Linha 26:
</math>  
</math>  
</center>  
</center>  
Se o domínio de <math>s(t)</math> é limitado, digamos <math>t \in [0, T]</math>, então uma base infinita com cardinalidade enumerável (em contraste com a base anterior, que possui cardinalidade não enumerável) é suficiente para representar <math>s(t)</math>, uma possível base ortonormal é a seguinte: <math>\bigg\{ \sqrt{\frac{1}{T}}, \sqrt{\frac{2}{T}}sen(\omega_n t), \sqrt{\frac{2}{T}}cos(\omega_n t) \bigg\}_{n \in \mathbb{N^*}}</math>, em que <math> \omega_n = \frac{n\pi}{T} </math>. Dessa forma, a representação e coordenadas de <math>s(t)</math> ficam
<center>
<math>
\begin{aligned}
s(t) &= a_0\sqrt{\frac{1}{T}} + \sqrt{\frac{2}{T}}\sum_{n=1}^{\infty}\bigg[a_n cos(\omega_n t) + b_n sen(\omega_n) \bigg] \\
a_0 &= \int_0^T s(t) \sqrt{\frac{1}{T}}dt \\
a_n &= \int_0^T s(t) \sqrt{\frac{2}{T}}cos(\omega_n t)dt \\
b_n &= \int_0^T s(t) \sqrt{\frac{2}{T}}sen(\omega_n t)dt \\
\end{aligned}
</math>
</center>
É impossível falar sobre bases enumeráveis de um sub-espaço de <math> \mathbb{R}^{\mathbb{R}} </math> sem representar esse canhão matemático com uma animação. Abaixo segue uma animação que calcula as primeiras <math>N</math> coordenadas (<math>a_0,\dots,a_N</math> e <math>b_1,\dots,b_N</math>) de um sinal qualquer e sobrepõem a série obtida incrementando <math>N</math> até as duas curvas serem indistinguíveis a olho nu.


Agora, considerando uma corda vibrante, o nosso sinal sonoro provém da vibração de um ponto específico da corda, digamos em <math>x=x_o</math>, então a função que representa esse sinal é <math>y(x_o, t)</math>
Agora, considerando uma corda vibrante, o nosso sinal sonoro provém da vibração de um ponto específico da corda, digamos em <math>x=x_o</math>, então a função que representa esse sinal é <math>y(x_o, t)</math>

Edição das 13h40min de 20 de abril de 2024

A equação da onda

Método FTCS

Sobre estabilidade

Análise espectral

Uma possível forma para quantitativamente analisar o som gerado por uma corda vibrante é estudar as frequências que compõem o seu movimento, técnica essa chamada de análise espectral. Antes de prosseguirmos vamos recapitular alguns resultados da álgebra linear

Supremacia da álgebra linear

O seguinte conjunto ={f|f:} é o espaço de funções reais de uma variável. Esse conjunto é um espaço vetorial, logo podemos utilizar toda a artilharia da álgebra linear, em especial, estamos interessados no sub-espaço gerado pela base B={sen(ωt)/π,cos(ωt)/π}ω+ [1], pois elementos de B, interpretados como sinais sonoros, representam um frequência pura de valor f=ω/(2π). Dessa forma, um sinal arbitrário s(t) pode ser escrito em termos das frequências puras que o formam

s(t)=1π0[a(ω)cos(ωt)+b(ω)sen(ωt]dω

E podemos extrair suas coordenadas (a(ω) e b(ω)), fazendo o produto escalar com os elementos da base

a(ω)=s(t)1πcos(ωt)dtb(ω)=s(t)1πsen(ωt)dt

Se o domínio de s(t) é limitado, digamos t[0,T], então uma base infinita com cardinalidade enumerável (em contraste com a base anterior, que possui cardinalidade não enumerável) é suficiente para representar s(t), uma possível base ortonormal é a seguinte: {1T,2Tsen(ωnt),2Tcos(ωnt)}n*, em que ωn=nπT. Dessa forma, a representação e coordenadas de s(t) ficam

s(t)=a01T+2Tn=1[ancos(ωnt)+bnsen(ωn)]a0=0Ts(t)1Tdtan=0Ts(t)2Tcos(ωnt)dtbn=0Ts(t)2Tsen(ωnt)dt

É impossível falar sobre bases enumeráveis de um sub-espaço de sem representar esse canhão matemático com uma animação. Abaixo segue uma animação que calcula as primeiras N coordenadas (a0,,aN e b1,,bN) de um sinal qualquer e sobrepõem a série obtida incrementando N até as duas curvas serem indistinguíveis a olho nu.

Agora, considerando uma corda vibrante, o nosso sinal sonoro provém da vibração de um ponto específico da corda, digamos em x=xo, então a função que representa esse sinal é y(xo,t)

Condição inicial para uma corda de violão

Notas

  1. A constante 1/π está presente por questão de normalização. Esse caso pode parecer um pouco estranho, dado que não é possível normalizar os cossenos e senos, pois sua integral em todo a reta não é definida, mas o que se deseja é a seguinte propriedade Aωcos(ωt)Aωcos(ωt)dt=δ(ωω) que é safisfeita quando Aω=Aω=1/π,ω,ω