MATLAB 코딩
ppt의 예시코딩
코딩 실습결과
위에서 왼쪽 예시코딩을 보면 A=inv(L)*R; B=inv(L)*D; U=V.*D;있고 for문 안의 식마다 B*U인데 B=inv(L)*U이므로 U를 U=V.*D;로 만들어 준 다음에 B=inv(L)*U;로 하고, for문 안의 식마다 있는 U를 지워줘야 올바로 된 값을 구할 수 있었다.
출력된 plot
ppt 예시코딩의 plot
위의
레포트 manual에 나와 있는 정보를 사용하였다. 실험 매뉴얼의 Fin 제원 정보를 통해서도 확인 가능하다.
온도 실험에 사용한 Fin은 Pure Copper로 이루어진 Thin Rectangular Fin으로써 2차원 형상을 가지고 있다. 이번 실험에서는 Fin의 가열 과정에서 Steady State 상태에 도달한 이후에 측정을 실시하므로 실험 전
1. Plot the 1-D temperature profile with analytical solution
(temperature vs fin length)
1) Analytical Solution
Fin의 미소면적에 대해 대류에 의한 열의 손실을 고려하여 열전달식을 세워 보면,
---------- (1)
한편 미소변화량은 다음과 같이 나타낼 수 있다.
---------- (2)
(1), (2)에서
---------- (3)
Four
(4) FDM으로 얻은 data의 Temperature profile
clear all
close all
T_fin = 28.8;
T_b = 34.8;
T_inf = 19;
k = 401;
h = 3.4752;
dx = 0.002;
d = 0.002;
A = eye(3750);
B = zeros(3750,1);
C = zeros(3750,1);
T = zeros(150,50);
for i=1:1:25
B(i,1) = T_fin;
end
for i=3726:1:3750
B(i,1) = T_b;
end
for i=26:1:3725
A(i,i-25) = d;
A(i,i-1) = d;
A(i,i) = -4*d
Matlab 으로 표를 그리는 명령어를 배웠다.
Bar(x,y,’k’) 지정하지 않은 좌표의 위치를 비워 두고 나타낸다. 하지만, Bar(y)는 지정한 좌표는 고려 하지 않고 순서대로 막대 그래프로 나타낸다.
Bar 이외에도 수평방향의 그래프를 위한 barh, 히스토그램을 위한 hist, 이산데이터 그래프를 위한stem,등등 다