Física Computacional - Contribuições do usuário [pt-br]

Histogramas

2011-09-19T18:35:13Z

Tekkito: Criou página com 'Os agentes estão distribuidos numa rede bi-dimensional e trocam capital com os segundos vizinhos através da regra do mínimo. A animação mostra a evolução da distribuição...'

Os agentes estão distribuidos numa rede bi-dimensional e trocam capital com os segundos vizinhos através da regra do mínimo. A animação mostra a evolução da distribuição de capital até a formação de um "gap", separando claramente os agentes em duas classe econômicas distintas.

[[Image:histograma1.gif]]

Cálculo do valor inicial do índice Gini (Gaspar)

2011-09-19T18:34:39Z

Tekkito: Criou página com 'Em muitas simulações, a condição inicial do sistema de agentes é uma distribuição uniforme de riquezas, com valores compreendidos entre 0 e 1. O valor do índice Gini ness...'

Em muitas simulações, a condição inicial do sistema de agentes é uma distribuição uniforme de riquezas, com valores compreendidos entre 0 e 1. O valor do índice Gini nessas condições é calculado abaixo.

Iniciamos com a definição do índice Gini. Para um sistema composto por <math>N</math> agentes, temos

<center>
<math>
G \equiv {1 \over (N-1) \sum_{i=1}^{N} w_i} \sum_{j=1}^{N-1} \sum_{k=j+1}^N \left| w_j - w_k \right|
</math>
</center>

Queremos calcular o valor médio de <math>G</math> quando a distribuição for uniforme, ou seja,

<center>
<math>
\begin{matrix}

\left \langle G \right \rangle &=& \left \langle {1 \over (N-1) \sum_{i=1}^{N} w_i} \sum_{j=1}^{N-1} \sum_{k=j+1}^N \left| w_j - w_k \right| \right \rangle \\
& & \\
& = & {1 \over (N-1) \left \langle \sum_{i=1}^{N} w_i \right \rangle} \left \langle \sum_{j=1}^{N-1} \sum_{k=j+1}^N \left| w_j - w_k \right| \right \rangle \\
& & \\
& = & {1 \over (N-1) \sum_{i=1}^{N} \left \langle w_i \right \rangle} \sum_{j=1}^{N-1} \sum_{k=j+1}^N \left \langle \left| w_j - w_k \right| \right \rangle \\

\end{matrix}
</math>
</center>

Temos, portanto, duas médias para calcular. A primeira é a riqueza média dos agentes, <math>\left \langle w_i \right \rangle</math>, que pode ser intuída pelo fato de que a distribuição inicial de riqueza é uniforme entre 0 e 1. Com isso, é claro que

<center>
<math>\left \langle w_i \right \rangle = {1 \over 2}</math>
</center>

A outra, mais trabalhosa, é a média sobre a diferença de riqueza entre todos os pares possíveis de agentes: <math>\left \langle \left| w_j - w_k \right| \right \rangle</math>.

Lembramos primeiro da definição de valor médio de uma variável aleatória contínua <math>x \;</math> sobre o intervalo <math>\left \lbrack 0, 1 \right \rbrack</math>:

<center>
<math>\left \langle x \right \rangle = \int \limits_{0}^{1} x \, P(x) \, dx</math>
</center>

sendo <math>P(x) \;</math> a função densidade de probabilidade associada à variável <math>x \;</math> sobre o intervalo. Analogamente,

<center>
<math>\left \langle x^2 \right \rangle = \int \limits_{0}^{1} x^2 \, P(x) \, dx</math>
</center>

Se <math>y \;</math> também for uma variável aleatória contínua, teremos

<center>
<math>
\left \langle xy \right \rangle = \int \limits_{0}^{1} \int \limits_{0}^{1} x \, P(x) \, y \, P(y) \, dx \,dy
</math>
</center>

Agora, temos que calcular <math>\left \langle \left| w_j - w_k \right| \right \rangle</math>, o que implica resolver

<center>
<math>
\left \langle \left| y - x \right| \right \rangle = \int \limits_{0}^{1} \int \limits_{0}^{1} \left| y - x \right| \, dx \,dy
</math>
</center>

Note que como a distribuição inicial de riqueza é uniforme e o intervalo de riqueza é <math>\left \lbrack 0, 1 \right \rbrack</math>, tanto <math>P(x) \,</math> como <math>P(y) \,</math> equivalem à unidade. O problema é que não podemos integrar diretamente o módulo no integrando. Temos que considerar os dois casos possíveis:

<center>
<math>
\left| y - x \right| =
\begin{cases}
y - x, \; \; \; y > x \\
x - y, \; \; \; y < x
\end{cases}
</math>
</center>

A integral dupla pode ser resolvida em duas partes, uma para cada caso. Consideremos então o caso <math>y > x</math>. Se a integral em <math>y</math> "varrer" todo o intervalo, a integral em <math>x</math> só poderá varrer o pedaço do intervalo entre 0 e <math>y</math>, pois somente assim satisfaremos <math>y > x</math> sempre. Portanto, o primeiro caso conduz a

<center>
<math>
\int \limits_{0}^{1} dy \int \limits_{0}^{y} \left( y - x \right) \, dx = {1 \over 6}
</math>
</center>

Usando um raciocínio análogo, reparamos que para o caso <math>x > y</math>, a integral em <math>x</math> deverá varrer valores suficientemente grandes, que sempre sejam superiores a <math>y</math>, conduzindo a

<center>
<math>
\int \limits_{0}^{1} dy \int \limits_{y}^{1} \left( x - y \right) \, dx = {1 \over 6}
</math>
</center>

Com isso,

<center>
<math>
\left \langle \left| y - x \right| \right \rangle = {1 \over 6} + {1 \over 6} = {1 \over 3}
</math>
</center>

Substituindo os resultados obtidos na expressão de <math>\left \langle G \right \rangle \;</math>, temos

<center>
<math>
\begin{matrix}

\left \langle G \right \rangle & = & {1 \over (N-1) \sum_{i=1}^{N} 1/2} \sum_{j=1}^{N-1} \sum_{k=j+1}^N 1/3 \\
& & \\
& = & {2 \over 3(N-1)N} \sum_{j=1}^{N-1} (N-j) \\
& & \\
& = & {2 \over 3(N-1)N} \left( N \sum_{j=1}^{N-1} - \sum_{j=1}^{N-1} j \right) \\
& & \\
& = & {2 \over 3(N-1)N} \left( N \left( N-1 \right) - {N \left( N-1 \right) \over 2} \right) \\
& & \\
& = & {2 \over 3(N-1)N} \left( {N \left( N-1 \right) \over 2} \right) \\
& & \\
& = & 1/3 \\

\end{matrix}
</math>
</center>

Note que usamos a "famosa fórmula" para a soma das bolinhas empilhadas em forma de triângulo:

<center>
<math>\sum_{i = 1}^{n} i = {n \left( n + 1 \right) \over 2}</math>
</center>

Modelo de gráfico em Gnuplot

2011-09-19T18:34:14Z

Tekkito: Criou página com 'Esta seção é dedicada à elaboração de um padrão para os gráficos produzidos em Gnuplot e "impressos" em formato EPS. Uma série de parâmetros podem ser ajustados, optimi...'

Esta seção é dedicada à elaboração de um padrão para os gráficos produzidos em Gnuplot e "impressos" em formato EPS. Uma série de parâmetros podem ser ajustados, optimizado a aparência e a legibilidade do gráfico. Atenção entretanto: ''o gráfico mostrado na janela do Gnuplot geralmente '''não corresponde exatamente''' à sua versão impressa em EPS!'' Veremos ocorrências disso mais adiante.

== <tt>linetype</tt>, <tt>linewidth</tt> & <tt>set size x,y</tt> ==

Usaremos como exemplo as funções <math>cos(x)</math>, <math>sin(x)</math>, <math>-cos(x)</math>, <math>-sin(x)</math>. Um plot das quatro funções no mesmo gráfico é obtido com:

set size 1,1
set grid
set yrange [-2:2]

plot cos(x) w l lt 1 lw 5 t "This is line type 1, thickness 5"
replot sin(x) w l lt 2 lw 4 t "This is line type 2, thickness 4"
replot -cos(x) w l lt 3 lw 3 t "This is line type 3, thickness 3"
replot -sin(x) w l lt 4 lw 2 t "This is line type 4, thickness 2"

Repare que:

'''<tt>w l</tt>''' = '''<tt>with lines</tt>'''

'''<tt>lt</tt>''' = '''<tt>linetype</tt>'''

'''<tt>lw</tt>''' = '''<tt>linewidth</tt>'''

'''<tt>set size 1,1</tt>''' = tamanho do gráfico em relação à "caixa original"

O gráfico gerado com isso é o seguinte:

<gnuplot>
set size 1,1
set grid
set yrange [-2:2]
plot cos(x) w l lt 1 lw 5 t "This is line type 1, thickness 5", sin(x) w l lt 2 lw 4 t "This is line type 2, thickness 4", -cos(x) w l lt 3 lw 3 t "This is line type 3, thickness 3", -sin(x) w l lt 4 lw 2 t "This is line type 4, thickness 2"
</gnuplot>

Vemos que '''<tt>linetype</tt>''' modifica a cor das curvas em questão, e que '''<tt>linewidth</tt>''' mexe com a espessura das curvas. Vejamos o que acontece se alterarmos '''<tt>set size 1,1</tt>''' para '''<tt>set size .75,.5</tt>''' :

<gnuplot>
set size .75,.5
set grid
set yrange [-2:2]
plot cos(x) w l lt 1 lw 5 t "This is line type 1, thickness 5", sin(x) w l lt 2 lw 4 t "This is line type 2, thickness 4", -cos(x) w l lt 3 lw 3 t "This is line type 3, thickness 3", -sin(x) w l lt 4 lw 2 t "This is line type 4, thickness 2"
</gnuplot>

Vemos que o a largura do gráfico diminuiu a 3/4 da unidade, e que a sua altura é agora 1/2 do tamanho original. Alterar esse parâmetro é um truque bastante útil quando se quer que o tamanho da fonte na legenda pareça maior em relação ao gráfico como um todo, sem ter que mexer diretamente no tamanho da fonte. Se fizermos '''<tt>set size 2,2</tt>''', certamente a fonte da legenda vai parecer muito pequena em relação ao todo! É o que vemos abaixo:

<gnuplot>
set size 2,2
set grid
set yrange [-2:2]
plot cos(x) w l lt 1 lw 5 t "This is line type 1, thickness 5", sin(x) w l lt 2 lw 4 t "This is line type 2, thickness 4", -cos(x) w l lt 3 lw 3 t "This is line type 3, thickness 3", -sin(x) w l lt 4 lw 2 t "This is line type 4, thickness 2"
</gnuplot>

== Imprimindo em EPS ==

O terminal Postscript do Gnuplot permite imprimir gráficos que contenham os detalhes que vemos nos gráficos acima, como diferença de cor e de espessura das curvas. Na verdade, o terminal PNG também permite visualizar essas características: de fato, os gráficos Gnuplot que são mostrados na Wikipédia são imagens PNG criadas com esse terminal. Entretanto, o terminal PNG perde uma informação importante nas curvas do gráfico: a ''estrutura'' da curva, ou seja, curva cheia, pontilhada, tracejada etc. O terminal Postscript permite caracterizar as curvas também pela sua estrutura. Para tanto, acessamos o terminal Postscript com as opções:

set terminal postscript enhanced dash color

Com isso, criamos gráficos que são "imunes" à impressões em branco e preto, no sentido de que se um periódico for totalmente impresso em branco e preto, sempre será possível identificar curvas diferentes nos gráficos pela sua estrutura, já que a informação da cor foi perdida. É claro que se o periódico for impresso a cores, ganhamos também a informação da cor para diferenciarmos as curvas além da estrutura.

Abaixo, temos um "screenshot" de uma imagem impressa em EPS com as opções enunciadas acima para o terminal Postscript. As demais características do gráfico são:

set size .5,.5
set grid
set yrange [-1.5:2]

plot cos(x) w l lt 1 lw 1 t "wave 1"
replot .4*cos(x+2.5) w l lt 2 lw 1 t "wave 2"
replot cos(x) + .4*cos(x+2.5) w l lt 3 lw 4 t "wave 1 + wave 2"

Repare no tamanho da fonte, que ficou maior já que diminuímos a caixa do gráfico (o retângulo tracejado) com '''<tt>set size .5,.5</tt>'''.

<gallery>
Imagem:Gaspar_Plotmodel.png|EPS immunity
</gallery>

== Modelo para artigo ==
Após seleção de figuras o ganhador é este:

[[Imagem:Modelo_paper.png]]

Ainda temos que modificar:
# Tirar Titulo (f=0.5)
# Eixo X a cada 0.2
# Legenda: wealth(0) = 0-10% (por exemplo)
# Tirar grade

Atualização do gráfico:

[[Imagem:Gaspar_Test_graph01.png]]

È possível fazer upload de imagens em EPS? ''--Gaspar''

Agora sim, ai vai um exemplo --[[Usuário:Sebas|Sebas]] 17:22, 29 Maio 2007 (BRT)

[[Image:Modelo_paper.eps]]

Porem, sempre ha, não vai diretamente para a figura. Não sei porque.
Mas, pelo menos podemos carregar eps.

Novo gráfico em png:
[[Imagem:Gaspar_06-MIN-10k-f10-risk2d.png|f=0.10, 10k passos]]

O original foi feito com set size = 0.5, 0.5
Há todo um espaço em branco por fora que não consigo remover via algum método automático. Abaixo, o link para o original em EPS:

[[Imagem:Gaspar_06-MIN-10k-f10-risk2d.eps]]

Note como as barras "colapsaram" para um mesmo valor após tempo o suficiente.

Distribuição de riquezas em função da aversão ao risco e capital inicial

2011-09-19T18:33:33Z

Tekkito: Criou página com '== 1. Introdução == Queremos acompanhar a riqueza <math>w_i</math> de alguns agentes de um sistema constituído de 300 agentes. Nessas simulações, o "saving" de cada agente ...'

== 1. Introdução ==

Queremos acompanhar a riqueza <math>w_i</math> de alguns agentes de um sistema constituído de 300 agentes. Nessas simulações, o "saving" de cada agente é mantido fixo ao longo do tempo ("quenched"), embora seja determinado no início de modo aleatório e unifome entre 0 e 1. Os resultados aqui apresentados se referem a amostras isoladas, sem calcular valores médios de qualquer tipo.

== 2. Método ==

Como a distribuição do "saving" é uniforme, acompanharemos a evolução de <math>w_i</math> de uma amostra de 20 agentes quaisquer do sistema (300 agentes), escolhidos no começo (por exemplo, acompanhamos <math>w_i</math> dos agentes de rótulo 10, 20, 30, ... , 190 e 200).
Em seguida, de posse do <math>w_i(t)</math> de cada um dos 20 agentes escolhidos, traçamos os respectivos gráficos de evolução de <math>w_i</math>.

Obs: notemos que o que importa é o valor do "saving" do agente, não o seu rótulo numérico.

== 3. Gráficos e discussão ==

Temos abaixo 3 gráficos para 3 valores do parâmetro <math>f</math>:

<gallery>
Imagem:Gaspar 2006-08-27-cashevol-f50-a-log-MINIMUM.jpg|f = 0.50
Imagem:Gaspar 2006-08-27-cashevol-f30-a-MINIMUM.jpg|f = 0.30
Imagem:Gaspar 2006-08-27-cashevol-f10-a-MINIMUM.jpg|f = 0.10
</gallery>

Para <math>f = 0.50</math>, vemos que o agente radical, com <math>1-saving=0.97</math>, apresenta uma flutuação muito maior do que os agentes mais conservadores. Entretanto, mesmo com uma flutuação variando cerca de 5 ordens de grandeza, o agente radical permanece no jogo mesmo após <math>10^4</math> passos de sistema (que é o tempo da simulação). Isso pode ser compreendido ao considerarmos a definição de probabilidade de vitória do agente mais pobre (apresentada em seções anteriores), que tende a 1 para <math>f=0.50</math> e para diferenças de riqueza muito grandes (como é o caso, nesse gráfico, para o agente radical). Portanto, enquanto a diferença e riqueza foi muito grande, o agente mais pobre praticamente ganha todas as apostas, recuperando sua riqueza em poucas transações.

Quando diminuímos <math>f</math>, notamos que há um aumento significativo na amplitude da flutuação dos agentes considerados. Nessas amostras, o que constatamos para valores menores de <math>f</math> é que os agentes radicais acabam sendo eliminados do jogo se esperarmos tempo o suficiente. Isso é bem visível para <math>f = 0.10</math>, em que apresentamos uma faixa mais completa de "savings". Notemos que parece haver um valor limitante de "saving", abaixo do qual agentes são eliminados, e acima do qual os agentes permanecem no jogo (agentes conservadores).

Isso deve estar relacionado com a própria regra do mínimo: agentes muito pobres não poderão, em geral, recuperar sua riqueza quando apostarem com agentes mais ricos que eles, pois há um limite (estipulado pela regra do mínimo) de riqueza que um agente poderá ganhar numa transação.

Variação de "w" na regra do mínimo (Gaspar)

2011-09-19T18:33:05Z

Investigações posteriores sobre autocorrelação (Gaspar)

2011-09-19T18:32:49Z

Tekkito: Criou página com '==Definição== Partimos da definição de autocorrelação de um sistema de agentes em relação ao estado do mesmo sistema em <math>t = \tau</math>: <center><math>corr(\tau,...'

==Definição==

Partimos da definição de autocorrelação de um sistema de agentes em relação ao estado do mesmo sistema em <math>t = \tau</math>:

<center><math>corr(\tau,t) = \frac{\sum_i^N w_i(\tau)w_i(\tau+t)}{\sum_i^N w_i(\tau)w_i(\tau)}</math></center>

sendo <math>N</math> o número de agentes e <math>w_i</math> a riqueza do <math>i</math>-ésimo agente. Também na forma na forma de produto escalar:

<center><math>corr(\tau,t) = \frac{\vec{W}(\tau) \cdot \vec{W}(\tau + t)}{\vec{W}(\tau) \cdot \vec{W}(\tau)}</math></center>

sendo <math>\vec{W}</math> o vetor que representa o sistema de agentes.

==Autocorrelação para distribuição uniforme entre 0 e 1==

Reconhecemos da definição acima que o denominador depende somente de <math>\tau</math>. Portanto, é esperado que ao longo da simulação, teremos para cada estado de referência <math>\vec{W}(\tau)</math> um denominador diferente, pois a distribuição de renda do sistema vai se alterando. Entretanto, o denominador será sempre o mesmo em <math>\tau = 0</math> pois trata-se do estado inicial do sistema, que é uma distribuição uniforme de riqueza entre 0 e 1. Pode-se mostrar que, para uma distribuição uniforme,

<center><math>\sum_{i=1}^N w_i^2 = \frac{N}{3}</math></center>

para <math>w_i = \lbrack 0,1 \rbrack</math> uniforme.

===Caso particular: quando um agente suga toda a riqueza===

Seja o <math>j</math>-ésimo agente o agente que vai concentrar toda a riqueza. Então,

<center>
<math>\lim_{t \rightarrow \infty} w_j(t) = W </math>
</center>

em que <math>W</math> representa a riqueza total do sistema. Como o sistema é conservativo, podemos usar o estado inicial do sistema para obter <math>W</math> de acordo com

<center>
<math>W = \sum_{i=1}^N w_{i}(0) = \frac{N}{2}</math>
</center>

A autocorrelação <math>t \rightarrow \infty</math> para <math>\tau = 0</math> vale, portanto,

<center>
<math>\lim_{t \rightarrow \infty} \left\langle corr(t) \right\rangle = \frac{3}{N} \lim_{t \rightarrow \infty} \sum_{i=1}^N \left\langle w_{i}(0) w_{i}(t) \right\rangle
</math>
</center>

Como a riqueza é concentrada em apenas um agente, sobrevive apenas a <math>j</math>-ésima parcela da soma, correspondente à riqueza do <math>j</math>-ésimo agente:

<center>
<math>\lim_{t \rightarrow \infty} \langle corr(t) \rangle=
\frac{3}{N} \lim_{t \rightarrow \infty} \langle w_{j}(0) w_{j}(t)\rangle =
\frac{3}{N} \langle w_{j}(0) \rangle \lim_{t \rightarrow \infty} w_{j}(t)
</math>
</center>

Mas
<center>
<math>\langle w_{j}(0) \rangle = 1/2</math>
</center>

de modo que
<center>
<math>
\lim_{t \rightarrow \infty} \langle corr(t) \rangle = \frac{3}{N} \frac{1}{2} \frac{N}{2} = \frac{3}{4}
</math>
</center>

==Denominador da autocorrelação em função de <math>\tau</math>==

Abaixo, gráficos do denominador em função de <math>\tau</math> e de <math>f</math> (parâmetro de assimetria) para um sistema de 300 agentes (média entre 100 rodadas, ou amostras) para as regras do mínimo e do perdedor. Observe como o denominador não é afetado por <math>f</math> para <math>\tau = 0</math>.

<gallery>
Image:Gaspar_2006-08-17-MINIMUM-corr_denom-tau(3).jpg|minimum
Image:Gaspar_2006-08-17-LOSER-corr_denom-tau(3).jpg|loser
</gallery>

=="Anomalias" da autocorrelação==

Abaixo, gráficos de <math>corr(f,\tau, t)</math> para um sistema de 300 agentes (média entre 100 rodadas, ou amostras), sendo um gráfico para cada <math>f</math>.
Para manter os gráficos limpos, eles são apresentados sem legendas.

'''Informações úteis:'''
<li>cada gráfico contem 41 curvas de autocorrelação;
<li>cada curva de autocorrelação é traçada em relação aos estados do sistema nos tempos <math>\tau = 0, 100, 200, 300, ... , 3900, 4000</math>;
<li>observe que, à medida que o sistema se aproxima do estado de equilíbrio, as curvas de autocorrelação se "empacotam" em torno de um certo valor;
<li>observe que, para a regra do perdedor, há um certo tempo durante o qual a autocorrelação é maior do que 1, mas que cai para um valor menor do que 1 ao se aproximar do estado de equilíbrio.
<li>observe que há uma curva de autocorrelação vermelha destacada de todas as outras: é a autocorrelação em relação ao estado inicial (<math>\tau = 0</math>), cujo valor assintótico visivelmente independe tanto de <math>f</math> como da regra em questão. Essa é uma anomalia ainda em estudo.

<gallery>
Image:Gaspar_2006-08-17-MINIMUM-corr-tau-f50.jpg|minimum, f = 0.50
Image:Gaspar_2006-08-17-MINIMUM-corr-tau-f40.jpg|minimum, f = 0.40
Image:Gaspar_2006-08-17-MINIMUM-corr-tau-f30.jpg|minimum, f = 0.30
Image:Gaspar_2006-08-17-MINIMUM-corr-tau-f20.jpg|minimum, f = 0.20
Image:Gaspar_2006-08-17-MINIMUM-corr-tau-f10.jpg|minimum, f = 0.10
Image:Gaspar_2006-08-17-LOSER-corr-tau-f50.jpg|loser, f = 0.50
Image:Gaspar_2006-08-17-LOSER-corr-tau-f40.jpg|loser, f = 0.40
Image:Gaspar_2006-08-17-LOSER-corr-tau-f30.jpg|loser, f = 0.30
Image:Gaspar_2006-08-17-LOSER-corr-tau-f20.jpg|loser, f = 0.20
Image:Gaspar_2006-08-17-LOSER-corr-tau-f10.jpg|loser, f = 0.10
</gallery>

=="Miolo" do programa utilizado para cálculo de autocorrelação==
Tendo em vista as anomalias apresentadas acima, apresentamos o pedaço do programa empregado para calcular a autocorrelação. Está sob análise pois, se há alguma real anomalia, deve haver um erro nessa parte do programa.

Breve legenda

B(agente) := vetor de $$$ dos agentes
B_TAU(agente,tau) := matriz com estado do sistema para cada tau
CORR_DENOM(tau) := vetor com denominadores da autocorrelação para
cada tau
MEAN_CORR(tempo,tau) := matriz com valores de autocorrelação para cada
tau e para cada tempo além de tau; no final, é
feita uma média entre todas as amostras para
obter melhor resolução

Cálculo do denominador no estado inicial:

CORR_DENOM(0) = DOT_PRODUCT(B,B)

DO j = 1, agents !!!salva o estado
B_TAU(j,0) = B(j) !!!inicial do
END DO !!!sistema

MEAN_CORR(0,0) = MEAN_CORR(0,0) + 1 !!!por def, a correlação
!!!no estado inicial = 1
...

Função que calcula autocorrelação:

corr = 0

DO j = 1, agents
corr = corr + B(j)*B_TAU(j,k)
END DO

corr = corr/CORR_DENOM(k)

Gini vs time

2011-09-19T18:32:06Z

Tekkito: Criou página com '= Tentativa de ajuste de curva = == 1. Resumo == Neste relatório, buscamos obter a forma funcional da curva de Gini gerada em uma simulação com as seguintes condições inici...'

= Tentativa de ajuste de curva =

== 1. Resumo ==
Neste relatório, buscamos obter a forma funcional da curva de Gini gerada em uma simulação com as seguintes condições iniciais:

<li>regra para transações: '''regra do perdedor'''
<li>número de agentes: '''100'''
<li>número de rodadas: '''9999'''
<li>número de passos de sistema em cada rodada: '''1010'''
<li>parâmetro(s) de assimetria <math>f</math> empregados: '''0.50'''

== 2. Tentativa de ajuste exponencial ==
Vemos abaixo a curva de Gini obtida de acordo com os parâmetros indicados logo acima.

[[Imagem:Gaspar_h1.jpg|Curva de Gini]]

Supomos que a curva de Gini, <math>G(t),</math> tenha dependência exponencial:

<center>
<math>G(t) = \alpha e^{\beta t} \;</math>
</center>

em que <math>\alpha</math> e <math>\beta</math> são constantes a determinar. Entretanto, reconhecemos que não se trata de uma exponencial decrescente tendendo a zero. O gráfico sugere que a forma da exponencial procurada deve conter uma constante aditiva <math>\gamma</math> para compensar a diferença entre o seu valor assintótico real e zero:

<center>
<math>G(t) = \alpha e^{\beta t} + \gamma\;</math>
</center>

Fazendo <math>H(t) \equiv \gamma - G(t)</math>, escrevemos

<center>
<math>H(t) = -\alpha e^{\beta t}\;</math>
</center>

Tomando o logaritmo natural lado a lado, obtemos

<center>
<math>\ln{H(t)} = \ln{-\alpha} + \beta t \;</math>
</center>

ou seja,

<center>
<math>\ln H(t) \propto t \; </math>
</center>

(Notemos que <math>\ln(-\alpha)</math> não é um problema, visto que <math>\alpha</math> ainda pode ser negativo).

A constante <math>\gamma</math> é o valor assintótico da curva de Gini mostrada no gráfico, e seu valor pode ser estimado por inspeção do gráfico (com mais resolução do que o mostrado aqui):

<center>
<math>\gamma \approx 0.519</math>
</center>

Tendo um valor aproximado para <math>\gamma \approx 0.519</math>, traçamos o gráfico de <math>H(t) \times t</math>:

[[Imagem:Gaspar_h5.jpg]]

Para verificarmos se há dependência exponencial, <math>\ln H(t)</math> deve ser proporcional a <math>t</math>, como mostrado acima. Portanto, tomamos a escala logarítmica

Comparação entre regra do mínimo e regra do perdedor (Gaspar)

2011-09-19T18:10:12Z

Tekkito: Criou página com '== 1. Regra do mínimo e regra do perdedor == Nas miniaturas, temos gráficos da relação entre os parâmetros investigados (o índice Gini, o índice <math>\theta</math> e a a...'

== 1. Regra do mínimo e regra do perdedor ==

Nas miniaturas, temos gráficos da relação entre os parâmetros investigados (o índice Gini, o índice <math>\theta</math> e a autocorrelação) com o parâmetro de assimetria <math>f</math> em duas variantes: pela a regra do mínimo e pela regra do perdedor. Na primeira, a quantidade de riquezas envolvida em cada transação é "limitada": os agentes põem em jogo uma fração da sua riqueza de acordo com o seu risco, e o vencedor da aposta recebe do perdedor uma quantia igual à menor das apostas, ou seja, o menor entre

<center>
<math>(\beta_i w_i, \beta_j w_j) \; </math>
</center>

sendo <math>\beta</math> o risco de um agente, e <math>w</math> a sua riqueza.

Na regra do perdedor, o agente vencedor da transação recebe tudo o que o agente perdedor apostou. Em termos pouco técnicos, com a regra do perdedor, um agente sempre tem a chance de "ganhar uma bolada", por mais pobre que ele seja. Na regra do mínimo, se um agente tem muito pouca riqueza, será muito difícil para ele vencer consecutivamente várias transações a ponto de poder acumular uma quantidade significativamente maior, visto que a quantidade de riqueza que ele pode adquirir em uma transação é limitada pela pela quantia que ele puder apostar.

==== Probabilidade de vitória do agente mais pobre ====

A probabilidade do agente mais pobre vencer uma transação é definida pela fórmula

<center>
<math>p(f,w_i,w_j) \equiv \frac{1}{2} + f \frac{|w_i-w_j|}{w_i+w_j}</math>
</center>

em que <math>w_i</math> e <math>w_j</math> são as riquezas dos <math>i</math>- e <math>j</math>-ésimos agentes, e <math>f</math> é denominado de ''parâmetro de assimetria,'' que pode assumir valores entre 0 e 0,5. O parâmetro de assimetria procura simular "ajuda externa" (''e. g.,'' programa de incentivo por parte do governo) aos indivíduos de baixa renda, visto que aumenta a probabilidade desses indivíduos vencerem uma transação. Além disso, vemos também que a diferença de riqueza entre os agentes da transação contribui com o favorecimento dos agentes mais pobres.

== 2. Não-estabilização de algumas curvas ==

Vale lembrar que o valor de cada um dos parâmetros é obtido em tempos suficientemente grandes da simulação: quanto mais tempo se passar, mais próximos estarão os parâmetros dos seus respectivos valores assintóticos (estes vindo a ser os valores procurados de fato). Entretanto, para parâmetros de assimetria pequenos (<math>f \lesssim 0.15</math>), o tempo necessário para que um parâmetro atinja um valor razoavelmente próximo do seu valor assintótico excede o tempo total da simulação. Como conseqüência, a discrepância entre o valor assintótico (que não foi atingido) e o valor medido para o parâmetro é apreciavelmente maior. Esse fato ocorre nos seguintes gráficos:

<li> <math>G(t) \times t</math> (regra do mínimo)
<li> <math>\theta(t) \times t</math> (regra do mínimo)
<li> <math>corr(t) \times t</math> com f = 0.20, 0.14, 0.08 (regra do mínimo)
----
As definições dos parâmetros são:

 <li><math>G = \frac{\sum_{i,j \atop i \ne j}^N |\omega_i - \omega_j|}{(N-1)\sum_i^N \omega_i}</math>

 <li> <math>\theta = \log \left(\langle \omega^2 \rangle - \langle \omega \rangle^2 \right) - \frac{2}{N}\sum_{j=1}^N \log \left( \omega_j + 10^{-20} \right)</math>

 <li> <math>corr(\tau,t) = \frac{\sum_i^N \omega_i(\tau)\omega_i(\tau+t)}{\sum_i^N \omega_i(\tau)\omega_i(\tau)}</math>

 sendo <math>N</math> o número de agentes do sistema, <math>w_i</math> a riqueza do i-ésimo agente, e <math>\tau</math> um certo tempo característico em relação ao qual será medida a autocorrelação do sistema, ou seja, medimos a autocorrelação em um certo tempo <math>t > \tau</math> em relação ao estado do sistema em <math>t = \tau</math>.

 Observação: a definição original de <math>\theta</math> apresenta uma média geométrica entre a riqueza dos agentes. A definição que apresentamos acima é na verdade o logaritmo decimal do parâmetro <math>\theta</math> original (é uma adaptação) por motivos de cálculo computacional. O cálculo de um produtório envolvendo fatores menores do que 1 (caso geral da riqueza dos agentes) conduz a um expoente cada vez menor, geralmente excedendo a resolução numérica do computador (underflow). Entretanto, consideramos essa adaptação válida pelo fato de buscarmos sempre o regime de estabilidade do parâmetro que, em princípio, não deve mudar com a transformação logarítmica.
----
Gráficos dos parâmetros em função do tempo para diversos <math>f</math>:
<gallery>
Imagem:Gaspar_2006-07-10-ginis-VS-t-MINIMUM.jpg|gini x t (minimum)
Imagem:Gaspar_2006-07-10-thetas-VS-t-MINIMUM.jpg|theta x t (minimum)
</gallery>
<gallery>
Imagem:Gaspar_2006-07-10-MINIMUM-corr-VS-t-f20.jpg|corr x t; f = 0.20(minimum)
Imagem:Gaspar_2006-07-10-MINIMUM-corr-VS-t-f14.jpg|corr x t; f = 0.14(minimum)
Imagem:Gaspar_2006-07-10-MINIMUM-corr-VS-t-f08.jpg|corr x t; f = 0.08 minimum)
</gallery>
----
Gráficos dos parâmetros medidos como função de <math>f</math> após o tempo de estabilização*:
<gallery>
Imagem:Gaspar_2006-07-10-GINI-loser-VS-minimum.jpg|gini x f
Imagem:Gaspar_2006-07-10-THETA-loser-VS-minimum.jpg|theta x f
Imagem:Gaspar_2006-07-10-CORR-loser-VS-MINIMUM.jpg|corr x f
</gallery>
----
Histogramas finais de distribuição de renda para regras do mínimo e do perdedor. Observe que, na regra do mínimo, as declividades na região de alta renda parecem não permanecer no valor -2 para nenhum f, opostamente à regra do perdedor (em que as declividades têm valores mais bem definidos):
<gallery>
Imagem:Gaspar_2006-09-03-histos-MINIMUM-gnulog.png|histo(minimum)
Imagem:Gaspar_2006-09-03-histos-LOSER-gnulog2.png|histo(loser)
</gallery>
----
Primeiras comparações entre regras com número de agentes diferentes:
<gallery>
Imagem:Gaspar_Zzzxxx-LOSER-MINIMUM-gini-f-1.jpg|gini x f
Imagem:Gaspar_Zzzxxx-LOSER-MINIMUM-theta-f-1.jpg|theta x f
Imagem:Gaspar_Zzzxxx-LOSER-MINIMUM-corr-f-1.jpg|corr x f
</gallery>

Figuras para o artigo

2011-09-19T18:09:11Z

Tekkito: Criou página com 'Rascunho de figuras para o artigo, algumas podem estar repetidas (estão na entrada principal de Tobias) == Regiões == <math>p_{ex}</math> probabilidade de favorecer o mais pob...'

Rascunho de figuras para o artigo, algumas podem estar repetidas (estão na entrada principal de Tobias)

== Regiões ==
<math>p_{ex}</math> probabilidade de favorecer o mais pobre, <math>\beta</math>: fração de saving o aversão ao risco.

=== <math>p_{ex}=0.9; \beta=0.9</math> ===

[[Imagem:p09-b09-cf95.jpg|444px|thumb|left|C=0.95]]
[[Imagem:p09-b09-cf05.jpg|444px|thumb|right|C=0.50]]
 

=== <math>p_{ex}=0.9; \beta=0.5</math> ===

[[Imagem:p09-b05-cf95.jpg|444px|thumb|left|C=0.95]]
[[Imagem:p09-b05-cf05.jpg|444px|thumb|right|C=0.50]]
 

=== <math>p_{ex}=0.7; \beta=0.9</math> ===

[[Imagem:p07-b09-cf95.jpg|444px|thumb|left|C=0.95]]
[[Imagem:p07-b09-cf05.jpg|444px|thumb|right|C=0.50]]
 

=== <math>p_{ex}=0.7; \beta=0.5</math> ===

[[Imagem:p07-b05-cf95.jpg|444px|thumb|left|C=0.95]]
[[Imagem:p07-b05-cf05.jpg|444px|thumb|right|C=0.50]]
 

=== <math>p_{ex}=0.5; \beta=0.9</math> ===

[[Imagem:p05-b09-cf95.jpg|444px|thumb|left|C=0.95]]
[[Imagem:p05-b09-cf05.jpg|444px|thumb|right|C=0.50]]
 

=== <math>p_{ex}=0.5; \beta=0.5</math> ===

[[Imagem:p05-b05-cf95.jpg|444px|thumb|left|C=0.95]]
[[Imagem:p05-b05-cf05.jpg|444px|thumb|right|C=0.50]]

Distribuição de Renda e Cooperação (Tobias)

2011-09-19T18:08:30Z

Tekkito: Criou página com '== Modelo == Um conjunto de <math>n</math> agentes descritos por 3 parâmetros trocam recursos (<math>w</math>), dois a dois, pela dinâmica de Monte Carlo. Os agentes podem ser ...'

== Modelo ==
Um conjunto de <math>n</math> agentes descritos por 3 parâmetros trocam recursos (<math>w</math>), dois a dois, pela dinâmica de Monte Carlo. Os agentes podem ser de dois tipos distintos descritos pelo dilema do prisioneiro, cooperadores (C) ou trapaceiros (D). Sendo assim podem haver 3 tipos de encontros, CC, CD (DC) e DD e cada um tem sua regra de troca de recursos.

'''Parâmetros:'''

Aversão ao risco (<math>b</math>) - os agentes não colocam todo o seu recurso em jogo, ao invés disso eles colocam somente quantias menores ou iquais a <math>w(1-b)</math>.

Probabilidade de favorecimento do pobre (<math>p</math>) - em um confronto CC esse parâmetro determina a probabilidade do agente mais pobre (menor <math>w</math>) ficar com os recursos postos em jogo.

Concentração de cooperados (<math>Cf</math>) - o sistema tem <math>nCf</math> agentes cooperadores e <math>n(1-Cf)</math> trapaceiros.

A quantidade de recursos posta em jogo é determinada pelo agente mais pobre, ou seja <math>w_p(1-b)</math> onde <math>w_p</math> é a riqueza do agente com menos recursos.

'''Regras de troca de recursos:'''

CC - Os dois agentes cooperadores promovem uma troca limpa, onde quem tem menor <math>w</math> tem probabilidade <math>p</math> de ganhar <math>w_p(1-b)</math> do agente mais rico, e este, probabilidade <math>(1-b)</math> de ganhar o mesmo <math>w_p(1-b)</math> do agente mais pobre.

CD (DC) - O agente trapaceiro leva todo recurso posto em jogo.

DD - Os dois agentes perdem os recursos postos em jogo e estes são imediatamente redistribuidos igualmente entre todos.

'''Em cada turno são feitos os seguintes passos:'''

-Dois agentes são escolhidos ao acaso;

-É determinado qual agente tem mais e menos recursos e o valor posto em jogo;

-De acordo com o comportamento dos escolhidos, é realizada a troca dos recursos;

-No caso de DD, os recursos perdidos são redistribuidos.

O número de turnos que o conjunto de agentes vai ser submetido é o suficiente para que a distribuição de riqueza estabilize sendo que inicialmente esta distribuição é homogênea entre 0 e 1.

Obs.: A regra de troca de recursos para encontros do tipo CC é, propositadamente, igual a um artigo anterior.

== Resultados ==

Inicialmente é interessante comparar a distribuição de riqueza entre este modelo, que inseriu agentes trapaceiros, e o modelo do artigo que o originou.

Distribuições de renda para vários valores de saving (<math>b</math>) e <math>p</math> para 3 valores de <math>Cf</math>. Observe que <math>Cf=1</math> reproduz o modelo utilizado no outro artigo.

<gallery>
Imagem:Regioes_bXp.jpg|Regiões estudas no artigo anterior
Imagem:dist_C1.jpg|Cf=1
Imagem:dist_C.75.jpg|Cf=0.75
Imagem:dist_C.5.jpg|Cf=0.5
</gallery>

----

'''Explorando a região <math>p<=0.5</math> e <math>Cf=0.95</math>'''

Os gráficos estão em escala logarítmica para facilitar a visualização.
<gallery>
Imagem:Mod7-p0.5-b0.9-cf0.95-rev1.jpg|p=0.5 b=0.9
Imagem:Mod7-p0.5-b0.5-cf0.95-rev1.jpg|p=0.5 b=0.5
Imagem:Mod7-p0.5-b0.1-cf0.95-rev1.jpg|p=0.5 b=0.1
</gallery>
<gallery>
Imagem:Mod7-p0.3-b0.9-cf0.95-rev1.jpg|p=0.3 b=0.9
Imagem:Mod7-p0.3-b0.5-cf0.95-rev1.jpg|p=0.3 b=0.5
Imagem:Mod7-p0.3-b0.1-cf0.95-rev1.jpg|p=0.3 b=0.1
</gallery>

'''Região <math>p>0.5</math> e <math>Cf=0.75</math>'''

Os gráficos estão em escala logarítmica.

<gallery>
Imagem:Mod7-p0.9-b0.9-cf0.75.jpg|p=0.9 b=0.9
Imagem:Mod7-p0.9-b0.7-cf0.75.jpg|p=0.9 b=0.7
Imagem:Mod7-p0.9-b0.5-cf0.75.jpg|p=0.9 b=0.5
Imagem:Mod7-p0.9-b0.3-cf0.75.jpg|p=0.9 b=0.3
Imagem:Mod7-p0.9-b0.1-cf0.75.jpg|p=0.9 b=0.1
Imagem:Mod7-p0.7-b0.9-cf0.75.jpg|p=0.7 b=0.9
Imagem:Mod7-p0.7-b0.7-cf0.75.jpg|p=0.7 b=0.7
Imagem:Mod7-p0.7-b0.5-cf0.75.jpg|p=0.7 b=0.5
Imagem:Mod7-p0.7-b0.3-cf0.75.jpg|p=0.7 b=0.3
Imagem:Mod7-p0.7-b0.1-cf0.75.jpg|p=0.7 b=0.1
</gallery>

----

'''Uma discussão sobre os comportamentos e a riqueza'''

Por simplicidade de notação constrói-se um espaço tridimencional onde cada coordenada é dada por um parâmetro e o vetor <math>\rho=(p,b,Cf)</math> denota uma dada combinação destes parâmetros.

Observe as tablelas abaixo:

{| border="1" cellpadding="2" cellspacing="0" align="center"
|-
! <math>\rho=(0.5,0.5,0.95)</math>
! style="background:#efefef;" | <math>w<0.5</math>
! style="background:#efefef;" | <math>w>0.5</math>
! style="background:#efefef;" | total
|-
! style="background:#efefef;" | C
| 0.8642154
| 0.0937847
| 0.9580002
|-
! style="background:#efefef;" | D
| 0.0000603
| 0.0419394
| 0.0419997
|-
! style="background:#efefef;" | total
| 0.8642758
| 0.1357242
| 1
|}

Normalizando as probabilidades para cada tipo de comportamento temos a chance de um agente ficar acima ou abaixo da linha de pobreza (<math>w=0.5</math>) para um determinado comportamento:

{| border="1" cellpadding="2" cellspacing="0" align="center"
|-
!
! style="background:#efefef;" | Pobre
! style="background:#efefef;" | Rico
|-
! style="background:#efefef;" | C
| 0.9021037
| 0.0978963
|-
! style="background:#efefef;" | D
| 0.0014357
| 0.9985643
|}

Comparando com <math>\rho=(0.5,0.5,1)</math> esta situação apresenta uma distribuição, mesmo entre os cooperadores, mais justa, o que é contra-intuitivo tendo em vista que adicionamos alguns trapaceiros a um grupo de cooperadores.

== Pequena alteração no modelo ==
Embora, como visto acima, a incerção de agentes trapaceiros torne mais justa a distribuição de riquezas para uma determinada região de parâmetros, esta mesma traz um incoveniênte: para todo o espaço dos parâmetros os agentes D tem melhores resultados, então, não haveria dilema algum e todos os agentes, se tiverem escolha, optariam por ser trapaceiros. Uma possível solução para este problema seria em vez de redistribuir os recursos perdidos pelos trapaceiros igualmente entre todos agentes, beneficiar agora, somente os cooperadores.

[[image:Mod7-p0.5-b0.5-cf0.30-mod1.jpg|300px|right|thumb|p=0.5 b=0.5 Cf=0.3]]

{| border="1" cellpadding="2" cellspacing="0" align="left"
|-
! <math>\rho=(0.5,0.5,0.3)</math>
! style="background:#efefef;" | <math>w<0.5</math>
! style="background:#efefef;" | <math>w>0.5</math>
! style="background:#efefef;" | total
|-
! style="background:#efefef;" | C
| 0.1616611
| 0.1283396
| 0.2900007
|-
! style="background:#efefef;" | D
| 0.6081919
| 0.1018074
| 0.7099993
|-
! style="background:#efefef;" | total
| 0.7698530
| 0.2301469
| 0.9999999
|}

{| border="1" cellpadding="2" cellspacing="0" align="center"
|-
!
! style="background:#efefef;" | Pobre
! style="background:#efefef;" | Rico
|-
! style="background:#efefef;" | C
| 0.5574507
| 0.4425493
|-
! style="background:#efefef;" | D
| 0.8566092
| 0.1433908
|}

Com essa alteração cada comportamento tem regiões onde é a melhor opção, em geral, quanto menor o número de agentes com o mesmo comportamento que um determinado indivíduo, melhor para este.

Para valores de <math>Cf</math> próximos de 1 os resultados não sofrem alterações significativas com esta alteração.

Artigos Relevantes

2011-09-19T18:07:47Z

Tekkito: Criou página com 'Artigo referente ao trabalho produzido por Romero, Lacasta e Sancho, em que é simulada a difusão de dímeros em superfícies compostas por átomos de Silicio. [http://scitation...'

Artigo referente ao trabalho produzido por Romero, Lacasta e Sancho, em que é simulada a difusão de dímeros em superfícies compostas por átomos de Silicio.
[http://scitation.aip.org/getabs/servlet/GetabsServlet?prog=normal&id=PLEEE8000069000005051105000001&idtype=cvips&gifs=yes]

Coolaboración con Miguel y Ricardo

2011-09-19T18:07:28Z

Tekkito: Criou página com ''''Primeros resultados''' * Efecto de la rampa de fuerza ( <math>\Delta f</math> ) 400px 1. Esto significa que la fricción medida depende de la velocida...'

'''Primeros resultados'''

* Efecto de la rampa de fuerza ( <math>\Delta f</math> )
[[Imagem:deltaf.png|400px]]

1. Esto significa que la fricción medida depende de la velocidad de variación de la fuerza.
2. Y que hasta ahora esa variación es demasiado rápida (el sistema no tiene tiempo para relajar).
3. Por otra parte es evidente para este caso que la fricción es muy chica (y va ser menor aún).
Así que tenemos que aumentar la interacción Adsorbato-Substrato.

----
'''Primer resultado potable'''

* Aumentando la interacción absorbato-substrato para <math>\epsilon_{A-S}=0.4</math> encontré el primer resultado de fricción
puramente "fonónica"

[[Imagem:eps0.4_fxv_1.png|500px]]

la fricción resultante es <math>f \approx 0.005 v</math>, siempre que <math>f < 0.0008</math>

Coolaboration with Claudio

2011-09-19T18:07:11Z

Tekkito: Criou página com '=== Role of anharmonicity and asimmetry of dimer potential === 1) Study anharmonic asimmetric potential: V(x)=0.5*k*x**2-beta*x**3/3 -Understand why anharmonicity destroys hyst...'

=== Role of anharmonicity and asimmetry of dimer potential ===

1) Study anharmonic asimmetric potential: V(x)=0.5*k*x**2-beta*x**3/3

-Understand why anharmonicity destroys hysteresis.
-Is there a resonance?
-Negative differential mobility near resonance?

Find parameters for Lennard-Jones interaction that mimicks the anharmonic
asimmetric potential

2) Study anharmonic simmetric potential
V(x)=0.5*k*x2+beta*|x3|/3

3) Double well potential
V(x)=alpha*x4-beta*x2

Transition between different states due to resonance

=== Excitation of internal vibrations (up to dissociation) by moving the substrate ===

1) 1D dimer

Consider "oscillating" substrate potential of the form
V_sub(x,w)=U*cos(2pi/a*(x+x0*sin(wt)))

Is there a resonance frequency w for which dimer internal motion is excited?

2) 2D cluster of nanoparticles on a substrate (connection with experiment of
Spasova: Co nanoparticles on metal substrate)

Idea: experimentally the nanoparticles form a triangular arrangement on the
substrate, but it is desirable to have a uniform coverage.
We could probably show that, by exciting a resonance in the cluster, the
nanoparticles could dissociate and be pushed away from each other, thus
balancing their tendency to clusterization.

At first we can study a model system of LJ particles and see what the
equilibrium configuration is and what happens if we move the substrate.

Potenciais de interação

2011-09-19T18:06:56Z

Tekkito: Criou página com 'Sobre os potencias de interacão nos modelos de atrito microscopico: Tipicamente entre as particulas deslizantes se usa potencias de interacao de Lennard-Jones ou harmonicos e o ...'

Sobre os potencias de interacão nos modelos de atrito microscopico:
Tipicamente entre as particulas deslizantes se usa potencias de interacao de Lennard-Jones ou harmonicos
e o sustrato é tipicamente é um potencial senoidal

'''Potencial Lennard-Jones'''

<math>V(r) = \epsilon\left[\left(\frac{\sigma}{r}\right)^{12} -
2\left(\frac{\sigma}{r}\right)^{6} \right]\,</math>

'''Derivada primeira'''

<math>\frac{dV(r)}{dr} = -\frac{12\epsilon}{\sigma}\left[\left(\frac{\sigma}{r}\right)^{13} -
\left(\frac{\sigma}{r}\right)^{7} \right]\,</math>

'''Derivada Segunda'''

<math>\frac{d^2V(r)}{dr^2} = \frac{12\epsilon}{\sigma^2}\left[13\left(\frac{\sigma}{r}\right)^{14} -
7\left(\frac{\sigma}{r}\right)^{8} \right]\,</math>

Esta última avaliada na posição de equilibrio <math>r=\sigma</math> é a constante elástica efetiva (''k'' da mola):

<math>k = \frac{72\epsilon}{\sigma^2}\,</math>

Numa cadeia 1D com interações de primeiros vecinhos a frequência maxima de vibração é

<math>\omega = \sqrt{2k/m} = 12\sqrt{\frac{\epsilon}{m\sigma^2}}</math>

O período de oscilação é então <math>\tau \approx 0.5 \sigma \sqrt{\frac{m}{\epsilon}}</math>.

'''Seguir'''
com relação de energias adsorbato-adsorbato/adsorbato-substrato, etc

Atrito de dímero em 2D

2011-09-19T18:06:28Z

Tekkito: Criou página com ''''Comportamento de um dímero num potencial periódico bidimensional com dissipação e força externa''' --''Trabalho de Mestrado de Italo Neide'' Objetivo desse tópico é ex...'

'''Comportamento de um dímero num potencial periódico bidimensional com dissipação e força externa'''
--''Trabalho de Mestrado de Italo Neide''

Objetivo desse tópico é expressar os resultados obtidos decorrentes do trabalho sobre o comportamento do dímero.
A seguir diversos passos que foram desenvolvidos e que estão sendo aperfeiçoados para o estudo do conteúdo proposto pelo tópico.

PARTÍCULA LIVRE EM POTENCIAL PERIÓDICO :

[[Imagem:imagem1.png]]

PARTÍCULA AMORTECIDA EM POTENCIAL PERIÓDICO :

A força relativa ao potencial usado aqui é da forma:

<center><math> F = -\frac{2{\pi}u_0}{l}\sin{(\frac{2{\pi}x}{l})} </math></center>

Aqui temos dois casos característicos para a partícula, ou ela não tem velocidade inicial suficiente para superar a primeira bareira do potencial (agora fica possível tal situação pois ela está com um termo dissipativo a mais na sua força resultante) ou ela vence a primeira e ou mais barreiras até que "estaciona" em algum ponto de mínimo do potencial periódico.

[[Imagem:imagem2.png]]

[[Imagem:imagem3.png]]

Agora os trabalhos tomam rumo com o dímero.

DÍMERO VIBRANDO EM SUA FREQUÊNCIA NATURAL DE VIBRAÇÂO

As simulações estão sendo feitas usando o modelo de um dímero, duas partículas presas por uma mola.As massas das partículas usadas nas simulações sempre serão normalizadas para 1Kg.A força é do tipo:

<center><math> F = k(x(1) -x(2) - B)\;</math></center>

Frequência natural de vibração de um dímero:

Sabe-se que a frequência de vibração w, de um bloco com massa m, preso por uma mola com constante elástica k, e a mola com sua outra extremidade presa a um ponto fixo é dada por:

<center><math> \omega = \sqrt{\frac{k}{m}}</math></center>

E que para qualquer movimento harmônico:

<center><math> \omega = \frac{2\pi}{T}</math></center>

Como temos um dímero, a relação pode ser interpretada da seguinte maneira, imagine que no sistema bloco mola temos a mola presa em outro sistema bloco mola, e não em um ponto fixo, resultando em um dímero.Isso leva ao fato de que o bloco agora sente o dobro da força que sentia antes, como a única influência sobre o bloco é a mola, pode-se concluir que agora k tenha duplicado, portanto temos que:

<center><math> \omega_{particula no dimero} = \sqrt{\frac{2k}{m}}</math></center>

[[Imagem:imagem4.png]]

<math> T \approx 0.44</math>

<math> \omega \approx 14.27</math>

<math> \omega \approx \sqrt{2k}</math>

<math> \omega \approx 14.24</math>

DÍMERO AMORTECIDO EM POTÊNCIAL PERIÓDICO

Abaixo alguns exemplos das simulações feitas.Para todas o parâmetro de rede l = 1.0

[[Imagem:imagem5.png]]

[[Imagem:imagem6.png]]

Próximo passo foi reproduzir o artigo feito por Sebastian et al. Aplicando uma força no dímero, fazendo essa força variar e analisando seu comportamento de acordo com os valores em que sua velocidade fica estável chega-se a um resultado bastante interessante.Abaixo o gráfico de um caso possível:

[[Imagem:imagem8.png]]

O que temos aqui é uma força aplicada no dímero, fazendo com que o dímero nunca pare, mesmo com o amortecimento, pois esse amortecimento está relacionado com a velocidade da partícula, ou seja com o passar do tempo a velocidade do dímero estabiliza em um valor e é esse o valor plotado no gráfico acima.
No inicio a força aplicada é pequena ainda para fazer o dímero se mover, pois ele esta preso no potencial do substrato.Conforme a força vai aumentando ele começa a se mover, até chegar a uma velocidade ressonante com a velocidade de vibração própia do dimero, portanto nesse ponto a força aplicada começa a parar de 'empurrar' o dímero.Essa força externa começa a contribuir para a vibração do dímero.Ela vai aumentando até chegar num ponto em que a força volta a ser transferida para a velocidade de translação do dímero, e portanto temos um salto na velocidade do dímero, e dae por diante ele entra no regime linear.

DÍMERO EM DUAS DIMENSÕES

Abaixo o primeiro gráfico obtido para as velocidades do centro de massa, de vibração e de rotação do dímero.

[[Imagem:imagem9.png]]

O ESTUDO EM DUAS DIMENSÕES ESTÁ DIVIDIDO EM DUAS PARTES : 1 - ANALISE DE SEU TEMPO DE PARADA, 2 - ANALISE COM FORÇA APLICADA.O ITEM 2 ESTÁ SUB-DIVIDIDO EM OUTRAS DUAS PARTES : 2a - INFLUENCIA DA ORIENTAÇÃO, 2b - FORÇA EM FASE COM GIRO

1 - ANALISE DE SEU TEMPO DE PARADA

Para compreender esse novo sistema sem refências bibliográficas foi desenvolvido um estudo detalhado de sua dinâmica.Abaixo alguns gráficos desta parte do trabalho:

Neste caso (casoa)temos o dímero posicionado com seu eixo 'em cima' de dois qualquer máximos de potencial da rede com menor distância(ou dois átomos mais próximos):

[[Imagem:imagem22.png]]

O dímero começa com uma velocidade inicial e obtem-se o tempo de parada para cada uma de suas componentes da velocidade: velocidade translacional, velocidade vibracional e velocidade rotacional.
Foi definido que quando o dímero começa a voltar contra o sentido em que ele inicialmente se movimentava, sua velocidade transalacional cessa.O mesmo vale para a rotacional, quando ele começa a girar no sentido inverso ao que ele girava inicialmente é considerado que ele não gira mais.Esta definido que a velocidade vibracional cessa quando alguns valores muito baixos de oscilação em relação aos de translação e ou rotação são alcançados pelos átomos do dímero.Convém notar aqui que não existe como seguir a RISCA esses parâmetros de parada, cada caso analisado é um caso diferente de dinâmica e são usados algoritoms diferentes para determinar a parada de cada velocidade.A cada novo tipo de dinâmica analisado surgiam mais possibiladades de ocorrências que fugiam do algoritmo anterior e então era necessário reformula-lo totalmente, adaptando para esse novo tipo de dinâmica.
Os pontos nesse gráfico são referentes a várias posições iniciais diferentes do dímero. No eixo y temos o ângulo de orientação do dímero, definido como o ângulo que o eixo do dímero faz com os dois átmos inicias da simulação(auto explicativo na figura anterior).Ou seja , temos ele variando de 0 a <math>\pi</math>.No eixo x temos o tempo de parada de cada velocidade.

[[Imagem:imagem10.png]]

Neste caso (casoa) percebe-se uma simetria no eixo y , e temos dois maximos para tempo de parada da velocidade de oscilação.É esperado que com velocidades iniciais mais altas a orientação do dímero em relação ao substrato não exerça influência na dinâmica do dímero e portanto, seus tempos de parada seja iguais, o que é observado para a figura seguinte, com velocidade 100 vezes maior que a figura do caso anterior.Convém ressaltar que a velocidade inicial é aplicada sempre no mesmo sentido, e nos dois casos, no sentido como a figura inicial propõe (figura auto explicativa).Percebe-se que o dímero não ganha rotação apenas com aplicação de velocidades desse tipo.

[[Imagem:imagem11.png]]

Vale ressaltar que casoa e casob estão com comensuração = 1.
Abaixo os mesmos casos com comensuração = 0.5

[[Imagem:imagem12.png]]

[[Imagem:imagem13.png]]

Sem comensuração observa-se que as velocidades não tomam muitas características que dependem da orientação.

Agora outra dinâmica, abaixo figura autoexplicativa:

[[Imagem:imagem23.png]]

Aqui nessa dinâmica o dímero é lançado com velocidade igual a da sua orientação.

Após simulação:

[[Imagem:imagem14.png]]

Agora outra dinâmica, abaixo figura auto-explicativa:

[[Imagem:imagem24.png]]

Aqui o dímero é lançado com uma velocidade de giro mais uma velocidade de translação.

[[Imagem:imagem15.png]]

Sem comensuração (C = 0.5):

[[Imagem:imagem16.png]]

Aqui nesse ponto ainda não é possível observar que a adição de um giro inicial faça alguma diferença na dinâmica do dímero.Ponto importante é saber se agora em duas dimensões a possibilidade do dímero rotar tem alguma influência decisiva na sua dinâmica, diferente do que acontecia em uma dimensão.
Primeiro ponto lógico a ser analisado é o caso acima para velocidades de giro variando, com velocidade de centro de massa fixas, abaixo duas simulações, uma com vcm : 15 e outra com vcm : 100

[[Imagem:imagem17.png]]

[[Imagem:imagem18.png]]

Perecebe-se que conforme vcm vai ficando maior que a velocidade de rotação ela vai se tornando uma linha 'comportada'.Alguns pontos a primeira vista parecem não estar certos, mas eles foram revisados e vistos um por um.
Agora é o ponto decisivo para ver se esse novo grau de liberdade, a rotação, tem alguma influência.O passo decisivo consiste em analisar um gráfico em que toma os tempos de parada para o dímero unidimensional com velocidades iniciais variando.Compara-lo com o bidimensional (que ficou idêntico), e compara-lo novamente com o bidimensional adicionado de velocidades de rotação com varios valores.Abaixo o gráfico obtido para diversas simulações(comparações foram feitas em escalas logaritmas para melhor verificação):

[[Imagem:imagem20.png]]

[[Imagem:imagem21.png]]

[[Imagem:imagem19.png]]

Aqui perecebe-se que para velocidades tangenciais(ou rotacionais, ou vt) suficientemente altas temos um tempo de parada diferente ao comparado sem velocidade rotacional inicial.Ficou provado então que existe influências decisivas diferentes comparadas ao caso unidimensional em relação ao giro.Neste ponto pareceu ser mais satisfatorio mudar o rumo do trabalho para estudar com forças aplicadas.

2 - ANALISE COM FORÇA APLICADA

2a - INFLUENCIA DA ORIENTAÇÃO

Nesse tópico foi pesquisado a influência de forças aplicadas no dímero, e como a comensuração afeta essa dinâmica.
As simulações foram feitas colocando o dímero com orientação variando de 0º a 90º, com uma força aplicada.Foi plotado os gráficos de velocidade final do dímero(velocidade estacionária) por orientação.A caomensuração está definida como b/l, em que b é a distância de equilíbrio do dímero e l é a distância entres dois máximos do potencial da rede, ou seja, o parâmetro de rede para uma célula simples.As simulações foram feitas para comensurações igual a 0.5, 1.0, 1.25 e 1.5.
Uma prévia da figura para um regime de forças altas seria uma mesma velocidade final para qualquer orientação inicial do dímero.Relembrando da figura em que é plotado a força variando, temos que o valor de quando a força entra em fase com a energia de vibração do dímero é de aproximadamente 0.5, portanto espera-se que para esse valor de força exista algum efeito diferente do normal.

b/l = 1.0

[[Imagem:imagem25.png]]

Aqui observa-se que conforme a orientação do dímero varia ele assume dois estados.Ou ele se movimenta e ganha velocidade, ou permanece parado.Isso existe pelo fato de existirem melhores 'caminhos' para seu movimento, ou poços de potenciais relativos a sua dinâmica.A direção da força na figura acima é de 0º, abaixo gráficos para forças nas direções de 22,5º, 45º, 67,5º , e 90º respectivamente:

[[Imagem:imagem26.png]]

[[Imagem:imagem27.png]]

[[Imagem:imagem28.png]]

[[Imagem:imagem29.png]]

Fica evidente que a direção da força aplicada não tem influência sobre em qual ângulo o dímero vencerá a barreira de potencial e ira se movimentar ou quando ele não consegue mais se mover.Os ângulos limites para o incio do movimento são mantidos.
Abaixo agora uma analise sobre a vizinhança dessa força.Forças aplicadas são de 0,3 , 0,4 e 0,6 respectivamente, todas na direção 0º.

[[Imagem:imagem30.png]]

[[Imagem:imagem31.png]]

[[Imagem:imagem32.png]]

As figuras acima mostram que a amplitude da velocidade estacionária do dímero vai decaindo conforme a força aplicada é diminuida, mas preserva-se o aspecto 'montanha' visto para com força aplicada igual a 0.5.Na ultima figura com força igual a 0.6 já não existe mais o efeito 'montanha', pois a força ja é suficiente para vencer qualquer potencial do substrato independente da orientação do dímero.
Agora resto o estudo com comensurações diferentes.Estão a baixo respectivamente figuras e comentários para b/l = 0.5, 1.25 e 1.5.

b/l = 0.5

[[Imagem:imagem33.png]]

A força é aplicada na direção de 0º, logo temos que quando o dímero está próximo da orientação de 90º ele começa a sentir mais os máximos do potencial da rede, e portanto enfrenta mais dificuldade para entrar em movimento, camaremos isso de efeito de borda de força.A prova para esse argumento é que se a força fosse aplicada numa direção de 90º teriamos esse efeito logo na proximidade de uma orientação de 0º, que é o que acontece abaixo:

[[Imagem:imagem34.png]]

Agora diminuindo a intensidade da força para 0.3 observam-se duas coisas, uma esperada, que é a diminuição da amplitude da velocidade estacionaria, e um segundo ponto interessante, o ângulo limite aumenta para a força aplicada na direção de 90º, ou seja , conforme se diminui a força temos um aumento no ângulo limite.Abaixo o gráfico.

[[Imagem:imagem35.png]]

Portanto com essa comensuração observamos o efeito de o ângulo limite "caminhar" conforme a força varia nas proximidades de 0.5.

b/l = 1.25

[[Imagem:imagem36.png]]

Aqui temos caso semelhante ao de comensuração = 0.5.Para analisar se o caso dessa queda do ângulo limite ser igual ao caso do l/b = 0.5 foi simulado para uma força aplicada no sentido de 90º.

[[Imagem:imagem37.png]]

As dinâmica das duas comensurações são bem parecidas, agora a analise de uma força menor aplicada para ver se o ângulo limite 'caminha' também.

[[Imagem:imagem38.png]]

E portanto o que no inicio parecia ser igual ao caso b/l = 0.5 provou-se ser totalmente diferente, agora com aspecto de 'montanha'.Indo mais longe a analisando o que acontece com uma força de 0,35 aplicada, obtem-se a seguinte figura:

[[Imagem:imagem39.png]]

E assim por diante diminuindo a força a dinâmica permanece a mesma.Temos então com essa comensuração algo que no inicio tem o efeito de borda de força, depois assume a característica de 'montanha' e finalmente cessa todo e qualquer movimento para forças logo abaixo de 0.4.

b/l = 1.5

[[Imagem:imagem40.png]]

Só com esse gráfico fica dificíl comentar muitas coisas sólidas, abaixo com força aplicada no sentido de 90º

[[Imagem:imagem41.png]]

Aqui observa-se uma troca simetrica de maximos e mínimos mantendo-se aproximado o ângulo limite.Primeira dinâmica em que o dímero surge com mais de uma possibiladade de entrar em movimento, até agora para o intervalo de estudo (0 - pi/2) não havia se observado essa característica ainda.Variando a força temos:

[[Imagem:imagem42.png]]

Portanto dinâmica parecida com a comensuração de 1 e temos diminuição da amplitude da velocidade estacionaria conforme se diminui a força aplicada.Agora um ponto interessante, de acordo com o estudo das outras comensurações surge uma duvida, conforme ia se diminuindo a força dependendo do caso, o ângulo limite 'andava' b/l = 0.5 , não fazia nada b/l = 1 ou cessava completamente b/l = 1.25.O que deve acontecer nesse caso, já que temos duas 'montanhas'? a figura abaixo responde essa pergunta:

[[Imagem:imagem43.png]]

A primeira 'montanha' não existe mais e a segunda continua com o ângulo limite igual. Portanto, nesse caso temos algo muito parecido com o de comensuração igual a 1, só que com duas 'montanhas'.Ao fazer mais simulações, observa-se que a uniformidade é apenas mantida para número inteiros de b/l, sempre retornando o aspecto de 'montanha' até que a distância de equilíbrio (b) se torna tão grande que não se observa mais efeitos do substrato relativo a comensurações.

Moodle

2011-09-19T17:59:19Z

Tekkito: Criou página com '== Instalação e configuração de Moodle == '''Receita sugerida''' E recomendável baixar o pacote completo (tar.gz) do site oficial de Moodle, pelas seguintes razões: # Ass...'

== Instalação e configuração de Moodle ==

'''Receita sugerida'''
E recomendável baixar o pacote completo (tar.gz) do site oficial de Moodle,
pelas seguintes razões:

# Assim temos a última versão estável. No momento de escrever isto, era a 1.8..., em tanto que no repositório debian estava a 1.6
# Na 1.8 formulas, tex, etc formam parte dos módulos standar, nas anteriores devem ser instalados depois.
# A configuração via web é tranqüila. A partir do apt-get, debconf, ... não funcionou na última vez.

Para detalhes ir a:

[http://docs.moodle.org/en/Installing_Moodle#Download_and_copy_files_into_place Instalando Moodle]

De todas formas deixo registro da instalação anterior:

'''Receita anterior'''

# Instalar e configurar ''mysql'' (criar base mysql e usuario root)
# Baixar ''moodle'' e extras (preferentemente com ''apt'' ou ''synaptic)''
# Responder cuidadosamente as perguntas de configuração:
# Servidor web -> ''apache2''
# Base de dados -> ''mysql''
# Servidor de mysql -> 127.0.0.1 (=localhost, porem deu problema com firefox porque não tinha configurado o localhost)
# Usuario de mysql -> root
# Senha -> a senha dada ao criar o usuario root em mysql
# Usuario de base moodle -> root
# Senha -> mesma é melhor
# Editar ''/etc/moodle/config.php'' e configurar server name (wwwroot) davinci.if.ufrgs.br por exemplo
<pre>
<?
# This file has been generated by debconf
# You can find a commented config file in /usr/share/doc/moodle/

/usr/share/doc/moodle/

unset($CFG);

$CFG->dbtype = 'mysql';
$CFG->dbhost = '127.0.0.1';
$CFG->dbname = 'moodle';
$CFG->dbuser = 'root';
$CFG->dbpass = 'xxxxxx';
$CFG->prefix = 'mdl_';

$CFG->dbpersist = "false";

$CFG->wwwroot = 'http://davinci.if.ufrgs.br/moodle';
$CFG->dirroot = '/usr/share/moodle';
$CFG->dataroot = '/var/lib/moodle';
$CFG->directorypermissions = 0750;
$CFG->admin = 'admin';

$CFG->respectsessionsettings = true;

if (file_exists("$CFG->dirroot/lib/setup.php")) { // Do not edit
include_once("$CFG->dirroot/lib/setup.php");
} else {
if ($CFG->dirroot == dirname(__FILE__)) {
echo "Could not find this file: $CFG->dirroot/lib/setup.php";
echo "Are you sure all your files have been uploaded?";
} else {
echo "Error detected in config.php";
echo "Error in: \$CFG->dirroot = '$CFG->dirroot';";
echo "Try this: \$CFG->dirroot = '".dirname(__FILE__)."';";
}
die;
}
?>
</pre>
12. editar ''/etc/moodle/apache.conf'' e mudar '''/moodle''' para '''/usr/share/moodle/''' (se é ai que está instalado)
e habilitar para todo mundo ('''deny from nobody''')
<pre>
<DirectoryMatch /usr/share/moodle/>

Options +FollowSymLinks
AllowOverride None

order deny,allow
deny from nobody
# allow from 127.0.0.0/255.0.0.0
# allow from 143.54.109.67
# allow from all
</pre>

* [http://docs.moodle.org/en/Teacher_documentation#Editing_course_section Teacher Docs]

Backup de servidor wiki

2011-09-19T17:59:06Z

Tekkito: Criou página com 'Exemplo de comando para salvar toda a base de dados numa maquina Debian GNU/Linux rodando mysqldump: mysqldump --user=USERNAME --password=PASSWORD --single-transaction --all-da...'

Exemplo de comando para salvar toda a base de dados numa maquina Debian GNU/Linux rodando mysqldump:

mysqldump --user=USERNAME --password=PASSWORD --single-transaction --all-databases > BACKUPFILE.sql

Se for so da base wiki "wikidb"
mysqldump -u root -p wikidb > wikidb.sql

e para recuperar essa base "wikidb"
mysql -u root -p -B wikidb < wikidb.sql

Tirado de [http://meta.wikimedia.org/wiki/Database_dump_and_restore_examples Database dump and restore examples]
Para mais dados ir a [http://meta.wikimedia.org/wiki/How_to_move_a_MediaWiki_Database How to move a wiki to another server]

Fórmulas

2011-09-19T17:58:52Z

Tekkito: Criou página com 'Para escrever Fórmulas em Wikimedia se usa o markup: <pre><nowiki><math> </math></nowiki></pre> para delimitar conteudos com fórmulas, que seguem, sem todas as funcionalidades...'

Para escrever Fórmulas em Wikimedia se usa o markup:
<pre><nowiki><math> </math></nowiki></pre>
para delimitar conteudos com fórmulas, que seguem, sem todas as funcionalidades, a sintaxe do LaTex.

'''Nota:''' o markup acima pode ser chamado com o botão de edição <math>\sqrt n</math>

Exemplo:
<pre><nowiki><math>x^2 = \alpha </math></nowiki></pre>

resulta em:

<math>x^2=\alpha</math>

Em geral o wikimedia minimiza recursos para mostrar as fórmulas no navegador; assim, sempre que for posível, vai mostrar em html extendido como no exemplo acima (isso se traduze em fórmulas de baixa qualidade e/ou fontes pequenas).
Fórmulas sofisticadas são tranformadas em figuras, resultando em qualidade de livro texto.
Contudo, se queremos forzar essa qualidade em fórmulas simples e só colocar \, ao final da formula, como no exemplo abaixo:

<pre><nowiki><math>x^2 = \alpha\,</math></nowiki></pre>

e voilá:

<math>x^2 = \alpha\,</math>

[http://meta.wikimedia.org/wiki/Help:Formula#Pros_of_HTML mais detalhes]

Modificando a barra de navegação (navigation bar)

2011-09-19T17:58:37Z

Tekkito: Criou página com '* Fazer login como Administrador (Sysop, WikiSysop ou usuário com permissão de Admin) * Ir a URL/wiki/index.php/MediaWiki:Sidebar * Editar, modificar e salvar (é aconselhável...'

* Fazer login como Administrador (Sysop, WikiSysop ou usuário com permissão de Admin)
* Ir a URL/wiki/index.php/MediaWiki:Sidebar
* Editar, modificar e salvar (é aconselhável comentar os link desativados com
<pre>

</pre>

Mudando o titulo da página principal (mainpage)

2011-09-19T17:58:19Z

Tekkito: Criou página com '* Fazer login como Administrador (Sysop, WikiSysop ou usuário com permissão de Admin) * Ir a Página principal -> URL/wiki/index.php * Abrir aba "mover" * Mudar o nome da "Pág...'

* Fazer login como Administrador (Sysop, WikiSysop ou usuário com permissão de Admin)
* Ir a Página principal -> URL/wiki/index.php
* Abrir aba "mover"
* Mudar o nome da "Página principal"
* Ir para URL/wiki/index.php/MediaWiki:Mainpage
* Editar e mudar nome para o MESMO colocado no passo 4

Mysql

2011-09-19T17:57:59Z

Tekkito: Criou página com '== Entrando no ''mysql'' == mysql -u root -p el pede a senha de root que deve ter sido configurada depois da instalação == Dentro do ''mysql'' == mysql> status mostra o ...'

== Entrando no ''mysql'' ==

mysql -u root -p

el pede a senha de root que deve ter sido configurada depois da instalação

== Dentro do ''mysql'' ==

mysql> status

mostra o status da base

mysql> show databases;

para ver as bases definidas

mysql> use moodle;

para trabalhar com uma base (''moodle'' por exmplo)

Dentro da base, podemos ver a tabelas com:

mysql> show tables;

e checar uma dada com (mdl_scale por exemplo é uma tabela do ''moodle'')

mysql> check table mdl_scale;
mysql> repair table mdl_scale;

a segunda linha é para repará-la, caso esteja corrupta (o que impediria o BACKUP externo)

== Backup de uma base ==

Na linha de comando, ou seja sem entrar no ''mysql''

mysqldump -u root -p moodle > moodle.sql

faz o backup da base do ''moodle'' por exemplo

mysql -u root -p moodle < moodle.sql

recupera a base "moodle" do backup "moode.sql"

== Backup completo ==

Isto é para fazer backup de todas as bases gerenciadas pelo ''mysql''

mysqldump -u root -p --all-databases > alldb_data.sql

Distribuição Docente 2007

2011-09-19T17:57:38Z

Tekkito: Criou página com 'O gráfico abaixo mostra a distribuição de carga horária definida pelo Departamento de Física para o semestre 2007/1, incluindo a atuação na Pós-graduação em Física. O...'

O gráfico abaixo mostra a distribuição de carga horária definida pelo Departamento de Física para o semestre 2007/1,
incluindo a atuação na Pós-graduação em Física.

O eixo horizontal representa a carga horária semanal.
No caso de horário noturno, após 20:30, cada hora aula contabiliza como 1,5.
Cada ajuda de laboratório conta como 1 (uma) hora semanal ---em ''history'' se pode consultar a primeira versão dos gráficos com
as horas de laboratório sem re-normalizar.

O eixo vertical representa o número de professores com uma determinada carga.

O total de horas de aula é de 457 (incluindo pós de física e contando como 1 hora as ajudas de laboratório) e o total de professores é 64. Só foram considerados os professores do quadro permanente, com encargos para o semestre.

A carga horária semanal média resultante é de 7,14h e está indicada pela linha azul.

A esquerda os gráficos incluindo os 64 professores, a direita o professor com carga de 24h semanal (Estágio Supervisionado em Ensino de Física=20h) é omitido o
que da uma média de horas semanal de 6,87

[[Imagem:distridoc2007.png|400px|left|thumb]]
[[Imagem:distridoc2007B.png|400px|right|thumb]]

Qual será o critério por trás desta peculiar distribuição?
Será que ela está correlacionada com a produção acadêmica--número de artigos publicados por exemplo?
Abaixo mostramos a carga horária de cada professor no eixo vertical, e a correspondente produção acadêmica (período 2002-2006, obtido do web of science, cv lattes e SABI)
no eixo horizontal:

[[Imagem:horas_x_artigos.png|400px|left|thumb]]
[[Imagem:horas_x_artigosB.png|400px|right|thumb]]

O coeficiente de correlação linear (-0.26) indica que esse não foi o critério utilizado.

Os dados completos utilizados estão na tabela seguinte e foram obtidos da página do departamento: http://www.if.ufrgs.br/fis/distrib/distrib_20071.html

[[Media:lis_doc_2007-0.png]]
[[Media:lis_doc_2007-1.png]]

Agradeço a Marcia Barbosa pelos rápidos esclarecimentos e a Fernando Lang pelas sugestões na apresentação dos dados.

--[[Usuário:Sebas|Sebas]] 12:00, 23 Fevereiro 2007 (BRT)

<analytics uacct="UA-379257-2"></analytics>

Colocar legendas num video

2011-09-19T17:57:15Z

Tekkito: Criou página com 'Esta é uma extração de mínimo conteúdo para conseguir adicionar (hardcoded) legendas num vídeo A situação típica é ter um vídeo em formato avi o mpeg baixado ou extrai...'

Esta é uma extração de mínimo conteúdo para conseguir adicionar (hardcoded) legendas num vídeo
A situação típica é ter um vídeo em formato avi o mpeg baixado ou extraido do qual não
temos as legendas na língua de nosso interesse e queremos gravar um CD ou DVD para assistir na TV
com as legendas.

Receita:

* baixe a legenda na língua desejada (arquivo .srt ou .sub) e salve no mesmo diretório do vídeo (.avi ou .mpeg)
* verifique que o ''mencoder'' está instalado (ele é usado pelo ''mplayer'')
* verifique a sincronia das legendas passando o vídeo com mplayer (ele permite carregar o arquivo de legendas)
* se tudo da certo então faça a versão .avi com as legendas "hardcoded" assim

mencoder filme.avi -o filme_legendado.avi -sub legendas.sub -ovc lavc -oac mp3lame

* filme.avi: o filme original sem legendas
* filme_legendado.avi : nome atribuído ao arquivo de saída gravado com legendas
* legendas.sub: é o arquivo de legendas
* -ovc e -oac indicam os codificadores de vídeo e som respectivamente

Pronto o "filme_legendado.avi"esta listo para ser gravado no CD ou DVD (use o k3b)

'''Atenção:''' a linha de comando acima gera um ''mpg4'' com ''codec'' de windows; depois de gravado, o meu DVD player respondeu que o codec não era reconhecido. A linha abaixo, especificando codec divx (-ffourcc), sim funcionou

mencoder filme.avi -o filme_leg.avi -sub filme.sub -ovc lavc -oac mp3lame -ffourcc DX50
----
Obrigado a Sérgio pela dica

'''NOTA''': ''mencoder'' tem infinitas opções, o exemplo é o mínimo indispensável para conseguir o filme com legendas;
para mais detalhes (qualidade, formato de tela, trechos, mixing, ripping, etc, etc) procurar por ''mencoder'' em google
ou ir para os seguintes links:

[http://gentoo-wiki.com/TIP_MEncoder_Tips_and_Tricks Mencoder Tips & Tricks] 
[http://www.misterhowto.com/index.php?category=Computers&subcategory=Video Video HowTo]

Sincronizar estação

2011-09-19T17:56:14Z

Tekkito: Criou página com 'No quesito sincronização de diretórios entre maquinas (trabalho e casa) acho que o ''UNISON'' é o canal ---- A diferença do ''rsync'', bom para rodar via ''cron'' e manter '...'

No quesito sincronização de diretórios entre maquinas (trabalho e casa) acho
que o ''UNISON'' é o canal
----
A diferença do ''rsync'', bom para rodar via ''cron'' e manter ''backup'', 
''unison'' permite sincronização bi-direcional, 
com possibilidade de resolução de conflitos na mesma jogada.

Exemplo de uso para sincronização de diretório '''''dados/''''' entre maquina local (casa ou laptop) onde é lançado o ''unison'' 
e a maquina remota (chamemos ela de '''''pcwork''''' no dominio '''''uni.br''''')

unison dados/ ssh://pcwork.uni.br/dados/

Como se trata de sincronização nos dois sentidos, a ordem dos diretórios (local e remoto) é irrelevante. 
Isto elimina alguns problemas que o ''rsync'' tem. 
Neste último escolher o sentido errado pode acarretar a perda de todos os dados novos!

Também pode ser feito via túnel (''port forwardin''g) assim: 
Primeiro estabelecer o túnel (por exemplo na porta '''1111''') , 
através da maquina '''tunel.uni.br''' (nome imaginário de maquina aberta para fora do domínio)
ssh -C -N -L 1111:pcwork.uni.br:22 -f tunel.uni.br

Significado das opções: '''C''', compressão; '''N''', não envio de comandos, apenas o túnel, 
'''L''': especifica o conexão entre a porta 1111 (local) e a 22 (remota);''' f''', coloca o ''ssh'' em ''background''

Depois

unison dados/ ssh://localhost:1111/dados/

e ai responder as perguntas

----

=== ATENÇÂO ===
''unison'' deve estar instalado nas duas máquinas e deve ser a '''mesma versão''' em ambas

--[[Usuário:Sebas|Sebas]] 11:41, 11 Junho 2008 (BRT)

C cedilha no Ubuntu

2011-09-19T17:55:56Z

Tekkito: Criou página com 'A solução para que '+c apareça como 'ç', e não como um c acentuado é: 1. Edite o arquivo /usr/share/X11/locale/en_US.UTF-8/Compose e substitua todas as ocorrências (minú...'

A solução para que '+c apareça como 'ç', e não como um c acentuado é:

1. Edite o arquivo /usr/share/X11/locale/en_US.UTF-8/Compose
e substitua todas as ocorrências (minúsculas e maiúsculas) do c acentuado
pelo ç.

2. Edite o arquivo /usr/lib/gtk-2.0/2.10.0/immodule-files.d/libgtk2.0-0.immodules

e acrescente ":en" no final da linha que se refere ao cedilha.

Em outras distribuições, é possível que o arquivo seja /etc/gtk-2.0/gtk.immodules
ou /etc/gtk-2.0/i386-redhat-linux-gnu/gtk.immodules .

3. Reinicialize o servidor X (feche sua sessão, e no menu de login, escolha
esta opção).

Latex, emacs, UTF-8 e acentos

2011-09-19T17:55:35Z

Tekkito: Criou página com 'As vezes os acentos colocados no '''''emacs''''' usando ''dead-keys'' não são interpretados corretamente ao compilar com '''''latex''''' . Entendo que por uma incompatibilidade...'

As vezes os acentos colocados no '''''emacs''''' usando ''dead-keys'' não são interpretados
corretamente ao compilar com '''''latex''''' . Entendo que por uma incompatibilidade entre
as '''codificações do ASCII estendido''' nesta fase de transição do '''latin1 (t1 ou ISO-8859-1)''' para o '''UTF-8'''.

Isto acontece mesmo tendo o

\usepackage{t1enc}

no cabeçalho do documento (ou outros pacotes equivalentes, como o "ae").

A maneira de consertar isso é ou usando outros pacotes para latin1 (t1 ou ISO-8859-1) ou o pacote com opção de UTF-8, assim:

\usepackage[latin1]{inputenc}

ou

\usepackage[T1]{fontenc}

para '''latin1 (t1 ou ISO-8859-1)''' ou

\usepackage[utf8]{inputenc}

para codificação '''utf8'''

e o problema deve se resolver.

---

As vezes, porém, é preciso converter o arquivo
entre diferentes encodings, o que pode ser feito com o comando ''iconv''. Por
exemplo:

iconv -t ISO-8859-1 -f UTF-8 inputfile > outputfile

Há outras maneiras (retirado da lista TeX-br, colaboração de G. Weber):

De ISO para UTF:

recode -d l1..u8 arquivo.tex

De UTF para ISO:

recode -d u8..l1 arquivo.tex

Caso queira converter o arquivo para a acentuação do latex:

recode -d l1..tex arquivo.tex

Lembrete: o recode grava em cima do próprio arquivo. Para sua segurança, faça um back-up do arquivo antes de executar o comando, ou:

recode -d l1..tex < arquivo.tex

isso faz com que o resultado do recode seja mostrado direto na tela.

Pdf

2011-09-19T17:55:10Z

Tekkito: Criou página com '== Tratamento de arquivos pdf == * Para colar varios pdf num só: gs -dNOPAUSE -sDEVICE=pdfwrite -sOUTPUTFILE=merge.pdf -dBATCH 1.pdf 2.pdf 3.pdf 1.pdf, 2.pdf e 3.pdf ficar...'

== Tratamento de arquivos pdf ==

* Para colar varios pdf num só:

gs -dNOPAUSE -sDEVICE=pdfwrite -sOUTPUTFILE=merge.pdf -dBATCH 1.pdf 2.pdf 3.pdf

1.pdf, 2.pdf e 3.pdf ficarão nessa ordem no arquivo merge.pdf. 
Podemos usar o * também mas ai irão em ordem alfabetica, no caso do exemplo é igual

XFig

2011-09-19T17:54:57Z

Tekkito: Criou página com 'Como adicionar uma biblioteca (com molas por exemplo): a) Crie um diretório no local correto: mkdir /usr/share/xfig/Libraries/Springs b) Coloque seu arquivo .fig lá (no ...'

Como adicionar uma biblioteca (com molas por exemplo):

a) Crie um diretório no local correto:

mkdir /usr/share/xfig/Libraries/Springs

b) Coloque seu arquivo .fig lá (no exemplo, um arquivo [http://www.if.ufrgs.br/~arenzon/spring.fig] com
uma mola desenhada):

cp spring.fig /usr/share/xfig/Libraries/Springs

Agora deve aparecer um novo item no menu de bibliotecas do XFig,
chamado "Springs".

Thunderbird

2011-09-19T17:54:16Z

Tekkito: Criou página com 'A receita de como instalar dicionários da página de Thunderbird-Firefox (por ventura a mistura de infos FAQ dos dois programas so faz aumentar a confusão e o tempo para resolv...'

A receita de como instalar dicionários da página de Thunderbird-Firefox (por ventura a mistura de infos FAQ dos dois programas
so faz aumentar a confusão e o tempo para resolver os problemas, seria bom que separesem as infos de cada um).
Em fim, mas a receita que esta no [http://www.mozilla.org/support/thunderbird/ Help] ou [http://www.mozilla.org/support/thunderbird/faq FAQs]
deve ser para windows pois em linux não consegui faze-lo para todos os dicionários.

Mas consegui faze-lo "a mão" assim:

0) Fechar o ''thunderbird'' se estiver aberto

1) Baixar as extensões dos dicionários desejados: pode procurar na página de mozilla o com google pelos archivos:

spell-es-ES.xpi spell-fr-FR.xpi spell-pt-BR.xpi

que são archivos comprimidos contendo os dicionarios (*.aff e *.dic) de Espanhol, Francês e Portugues BR respectivamente,
e salva-los num diretorio vazio para melhor identifica-los. Ex: dic/

2) Descompactar os dicionários desejados, por exemplo:

unzip spell-pt-BR.xpi

o que deve resultar em
Archive: spell-pt-BR.xpi
inflating: install.js
inflating: pt-BR.aff
inflating: pt-BR.dic
inflating: README_pt_BR.txt

3) Finalmente como root ou com ''sudo'' mover os archivos .aff e .dic para o /usr/share/myspell/dicts/ (*)

sudo mv pt* /usr/share/myspell/dicts/

4) Lanzar o ''Thunder'' e verificar com o botão "spell" que o(s) novos dicionários estão disponíveis.
Assim consegui instalar o BR, ES e FR.

(*) Ignoro quan geral é a existencia do diretório /usr/share/myspell/dicts/
Se ele no existe no seu sistema procure por outros .aff ou .dic assim:

locate .aff

e prove colocar ali os .aff e .dic novos

Se o resultado de "locate ..." for nulo então instale com ''synaptic'' ou ''apt-get'' algum dicionário extra
e retorne ao item 3)

Thunderbird & firefox

2011-09-19T17:53:56Z

Tekkito: Criou página com '== Para o Firefox ser o navegador web do Thunderbird == * fechar o ''thunderbird'' (senão as modificacões serão perdidas): * editar o archivo: ''prefs.js'' Geralmente esse ...'

== Para o Firefox ser o navegador web do Thunderbird ==

* fechar o ''thunderbird'' (senão as modificacões serão perdidas):

* editar o archivo: ''prefs.js''

Geralmente esse archivo se encontra em .''thunderbird/zxfrs5rd.default/''

(ou ''.firefox-thunderbird/...''), onde "zxfrs5rd" é uma sequência única para cada usuário

* acrescentar ao ''prefs.js'' as seguintes duas linhas:
<pre>
user_pref("network.protocol-handler.app.http", "/usr/bin/firefox");
user_pref("network.protocol-handler.app.https", "/usr/bin/firefox");
</pre>
----

'''NOTA1:''' verificar se o firefox está em /usr/bin/firefox
conferir com "which firefox" e substituir pela resposta se for diferente de /usr/bin/firefox

'''NOTA2:''' ha um ''prefs.js'' que esta em ''.thunderbird/...'' mas não é o que tem de ser editado

----

Outras informações sobre configuração (link externo): http://www.geocities.com/pratiksolanki/

32 bits rodando em 64

2011-09-19T17:53:30Z

Tekkito: Criou página com 'Muitos programas úteis como openoffice (BrOffice) ou flash plugins para navegadores não tem versões estáveis para 64 bits. Os passos abaixo mostram como fazer isto em um AM...'

Muitos programas úteis como openoffice (BrOffice) ou flash plugins para navegadores não tem versões estáveis para 64 bits. Os passos abaixo mostram como fazer isto em um AMD64 com instalação [http://www.debian.org Debian].

Resumidamente, o método consiste em fazer uma instalação básica de 32 bits em um diretório (/var) da instalação de 64 bits e, quando os programas de 32 bits forem solicitadas, passa-se a usar (chroot) automaticamente
esse diretório como raiz para a execução do mesmo.

Os passos abaixo são uma tradução aproximada de [http://www.debian-administration.org/articles/356 artigo]
encontrado em uma página do administrador Debian.

'''1. Instalando o chroot'''

1.1) Criar (com o usuário root) o diretório vazio na raiz:

$ mkdir /var/etch-386-chroot

1.2) Instalar o sistema base Debian etch versão i386 nesse diretório (abaixo é uma única linha de comando):

$ debootstrap --arch i386 etch /var/etch-386-chroot http://debian.if.ufrgs.br/debian/

Escolha o espelho mais próximo... Isto deve demorar de 10 a 15 min.
Se tudo deu certo até aqui deve aparecer ao final a mensagem:

I: Base system installed successfully.

Se o comando debootstrap não existir, instale com:

$ apt-get install debootstrap

1.3) Coloque o caminho(path) de suas bibliotecas no /etc/ld.so.conf. Para isto edite esse arquivo de forma
que ele fique assim:

/usr/X11R6/lib
# chroot i386 system libs
/var/etch-386-chroot/lib
/var/etch-386-chroot/usr/lib
/var/etch-386-chroot/usr/X11R6/lib
/var/etch-386-chroot/usr/local/lib

1.4) Coloque o link para o seu linker de 32 bits no caminho
de /lib, mas preserve o link original, caso algo dê errado...

$ cd /lib
$ mv ld-linux.so.2 ld-linux.so.2.orig
$ ln -s /var/etch-386-chroot/lib/ld-linux.so.2 ld-linux.so.2

1.5) Rode ldconfig para atualizar o cache do seu linker.

$ ldconfig

1.6) Para rodar um programa dentro do chroot necessita-se que algumas partes do sistema 64 bits estejam dentro do chroot. Isto pode ser feito com um "bind mount". No exemplo abaixo, ligaremos o /tmp, o /home e o proc. Edite o arquivo /etc/fstab adicionando as linhas abaixo:

#ia32 chroot
/home /var/etch-386-chroot/home none bind 0 0
/tmp /var/etch-386-chroot/tmp none bind 0 0
proc /var/etch-386-chroot/proc proc defaults 0 0

1.7) Monte todos eles...
$ mount /var/etch-386-chroot/home
$ mount /var/etch-386-chroot/tmp
$ mount /var/etch-386-chroot/proc

1.8) Passe os usuários e grupos para o chroot:
$ cp /etc/passwd /var/etch-386-chroot/etc/passwd
$ cp /etc/shadow /var/etch-386-chroot/etc/shadow
$ cp /etc/group /var/etch-386-chroot/etc/group

1.9) Nesse ponto pode-se já rodar programas dentro do chroot. Troque o shell do diretório raiz para o do seu chroot:
$ chroot /var/etch-386-chroot

1.10) Experimente rodar algo simples como:
$ cat /etc/apt/sources.list

'''2. Instalar programas 32-bits no chroot'''

2.1) Ainda dentro do shell do chroot instale uns pacotes que precisaremos (o comando abaixo é de apenas uma linha)...

$ apt-get install libogg0 libogg-dev libvorbis0a libvorbis-dev vorbis-tools
libpng12-0 libpng12-dev liblogfile-rotate-perl libconfhelper-perl
debhelper fakeroot gcc g++ libgtk1.2 libgtk1.2-dev svgalibg1
svgalibg1-dev libsdl1.2-dev libaa1 libaa1-dev libxv1 libxv-dev

2.2) Além de outros que sempre precisamos:

$ apt-get install vim bzip2 jed xjed

2.3) Perl vai reclamar de não ter locale configurado... Adicione
export LC_ALL=C
ao arquivo /root/.bashrc e carregue-o:
source ~/.bashrc

2.4) Agora que chegamos até aqui, instale da maneira usual aqueles programas que só rodam em 32 bits, e.g., para instalar o firefox:

$ apt-get install firefox

Saia do chroot e não esqueça de desinstalar a versão em 64 bits.

$ exit
$ apt-get remove firefox

Teste para ver se o firefox 64 foi desinstalada:

$ firefox
bash: firefox: command not found

Verifique se a versão 32 bits funciona:

$ xhost +
$ chroot /var/etch-386-chroot
$ firefox

Instale os plugins necessários como de hábito.
O item 4, ao final desta, mostra detalhes para a
instalação do flashplayer.

'''3. Tornando transparente ao usuário'''

Até agora tudo funcionou porque estamos trabalhando como root e fazendo
manualmente a passagem para 32 bits via chroot. Para o usuário normal isto
tudo deve ser feito quando ele chama o programa do ambiente de 64 bits.

3.1) Saia do chroot e instale o pacote dchroot:

$ exit
$ apt-get install dchroot

3.2) Crie o arquivo /etc/dchroot.conf com o conteúdo abaixo:

#etchi386 chroot
etchi386 /var/etch-386-chroot

3.3) Agora um usuário normal deve ser capaz de rodar comandos no chroot:

$ dchroot -d firefox

3.4) Adicione um script para facilitar as coisas de forma que se possa chamar
os programas 32 bits transparentemente da linha de comando. No diretório
/usr/local/bin crie o script abaixo, nomeie-o do_chroot e troque as permissões
para 755:

#!/bin/sh
ARGS=""
for i in "$@" ; do
ARGS="$ARGS '$i'"
done
exec dchroot -d -q "`basename $0`" "$ARGS"

3.5) Último passo crie os links simbólicos dos programas para o script:

$ln -s do_chroot firefox

'''4.) Instalando o flashplayer'''

4.1) Crie e entre em um diretório temporário, baixe o pacote flashplayer da Macromedia e
abra-o.

$ mkdir tmp; cd tmp
$ wget http://fpdownload.macromedia.com/get/flashplayer/current/install_flash_player_7_linux.tar.gz
$ tar xvzf install_flash_player_7_linux.tar.gz

4.2) Como root, entre no diretório install_flash_player_7_linux (criado no desempacotamento) e copie
os arquivos flashplayer.xpt e libflashplayer.so para o diretório de plugins de instalação do firefox:

$ cp flashplayer.xpt /var/etch-386-chroot/usr/lib/firefox/plugins
$ cp libflashplayer.so /var/etch-386-chroot/usr/lib/firefox/plugins

4.3) Acerte as permissões:

$ cd /var/etch-386-chroot/usr/lib/firefox/plugins
$ chmod a+rw flashplayer.xpt
$ chmod a+rw libflashplayer.so

4.4) Provavelmente faltará a biblioteca libXmu6 na sua instalação i386; faça o chroot para a instalação 386
e instale a biblioteca:

$ chroot /var/etch-386-chroot/
$ apt-get update
$ apt-get install libXmu6

4.5) Saia do chroot, abra o firefox como usuário normal e digite na linha de endereços internet do navegador:

about:plugins

O Shockwave flash deve aparecer entre os plugins instalados.

Upgrade para Edgy

2011-09-19T17:52:17Z

Tekkito: Criou página com '=== ifort (intel f90) ERROR === > export: 27: Illegal option -n The problem is that edgy uses dash rather than bash. You can change the first line of the ifort and ifc script...'

=== ifort (intel f90) ERROR ===

> export: 27: Illegal option -n

The problem is that edgy uses dash rather than bash. You can change the first line of the ifort and ifc scripts from

#!/bin/sh
to

#!/bin/bash

and all should work.

Kubuntu

2011-09-19T17:51:57Z

Tekkito: Criou página com 'kubuntu ou ubuntu são instalados sem configurar a senha de "root" Para faze-lo a mão: $ sudo passwd root $ Password: senha_de_usuario $ Enter new UNIX password: senha_de_r...'

kubuntu ou ubuntu são instalados sem configurar a senha de "root"

Para faze-lo a mão:

$ sudo passwd root
$ Password: senha_de_usuario
$ Enter new UNIX password: senha_de_root
$ Retype new UNIX password: senha_de_root

Sources.list

2011-09-19T17:50:32Z

Tekkito: Criou página com '''sources.list'' é um archivo que contem os endereços ''http'' dos repositorios de pacotes de distribuições '''Debian''', o que se aplica também para '''Ubuntu''' e '''Kubun...'

''sources.list'' é um archivo que contem os endereços ''http'' dos repositorios de pacotes de distribuições '''Debian''',
o que se aplica também para '''Ubuntu''' e '''Kubuntu'''.

Este archivo reside em ''/etc/apt'' pois fornece essa informação ao ''apt-get'' que é a interfase de instalação de programas. 
Também é usado pelo ''synaptic'' que é uma interfase gráfica amigável entre o ''apt-get'' (que funciona por linha de comando)
e o usuário.

Uma pergunta que surge sempre é como escolher esses repositorios. 
Neste sitio é posível configurar um conjunto de repositorios ao gosto de cada um:

[http://www.ubuntu-nl.org/source-o-matic/ Source-O-Matic]

O Ubuntu Source-O-Matic faz a geração automática do arquivo ''/etc/apt/source.list'' tomando por base a sua localização geográfica, a versão do ubuntu que você usa e a arquitetura de seu computador.

Além dos repositórios padrão, o aplicativo oferece também a possibilidae de incluir repositórios contendo um grande número de alternativas interessantes (Wine, Opera, etc).

Reinstalando todos os pacotes de instalação antiga

2011-09-19T17:50:03Z

Tekkito: Criou página com '== Backup list of installed software == Debian Linux If you are using Debian Linux use dpkg command to list installed software: $ dpkg --get-selections Store list of installe...'

== Backup list of installed software ==

Debian Linux
If you are using Debian Linux use dpkg command to list installed software:

$ dpkg --get-selections

Store list of installed software to a file called /backup/installed-software.log

$ dpkg --get-selections > /backup/installed-software.log

Now you have a list of installed software. After installing base system you can immediately install all software.

Debian Linux
Debian Linux makes your life easy. All you have to do is type following two commands:

# dpkg --set-selections < /backup/installed-software.log
Now your list is imported use dselect or other tools to install the package.

# dselect

Select ‘i‘ for install the software.

Instalação de linux basica pela rede

2011-09-19T17:49:43Z

Tekkito: Criou página com 'Abaixo está exemplificado um processo de instalação do sistema operacional (S.O) pela conexão de rede. Começamos usando o [http://pcleon/pub/utils/tomsrtbt-2.0.103.tar.gz ...'

Abaixo está exemplificado um processo de instalação do sistema operacional (S.O) pela conexão de rede.

Começamos usando o [http://pcleon/pub/utils/tomsrtbt-2.0.103.tar.gz Tomsrtbt], que é um S.O muito básico, o qual é instalado em um disquete, no qual mais tarde é feito o boot.
O donwnload do Tomsrtbt pode ser feito clicando-se no link acima. Para fazer sua instalação, basta descompactar o arquivo.
&tar -xvvzf tomsrtbt-2.0.103.tar.gz -> descompacta o arquivo
Feito isto, precizamos fazer a instalação do sistema em um disquete, para isso entre no diretório descompactado do toms
$cd tomsrtbt-2.0.103
Uma vez no diretório, basta executar o seguinte comando para que a instalação tenha início:
$./install.s
Observe que não é preciso ser feita a montagem do disquete, simplismente basta coloca-lo o drive que o programa de instalação se encarregará de todo resto.
Feito isto está concluído nosso disquete de boot.

Lembre-se, que poderá utilizar qualquer meio para fazer a instlação via rede, por algum sistama que você já tenha ou por um live CD, ou outro disquete como este do toms. O Toms é só um exemplo.

Agora sim passamos para a parte da instalão via rede.

De o boot da máquina no seu sistema (live CD, disquete), que neste casa estamos utilizando um disquete com o Toms.
Se for usar o toms, ele vai pedir para que faça um login
usuário: root
senha: xxxx
Feito o boot no toms, passamos para a configuração da máquina na qual será feita a instalação.

1º) Configurção do pc na rede
$ifconfig eth0 netmask yyy.yyy.yyy.yyy xxx.xxx.xxx.xxx
Na parte dos "y" deve ser colocada a máscara da rede, ex.: 255.255.255.0
Na parte dos "x" deve ser colocado o endereço de ip do pc, ex.: 101.90.68.01

2º) Preparando o HD, criando uma tabela de partições.
$fdisk /dev/hda
Vai abrir um software. Verifique se já há alguma partição em seu HD com o comando
$p
Se já hover, é recomendado que delete as partições existentes. "Mas somente se este for o único sistema a ser instalado, caso haja outro sistema instalado e que desejas manter, mantenha a(s)partição(ões) daquele sistema."
Para deletar uma partição basta digitar o comando
$d
e então o número da partição que desejas remover, ex:
$1
3º) Criar as novas partições

3º-a)O Swap
$n
$p
$1
$opção default
$+256M ->o tamanho do swap fica a seu critério
neste caso estamos criando um swap de 256MB, o recomendado é para o swap é o quanto você tem de memória RAM mais 50% do tamanho da RAM, ex.: 128MB de RAM, swap = 128 + 128 * 50% => swap = 192MB

3º-b)A partição onde vai ficar o sistema.
$n
$p
$2
$opção default
$opção default -> default aqui = utiliza o restante do espaço do disco
Não estou dando muitas explicações nesta parte, pois o software (que está no disco do toms) de particionamento as fornece.
3º-c)Concluindo o particionamento e escrevendo as informações no disco.
$w
Pronto, o software irá preparar o disco e será fechado automaticamente.

4º)Concluindo a preparação do HD.

4º-a)"Escrevendo" o swap no disco.
$mkswap /dev/hda1
4º-b)"Escrevendo" a partição de dados no disco.
$mke2fs -j /dev/hda2
Obs.:/dev/hdaX, o X representa o número da partição, que é aquele q você estabeleceu para a partição quando preparou seu disco com o "fdisk"

5º)Criando o diretório /root e montando a partição de dados nela.
$mkdir /root
$mount /dev/hda2 /root
6º)Montando a unidade da rede onde encontra-se a imagem de instalação, em nosso caso está em:
143.54.2.244:/linux/imagem
$mount 143.54.2.244:/linux/imagem /mnt
Vão aparecer algumas mensagens de erro, mas não se preocupe, pois é normal.

7º)Fazendo a cópia da imagem da rede para a sua máquina
$pax -rw /mnt /root
Agora a única coisa que você tem para fazer é esperar a cópia terminar, você pode acompanhar o processo de cópia verificando o quanto do espaço de seu disco está sendo usado. Para isso basta abrir uma nova "janela", precione Ctrl+Alt+F2 para faze-lo. Novamente será solicitado fazer o login,
usuário: root
senha: xxxx

Para ver o quanto do espaço de seu disco está sendo utilizado, basta digitar o seguinte comando
$df
e pronto, assim você poderá acompanhar como está o processo de cópia da imagem.

Para retornar a janela inicial onde você está fazendo a intalação basita precionar Ctrl+Alt+F1.
Assim que a cópia tiver sido finalizada lhe será devolvida a prompt de comando
$#
8º)Com a cópia finalizada, vamos para a parte da configuração do sistema que será instalado.
$mv /root/mnt/* /root
$chroot /root
9º)Configuração de rede:
Obs.:Para salvar as alterçoes que serão feitas com o programa "jed", precione as seguinte seqüencia de teclas: Ctrl+X -> Ctrl+S
E para sair do jed digite a seguinte seqüencia de teclas: Ctrl+X -> Ctrl+C
$jed /etc/network/interfaces
Aqui você deverá colocar o endereço de IP de sua máquina (seu pc)
$jed /etc/hostname
Aqui você deve colocar no nome de sua máquina (pc)
10ºConfiguração de teclado, mouse e monitor (placa de vídeo)
$jed /etc/X11/XF86Config-4
Aqui você poderá configurar o sistema de acordo com o teclado que irá utilizar, sua placa de vídeo e o tipo de mouse que iras utilizar, se for PS2 ou Serial
Para um mouse serial vai ficar assim:
Identifier "Configured Mouse"
Driver "mouse"
Option "CorePointer"
Option "Device" "/dev/ttyS0"
Option "Protocol" "microsoft"
Option "Emulate3Buttons" "true"
Option "ZAxisMapping" "4 5"
Onde o "0" no /dev/ttyS0 é o número da porta onde o mouse está conectado. 0 = porta 1.
E para um mouse PS2 vai ficar assim:
Section "InputDevice"
Identifier "Configured Mouse"
Driver "mouse"
Option "CorePointer"
Option "Device" "/dev/ttyS0"
Option "Protocol" "microsoft"
Option "Emulate3Buttons" "true"
Option "ZAxisMapping" "4 5"
É imprescindível que você salve as alterções em todas configurações de rede (IP e nome do PC), de teclado, mouse, placa de vídeo.

Feito toda esta parte de configuração só resta é executar o comando de instalação.
$/sbin/grub-install /dev/hda
Se ocorrer algum erro o sistema irá lhe informar após o término da execução deste comando.

Para concluir, basta sair do chroot, para isso:
$exit
e reiniciar sua máquina:
$reboot

E pronto, agora você não precisa mais do disquete do toms, pode retirá-lo do drive.
Seu PC já será iniciado com o novo sistema.

Para uma reconfiguração de teclado após concluída a intalação poderá ser feita através do seguinte comando, mas para isso, vc precisa ir para uma interface de comando, pode ser precionando Ctrl+Alt+F1 então digite o comando:
$dpkg-reconfigure console-data
Vai aparecer um ambiente gráfico no qual você poderá escolher melhor o seu tipo de teclado. Finalize precionando as teclas Ctrl+Alt+Backspace.
Para retornar para o ambiente gráfico basta precionar: Ctrl+Alt+F7.

Pronto, agora sim, definitivamente está conluído o processo de intalação e configuração.

Linux & vaio sz

2011-09-19T17:49:23Z

Tekkito: Criou página com 'Atualmente (Ubuntu hardy - kernel 2.6.24-21) quase tudo funciona de partida. Por isso retiro a parte que se refere ao som pois isso não é mais problema. As únicas coisas que a...'

Atualmente (Ubuntu hardy - kernel 2.6.24-21)
quase tudo funciona de partida. Por isso retiro a parte que se refere ao som pois
isso não é mais problema. As únicas coisas que ainda não funcionam de partida são
o "leitor de digitais" ou (pouca o nenhuma utilidade) e a "webcam"

Então retiro a parte de som e coloco a de webcam. Se alguém quiser ainda ver a versão antiga desta entrada
ir no history.

== Webcam ==

A dica abaixo encontrei no foro de Ubuntu [http://ubuntuforums.org] e pode ser útil não só para
os Vaios como também para outros notebooks cuja webcam use drivers da Ricoh
Vou detalhar abaixo pois depois de resolvido e com a webcam funcionando não encontrei mais o
item no foro.

* Acrescentar a linha no /etc/apt/sources.list

deb http://download.tuxfamily.org/arakhne/ubuntu hardy-arakhne universe</nowiki>

substituindo hardy pela distribuição correspondente

* Baixar a chave publica e acrescentá-la ao apt key list (para nao ter warning durante update de pacotes)

wget -q http://download.tuxfamily.org/arakhne/public.key -O- | sudo apt-key add

caso alguma coisa mude com o tempo ir a [http://www.arakhne.org/download.html]

* Atualizar lista de pacotes

sudo apt-get update

* procurar o driver da camera

sudo apt-cache search ricoh

* eu instalei o ricoh-webcam-r5u870:

sudo apt-get install ricoh-webcam-r5u870

* pronto a camera deve estar funcionando

----
Os últimos dois links estão desatualizados

* http://www.linlap.com/wiki/Sony+VAIO+VGN-SZ+Series

* http://ariel.vardi.free.fr/ariel//vaiosz.html
* http://www.clip.dia.fi.upm.es/~cochoa/wiki/index.php/Installing_Ubuntu_On_a_Vaio

Impressoras multifuncionais

2011-09-19T17:49:08Z

Tekkito: Criou página com 'Como instalar/configurar uma Epson CX7700 no Debian (deve ser parecido para outros modelos da Epson e outras distribuições, embora para Ubuntu, a partir da 7.10, não seja mais...'

Como instalar/configurar uma Epson CX7700 no Debian
(deve ser parecido para outros modelos da Epson e outras distribuições,
embora para Ubuntu, a partir da 7.10, não seja mais necessário
instalar a mão, os pacotes necessários já estão disponíveis)

1) Instale os pacotes referentes ao CUPS (impressora) e
SANE (scanner), usando o synaptic, por exemplo (recomendando).

2) Baixe os drivers referentes ao seu modelo de

http://www.avasys.jp/english/linux_e/dl_spc.html

3) Descompacte o arquivo e vá para o diretório recém criado.
Procure onde estao as bibliotecas do SANE (libsane) instaladas
no seu sistema. Se foram instaladas via apt-get/synaptic, devem
estar em /usr/lib. O pacote acima, por default, instala tudo em
/usr/local, então precisamos ensiná-lo:

./configure --prefix=/usr

4) Como root:

make
make install

5) O comando

sane-find-scanner

deve encontrar agora o scanner, algo como:

found USB scanner (vendor=0x04b8 [EPSON], product=0x081f [USB2.0 MFP(Hi-Speed)]) at libusb:005:006

Ou também através do comando

lsusb

que retorna algo como:

Bus 005 Device 006: ID 04b8:081f Seiko Epson Corp.

6) Tente rodar o executável. Deve funcionar somente como root, por enquanto:

iscan

7) Não tenho certeza se os passos abaixo são realmente necessários, ou se o
item 6 acima já deu resultado, mas há dois arquivos em /etc/sane.d que parecem
configurar o scanner:

/etc/sane.d/epkowa.conf
/etc/sane.d/epson.conf

e descomentamos as seguintes linhas:

usb
usb 0x04b8 0x081f

onde os argumentos da segunda linha foram obtidos na saída do comando sane-find-scanner (ver acima).

8) Finalmente, devemos dar permissão aos usuários para utilizarem o scanner. Há uma dica
em

http://www.sane-project.org/man/sane-usb.5.html

The permissions for the device files used by libusb must be adjusted
for user access. Otherwise only root can use SANE devices. For Linux,
the devices are located in /proc/bus/usb/ or in /dev/bus/usb, if you
use udev. There are directories named e.g. "001" (the bus name) con-
taining files "001", "002" etc. (the device files). The right device
files can be found out by running scanimage -L as root. Setting permis-
sions with "chmod" is not permanent, however. They will be reset after
reboot or replugging the scanner.
Usually udev or for older distributions the hotplug utilities are used,
which support dynamic setting of access permissions. SANE comes with
udev and hotplug scripts in the directory tools/udev and tools/hotplug.
They can be used for setting permissions, see
/usr/local/doc/sane-1.0.18/README.linux, tools/README and the README in
the tools/hotplug directory for more details.

Por exemplo, no meu caso, fiz:

chmod a+w /dev/bus/usb/005/*

Funciona, mas como diz no texto acima, não é uma modificação permanente. Uma
solução permanente é acrescentar o nome do usuário no grupo scanner (pode
ser feita via kuser).

Touchpad (Synaptic/ALPS)

2011-09-19T17:48:52Z

Tekkito: Criou página com 'Como evitar que o cursor pule descontroladamente ao editar texto por um arrastre involuntario da mão sobre o wild touchpad das notebooks? Como qualquer dispositivo que funcio...'

Como evitar que o cursor pule descontroladamente ao editar texto por um arrastre involuntario da mão sobre o
wild touchpad das notebooks?

Como qualquer dispositivo que funcione (mal ou bem), isto é que o driver foi achado e está vivo (o dispositivo,
seja mouse, tela ou touchpad=mouse de laptop) a configuracão é controlada no archivo "xorg.conf"

Para ver o estado dele:

cat /etc/X11/xorg.conf

Caso não esteja em /etc/X11 procure ele com

locate xorg.conf

A dois tipos de entrada nesse archivo: tipo de dispositivo com opcoes (colocadas em qualquer ordem, é irrelevante)

A secao que nos interesa agora é

Section "InputDevice"
Identifier "Synaptics Touchpad"
Driver "synaptics"
Option "SendCoreEvents" "true"
Option "Device" "/dev/psaux"
Option "Protocol" "auto-dev"
Option "HorizScrollDelta" "0"
EndSection

assim normalmante fica instalado o Touchpad Synaptic, mas com funcionalidades mínimas (ou pior, ruins)

e o bloco que carrega os dispositivos no fim.

Section "ServerLayout"
Identifier "Default Layout"
Screen "Default Screen"
InputDevice "Generic Keyboard"
InputDevice "Configured Mouse"
InputDevice "stylus" "SendCoreEvents"
InputDevice "cursor" "SendCoreEvents"
InputDevice "eraser" "SendCoreEvents"
InputDevice "Synaptic Touchpad"
EndSection

----

* Faca copia de seguranca do xorg.conf ativo assim

sudo cp /etc/X11/xorg.conf /etc/X11/xorg.conf.bck

* Agora podemos editar tranquilamente ele (se der zebra substituimos pelo .bck)

Acrescentando as definicoes que seguem (com opcoes novas)

Section "InputDevice"
Identifier "Synaptics Touchpad"
Driver "synaptics"
Option "SendCoreEvents" "true"
Option "Device" "/dev/psaux"
Option "Protocol" "auto-dev"
Option "HorizScrollDelta" "0"
# New options !!!
Option "LeftEdge" "120"
Option "RightEdge" "830"
Option "TopEdge" "120"
Option "BottomEdge" "650"
Option "FingerLow" "14"
Option "FingerHigh" "15"
Option "MaxTapTime" "180"
Option "MaxTapMove" "110"
Option "ClickTime" "10"
Option "EmulateMidButtonTime" "75"
Option "VertScrollDelta" "10"
Option "MinSpeed" "0.45"
Option "MaxSpeed" "0.75"
Option "AccelFactor" "0.020"
Option "EdgeMotionMinSpeed" "200"
Option "EdgeMotionMaxSpeed" "200"
Option "UpDownScrolling" "1"
Option "CircularScrolling" "0"
Option "CircScrollDelta" "0.1"
Option "CircScrollTrigger" "2"
Option "SHMConfig" "true"
EndSection

* Salvar

* Relançar o X

teremos por exemplo que o touchpad detecta quando o teclado esta em uso e desativa o arrastro de mao

* Nota: eu não sei se todas essa opções são realmente necessarias. Para ter certeza haveria que desativar/ativar
trocar valor uma por uma e analisar o resultado (o que é muito chato pois precisa relançar o X a cada mudença)

Mas a que achei chave para controlar o touchpad é

Option "ClickTime" "10"

Aperentemente esse valor (10) é o tempo que o teclado deve estar sem atividade para o touchpad
voltar a controlar a posição do cursor

* Nota2: Nessa busca encontrei um site que mostra um pouco os detalhes do que significa cada opção do synaptics, o link é [http://linux.die.net/man/5/synaptics Lista de comandos]
----

Mais informações:
[http://ubuntu.wordpress.com/2005/11/15/fixing-my-alps-touchpad-with-the-synaptics-driver/ ALPS Touchpad]

Dicas para apresentações

2011-09-19T17:46:08Z

Tekkito: Criou página com 'Algumas dicas para os alunos de IC (mas que podem ser uteis até para os de doutorado) sobre apresentacões em geral: ---- [http://citnews.unl.edu/presentations/graphs.html Para...'

Algumas dicas para os alunos de IC (mas que podem ser uteis
até para os de doutorado) sobre apresentacões em geral:

----
[http://citnews.unl.edu/presentations/graphs.html Para ter em conta ao fazer gráficos]

Applets

2011-09-19T17:45:53Z

Tekkito: Criou página com '== Applets de Física == * [http://www.surendranath.org/Applets/Heat/MolMotion/MolMotionApplet.html Gas Ideal] * [http://comp.uark.edu/~jgeabana/mol_dyn/KinThI.html Kinetic The...'

== Applets de Física ==

* [http://www.surendranath.org/Applets/Heat/MolMotion/MolMotionApplet.html Gas Ideal]

* [http://comp.uark.edu/~jgeabana/mol_dyn/KinThI.html Kinetic Theory]

* [http://www.phy.ntnu.edu.tw/ntnujava/ ntnujava]

Lei de Bragg

2011-09-19T17:45:10Z

Tekkito: Criou página com '     Depois da descoberta de que os raios catódicos eram elétrons, feita por J.J. Thomson faltava ainda descobrir a natureza dos raios X. Helmhotz prev...'

Depois da descoberta de que os raios catódicos eram elétrons, feita por J.J. Thomson faltava ainda descobrir a natureza dos raios X. Helmhotz previa que os raios X eram ondas eletromagnéticas.

     Para confirmar a hipótese de que os raios X eram [[ondas]] eletromagnéticas foi necessário verificar se feixes de raio X apresentavam interferência. Se essa hipótese se confirmasse os raios X seriam considerados como ondas, pois interferência é um fenômeno característico de ondas. Porém a interferência não era observada.

     A não visualização de interferência estava relacionada com o comprimento de onda, conforme Max von Laue. Ele propôs que o comprimento de onda dos raios X era tão pequeno que tornava impossível construir fendas que estivessem juntas o suficiente.

     [[von Laue]] utilizou então, cristais em suas experiência, uma vez que estes possuem planos regulares bem próximos, caso a distância entre estes planos fosse comparável com o comprimento de onda dos raios X seria possível difratar feixes através destes planos.

     von Laue, W. Friedich e P. Knipping realizaram um experimento onde incidiam raios X em cristais e conseguiram observar manchas em uma placa fotográfica que indicavam difração. Porém não conseguiam explicar por que estas manchas apareciam em direções específicas.

     [[William Henry Bragg]] e [[William Lawrence Bragg]] trabalharam em cima de uma teoria para explicar tais manchas, esta teoria lhes concedeu o prêmio Nobel de Física em 1915.

     Quando um feixe de raios X incide sobre um cristal, os planos de átomos deste cristal funcionarão como espelhos e farão com que o feixe seja refletido. Esta reflexão ocorre porque o espaçamento entre os planos atômicos de um cristal é da ordem do comprimento de onda do raio X.

     Quando dois raios incidem em planos diferentes do cristal um deles percorrerá um caminho um pouco maior. Esta diferença de caminho caracterizará a diferença de fase entre os raios difratados, que está diretamente relacionada com o tipo de interferência que ocorrerá entre os raios.

     Como se sabe sempre que a diferença de fase entre duas ondas for um número inteiro de comprimentos de onda (nλ, para n= 0, 1, 2, 3, etc), estas irão interferir construtivamente, isto é, suas amplitudes se somam. E quando a diferença de fase entre duas ondas for um número inteiro mais meio comprimento de onda [(n+1/2) λ, para n= 1, 2, 3, etc], estas irão interferir destrutivamente, isto é, suas amplitudes se subtraem.
Observando a figura podemos obter uma equação que descreve o tipo de interferência gerado a partir da diferença de caminho percorrido pelos raios.Esta é a lei Bragg:

::<math>2 d \sin\theta = n\lambda\,</math>

Natureza dos Raios-X

2011-09-19T17:44:53Z

Tekkito: Criou página com '     A descoberta dos raios-X gerou uma grande questão a ser esclarecida. O que eram os raios-X? As primeiras análises de Röntgen indicavam que os...'

     A descoberta dos raios-X gerou uma grande questão a ser esclarecida. O que eram os raios-X? As primeiras análises de [[Röntgen]] indicavam que os raios-X poderiam ser algo parecido com o ultravioleta. Após verificar outras propriedades, Röntgen conclui sugerindo que os raios-X poderiam ser vibrações eletromagnéticas longitudinais no éter.

     Em 1912, Max [[von Laue]] propôs que se os cristais fossem compostos por átomos regularmente espaçados e os raios-X fossem efetivamente ondas eletromagnéticas com comprimento de onda da ordem de distâncias interatômicas; então, os cristais funcionariam como centros espalhadores para os raios-X e seria possível a difração de raios-X. A suposição de Laue estava correta e o fenômeno de difração de raios-X foi descoberto corroborando o caráter ondulatório dos raios-X. Devido aos seus trabalhos, Laue ganhou o prêmio Nobel de Física em 1914.

     Em 1918, [[Arthur Holly Compton]] inicia seus estudos sobre o espalhamento de raios-X. Como resultado de suas pesquisas, em 1923, verifica que quando incidia raios-X de um determinado comprimento de onda sobre um alvo de grafite os raios-X espalhados apresentavam dois picos de intensidade; um do mesmo comprimento de onda do feixe incidente e outro, que variava conforme o ângulo do espalhamento, de comprimento de onda maior.

     Esse fato não poderia ser explicado através das características ondulatórias dos raios-X; logo, a interpretação de Compton propunha uma natureza corpuscular dos raios-X. Propôs que o feixe de raios-X era composto por fótons que colidiam com elétrons livres do alvo e eram espalhados devido a colisão. Ao colidirem cediam energia para o elétron resultando em um feixe espalhado de menor freqüência e conseqüentemente de maior comprimento de onda. Essa variação no comprimento de onda do feixe espalhado ficou conhecida como [[Efeito Compton]] e foi posteriormente comprovado por Charles Thomson Rees Wilson com quem dividiu o prêmio Nobel de 1927.

     Através dessas descobertas ficou evidenciada a natureza dualística dos raios-X, ou seja, de possuir comportamento ondulatório e corpuscular. Trata-se de uma onda eletromagnética situada no espectro eletromagnético entre os raios gama e o ultravioleta e tendo o fóton como partícula elementar.

Produção e emissão de Raios-X

2011-09-19T17:44:36Z

Tekkito: Criou página com '     Basicamente, os raios-X são produzidos pela colisão de elétrons, com energia cinética da ordem de keV, contra o anodo em um tubo de raios-X. &n...'

Basicamente, os raios-X são produzidos pela colisão de elétrons, com energia cinética da ordem de keV, contra o anodo em um tubo de raios-X.

     Elétrons são arrancados do cátodo e acelerados contra o anodo, geralmente de tungstênio, devido uma alta diferença de potencial estabelecida entre ambos. Com a desaceleração sofrida no impacto com o anodo, raios-X são emitidos pelo alvo.

[[Image:XraySpectra.gif|thumb| Figura I]]

     A emissão de raios-X pode ser analisada através de seu espectro (Figura I), que consiste em um gráfico contínuo de intensidade de emissão por comprimento de onda com picos característicos.

     A parte contínua refere-se ao tipo de emissão chamado ''Bremsstrahlung'' (do alemão; brems = frenagem, strahlung = radiação) e é devido aos elétrons que são desacelerados pela interação e colisões com os átomos do alvo. À medida que o elétron penetra no alvo, vai perdendo energia e emitindo raios-X em diversos comprimentos de onda, aos quais, as intensidades em que cada comprimento de onda é emitido são explicitadas no espectro. Analisando o espectro verifica-se um comprimento de onda mínimo para a emissão que ocorre quando o elétron incidente perde toda sua energia cinética em uma única colisão, e emitindo um fóton com energia correspondente. Logo, utilizando o como referência a equação do efeito fotoelétrico de Einstein (<math>E = hf\;</math>) temos que: <math> K = \frac{hc}{\lambda_{m\acute{i}n}}</math>; onde podemos escrever a energia cinética, adquirida devido a diferença potencial do tubo de raios-X, como (<math>K = eV\;</math>) e isolando o comprimento de onda mínimo obtendo:

<math> \lambda_{m\acute{i}n} = \frac{hc}{eV}</math>

     Os picos no gráfico devem a outro tipo de emissão chamada ''Emissão Característica'', pois a localização dos picos depende do material do anodo. Quando o elétron do feixe provindo do catodo penetra no anodo ele pode colidir com um elétron de uma subcamada mais interna fornecendo energia suficiente para expeli-lo do átomo, deixando o átomo excitado. Assim, para completar novamente a camada mais interna, um elétron de camada mais externa salta para a camada mais interna emitindo um fóton com energia igual a diferença de energia entre as camadas envolvidas. Esse processo de desexcitação pode ocorrer em várias etapas e com transições entre diferentes níveis de energia emitindo mais de um fóton representando diferentes picos no espectro de raios-X.

História dos Raios-X

2011-09-19T17:44:12Z

Tekkito: Criou página com '==A Busca ==      No ano de 1938, Michael Faraday realizou inúmeros experimentos com descargas elétricas em gases rarefeitos, associando seu nome à d...'

==A Busca ==

     No ano de 1938, Michael Faraday realizou inúmeros experimentos com descargas elétricas em gases rarefeitos, associando seu nome à descoberta dos raios catódicos. Porém, devido a dificuldades técnicas para a produção de vácuo de boa qualidade, as pesquisas só foram impulsionadas vinte anos depois.

     Por volta de 1958, Johann Hittorf (1824 – 1914) com seu mestre Julius Plücker (1801-1868) observou tubos com raios de grande energia que eram expelidos de um eletrodo negativo que, ao colidir com o vidro produziam uma certa fluorescência.

     Em 1876, Eugene Goldstein (1850-1931) interpretou que esses raios eram ondas no éter e denominou-os raios catódicos. Para o físico William Crookes, os raios catódicos eram moléculas carregadas aos quais constituíam o quarto estado da matéria. Em 1897, J.J. Thomson desvenda o mistério e demonstra que os raios catódicos são, na verdade, elétrons.

     Ao longo desses 40 anos, muitos pesquisadores estiveram muito próximos à descoberta dos raios-X.
Num artigo publicado em 1880, Goldstein menciona que uma tela fluorescente podia ser excitada, mesmo quando protegida dos raios catódicos.
Thomson também chegou perto; um ano antes da descoberta do raio-X, relatou que havia observado fosforescência em peças de vidros colocadas a centímetros de distância do tubo de vácuo.

     Crookes chegou a queixar-se da fábrica de insumos fotográficos Ilford, por lhe enviar papéis “velados”. Esses papéis, que protegidos da luz eram colocados próximos aos tubos de raios catódicos, velavam devido aos raios-X ali produzidos.

     Heinrich Hertz, em 1892, verificou que um pedaço de vidro de urânio envolto por finas folhas de alumínio e colocado dentro do [[tubo de raios catódicos]] ficava luminoso devido a estes atravessarem o alumínio. Na época Hertz estava muito ocupado para pesquisar mais sobre o fenômeno e deu permissão a seu aluno Philipp Lenard para aprofundar-se mais nesses estudos. Assim, Lenard construiu uma versão própria de um tubo de Crookes que continha uma pequena janela de alumínio por onde os raios catódicos podiam transpassar para o exterior do tubo permitindo o estudo da interação dos raios com diversos materiais fosforescentes.

     Após suas primeiras observações Lenard publicou que os raios sensibilizavam uma chapa metálica; um disco de alumínio previamente eletrizado era descarregado quando atingido pelos raios mesmo quando colocado a alguns metros do tubo, à distância muito maior do que o alcance dos raios catódicos no ar que é de aproximadamente poucos centímetros; quando a mão era colocada no trajeto dos raios o efeito de descarga elétrica desaparecia; e que alguns raios eram defletidos continuamente por um campo magnético em diferentes intensidades e outros nem chegavam a defletir.

     Acredita-se que os raios estudados por Lenard eram constituídos de raios catódicos (elétrons) e de raios-X, porém ele acreditava ser apenas raios catódicos. Caso Lenard tivesse utilizado uma janela de alumínio espessa suficiente para que os elétrons fossem bloqueados, observaria apenas um feixe de raios-X.

     Pode-se dizer que Lenard foi, entre todos os pesquisadores, o que mais se aproximou da descoberta de [[Röntgen]]. A descoberta dos raios-X estava muito próxima de ser atingida faltava apenas alguém com uma visão privilegiada para detectá-los.

== A Descoberta ==

     Os raios-X foram descobertos em 1895 por Wilhlem Conrad Röntgen, que recebeu o primeiro Prêmio Nobel de Física em 1901.

     Na noite do dia 8 de novembro de 1895, Röntgen estava em seu laboratório pesquisando sobre os raios catódicos. Observou que quando ionizava um gás, a extrema baixa pressão, contido no interior do tubo de raios catódicos com a ajuda de uma bobina de Ruhmkorff, uma tela coberta de platinocianeto de bário fluorescia.

     Röntgen tentou bloquear, em vão, essa luminescência de vários modos, cobrindo o tubo com cartolina preta e colocando os mais variados materiais com diferentes espessuras na região entre o tubo e a tela. Porém, posicionando um pedaço de chumbo perto da tela verificou sua sombra e o mais intrigante, observou a uma imagem escura no formato dos ossos de sua mão.

     Röntgen interpretou esse acontecimento como um novo tipo de raios que eram originados a partir da colisão dos raios catódicos com a parede do tubo. Desconhecendo a natureza de tais raios denominou-os raios “X”.

     Seis semanas depois da descoberta, convidou sua mulher Bertha para ir ao laboratório e tirar uma radiografia dos ossos de sua mão esquerda com seu anel de casamento claramente visível, gerando uma das mais famosas imagens na história da fotografia.

     Em menos de oitos semanas Röntgen já havia descoberto a maioria das propriedades fundamentais desses raios e em 1897 anos abandonou os estudos nessa área para preocupar-se em estudos de seu maior interesse. Em 28 de dezembro 1896 envia ao presidente da Sociedade de Física e Medicina de Würzburg o primeiro de um total de três trabalhos sobre o assunto, o manuscrito “Sobre um novo tipo de raios” ao qual considerou uma comunicação preliminar. O segundo trabalho com o mesmo nome do anterior foi enviado em 9 de março de 1896 e o terceiro e último trabalho é datado de 10 de março de 1897.

     A descoberta de Röntgen teve imediata aplicação na área da medicina sendo as primeiras imagens de ossos fraturados feitas em janeiro de 1896 antes mesmo da descoberta da natureza desses raios.

Raio X

2011-09-19T17:43:56Z

Tekkito: Criou página com '==História dos Raios-X== ==Produção e emissão de Raios-X== ==Natureza dos Raios-X== ==Lei de Bragg== ==Emissão e absorção de raios X== *Origem da emis...'

==[[História dos Raios-X]]==

==[[Produção e emissão de Raios-X]]==

==[[Natureza dos Raios-X]]==

==[[Lei de Bragg]]==

==Emissão e absorção de raios X==

*Origem da emissão;

*Espectroscopia de Raios-X;

==Utilização dos raios X==

*Absorção de raios X;

*Medicina;

*Raios-X;

*Tomografia;

*Metalurgia;

*Difração de Raios-X;

*Ciência dos Materiais;

*Espectroscopia de Raios-X;

*Análise de composição de materiais;

*PIXE, EIXE;

==Segurança com raios X==

Ondas

2011-09-19T17:43:04Z

Tekkito: Criou página com 'Os fenômenos ondulatórios têm um papel destacado na natureza. São inúmeros os fenômenos da natureza que podem ou necessitam ser descritos por modelos ondulatórios. [[Im...'

Os fenômenos ondulatórios têm um papel destacado na natureza. São inúmeros os fenômenos da natureza que podem ou necessitam ser descritos por modelos ondulatórios.

[[Image:Scheila6.jpg|thumb|Ondas Bidimensionais|300px|right|Fig.1: Ondas na água]]
Podemos observar facilmente ondas nas superfícies da água em um copo, na piscina ou no mar. Podemos até surfar sobre algumas destas ondas. Também nos assustam as descrições de um tsunami ou os tremores devido à onda gerada por um terremoto. Algumas categorias de fenômenos ondulatórios foram e são essenciais para a evolução e sustentação da vida sobre a terra e estão presentes na maioria das modernas tecnologias. Entre esses se destacam todos aqueles que se relacionam com a transmissão de informação e/ou de energia, como o som, a luz e mais genericamente, todas as formas de ondas eletromagnéticas.

[[Imagem:Peoplewave.gif|thumb| "La Ola"|300px|left| La Ola - "Animation courtesy of Dr. Dan Russell" http://www.kettering.edu/~drussell/Demos/waves-intro/waves-intro.html, Kettering University]]
"Onda" é uma perturbação do meio que se propaga espacialmente por algum tempo. Esta propagação deverá, portanto ser descrita por uma função de variáveis espaciais e temporal. Numa onda mecânica as partículas do meio oscilam em torno de sua posição equilíbrio quando a perturbação passa por elas. O que se propaga no espaço é apenas a perturbação do meio, não havendo transporte de matéria, apenas a energia é propagada pela onda. Observamos isso quando deixamos cair uma pedra numa lagoa parada. Veremos uma perturbação que se propaga radialmente sobre a superfície da água. Após a passagem da perturbação a água volta para a sua condição de repouso. Num estádio de futebol a torcida faz a "ola", onde cada pessoa reproduz em seqüência o movimento do colega e retorna ao estado que estava antes. Se provocarmos um pulso numa corda, também observaremos a perturbação passar e a corda voltar ao estado inicial. Então, durante a passagem da onda, cada ponto oscila em torno de sua posição de equilíbrio.

Há algumas '''[[Propriedades e Características|propriedades e características]]''' que precisamos conhecer a respeito de ondas para iniciar o entendimento desses fenômenos.

A principal assinatura dos sistemas ondulatórios resulta da propriedade da '''[[superposição]]''' das ondas, da qual derivam duas de suas propriedades mais características, a '''[[interferência]]''' e a '''[[difração]]'''. É também baseado na superposição que o '''[[Princípio de Huygens|princípio de Huygens]]''' explica satisfatoriamente o comportamento das ondas constituindo-se como base das modernas teorias da difração.

A experiência da '''[[Fenda Dupla|fenda dupla]]''' de Young é fundamental para a compreensão das idéias da física moderna como, por exemplo, o comportamento dual da luz. Assim como o estudo das fendas múltiplas, as '''[[redes de difração]]''', que possibilitam o conhecimento a cerca do espectro da luz.

Por meio da interferência e da difração é possível entender o funcionamento dos instrumentos ópticos como lunetas, telescópios e máquinas fotográficas e digitais (o '''[[poder de resolução]]''' destes instrumentos, devido à difração); a continuidade nas imagens de TV, que resulta também da difração; as cores que são geradas por interferência em bolhas de sabão, asas de borboleta, película de óleo; entre outros fenômenos naturais e tecnológicos presentes no nosso quotidiano.

=='''''[[Faça em aula!]]'''''==
Aqui você encontra diversas demonstrações para tornar a aula mais interativa.

=='''''[[Tem que resolver!]]'''''==
Testes seus conhecimentos resolvendo os exercícios.

=='''''[[Para visitar!]]'''''==
Sugestão de sitios interessantes com textos, imagens e aplicativos.

=='''''Textos complementares.'''''==
# [[Matemática da Onda]]
# [[Matemática da Interferência]]
# [[Matemática da Difração]]
# [[Laser]]
# [[A Fisica do Surf]]
# [[Difração de Raio X]]

=='''''[[Referências]]'''''==
Para desenvolver este material consultamos...

=='''''Sumário'''''==
# [[Propriedades e Características]]
# [[Superposição]]
# [[Difração]]
##[[Fendas simples]]
##[[Fios]]
##[[Orifícios]]
##[[poder de resolução|Poder de resolução]]
##[[Formação das imagens]]
##[[Difração nas Sombras]]
#[[Interferência]]
##[[Interferência de filmes finos]]
##[[Interferômetro de Michelson]]
##[[Interferômetro de Mach Zender]]
# [[Princípio de Huygens]]
# [[Fenda Dupla]]
# [[Coerência]]
# [[Redes de difração]]
# [[Espectroscópio]]

=='''''[[Textos de Apoio ao Ensino de Física]]'''''==

Redes

2011-09-19T17:42:08Z

Tekkito: Criou página com '== Redes Aleatórias == <math>G_{n,m}</math>: conjunto dos grafos de <math>n</math> vértices e <math>m</math> arestas (se escolhe aleatoriamente um par de vértices e se conec...'

== Redes Aleatórias ==

<math>G_{n,m}</math>: conjunto dos grafos de <math>n</math> vértices e <math>m</math> arestas (se escolhe aleatoriamente um par de vértices e se conectam ate completar <math>m</math>)
versão microcanônica (sem flutuação). Existe probabilidade de mais de uma aresta entre vértices, porem probabilidade é <math>m/n^2</math>

O relação de arestas/vértices é <math>m/n</math>

----

<math>G_{n,p}</math>: conjunto dos grafos de <math>n</math> vértices onde cada aresta (das <math>n(n-1)/2</math> possíveis) é estabelecida com probabilidade <math>p</math> - versão canônica (com flutuação). Nesta versão não é possível definir duas vezes a mesma aresta.

O número total esperado de arestas é portanto <math>m \approx pn(n-1)/2</math>, que para <math>n</math> grande resulta ser <math>m = pn^2/2</math>. Assim o a relação arestas/vértices é <math>m/n = pn/2</math>

O número de médio de arestas por vértice é <math>z = <k> = pn</math>

<math>G_{n,p}</math> tem distribuição de grau binomial: <math>p_k = \frac{n}{k}p^k(1-p)^{n-k}</math>

No limite <math>n</math> grande é a distribuição de Poisson:

<math>p_k = \frac{z^k e^{-z}}{k!}</math>

O resultado mais notável dos trabalhos de Erdös e Rényi é que quando <math>n \rightarrow \infty</math> para <math>z = 1</math> ha uma transição onde aparece uma componente gigante. <math>z=1 \Rightarrow pn=1 \Rightarrow p = 1/n \Rightarrow m=n/2
</math>

== Links ==
[http://ccl.northwestern.edu/netlogo/ Netlog]

--[[Usuário:Sebas|Sebas]] 15:38, 25 Maio 2008 (BRT)
----

[[FIP20706 Tópicos em Física Estatística: redes e dinâmica de epidemias|voltar]]

FIP20706 Tópicos em Física Estatística: redes e dinâmica de epidemias

2011-09-19T17:41:52Z

Tekkito: Criou página com '== Epidemias == === Modelo SIR (suscetível-infectado-removido) === A equação de partida do modelo expressa a taxa de conversão de sadio em infectado como resultado do conta...'

== Epidemias ==

=== Modelo SIR (suscetível-infectado-removido) ===

A equação de partida do modelo expressa a taxa de conversão de sadio em infectado como resultado do contato com estes:

:<math>\frac{dS}{dt} = -\alpha I S</math>

onde '''S''' e '''I''' são o número de suscetíveis e infectados respectivamente, que somados aos removidos '''R,''' conformam o total de indivíduos da população''' N'''.
Ou seja

:<math>N = S + I + R</math>

Neste primeiro modelo a quantidade '''N''' é conservada, ou seja não a entrada nem saída de indivíduos, nem mortes (ao não ser pela própria doença que são os R) nem nascimentos. Se bem limitado e o marco para começar os estudos de epidemia além disso pode se considera que se a doença se desenvolve rápido as taxas de nascença e mortalidade ou migração podem ser omitidas.

Voltando à primeira equação ela explicita o fato dos suscetíveis saírem desse estado numa taxa proporcional ao número de indivíduos infectados imersos na população. Isto é razoável, se há dois infectados espalhados, a taxa de infecção deve dobrar. Mas se '''S''' dobra (o que se '''I''' é proximo de 1 significa que '''N''' dobra) não significa que o infectado vai conseguir infectar com o dobro de rapidez. Isto significa que <math>\alpha</math> na primeira equação depende de '''N''' inversamente.

De outra forma:
:<math>\alpha = \beta/ N</math>, onde <math>\beta</math> sim é o parâmetro de infecciosidade que só depende da doença.

Rescrevemos a primeira equação usando este parâmetro assim:

:<math>\frac{dS}{dt} = -\beta I \frac{S}{N}</math>

=== [[Redes]] ===

Histogramas e Densidade de Probabilidade

2011-09-19T17:40:31Z

Tekkito: Criou página com 'Geradores aleatórios tem amplo uso em simulações físicas, em particular nas simulações da área de Mecânica Estatística onde são usados para imitar o efeito da temperatu...'

Geradores aleatórios tem amplo uso em simulações físicas, em particular
nas simulações da área de Mecânica Estatística onde
são usados para imitar o efeito da temperatura. Como mencionado nas
seções anteriores, as linguagens de programação usualmente apresentam
um gerador de números aleatórios (congruente) intrínseco, supostamente
uniforme, ou seja, um gerador que produz com igual probabilidade os números
aleatórios dentro do intervalo de definição do gerador. (De fato, há uma segunda
condição importante para um bom gerador que é mais difícil de ser obtida
em geradores congruentes: '''a descorrelação''', mas não vamos tratar dela agora.)

A primeira condição pode ser testada na prática fazendo-se um histograma
dos números aleatórios gerados. Este é o objetivo inicial desta seção.

1) Faça um programa que gere <math>N</math> números aleatórios inteiros
no intervalo <math>0<x<L</math>. Divida esse intervalo em <math>M</math> partes
e conte o número de números aleatórios, <math>h(i)</math>, gerados em cada divisão,<math> i</math>
deste intervalo. Use, por exemplo, <math>L=100</math> e teste diferentes valores
de <math>N=100, 1000, 10000</math>. Faça um gráfico (histograma) de <math>h(i)\times i</math>.

2) Calcule o valor médio e o desvio quadrático médio de <math>h</math> para os três casos
anteriores. Faça um gráfico do desvio quadrático médio como função de <math>L</math>.

3) Calcule o valor médio e o desvio quadrático médio dos números aleatórios gerados. Note
que pode-se calcular esses valores de duas formas:
\bar

'''Densidade de Probabilidade'''

Os resultados do item 1 do exercício acima sugerem a forma de calcular a probabilidade <math>P(x)</math>
de que um número aleatório <math>x</math> seja gerado dentro do intervalo <math>0\leq x \leq i</math>:

<math>P(x_i)=\frac{1}N{\sum_{j=1}^{j<i+1} h(x_j)}</math>.

Pode-se generalizar o cálculo de P(x) para um intervalo <math>\Delta x</math> qualquer usando
o conceito de densidade de probabilidade:

<math> P(x)=\int_0^x\rho(x)dx</math>

assim, <math>\rho(x)</math>, é a densidade de probabilidade e <math>\rho(x)dx</math> é a probabilidade que <math>x</math> seja sorteado no intervalo <math>[x,x+dx]</math>.

A probabilidade de que um número aleatório seja sorteado dentro do intervalo de definição do gerador
deve ser unitária, portanto,

<math>P(L)=\int_0^L \rho(x)dx =1</math>.

que é o que se chama de normalização.
Um gerador uniforme deve apresentar os números aleatórios com igual probabilidade
em qualquer região de seu intervalo de definição, ou seja, a densidade de probabilidade
deve ser constante.

<math>\rho(x)=A \;\;\;\;\;\; \{0<x<L\}</math>

Podemos encontrar <math>A</math> normalizando-se a expressão acima:

<math>\int_0^L \rho(x)dx = \int_0^L A dx = 1 \;\;\;\;\ \Rightarrow A=\frac{1}{L}</math>.

'''Valores Médios de Distribuições Aleatórias'''

O cálculo dos valores médios de variáveis aleatórias ou potências dessas variáveis é, por
definição,

<math>\bar{x^k} = \frac{1}N\sum_{i=1}^{i=N} x^k</math>,

mas, lembrando o histograma que fizemos no início desta seção, podemos equivalentemente usar a contagem
de números que caem em cada intervalo multiplicada por <math>x^k</math> e somar sobre os <math>M</math> valores
em que partimos o intervalo <math>L</math>.

<math>\bar{x^k} = \frac{1}N\sum_{i=1}^{M} h(i) x(i)^k</math>,

lembre que <math>\frac{h(i)}N</math> é a probabilidade que um número aleatório caia no intervalo entre
i e i+1. Podemos estender esse cálculo para limite em que o intervalo passa a ser infinitesimal <math>dx</math>. A
probabilidade de que um número aleatório aí caia é dada por <math>\rho(x)dx</math> e a soma passa a uma integral,

<math>\bar{x^k}=\int_0^L \rho(x) x^k dx</math>.

O valores médios, <math>\bar{x^k} </math>, são chamados de '''momentos''' da distribuição e o conjunto desse
valores (para <math>1<k<\infty</math>) define-a univocamente.

'''Exercícios'''

1- Mostre que o k-ésimo momento de uma distribuição uniforme definida em um intervalo L é dada por

<math>\bar{x^k} =\frac{L^k}{k+1}</math>.

2- Faça um programa que calcule numericamente o resultado obtido no ítem anterior.