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

Bubble Sort em Scilab

WhiteHawk

Função Scilab que ordena um vetor pelo método Bubble Sort.

Por: White Hawk

PJEOffice - Baixa automaticamente última versão do CNJ (Conselho Nacional de Justi&cce

mamp

Repositorio Nacional do Conselho de Justiça criou uma nova versao PjeOffice para Processo eletronico PJE. Link oficial do repositório dá muitas instrucões compiladas na internet sem muita efetvidade. Mas este script é funcional, simples e testado. Para rodar script com código abaixo. Detalhe importante: para evocar, o administrador não foi cuidadoso com nomes dos arquivos e mudou ... vou fazer artigo para simplificar uso do PjeOffice, mas para chamar Pje precisa fazer o seguinte: no Diretório que vc rodou script cd pje-office e digite ./pjeOffice.sh Sim o administrador do CNJ mudou letrinhas maiúsculas do diretório e do programa intencionalmente, mas funcina.

Por: Marcos Moraes

Octave - Sistemas lineares por Jacobi

danielms

O método de Jacobi trata-se dum algoritmo para determinar a solução de um sistema de equações lineares com os maiores valores absolutos em cada linha e coluna dominados pelo elemento da sua diagonal. Trata-se duma versão simplificada do algoritmo de valores próprios de Jacobi.

Por: Daniel Moreira dos Santos

Octave - Calcular raiz pelo método de Newton

danielms

Em análise numérica, o método de Newton (ou método de Newton-Raphson) tem o objetivo de estimar as raízes de uma função. Para isso, toma-se um ponto qualquer da função, calcula-se a equação da tangente (derivada) da função nesse ponto, calcula-se o intercepto da tangente ao eixo das abcissas, calcula-se o valor da função nesse ponto, e repete-se o processo, que deve tender a uma das raízes da função rapidamente, ou não tender a nada, deixando isso claro logo.

Por: Daniel Moreira dos Santos