How to correct the regenerated plot of given analytical solution of ODE?
2 次查看(过去 30 天)
显示 更早的评论
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/600145/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/600150/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/600155/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/600160/image.png)
The analytical solution for this ODE is given belown and written as MATLAB in code routine:
for ![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/600150/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/600150/image.png)
t = -log(1-a)/((K1+K2)*(B-1)) + log(1-(a/B))/((K1+(K2*B))*(B-1)) + (K2*log(1+((a*K2)/K1)))/((K1+(K2*B))*(K1+K2));
and for ![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/600160/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/600160/image.png)
t = -(1/K3)*(log((1-a)/0.7))+tc;
where,
tc = -log(0.7)/((K1+K2)*(B-1)) + log(1-(0.3/B))/((K1+(K2*B))*(B-1)) + (K2*log(1+((0.3*K2)/K1)))/((K1+(K2*B))*(K1+K2));
other constants used are:
R = 8.314;
T = 450;
B = 0.4:0.54;
A1 = 2.101e9;
A2 = -2.014e9;
A3 = 1.960e9;
delE1 = 8.07e4;
delE2 = 778e4;
delE3 = 5.66e4;
K1 = A1*exp(-delE1/(R*T));
K2 = A2*exp(-delE2/(R*T));
K3 = A3*exp(-delE1/(R*T));
k = K1/K2;
I regenerated the plot for
from the the given solution:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/600175/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/600180/image.png)
But the actual plot of the given solution is
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/600185/image.png)
where I'm doing WRONG?
Please guide me.
Thanks
0 个评论
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!