아래는 c / c + + 프로그램과 소스 코드를 계산합니다 ybus 및 전기 zbus 네트워크의 소요 기간은 어떤 크기합니다.의 입력 시스템이 적혀 "rx.txt"를 다음과 같은 열 형식; '버스', '이 버스를 ','연구 ', 그리고 점령한'합니다. '보낸 사람 버스'와 '버스'연속적인 양의 정수 여야합니다 제로부터 시작합니다. 버스를 0으로 수락 또는 스윙 버스로 추측합니다. 광고는 절차가 복잡 헤더 파일을 사용 매트릭스 작업을 지원합니다. 계산은 간단하게 zbus 반전시킬 ybus합니다.

//-----------------------------
# pragma hdrstop
//----------------------------
# pragma argsused
# 포함 <stdio.h>
# 포함 <conio.h>
# 포함 "matrix.h"
# 포함 <complex.h>
를 사용하여 표준 : : 복잡한;
네임 스페이스를 사용하여 수학;
형식 복잡 복잡;
형식 매트릭스 매트릭스;
무효 메인 (무효)
(
/ / 선언할 변수를
int cnt, 전, 일본, matsize = 1;
이중 tmp3, tmp4;
복잡한 tmp6;
매트릭스 그리고, ybus, zbus;
파일을 * 안에;
그 setsize (matsize, matsize); / / 임시 매트릭스 크기를 설정하는
/ / 가져와 분기 데이터를
이 = fopen ( "rx.txt", "연구"); / / 열려있는 파일
만일 (== 널)
(printf ( "\ nrx.txt를 찾을 수없습니다");
getch ();
)
cnt = 일본 = 나는 = 0;
을 (;;)
(경우 (fscanf ( "% 패", & 전) == 못하게 파일이 끝났습니다)
(휴식 시간;)
만약 (나는> matsize)
(matsize = 나는;)
fscanf ( "% 패", &日);
만일 (日> matsize)
(matsize = 일본;)
그 setsize (matsize +1, matsize +1);
fscanf ( "% 만약", & tmp3);
fscanf ( "% 만약", & tmp4);
tmp6 = 컴플렉스 (tmp3, tmp4);
그리고 (전, 일본) = 1.0/tmp6;
cnt + +;
)
fclose (인치);
/ / 끝에 도착 지점 데이터를
/ / 만들기 ybus
ybus.setsize (matsize +1, matsize +1);
/ / 대각선 항목을
을 (나는 = 0; 나는 <= matsize; 나는 + +)
(에 대한 (일본 = 0; 일본 <= matsize; 일본 + +)
(ybus (난) = ybus (난) + 그리고 (전, 일본);)
)
을 (日= 0; 일본 <= matsize; 일본 + +)
(위해 (나는 = 0; 나는 <= matsize; 나는 + +)
(ybus (일본, 일본) = ybus (日,日) + 그리고 (전, 일본);)
)
/ / 끝
/ / 오프 대각선 항목을
cnt = 0;
에 대한 (일본 = cnt; 일본 <= matsize; 일본 + +)
(위해 (나는 = 0; 나는 <= matsize; 나는 + +)
(만일 (日! = 나는)
(ybus (전, 일본) = ybus (日, 나) =- 1.0 * 그리고 (전, 일본);)
)
cnt + +;
)
/ / 끝
/ / 엔드를 생성 및 버스
/ / 버스를 0이 매트릭스에 포함되어 있지 않은 ybus
/ / 때문에 버스에 대한 참조
에 대한 (일본 = 1; 일본 <= matsize; 일본 + +)
(위해 (나는 = 1; 나는 <= matsize; 나는 + +)
(ybus (- 1, 일본 - 1) = ybus (전, 일본);
)
)
ybus.setsize (matsize, matsize);
zbus.setsize (matsize, matsize);
printf ( "\ n \ nybus 매트릭스 \ n");
법원 <<ybus <<endl;
zbus =! ybus; / / ybus의 반전
printf ( "\ n \ nzbus 매트릭스 \ n");
법원 <<zbus <<endl;
getch ();
)