Nothing Special   »   [go: up one dir, main page]

Primerofactir de Idealidad

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 5

TECNOLÓGICO NACIONAL DE MÉXICO

INSTITUTO TECNOLÓGICO DE OAXACA

DEPARTAMENTO DE INGENIERÍA ELÉCTRICA

OSTENTACION DE UN SCRIBD EN OCTAVE PARA EL ANALISIS


DEL COMPORTAMIENTO DE UNA CELULA FOTOVOLTAICA.
(FACTOR DE IDEALIDAD).

PRESENTA:
VICENTE ZARATE MILLAGUI BONIFACIO

N° CONTROL:
20160390

TECNOLOGIA DE LA CELDA FOTOVOLTAICA


7EB

A CARGO DE:
ING. PEREZ ARREORTUA NOÉ

LUBES 11 DE DICIEMBRE DE 2023.


Esta es una breve descripción de lo que puede realizar el siguiente script
elaborado en Octave, su objetivo está basado principalmente en el análisis de la
curva corriente-voltaje de las células solares. considerando y dando la primordial
relevancia a los parámetros eléctricos. Empleando el método de Newton-
Raphson, método numérico para poder resolver ecuaciones no lineales, también
utiliza un ajuste polinomial que en seguida se describe.

Este código a un inicio está definido para los parámetros eléctricos de la célula
1M156B5 de una marca nombrada “MOTECH”, referido con un código de
eficiencia de 180.

Definir parámetros
Para el análisis de la curva corriente voltaje, es necesario considerar y definir
nuestras constantes, es el ejemplo, de la corriente a cortocircuito, Vms, Vocs

%parametros electricos de la celula (180)

Vms=0.532;

Vocs=0.629;

Ims=8.31;

Iscs=8.79;

k=1.3806e-23;

Ts=25;

q=1.602e-19;

Cálculo de coeficientes
para avanzar en el cálculo de nuestra gráfica, calculamos nuestras C´s,

utilización de un método numérico para resolver ecuaciones no lineales

Para el caso, se empleó el método de newton Raphson, así para encontrar la


raíz de la ecuación F(x)=0, mismo resultado que fue almacenado en x1

% Expresiones definidas

Tcs=Ts+273;

Vtcs=(k*Tcs)/q;
C1=Vms/Vtcs;

C2=Vocs/Vtcs;

C3=1-(Ims/Iscs);

Valores necesarios
primero se calculó el factor de corrección voltaje nombrado (An), así también la
corriente de saturación inversa. Para ello definimos un rango para graficar.

%Valor inicial para graficar

%An = x0

x0 = [.5:.001:.8];

Iniciamos definiendo nuestra función.

con un valor inicial, para este caso: 0.77, con una tolerancia de t = 0.000000001.

así también ordenamos que limite el rango de iteraciones a I=1000


%Definimos la función

f0=(exp(C1 * x0)) - C3 *(exp(C2 * x0)) - 1 + C3;

%Valor inicial

xi =0.77;

%Tolerancia

t = 0.000000001;

%Definimos el numero maximo de iteraciones

I=1000

Con el comando Plot graficamos con el título “Factor de idealidad” tomando a x0 y f0

%Graficamos

plot(x0,f0,'g','linewidth',.8)

title("Función del factor de idealidad")

En la Figura 1 podemos ver el comportamiento de la función


Figura 1. Función del factor de idealidad.
Ahora bien, aplicando la derivada a nuestra función tenemos que:

%Aplicando la función

f=(exp(C1 * xi)) - C3 *(exp(C2 * xi)) - 1 + C3;

%Derivada de la función

df= C1 * (exp(C1 * xi)) - C3 * C2 * (exp(C2 * xi));

y esto a su vez para i: a 1000 iteraciones

for i = 1 : I

f=(exp(C1 * xi)) - C3 *(exp(C2 * xi)) - 1 + C3;

df= C1 * (exp(C1 * xi)) - C3 * C2 * (exp(C2 * xi));

Aplicando el método de Newton-Raphson para encontrar la raíz función real


derivable, La idea básica del método de Newton-Raphson es utilizar la derivada
de la función para aproximarse a la raíz. Se parte de una suposición inicial para
la raíz y luego se utiliza la derivada para calcular una corrección que se aplica a
la suposición inicial para obtener una mejor aproximación. Este proceso se repite
hasta que se alcanza la convergencia, es decir, hasta que la diferencia entre dos
aproximaciones consecutivas es menor que un valor de tolerancia
preestablecido.
%Aplicamos el metodo de Newton-Raphson

x = xi - f/df;

dist = abs (x-xi);

partiendo del método de Newton tenemos entonces que si


if dist < t

disp('El sistema converge')

%El break se utiliza para salir de la evaluación

break

end

if f < t

disp('El sistema Converge')

end

xi = x;

end

dado que nuestro sistema converge imprimimos nuestro factor de corrección de


voltaje
An = 1/xi

También podría gustarte