MATLAB拟合Aa，aa，AA遗传比例，证明AA，aa，Aa在遗传过程中比例不变，写出程序代码

MATLAB是一种强大的数值计算工具，常用于数据处理和统计分析。要通过它来模拟和验证AA、aa和Aa三种基因型的比例是否在遗传过程中保持不变，我们可以创建一个简单的概率模型，并利用迭代法或蒙特卡洛模拟来进行实验。首先，假设我们有一个已知的初始比例，比如AA占25%，aa占50%，Aa占25%。

% 定义初始遗传比例
initial_ratio = [0.25, 0.5, 0.25]; % AA, aa, Aa

% 设定迭代次数和随机种子以保证结果复现
num_iterations = 1000;
rand('seed', 1); % 设置随机数种子

% 模拟遗传过程
for iteration = 1:num_iterations
    % 随机选择亲代基因
    parent1 = randi([1 3], 1, 2); % 父亲的两个基因
    parent2 = randi([1 3], 1, 2); % 母亲的两个基因
    
    % 计算后代的基因组合并更新比例
    offspring = [parent1(1), parent2(1)]; % 子女从父母各继承一个
    if parent1(2) == parent2(2)
        offspring(end+1) = parent1(2); % 如果双亲同型，则子女也为该基因型
    else
        offspring(end+1) = randi([1 3]); % 否则，随机获得第三个基因
    end
    
    % 更新比例，注意加上新出生个体
    initial_ratio = (initial_ratio * (iteration - 1)) + (ones(size(initial_ratio)) / 4);
    
    % 判断当前比例是否稳定
    is_stable = norm(diff(initial_ratio, 1), 'fro') < 1e-6; % 比较连续两次的比例变化
    if is_stable
        break; % 如果比例趋于稳定，提前结束循环
    end
end

% 输出最后的遗传比例
disp("经过迭代后的遗传比例: ");
disp(initial_ratio);

% 相关问题: