Fazendo um ajuste não linear em dados experimentais - FORTRAN 90

iagolira

Olá pessoal, como sou amante do Fortran, resolvi criar um programa que faz um ajuste não linear em dados experimentais. - Eu sei, já existe programas para tais! A grande utilidade é quando usa-se muitos parâmetros a serem determinados, o que é vantajoso em relação aos demais. No programa existe a função PLOT, onde nesta dá-se a entrada da função a fazer o ajuste. Exemplo: FUNCTION PLOT(X,A,Qp) IMPLICIT NONE INTEGER, PARAMETER :: qpl = selected_real_kind(15, 16) INTEGER :: Qp real(KIND=qpl) :: X, A(Qp), PLOT PLOT=1.0_qpl+A(1)*exp(-A(2)*X)*cos(A(3)*X)-& &A(4)*exp(-A(5)*X)*cos(A(6)*X)-& &A(7)*exp(-A(8)*X)*cos(A(9)*X)+& &A(10)*exp(A(11)*X)*cos(A(12)*X)-& &A(13)*exp(-A(14)*X)*cos(A(15)*X) END FUNCTION PLOT Percebe-se que a função PLOT têm 16 parâmetros a serem determinados, então percebe-se que é fácil entrar com os valores. COMPILANDO No terminal digite: gfortran Fit.Date.f90 -o FitDate.x -O3 O "-O3" é opcional, pois é um parâmetro de otimização. EXECUTANDO Ainda no terminal, digite: ./FitDate.x Então, aparecerá uma tela pedidos os arquivo que contém os dados a serem analisados, a quantidade de parâmetros e o erro que você quer cometer. Quanto menor o erro, mais demorado. No programa existe uma variável chamada de "tol" (PARAMETER(tol=0.000000000001)), esta é a precisão do cálculo, então ajuste para suas necessidades.

Por: Iago Lira

Método das Cordas em SCILAB

arieldll

Método das Cordas/Contraposição em SCILAB.

Por: Ariel Galante Dalla Costa

Método de Jacobi em SCILAB

arieldll

Solução computacional para sistemas de equações lineares.

Por: Ariel Galante Dalla Costa

Matriz de Hilbert e resolução de sistemas lineares

danielms

Uma matriz de Hilbert é uma matriz quadrada com a seguinte forma: H_{ij} = \frac{1}{i+j-1} Aqui vamos mostrar como resolver um sistema linear HX=b utilizando o octave, onde H é uma matriz de Hilbert.

Por: Daniel Moreira dos Santos

Controle de maior e menor de idade em Ruby

andrezc

Controle de maior/menor de idade para entrar em uma festa, escrito em Ruby. Faz uso de if / else.

Por: André