动态规划实现rna二级结构预测的标记函数
时间: 2023-12-25 20:03:24 浏览: 124
rna二级构预测
RNA二级结构预测中,标记函数是指对于RNA序列中每个位置i和j,标记函数M(i,j)代表序列中i到j位置是否匹配。动态规划实现RNA二级结构预测的标记函数的过程如下:
1. 构建一个二维数组M,用于存储标记函数的结果。数组的大小为RNA序列的长度,即M[i][j]表示序列中i到j位置是否匹配。
2. 初始化数组M。对于任意的i=j,标记函数M(i,j)为0。
3. 根据以下递推式计算标记函数M(i,j):
M(i,j) = max{M(i+1,j), M(i,j-1), M(i+1,j-1) + P(i,j)}
其中,P(i,j)表示序列中i和j位置的配对得分,即i和j配对的概率。
4. 递推计算完M后,标记函数M(1,n)即为RNA序列的最优二级结构预测结果。
需要注意的是,标记函数M(i,j)的计算需要考虑到i和j位置是否配对的情况,因此需要结合RNA序列的配对得分进行计算。同时,为了使得标记函数的计算更加高效,可以采用动态规划的方法,将计算过程分解为子问题,避免重复计算。
阅读全文