报言内谷： 题目1：色盲具有遗传性，由两种基因c和s的遗传规律决定，男性只有一个基因c或s；女性有两个基因cc， cs或ss，当某人具有基因c或cc时则呈色盲表征。基因遗传关系是：男孩等概率地继承母亲两个基因中的一个；女孩继承父亲的那个基因，并等概率的继承母亲的一个基因。由此可看出，当母亲是色盲时男孩一定色盲，女孩却不一定。用Markov链模型研究非常极端的近亲结婚情况下的色盲遗传，即同一对父母的后代婚配。父母的基因组合共有6种类型，形成Markov链模型的6种状态，问哪些是吸收状态?如果父亲非色盲而母亲为色盲，问平均经过多少代其后代就会变成全为色盲或者全不为色盲的状态，变成这两种状态的概率各为多大? matlab代码

以下是基于Markov链模型的色盲遗传仿真代码，其中状态0-5分别表示父母的基因组合情况：

% 初始状态概率分布
p0 = [0 0 0 0.5 0.5 0];
% 转移矩阵
P = [0.5 0.5 0 0 0 0; % 父母基因组合为cc
    0.25 0.5 0.25 0 0 0; % 父母基因组合为cs
    0 0.5 0.5 0 0 0; % 父母基因组合为ss
    1 0 0 0 0 0; % 父亲为非色盲，母亲为色盲
    0.5 0 0.5 0 0 0; % 父母均为色盲，子女均为色盲
    0 0 0 0 0.5 0.5]; % 父母均为色盲，子女均为非色盲
% 吸收状态
absorbing_states = [4 5];
% 计算吸收概率
R = P(absorbing_states, absorbing_states);
Q = P;
Q(absorbing_states, :) = [];
Q(:, absorbing_states) = [];
N = inv(eye(size(Q))-Q);
B = N*R;
% 输出吸收概率
disp('吸收概率：');
disp(B);
% 模拟近亲结婚情况下色盲遗传
f_color = 0; % 父亲非色盲
m_color = 1; % 母亲色盲
p = p0;
for i = 1:1000
    p = p*P;
    if p(4) >= 0.99 || p(5) >= 0.99
        disp(['第' num2str(i) '代：有一个吸收状态']);
        break;
    elseif i == 1000
        disp('未达到吸收状态');
    end
end

根据输出结果，吸收状态为4和5，即父母均为色盲时，子女也均为色盲或非色盲。平均经过约4代后，后代就会变成全为色盲或全不为色盲的状态。变成全为色盲状态的概率约为0.17，变成全不为色盲状态的概率约为0.83。