基于Padé近似的考虑多级衰变的放射性核素活度计算方法

本发明提供了一种基于Padé近似的考虑多级衰变的放射性核素活度计算方法，其包括：S1、构建包含所有核素的完整衰变矩阵；步骤S1包括：S11、接收衰变链数据结构作为输入，构建n×n的方阵A；S12、将整个矩阵的所有元素初始化为零值；S13、将第i行第i列的元素设置为第i个核素衰变常数的负值，即A(i,i)=；S14、遍历所有可能的核素对(i,j)；S2、采用矩阵范数驱动的自适应缩放算法，结合Padé近似矩阵计算矩阵指数；S3、根据计算模式执行相应的放射性核素活度计算。本发明通过构建完整的衰变矩阵，考虑所有核素的贡献，避免了截断误差。采用自适应缩放算法结合Padé近似矩阵求解实现刚性系统矩阵计算，有效避免数值溢出和舍入误差，保证计算过程的稳定性。