Commit a7194b80 authored by Angioni's avatar Angioni

minor updates

parent f967d1a4
......@@ -123,6 +123,7 @@ if strcmp(GridData.type_of_model,'single_phase')==1
std_err_Iphase_mat = std_err_polar_mat_IRI(3:2:2*GridData.Lines_num+1);
std_err_Ireal_mat = GridData.base_current*std_err_rectangular_mat_IRI(2:2:2*GridData.Lines_num+1);
std_err_Iimag_mat = GridData.base_current*std_err_rectangular_mat_IRI(3:2:2*GridData.Lines_num+1);
std_err_Iri = GridData.base_current*std_err_rectangular_mat_IRI;
else
ref = 11*GridData.DM.NGF+10*GridData.DM.NGS+2*GridData.DM.Nload;
......
......@@ -26,7 +26,7 @@ function [ rotI ] = calc_rotI( GridData,PowerData,use_pu)
rm_column = 1;
if use_pu == 0
PowerDataVmagn = PowerData.Vmagn * GridData.base_voltage;
PowerDataImagn = PowerData.Imagn* (GridData.base_power / GridData.base_voltage);
PowerDataImagn = PowerData.Imagn* GridData.base_current;
else
PowerDataVmagn = PowerData.Vmagn ;
PowerDataImagn = PowerData.Imagn ;
......
......@@ -100,8 +100,8 @@ if strcmp(GridData.type_of_model,'single_phase') == 1
x = GridData.topology(3,b);
Vnode(x,1) = Vnode(x,1) + 1;
Volts(x,1) = Volts(GridData.topology(2,b),1) - ( GridData.R1(1,b)+1i*GridData.X1(1,b))*Amps(b,1);
Vmagn_status(:,x) = abs(Volts(:,x));
Vph_status (:,x) = (phase(Volts(:,x)));
Vmagn_status(x,1) = abs(Volts(x,1));
Vph_status(x,1) = (phase(Volts(x,1)));
end
end
end
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment