Mudanças entre as edições de "Análise espectral e estabilidade"

De WikiLICC
Ir para: navegação, pesquisa
m (Condição de Neumann em todos os lados)
m (O problema)
Linha 8: Linha 8:
  
 
===Condição de Neumann em todos os lados===
 
===Condição de Neumann em todos os lados===
O espectro é real com autovalores entre &sigma C (-1,1] incluindo 0 e 1.
+
O espectro é real com espectro σ⊂(-1,1] incluindo 0 e 1.
  
 
===Condição de Neumann em todos os lados fixando um ponto, P(1,2)=1 ===
 
===Condição de Neumann em todos os lados fixando um ponto, P(1,2)=1 ===
O espectro tem o maior autovalor λ=1.6537 e os outros autovalores próximos do eixo real (talvez devido a erros de ponto flutuante) entre (-1,1) incluindo 0.
+
O espectro tem autovalores próximos do eixo real (talvez devido a erros de ponto flutuante) com σ⊂ 1.6537 ∪ (-1,1) incluindo 0.
  
 
===Condição de Neumann em 3 lados fixando P=1 na entrada (ou na saída)===
 
===Condição de Neumann em 3 lados fixando P=1 na entrada (ou na saída)===
O espectro tem o maior autovalor λ=11.23 e os outros autovalores próximos do eixo real (com um pequeno circulo de valores complexos) entre (-1,1) incluindo 0.
+
O espectro tem a maioria dos autovalores próximos do eixo real com alguns autovalores complexos (um pequeno circulo de valores complexos dentro do círculo unitário) com σ ⊂ 11.23 ∪ (-1,1) incluindo 0.
 
===Teste 4===
 
===Teste 4===
  
 
===Teste 5===
 
===Teste 5===

Edição das 23h11min de 21 de julho de 2009

A análise dos autovalores de uma matriz de iteração pode ser usada para estudar a estabilidade de um método iterativo.

Vamos relatar um estudo para um problema específico.

O problema

Queremos aproximar a solução da equação de Navier Stokes em um duto. Para isso devemos resolver a cada passo de tempo uma equação de Poisson como <math>\frac{\partial^2 p}{\partial x^2}+\frac{\partial^2 p}{\partial y^2}= f(u)\,\!</math>

Condição de Neumann em todos os lados

O espectro é real com espectro σ⊂(-1,1] incluindo 0 e 1.

Condição de Neumann em todos os lados fixando um ponto, P(1,2)=1

O espectro tem autovalores próximos do eixo real (talvez devido a erros de ponto flutuante) com σ⊂ 1.6537 ∪ (-1,1) incluindo 0.

Condição de Neumann em 3 lados fixando P=1 na entrada (ou na saída)

O espectro tem a maioria dos autovalores próximos do eixo real com alguns autovalores complexos (um pequeno circulo de valores complexos dentro do círculo unitário) com σ ⊂ 11.23 ∪ (-1,1) incluindo 0.

Teste 4

Teste 5