Versão 1.0 Objetivo: Ajusta uma equação teórica a dados experimentais. Síntaxe: [CuAj PaOt]=ajuste(dados,função,starting,varsym), onde 'dados' são os dados experimentais, 'função' a função teórica, 'starting' o valor inicial de cada parâmetro livre (neste caso um vetor linha), e, finalmente, 'varsym' o mesmo vetor linha anterior, porém com as variáveis simbólicas. NOTA IMPORTANTE: a variavel independente da 'função' deve, obrigatóriamente, se chamar 'xt'. Como resultado, CuAj é uma curva discreta obtida a partir da 'função', com os valores obtidos da otimização executada. Faça plot(dados(:,1),dados(:,2),'o',CuAj(:,1),CuAj(:,2)) para ver a diferença. dados(:,1)=CuAj(:,1). Ainda, PaOt são os parâmetros otimizados. Exemplo: v=[0.1:0.1:5]';u=tanh(v)+diag(rand(length(v))/10); exp=[v u]; syms a b c xt real; teo=a*tanh(b*xt)+c; [CuAj PaOt]=ajuste(exp,teo,[2 2 2],[a b c]); Observações: 1. Esta rotina somente funciona com outra rotina: ajuste_rot.m 2. Caso a 'função' possua magneton de Bohr e constante de Boltzmann, estes já estão inseridos na rotina e devem se chamar 'mibor' e 'kb'. Valores utilizados:mibor=9.27e-21 (erg*Oe-1) e kb=1.38e-16 (erg*K-1) %%%%%%%%%%%%%%%%%% Rotina otimizada para MATLAB 7.1 Rotina elaborada por Mario Reis (marior@if.uff.br). Favor comunicar possíveis erros e/ou 'bugs'. Histórico de alterações: Versão 1.0 (Julho 2009) *ORIGINAL*