%%%%%%%%%%%%%%% 벡터와 행렬 문제 %%%%%%%%%%%%%%%%%%%%%%
% 김재헌, 오규남, 이주현, 이창진
A = [3 0 2; 0 3 -2; 2 -2 1] %행렬 A의 입력
B = [21; -12; 17] %행렬 B의 입력
I_A = inv(A) %행렬 A의 역행렬
D_A = det(A)
%inv(B) %행렬 B의 역행렬
x = inv(A)*B %Ax=B꼴에서 x를 구함
E_A = eig(A) %행렬 A의 고유치
[X, D] = eig(A) %고유벡터, 대각화
T_A = A' %전치행렬
A^7 %행렬 A의 7제곱
DIAG_A = diag(A) %행렬 A의 대각값을 구함
diag(diag(A)) %행렬 A의 대각행렬을 구함
jordan(A) %행렬 A의 jordan형을 구함
A = [2 1; 3 4]
poly(A) % det(A-λI)=0의 계수를 내림차순으로 .... 결과는 λ^2 - 6λ + 5 = 0
pause
% 황승희
A = [4 -2 5]
B = [3 1 -1]
Dot_Vec = dot(A,B) % 내적
Cross_Vec = cross(A,B) % 외적
%%%%%%%%%%%%%%%%%%%%%%% 선형 미분 방정식의 풀이 %%%%%%%%%%
% 박세환
y = dsolve('Dy = sin5x', 'y(1)=3') % y' = sin5x
ezplot(y)
pause
y = dsolve('Dy-y = sin(t)', 'y(0) = 10') % y' - y = sint, y(0) = 10
ezplot(y)
pause
dsolve('t^2*D2y-6*y = 0') % t^2 * y(t)'' - 6*y = 0일 때의 y(t)를 구함
pause
% 연립 미분 방정식-이창진
% [y1, y2] = dsolve('수식1', '수식2')
% y1(t)' - y1 + 2*y2 = 4*cos(t) - 2*sin(t)
% y2(t)' -3*y1 + 4*y2 = 5*cos(t) - 5*sin(t)
% y1(0)=1, y2(0)=2
[y1, y2] = dsolve('Dy1 - y1 + 2*y2 = 4*cos(t) - 2*sin(t)', 'Dy2 -3*y1 + 4*y2 = 5*cos(t) - 5*sin(t)', 'y1(0)=1, y2(0)=2')
% 황승희
y = int('sin(x)') % y' = sin(x) 구하기, 적분을 취함
ezplot(y)
pause
% 이성열, 김재헌, 정영락
X = []
X0 = [2; -1; 4]
A = [0 -6 -1; 6 2 -16; -5 20 -10]
for t=0:0.01:1
X = [X expm(t*A)*X0];
end
plot3(X(1,:), X(2,:), X(3,:), 'o')
grid on
pause
% 미분
sym x
diff x^3 %미분 - 오규남
syms x s % 기호 변수를 정의
y = 5*x^2 % 기호식 정의
diff(y,x) % y를 x에 관하여 미분, - 이주현, 김현주
int(y) % x에 대해 적분
int(y,x) % y를 x에 관하여 적분한다.
% 적분 - 오규남
sym y
int y^3
% 김재헌, 오규남, 이주현, 이창진
A = [3 0 2; 0 3 -2; 2 -2 1] %행렬 A의 입력
B = [21; -12; 17] %행렬 B의 입력
I_A = inv(A) %행렬 A의 역행렬
D_A = det(A)
%inv(B) %행렬 B의 역행렬
x = inv(A)*B %Ax=B꼴에서 x를 구함
E_A = eig(A) %행렬 A의 고유치
[X, D] = eig(A) %고유벡터, 대각화
T_A = A' %전치행렬
A^7 %행렬 A의 7제곱
DIAG_A = diag(A) %행렬 A의 대각값을 구함
diag(diag(A)) %행렬 A의 대각행렬을 구함
jordan(A) %행렬 A의 jordan형을 구함
A = [2 1; 3 4]
poly(A) % det(A-λI)=0의 계수를 내림차순으로 .... 결과는 λ^2 - 6λ + 5 = 0
pause
% 황승희
A = [4 -2 5]
B = [3 1 -1]
Dot_Vec = dot(A,B) % 내적
Cross_Vec = cross(A,B) % 외적
%%%%%%%%%%%%%%%%%%%%%%% 선형 미분 방정식의 풀이 %%%%%%%%%%
% 박세환
y = dsolve('Dy = sin5x', 'y(1)=3') % y' = sin5x
ezplot(y)
pause
y = dsolve('Dy-y = sin(t)', 'y(0) = 10') % y' - y = sint, y(0) = 10
ezplot(y)
pause
dsolve('t^2*D2y-6*y = 0') % t^2 * y(t)'' - 6*y = 0일 때의 y(t)를 구함
pause
% 연립 미분 방정식-이창진
% [y1, y2] = dsolve('수식1', '수식2')
% y1(t)' - y1 + 2*y2 = 4*cos(t) - 2*sin(t)
% y2(t)' -3*y1 + 4*y2 = 5*cos(t) - 5*sin(t)
% y1(0)=1, y2(0)=2
[y1, y2] = dsolve('Dy1 - y1 + 2*y2 = 4*cos(t) - 2*sin(t)', 'Dy2 -3*y1 + 4*y2 = 5*cos(t) - 5*sin(t)', 'y1(0)=1, y2(0)=2')
% 황승희
y = int('sin(x)') % y' = sin(x) 구하기, 적분을 취함
ezplot(y)
pause
% 이성열, 김재헌, 정영락
X = []
X0 = [2; -1; 4]
A = [0 -6 -1; 6 2 -16; -5 20 -10]
for t=0:0.01:1
X = [X expm(t*A)*X0];
end
plot3(X(1,:), X(2,:), X(3,:), 'o')
grid on
pause
% 미분
sym x
diff x^3 %미분 - 오규남
syms x s % 기호 변수를 정의
y = 5*x^2 % 기호식 정의
diff(y,x) % y를 x에 관하여 미분, - 이주현, 김현주
int(y) % x에 대해 적분
int(y,x) % y를 x에 관하여 적분한다.
% 적분 - 오규남
sym y
int y^3
