Matlab 編程繪制二元相圖

R=8.314472;

dHb=11000;

dHc=20080;

Tb=546;

Tc=903;

k=R*T;

t=T-273.15;

dGb=dHb*(1-T/Tb);

dGc=dHc*(1-T/Tc);

%%plot(dGb,t,dGc,t,k,T)

xbl=(1-exp(dGc./k))./(exp(dGb./k)-exp(dGc./k));

xbs=exp(dGb./k).*(1-exp(dGc./k))./(exp(dGb./k)-exp(dGc./k));

%%plot(dGb,t,dGc,t)

plot(xbl,T,'b-');

hold on;

plot(xbs,T,'r-');

xlabel('B物質濃度'),ylabel('T');

legend('液相濃度xbl','固相濃度xbs');