Logo Studenta

Ajuste_No_lineal

Esta es una vista previa del archivo. Inicie sesión para ver el archivo original

clc
clearvars
%Datos%
datos=[1 99;2 95;5 85;15 55;25 30;30 24;35 20;40 15];
x=datos(:,1);
y=datos(:,2);
%Ajuste por modelo% %sistema de ecuaciones lineales asociado al problema de encontrar la funcion exponecial (enunciado)%
A=[ones(size(x)),x]; %MISMA FUNCIÓN QUE PONES EN EL FUN%
b=y;
c=A\b; %Como el sistema NO es rectangular (sin determinante y sin inversa), MATLAB usa MINIMOS CUADRADOS para ajustar una funcion %
%Hacer funcion que simula un polyval%
fun=@(t) c(1)+c(2).*t; 
xx=linspace(min(x)-0.1,max (x)+0.1,250); %250 es un largo arbitrario%
hold on
%ploteo datos%
plot(x,y,'mo');
%ploteo curva%
plot(xx,fun(xx),'b--'); %FUN ES MI NUEVO POLYVAL%

Otros materiales