风电入网潮流计算14节点-VrPQ.m

风电入网潮流计算14节点-VrPQ.m
关于风电入网的matlab程序：风电节点无功对电压偏导求修正newton拉夫逊直角迭代时，对电压进行分解u=e j*f,把修正的风节点电压迭代量代入后潮流不收敛呀？是不是偶的风节点迭代加入错了？
if j1==14
                          %X3=2*P^2*X_deta^2*r2* j*f)- j*f)^3)/sqrt j*f)^4-4*P^2*X_deta^2))/ j*f)^2-sqrt j*f)^4-4*P^2*X_deta^2))^2)- j*f)- j*f)^3/sqrt j*f)^4-4*P^2*X_deta^2))*/;
                           %X4=j*X3;
                       %end;
%%%%%%%%%%节点为风电节点。