固体填料床中的非稳态传热模型，将流体视为一维牛顿塞流，不考虑固体填料床径向传热，可计算轴向温度分布以及时间变化下温度分布的matlab代码

固体填料床中的非稳态传热通常涉及到傅里叶定律，用于描述热传导过程。在这种简化的一维模型中，我们假设固体填料颗粒对流体的冷却作用沿轴向均匀，并忽略了径向的传热影响。这种简化情况下的数学模型常常基于欧拉-贝努利方程结合能量守恒定律。

对于MATLAB代码实现，这里提供一个大致的框架：

function [tempDist, timeSeries] = solid_bed_heat_transfer(Q_in, T_wall, L, k, dt, numSteps)
    % 定义常量和变量
    alpha = k / (rho*c); % 热扩散系数 (导热系数/k，密度/rho，比热容/c)
    T_bed = zeros(1, numSteps); % 温度数组
    T_bed(1) = T_wall; % 初始条件 - 管壁温度
    dx = L / (numSteps - 1); % 分段长度

    for i = 1:numSteps - 1
        % 一阶差分法
        T_bed(i+1) = T_bed(i) + alpha * Q_in * dt / dx - alpha * (T_bed(i+1) - T_bed(i)) * dt;
    end
    
    % 轴向温度分布
    tempDist = linspace(T_wall, T_bed(end), numSteps);
    
    % 时间序列
    timeSeries = linspace(0, numSteps*dt, numSteps);
end

在这个函数中，
- `Q_in` 是单位时间内通过填料床的热量流量；
- `dt` 是时间步长；
- `numSteps` 是模拟的时间步数。

请注意，这只是一个基础版本的代码，实际应用中可能需要考虑更多的因素，比如边界条件、动态传质等。此外，这个代码没有包含任何图形展示功能，如果需要可视化结果，可以添加相应的绘图代码。