修正文字显示错误 >> % 参数设置
mu = 100; % 原始电阻均值 (Ω)
sigma = 8; % 原始电阻标准差 (Ω)
N = 1e5; % 样本数
sigma2_R = sigma^2; % 原始电阻方差 = 64 Ω²
% 生成独立正态分布样本
R = normrnd(mu, sigma, [N, 4]);
% 计算各连接方式等效电阻
R1 = R(:,1) + R(:,2); % 串联
R2 = (R(:,1).*R(:,2)) ./ (R(:,1)+R(:,2)); % 并联
R3 = ((R(:,1)+R(:,2)) .* (R(:,3)+R(:,4))) ./ ((R(:,1)+R(:,2)) + (R(:,3)+R(:,4))); % 先串后并
R4 = (R(:,1).*R(:,2))./(R(:,1)+R(:,2)) + (R(:,3).*R(:,4))./(R(:,3)+R(:,4)); % 先并后串
% 绘制四个电阻的分布图(最终修正)
figure;
subplot(2,2,1);
histfit(R1, 50, 'normal');
title({...
'R1 (串联)', ...
sprintf('原始电阻: $\\mu=%.1f\\ \\Omega$, $\\sigma^2=%.1f\\ \\Omega^2$', mu, sigma2_R), ...
sprintf('理论均值 $=%.1f\\ \\Omega$, 仿真均值 $=%.1f\\ \\Omega$', 2*mu, mean(R1)), ...
sprintf('理论方差 $=%.1f\\ \\Omega^2$, 仿真方差 $=%.1f\\ \\Omega^2$', 2*sigma2_R, var(R1))...
}, 'Interpreter', 'latex', 'FontSize', 10);
xlabel('阻值 ($\\Omega$)', 'Interpreter', 'latex'); % 修正转义
ylabel('频数');
subplot(2,2,2);
histfit(R2, 50, 'normal');
title({...
'R2 (并联)', ...
sprintf('原始电阻: $\\mu=%.1f\\ \\Omega$, $\\sigma^2=%.1f\\ \\Omega^2$', mu, sigma2_R), ...
sprintf('理论均值 $\\approx%.1f\\ \\Omega$, 仿真均值 $=%.1f\\ \\Omega$', mu/2, mean(R2)), ...
sprintf('理论方差 $\\approx%.3f\\ \\Omega^2$, 仿真方差 $=%.3f\\ \\Omega^2$', sigma2_R/8, var(R2))...
}, 'Interpreter', 'latex', 'FontSize', 10);
xlabel('阻值 ($\\Omega$)', 'Interpreter', 'latex'); % 修正转义
ylabel('频数');
subplot(2,2,3);
histfit(R3, 50, 'normal');
title({...
'R3 (先串后并)', ...
sprintf('原始电阻: $\\mu=%.1f\\ \\Omega$, $\\sigma^2=%.1f\\ \\Omega^2$', mu, sigma2_R), ...
sprintf('理论均值 $=%.1f\\ \\Omega$, 仿真均值 $=%.1f\\ \\Omega$', mu, mean(R3)), ...
sprintf('理论方差 $\\approx%.1f\\ \\Omega^2$, 仿真方差 $=%.1f\\ \\Omega^2$', sigma2_R/4, var(R3))...
}, 'Interpreter', 'latex', 'FontSize', 10);
xlabel('阻值 ($\\Omega$)', 'Interpreter', 'latex'); % 修正转义
ylabel('频数');
subplot(2,2,4);
histfit(R4, 50, 'normal');
title({...
'R4 (先并后串)', ...
sprintf('