Estoy tratando de modelar una parte de la red eléctrica (sí, ya he revisado los foros de MATLAB y encontré la solución de flujo de carga), pero esto es más bien un problema de modelado para la estimación del estado.
Tengo un conjunto de ecuaciones diferenciales, no lineales, que representan una estructura de 5 buses (3 generadores y 2 cargas). Tengo un vector de estado que acumula todas las variables de mi estado de interés.
Intentando investigar en los foros, ¿cuál es la forma adecuada de modelar este sistema? ¿Es simplemente declarar un sistema de ecuaciones que luego puede ser resuelto por ODE45? Por favor, vea un código de ejemplo en la parte inferior
También se debe solicitar una subpregunta. ¿Es posible modelar un sistema no lineal en MATLAB? O tenemos que usar Simulink.
Otra subpregunta. ¿Estoy confundiendo el modelado del sistema con la resolución de las ecuaciones diferenciales? ¿Hay una diferencia? Realmente apreciaría alguna orientación con este poco de teoría.
Código de ejemplo (solo para una de las 12 ecuaciones):
function [dx_dt] = statevector (t,x)
%This function contains the systems of DE describe the state vector
%State Vector x(t) = (Xg(t), Xl(t))
%x(t)=[V1(t),ω1(t),θ2(t),V2(t),ω2(t),θ3(t),V3(t),ω3(t)},θ4(t),V4(t),θ5(t),V5(t}) %x(1) = V1
%x(2) = ω1
%x(3) = θ2
%x(4) = V2
%...and so on
%Generators are busses 1 - 3, Load are busses 4 - 5
%Define Constants J = 1.26; % Inertia Constant D = 2; % Damping Coefficient Tdo = 0.25; % Time Constant Xd = 1.05; Xdp = 0.185; Pmi = 1; % Mechanical Input Power Ef = 1; % Electromagnetic Field
dx_dt = zeros(12,1); %initializing column vector which will contain eqns
%do i need to transpose dx_DT?
dx_dt(1) = 4 * Ef - 4 * x(1) + 3.46 * (2*x(4)*sin(x(3)) + 4*x(4)*cos(x(3)) + x(10)*sin(x(9)) + 4*x(10)* cos(x(4)));