//?30節點系統潮流計算程序.cpp?:?Defines?the?entry?point?for?the?console?application.
//
#include?“stdafx.h“

//*************************************************************************************************************//
//IEEE-30節點系統潮流計算程序

#include
#include

#define?N?30??????????????//節點數
#define?n_PQ?24????????????//PQ節點數
#define?n_PV?5????????????//PV節點數
#define?n_br?41???????????//串聯支路數
#define?n_C?2?????????????//裝有并聯電容節點數
//#define?n_L?1

#define?PI?3.1415926
double?K1K2K3K4;??????????//變壓器變比
double?Yc1Yc2;
//double?Xmcr;

void?main（）
{
	//void?disp_matrix（double?*disp_pint?disp_mint?disp_n）;???????				//矩陣顯示函數

	//節點電壓賦初值，PV節點電壓幅值已知，相角置0；平衡節點電壓幅值和相角均已知；PQ節點電壓幅值設1，相角設0.
	//				??1????????2????????3????????4????????5????????6????????7????????8????????9????????10?
	/*double?Us[2*N]={1.000001.000001.000001.000001.000001.000001.000001.000001.000001.00000
					1.000001.000001.000001.000001.000001.000001.000001.000001.000001.00000
					1.000001.000001.000001.000001.033801.005801.023001.091301.088301.05000
					};?*/
	/*double?Us[2*N]={1.000001.000001.000001.000001.000001.000001.000001.000001.000001.00000
					1.000001.000001.000001.000001.000001.000001.000001.000001.000001.00000
					1.000001.000001.000001.000001.050001.050001.050001.050001.050001.05000
					};*/
	double?Us[2*N]={1.000001.000001.000001.000001.000001.000001.000001.000001.000001.00000
					1.000001.000001.000001.000001.000001.000001.000001.000001.000001.00000
					1.000001.000001.000001.000001.000001.000001.000001.000001.000001.05000
					};
????
	//PV節點發電機有功賦初值（除平衡節點外）
	//			????1??????2??????3??????4??????5??????6??????7??????8??????9??????10?
	double?Ps[N]={0.00000.00000.00000.00000.00000.00000.00000.00000.00000.0000
				??0.00000.00000.00000.00000.00000.00000.00000.00000.00000.0000
				??0.00000.00000.00000.00000.57560.24560.35000.17930.1691
				};?
	//發電機無功初值置0
	double?Qs[N]={0};???????????????????????????????????????????????????????????????????????????????????????????

	//各節點有功負荷賦值
	//			????1?????2?????3?????4?????5?????6?????7?????8?????9?????10?
	double?PL[N]={0.0240.0760.0000.2280.0000.0580.1120.0620.0820.035
				??0.0900.0320.0950.0220.1750.0000.0320.0870.0000.035
				??0.0000.0000.0240.1060.2170.9420.3000.0000.000
				};??
	
	//各節點無功負荷賦值
	//			????1?????2?????3?????4?????5?????6?????7?????8?????9?????10?
	double?QL[N]={0.0120.0160.0000.1090.0000.0200.0750.0160.0250.018
				??0.0580.0090.0340.0070.1120.0000.0160.0670.0000.023
				??0.0000.0000.0090.0190.1270.1900.3000.0000.000
				};????

	//double?K1=1.0155;K2=0.9629;K3=1.0129;K4=0.9581;??????????//變壓器當前變比
	//double?Yc1=0.19;Yc2=0.04;
	double?K1=1;K2=1;K3=1;K4=1;??????????//變壓器當前變