opts = spreadsheetImportOptions("NumVariables", 32);
% Specify sheet and range
[Link] = "Hoja1";
[Link] = "A2:AF65";
% Specify column names and types
[Link] = ["Ajustar", "MarcaDeTiempo", "O2", "CO2", "COppm",
"COnppmn", "CO", "NOppm", "NOnppmn", "NOppm1", "NOnppmn1", "NOppm2", "NOnppmn2",
"NOppm3", "SOppm", "HSppm", "HC", "FlueDraughtPa", "T1C", "T2C", "NetTempC",
"AmbientTempC", "TotalLosses", "DryLosses", "WetLosses", "COLosses",
"PoisonIndex", "COCORatio", "Efficiency", "ExcessAir", "Lambda",
"BarometricPressuremBar"];
[Link] = ["double", "datetime", "double", "double", "double",
"double", "double", "categorical", "categorical", "categorical", "categorical",
"categorical", "categorical", "categorical", "categorical", "categorical",
"double", "categorical", "categorical", "categorical", "categorical", "double",
"double", "double", "double", "double", "double", "double", "double", "double",
"double", "double"];
% Specify variable properties
opts = setvaropts(opts, ["NOppm", "NOnppmn", "NOppm1", "NOnppmn1", "NOppm2",
"NOnppmn2", "NOppm3", "SOppm", "HSppm", "FlueDraughtPa", "T1C", "T2C",
"NetTempC"], "EmptyFieldRule", "auto");
% Import the data
emision = readtable("C:\Users\lab\MATLAB Drive\[Link]", opts,
"UseExcel", false)
clear opts
emision(1,:)=[]
Se acortar los nombres de variables a trabajar para agregar en la figura 3.
CO = [Link];
CO2 = emision.CO2;
HC = [Link];
O2 = emision.O2;
lambda=[Link]
En este punto las graficas de CO no resultan funcionales para el analis, ya que
dentro de su matriz existe datos vacios que afectan al calculo.
figure,scatter([Link],[Link],5,"black") %cambiar parte y color
xlabel('dosado')
ylabel('CO (%vol')
grid on
En este apartado se procede en cada emision ha analizar eliminar de la matriz
datos NAN, de esta forma se completa los datos al final para mejorar el calculo
de CO.
% Fill missiong data
[emision2, missingIndices]= fillmissing(emision,"linear",...
"DataVariables",@isnumeric);
figure,scatter([Link],[Link],5,"black")
xlabel('dosado relativo')
ylabel('CO(%vol)')
grid on
lam = 0.6:0.05:1.4
CO_ajust = 30.5*lam.^2-59.07*lam+28.96;
figure,plot(lam,CO_ajust)
figure,scatter([Link],emision.CO2,5,"filled") %cambiar parte y color
xlabel('dosado')
ylabel('CO2 (%vol')
grid on
En este apartado se procede en cada emision ha analizar eliminar de la matriz
datos NAN, de esta forma se completa los datos al final para mejorar el calculo
en CO2.
% Fill missiong data
[emision2, missingIndices]= fillmissing(emision,"linear",...
"DataVariables",@isnumeric);
figure,scatter([Link],emision2.CO2,5,"black")
xlabel('dosado relativo')
ylabel('CO2(%vol)')
grid on
lam = 0.6:0.05:1.4
%CO2_ajust = -40*lam.^2+61.33*lam-12.22;
CO2_ajust = 121.2*lam.^3-416.9*lam.^2+442.8*lam-138.2;
figure,plot(lam,CO2_ajust)
figure,scatter([Link],[Link],5,"filled") %cambiar parte y color
xlabel('dosado')
ylabel('HC (ppm')
grid on
En este apartado se procede en cada emision ha analizar eliminar de la matriz
datos NAN, de esta forma se completa los datos al final para mejorar el calculo
en HC.
% Fill missiong data
[emision2, missingIndices]= fillmissing(emision,"linear",...
"DataVariables",@isnumeric);
figure,scatter([Link],[Link],5,"cyan")
xlabel('dosado relativo')
ylabel('HC (ppm)')
grid on
lam = 0.6:0.01:1.4
HC_ajust = 3.099*exp(04)*lam.^2-5.121*exp(04)*lam+2.749*exp(4);
plot(lam,HC_ajust)
%HC1_ajust = -3.567*exp(1)+16*lam.^3+8.91*exp(1)+16*lam.^2-
2.749*exp(1)+16*lam+1.998*exp(1)+16;
%figure,plot(lam,HC1_ajust)
%hold off
figure,scatter([Link],emision.O2,5,"filled") %cambiar parte y color
xlabel('dosado')
ylabel('O2 (%vol')
grid on
En este apartado se procede en cada emision ha analizar eliminar de la matriz
datos NAN, de esta forma se completa los datos al final para mejorar el calculo
en O2.
% Fill missiong data
[emision2, missingIndices]= fillmissing(emision,"linear",...
"DataVariables",@isnumeric);
figure,scatter([Link],emision2.O2,5,"cyan")
xlabel('dosado relativo')
ylabel('O2 (%vol)')
grid on
lam = 0.6:0.05:1.4
O2_ajust = 78.7*lam.^2-114*lam+41.8;
figure,plot(lam,O2_ajust)
Ventana Lambda
lam = 0.6:0.05:1.4
CO_ajust = 30.5*lam.^2-59.07*lam+28.96;
CO2_ajust = 121.2*lam.^3-416.9*lam.^2+442.8*lam-138.2;
HC_ajust = 3.099*exp(04)*lam.^2-5.121*exp(04)*lam+2.749*exp(4);
O2_ajust = 78.7*lam.^2-114*lam+41.8;
ax1 = nexttile;
xlabel(ax1,'\lambda')
yyaxis left
plot(lam,CO_ajust,'k'); hold on;
plot(lam,CO2_ajust,'g');
plot(lam,O2_ajust,'b');
ylabel(ax1,'vol(%)')
axis([0.6 1.4 0 100])
yyaxis right
plot(lam,HC_ajust,'m');
ylabel(ax1,'HC (ppm)')
legend('CO','CO2','HC','O2')
grid
