가우스 자이델 법 매트랩 프로그래밍

소개

가우스 자이델법으로 연립방정식 푸는 매트랩 코드

목차

◉ 문제 11.29 회로도
◉ 프로그램 (Gauss-Seidel 이용)
◉ 실행 결과
◉ 결과 분석

본문내용

회로에 전류법칙을 적용하면,



를 구하라.


◉ 프로그램 (Gauss-Seidel 이용)

%% << step.1 - 변수입력---여기서 원하는 입력값 넣어줌>>%% %%------------------------

R1=10; R2=20; R3=25; R4=4; R5=25; R6=8; %% 문제에 설정된 저항변수
V=100; I2=10; %% 문제에 설정된 전압원,전류원 변수
N = 10; %% 실행횟수
Es=0.5; %% 허용오차


%% << step.2 - 행렬표현 & 초기행렬지정>> %%-----------------------------------------------
A = [(R1+R2+R3) 0 -R3 ; 0 (R4+R5+R6) -R4 ; -R3 -R4 R3+R4 ];%% 임의로 열거된행렬 (문제의나온식과 같이 열거)
a = [ 0 (R4+R5+R6) -R4 ; -R3 -R4 R3+R4 ; (R1+R2+R3) 0 -R3 ]; %% 입력변수를 바꾸면 나중에 피벗과 scaling & initializing 될 행렬
adummy = [ 0 (R4+R5+R6) -R4 ; -R3 -R4 R3+R4 ; (R1+R2+R3) 0 -R3 ] ; %% 피벗화 할때 이용할 행렬

a1 = [(-R2*I2) 0 -R3 ; -R5*I2 (R4+R5+R6) -R4 ; V -R4 R3+R4]; %%행렬식 만들기위한 행렬
a2 = [ (R1+R2+R3) (-R2*I2) -R3 ; 0 -R5*I2 -R4 ; -R3 V R3+R4];%%Cramer법으로 참값 구하기위한 행렬
a3 = [ (R1+R2+R3) 0 (-R2*I2) ; 0 (R4+R5+R6) -R5*I2 ; -R3 -R4 V];

참고자료

· 없음