I am not sure I understand the issue here. Any help would be appreciated
%% Initial Variables
s1=0.27; s2=1.04; s3=0.716; s4=0; s5=0; Vref=1.04; Pref=0.716;
%% Define Diferential Equations
function dsdt= pv_system(t,s)
%parameters
Ki=20; Kp=4; Kip=10; Kiq=10; Tr=0.02; TGpv=0.01; Rq=0; Rp=0.05;
Vref=1.04; Qactual=0.27; Vt=1.04;Pref=0.716;DeltaOmega=0;Pactual=0.716;
Qcmd =0.27;
%State Variables
s1=s(1); s2=s(2); s3=(s); s4=s(4); s5=s(4);
%diferential Equations
ds1dt=Ki*(Vref-s2-Rq*Qactual);
ds2dt = (1 / Tr) * (Vt – s2);
ds3dt = (1 / TGpv) * (Pref – (DeltaOmega / Rp) – s3);
ds4dt = Kip * (s3 – Pactual);
ds5dt = Kiq * (Qcmd – Qactual);
dsdt = [ds1dt; ds2dt; ds3dt; ds4dt; ds5dt];
end
%% inital Vector
s0 = [s1;s2;s3;s4;s5];
%% Time span for simulation
tspan = [0 10]; % 0 to 10 seconds
%% Solve using ode45
[t, s] = ode45(@pv_system, tspan, s0);I am not sure I understand the issue here. Any help would be appreciated
%% Initial Variables
s1=0.27; s2=1.04; s3=0.716; s4=0; s5=0; Vref=1.04; Pref=0.716;
%% Define Diferential Equations
function dsdt= pv_system(t,s)
%parameters
Ki=20; Kp=4; Kip=10; Kiq=10; Tr=0.02; TGpv=0.01; Rq=0; Rp=0.05;
Vref=1.04; Qactual=0.27; Vt=1.04;Pref=0.716;DeltaOmega=0;Pactual=0.716;
Qcmd =0.27;
%State Variables
s1=s(1); s2=s(2); s3=(s); s4=s(4); s5=s(4);
%diferential Equations
ds1dt=Ki*(Vref-s2-Rq*Qactual);
ds2dt = (1 / Tr) * (Vt – s2);
ds3dt = (1 / TGpv) * (Pref – (DeltaOmega / Rp) – s3);
ds4dt = Kip * (s3 – Pactual);
ds5dt = Kiq * (Qcmd – Qactual);
dsdt = [ds1dt; ds2dt; ds3dt; ds4dt; ds5dt];
end
%% inital Vector
s0 = [s1;s2;s3;s4;s5];
%% Time span for simulation
tspan = [0 10]; % 0 to 10 seconds
%% Solve using ode45
[t, s] = ode45(@pv_system, tspan, s0); I am not sure I understand the issue here. Any help would be appreciated
%% Initial Variables
s1=0.27; s2=1.04; s3=0.716; s4=0; s5=0; Vref=1.04; Pref=0.716;
%% Define Diferential Equations
function dsdt= pv_system(t,s)
%parameters
Ki=20; Kp=4; Kip=10; Kiq=10; Tr=0.02; TGpv=0.01; Rq=0; Rp=0.05;
Vref=1.04; Qactual=0.27; Vt=1.04;Pref=0.716;DeltaOmega=0;Pactual=0.716;
Qcmd =0.27;
%State Variables
s1=s(1); s2=s(2); s3=(s); s4=s(4); s5=s(4);
%diferential Equations
ds1dt=Ki*(Vref-s2-Rq*Qactual);
ds2dt = (1 / Tr) * (Vt – s2);
ds3dt = (1 / TGpv) * (Pref – (DeltaOmega / Rp) – s3);
ds4dt = Kip * (s3 – Pactual);
ds5dt = Kiq * (Qcmd – Qactual);
dsdt = [ds1dt; ds2dt; ds3dt; ds4dt; ds5dt];
end
%% inital Vector
s0 = [s1;s2;s3;s4;s5];
%% Time span for simulation
tspan = [0 10]; % 0 to 10 seconds
%% Solve using ode45
[t, s] = ode45(@pv_system, tspan, s0); ode45 MATLAB Answers — New Questions

​