열역학 1-75.docx



열역학 Thermodynamics 5th Edition.

Fundamentals and Applications

-Yunus A. Cengel

-Michael A. Boles

-부준홍 김덕줄 김세웅 김수현 신세현 이교우 정우남 최경민 공역

McGraw-Hill



문제 1-75


공학용 소프트웨어를 이용하여 방정식의 해를 찾는다.

풀이: 가정 먼저 EES Limited Academic Version을 이용하여 방정식의 해를 구하며, 코드 내용은 아래와 같다.
(2*x^3)-(10*x^0.5)-(3*x)=-3


EES의 Calculate>Slove를 선택하거나 F2키, 또는 아래의 아이콘을 선택하여 방정식을 푼다.


아래와 같이 오류가 있을 경우


Options의 Variable Info를 선택하거나 F9키 또는 를 바로 눌러 변수 정보 창을 연다.


그 후 변수 x의 초기 추측값을 2로 수정하여 Calculate>Slove 또는 F2키, Slove 아이콘을 선택하여 방정식을 푼다.


계산된 해는 아래와 같다.


매트랩(Matlab)의 fzero() 함수를 이용하여 주어진 방정식의 해를 구한다.
fzero() 함수를 사용하기 전에 ezplot() 함수를 이용하여 방정식의 해가 정수 2 근처에 있음을 알 수 있다.
ezplot('(2*x^3)-(10*x^0.5)-(3*x)+3'), grid on

다음 fzero() 함수를 이용하여 방정식의 해를 구한다.
fzero('(2*x^3)-(10*x^0.5)-(3*x)+3',2)


아래와 같이 초기값 선택이 좋지 않으면 해를 제대로 계산되지 않을 수 있다.


매트랩(matlab)의 slove()를 이용하여 주어진 방정식의 양의 실근은 구한다.
syms s
fzero('(2*x^3)-(10*x^0.5)-(3*x)+3',[0.1 3])


slove() 함수를 이용하면 방정식의 모든 해를 계산할 수 있다.


Posted by 귀여운촌아
: