急!Matlab应用题这是一道Matlab在化工中的应用题,麻烦哪位大侠给修改一下源程序到可以在7.10.0版里运行（在源程序里纠错）.在下感激不尽! 1、问题描述（管式反应器的计算的温度分布）已知
来源：学生作业帮助网 编辑：六六作业网 时间：2024/11/23 00:10:59
急!Matlab应用题这是一道Matlab在化工中的应用题,麻烦哪位大侠给修改一下源程序到可以在7.10.0版里运行（在源程序里纠错）.在下感激不尽! 1、问题描述（管式反应器的计算的温度分布）已知
急!Matlab应用题
这是一道Matlab在化工中的应用题,麻烦哪位大侠给修改一下源程序到可以在7.10.0版里运行（在源程序里纠错）.在下感激不尽!
1、问题描述（管式反应器的计算的温度分布）
已知在管式反应器内进行液相反应A＝R＋S.该反应为吸热反应,所需热量由管外油浴供给,油温为340℃.实验确定反应温度与转化率的关系为：
其中,反应速率常数 , ,反应器壁温 若反应器入口温度 ,入口转比率 ,要求反应器出口转化率为 .试确定不同转化率下反应器的温度.
源程序：
function TFR
% 管式反应器的温度分布
clear
clc
global r k0 E tc
r=1.987;
k0=1.17e17;
E=44500;
t0=340;
tc=340;
xA0=0;
% dx=1/100;
% 调用ode45求解微分方程
[x,t]=ode45(@XTEqs,[0 0.999],[t0]);
% 通过spline函数插值获得0~1之间的函数值,间隔为0.1
xx = 0:0.1:1;
tt = spline(x,t,xx);
N=length(xx);
disp(sprintf('\n转化率 反应温度（K） '));
for i=1:N;
disp(sprintf('\n %5.1f,%8.2f%',xx(i),tt(i)));
end
% 对函数曲线和离散点画图
plot(xx,tt,'o',x,t,'b-');
xlabel('Fractional Conversion');
ylabel('Temperature (K)');
title('x~T');
grid;
% 定义微分方程
function dydt=XTEqs(x,t)
global r k0 E tc
k=k0*exp(-E/r/(t+273.15));
dydt=-65.0-15.58*(t-tc)/k/(1-x);
计算结果
执行程序：
在MATLAB的COMMAND窗口输出命令：>> TFR
急!Matlab应用题这是一道Matlab在化工中的应用题,麻烦哪位大侠给修改一下源程序到可以在7.10.0版里运行（在源程序里纠错）.在下感激不尽! 1、问题描述（管式反应器的计算的温度分布）已知
MATLAB是矩阵实验室（Matrix Laboratory）的简称,是美国MathWorks公司出品的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境,主要包括MATLAB和Simulink两大部分.