비선형방정식의 해 값을 간결하게 표현하려면 어떻게 해야하나요?

Question

0 个投票

안녕하세요

Symbolic Toolbox 를 통해 비선형방정식의 해를 구하려고 하는데요,

예를들어, 아래와 같이 코드를 작성한 후 계산하면 해 값으로 1/2 - (3355~~ 식으로 출력이 됩니다.

사실 이 값을 command에 입력하면 0.1558 로 간단하게 출력이 됩니다.

이런식으로 해를 간단한 형태로 출력되게 하려면 어떻게 해야할까요?

syms f y_A0 y_B0 y_C0 P T
for y_A0 = 0.5:0.1:0.9
y_B0 = 1-y_A0;
y_C0 = 0;
P = 2;
T = 423;
equation1 = (y_C0+f*y_A0)*(1-2*f*y_A0)^2/(y_A0-f*y_A0)/(y_B0-2*f*y_A0)^2/P^2 - exp(11454/T-28.36);
equation2 = y_A0-f*y_A0 > 0 ;
equation3 = y_B0-2*f*y_A0 > 0 ;
equ = [equation1 equation2 equation3 ];
solution = solve(equ,f)
end
solution =
 
1/2 - (33554432*9502337923034183^(1/2))/9502337923034183
 
 
solution =
 
5/9 - (2*3^(1/2)*9502337923034183^(1/2)*12035612713430087^(1/2)*sin(pi/3 - atan((528482304*28507013769102549^(1/2)*1917227268151034126266240444162^(1/2))/36083747825255125136499363251279)/3))/85521041307307647 - (2*9502337923034183^(1/2)*12035612713430087^(1/2)*cos(pi/3 - atan((528482304*28507013769102549^(1/2)*1917227268151034126266240444162^(1/2))/36083747825255125136499363251279)/3))/85521041307307647
 
 ...