A goal is a dream with a plan : Ajuste de curvas

INTERPOLACIÓN DE DATOS

Se ingresan los valores medidos de los aguacates durante 6 dias en matrices para su fácil manejo, y se busca un polinomio que se ajuste a la tendencia

clear all
close all
clc

Aguacate uno

MAT_1=[
    209.03  191.63  192.84  192.07  191.46  187.38 % Blanco
    133.36  105.41  112.92  110.96  110.76  105.38% Azul
    188.81  165.47  165.20  163.79  163.84  156.64% Verde
    187.38  164.81  165.75  166.02  164.39  160.65% Rojo
    20.69   18.57   18.70   17.77   17.44   17.02];% Negro

Agucate dos

MAT_2=[
    206.63  195.98  194.47  189.39  186.99  183.35 % Blanco
    130.77  109.38  110.63  104.20  103.98  96.64 % Azul
    183.95  166.27  161.75  154.13  153.42  147.95 % Verde
    185.91  174.38  173.50  167.63  162.85  159.29 % Rojo
    20.80   19.34   19.31   17.90   17.41  16.33]; % Negro

Aguacate tres

MAT_3=[
    208.31  190.52  199.68  195.32  190.11  185.55 % Blanco
    132.64  108.14  121.72  116.17  110.10  99.10  % Azul
    187.21  162.19  175.37  170.38  162.27  154.77 % Verde
    187.10  163.60  173.80  168.62  163.54  159.59 % Rojo
    20.94   19.01   19.92   18.32   17.63   16.97]; % Negro

Corrección día 3

MAT_1(:,3)=(MAT_1(:,2)+MAT_1(:,4))/2;
MAT_2(:,3)=(MAT_2(:,2)+MAT_2(:,4))/2;
MAT_3(:,3)=(MAT_3(:,2)+MAT_3(:,4))/2;

Porcentaje de cambio en verde

Ag_1=MAT_1(3,:)./(MAT_1(1,:).*MAT_1(5,:));
Ag_2=MAT_2(3,:)./(MAT_2(1,:).*MAT_2(5,:));
Ag_3=MAT_3(3,:)./(MAT_3(1,:).*MAT_3(5,:));
Ag_P=(Ag_1+Ag_2+Ag_3)/3;
dia=[1,2,3,4,5,6];
figure(1)
plot(dia,Ag_P,'--dg','Linewidth',2),title('Comportamiento verde','FontSize',12,'FontWeight','bold'),xlabel('Dias'),ylabel('Valores');

Gráfica de los tres agucates y su promedio

figure(2)
plot(dia,Ag_1,'--*b',...
     dia,Ag_2,'--or',...
     dia,Ag_3,'--pm',...
     dia,Ag_P,'--^g'),title('Cambio de color a traves de los días','FontSize',12,'FontWeight','bold'),xlabel('Dias'),ylabel('Valores');
legend('Aguacate uno ','Aguacate dos ','Aguacate tres ','Promedio','Location','BestOutside');

Polinomio y su ajuste

horas=1:0.1:6;

P son los coeficientes del polinomio de 5to grado S es el grado aproximadoo de error

[P,S]= polyfit(dia,Ag_P,5)

figure(3)
plot(dia,Ag_P,'--^b',...
    horas,polyval(P,horas),'--*g'),title('Ajuste del polinomio','FontSize',12,'FontWeight','bold'),xlabel('Dias'),ylabel('Valores');
legend('Promedio','Polinomio grado 5','Location','BestOutside');

P =

    0.0000   -0.0001    0.0008   -0.0034    0.0078    0.0380


S = 

        R: [6x6 double]
       df: 0
    normr: 6.3216e-17

A goal is a dream with a plan

martes, 15 de marzo de 2016

Ajuste de curvas