How to find the points where two lines cross?

Question

Will Jeter 2020-11-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/644928-how-to-find-the-points-where-two-lines-cross

编辑： Matt J 2020-11-11

I need help finding the intersection of the two functions from this graph. I would like to have markers and have the command window tell me what the tmeprature value is at each point

T = [310:1:450];

v = 1;

V = 10;

dHrxn = -80000;

UA = 3600;

kf0 = 1;

kc0 = 100;

T0 = 400;

ER = 20000;

CpA = 40;

kf = exp(-20000*(1./T-1/400));

kc = 100*exp((80000/1.987)*(1./T-1/400));

x = 1./(v./(V.*kf)+1+(1./kc));

R = 400*(T-310);

G = x*-dHrxn;

figure

hold all

plot(T,R)

plot(T,G)

legend('R(T)','G(T)')

xlabel('Temperature (K)')

ylabel('G(T) and R(T) (cal/mol)')

title('G(T) and R(T) vs. T')

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Matt J 2020-11-11

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/644928-how-to-find-the-points-where-two-lines-cross#answer_541733

编辑：Matt J 2020-11-11

在 MATLAB Online 中打开

T1=fzero(@GRdiff,[350,400])
T1 = 377.5263
T2=fzero(@GRdiff,[400,450])
T2 = 418.3883
function out=GRdiff(T)
    v = 1;
    V = 10;
    dHrxn = -80000;
    UA = 3600;
    kf0 = 1;
    kc0 = 100;
    T0 = 400;
    ER = 20000;
    CpA = 40;
    kf = exp(-20000*(1./T-1/400));
    kc = 100*exp((80000/1.987)*(1./T-1/400));
    x = 1./(v./(V.*kf)+1+(1./kc));
    R = 400*(T-310);
    G = x*-dHrxn;
    out=R-G;
end