并网风电场运行静态电压稳定性分析中PQ仿真模型创建与分析，给出matlab实现代码

创建PQ仿真模型：

% 定义风电场节点
Vwind = 1; % 风电场电压幅值，pu
theta_wind = 0; % 风电场电压相角，rad
Vwind_complex = Vwind*exp(1i*theta_wind); % 风电场电压复数形式
Pwind = 10; % 风电场有功出力，MW
Qwind = -5; % 风电场无功出力，MVAr

% 定义负荷节点
Vload = 1; % 负荷电压幅值，pu
theta_load = 0; % 负荷电压相角，rad
Vload_complex = Vload*exp(1i*theta_load); % 负荷电压复数形式
Pload = 20; % 负荷有功功率，MW
Qload = 10; % 负荷无功功率，MVAr

% 定义母线节点
Vbus = 1; % 母线电压幅值，pu
theta_bus = 0; % 母线电压相角，rad
Vbus_complex = Vbus*exp(1i*theta_bus); % 母线电压复数形式
Pbus = Pwind - Pload; % 母线有功功率，MW
Qbus = Qwind - Qload; % 母线无功功率，MVAr

% 创建PQ节点矩阵
Y = [1/Vwind_complex -1i*Qwind/Vwind_complex 0;
     -1i*Qwind/Vwind_complex 1/Vwind_complex 0;
     0 0 0];
S = [Pwind-Pload; Qwind-Qload; 0];
V = [Vwind_complex; Vload_complex; Vbus_complex];

% 计算电流
I = Y*V;

静态电压稳定性分析：

% 定义潮流潮流限制
Vmin = 0.95; % 最小电压幅值，pu
Vmax = 1.05; % 最大电压幅值，pu

% 进行潮流计算
iter = 0; % 迭代次数
max_iter = 100; % 最大迭代次数
tol = 0.01; % 收敛容差
while iter < max_iter
    % 更新电压幅值
    V(1) = max(min(abs(V(1)), Vmax), Vmin)*exp(1i*angle(V(1)));
    V(2) = max(min(abs(V(2)), Vmax), Vmin)*exp(1i*angle(V(2)));
    
    % 计算潮流误差
    S_cal = V.*conj(I)/1e6;
    S_error = S - S_cal;
    
    % 判断是否收敛
    if abs(S_error(1)) < tol && abs(S_error(2)) < tol && abs(S_error(3)) < tol
        break;
    end
    
    % 更新迭代次数
    iter = iter + 1;
    
    % 更新电压相角
    theta = angle(V);
    theta(1) = theta(1) + atan2(imag(I(1)), real(I(1)));
    theta(2) = theta(2) + atan2(imag(I(2)), real(I(2)));
    V(1) = abs(V(1))*exp(1i*theta(1));
    V(2) = abs(V(2))*exp(1i*theta(2)));
end