```matlab
clc;
clear;

num_simulations = 1e5;
threshold_male = 2250;
threshold_female = 1800;

inputs = struct();
inputs.male.knee_angle_mean = 36.8;
inputs.male.knee_angle_std = 15.7;
inputs.male.shank_angle_mean = -11.6;
inputs.male.shank_angle_std = 8.9;
inputs.male.posterior_grf_mean = 0.6;
inputs.male.posterior_grf_std = 0.29;
inputs.male.varus_valgus_mean = -0.003;
inputs.male.varus_valgus_std = 0.017;
inputs.male.internal_external_mean = -0.016;
inputs.male.internal_external_std = 0.003;
inputs.male.hamstring_force_mean = 190.5;
inputs.male.hamstring_force_std = 118.3;
inputs.male.gastrocnemius_force_mean = 59.0;
inputs.male.gastrocnemius_force_std = 57.7;

inputs.female.knee_angle_mean = 34.5;
inputs.female.knee_angle_std = 20.4;
inputs.female.shank_angle_mean = -7.2;
inputs.female.shank_angle_std = 13.2;
inputs.female.posterior_grf_mean = 0.57;
inputs.female.posterior_grf_std = 0.22;
inputs.female.varus_valgus_mean = -0.006;
inputs.female.varus_valgus_std = 0.016;
inputs.female.internal_external_mean = 0.0002;
inputs.female.internal_external_std = 0.003;
inputs.female.hamstring_force_mean = 156.3;
inputs.female.hamstring_force_std = 107.9;
inputs.female.gastrocnemius_force_mean = 63.1;
inputs.female.gastrocnemius_force_std = 48.8;

male_injury_count = 0;
female_injury_count = 0;

for i = 1:num_simulations
    male_knee_angle = normrnd(inputs.male.knee_angle_mean, inputs.male.knee_angle_std);
    male_shank_angle = normrnd(inputs.male.shank_angle_mean, inputs.male.shank_angle_std);
    male_posterior_grf = normrnd(inputs.male.posterior_grf_mean, inputs.male.posterior_grf_std);
    male_varus_valgus = normrnd(inputs.male.varus_valgus_mean, inputs.male.varus_valgus_std);
    male_internal_external = normrnd(inputs.male.internal_external_mean, inputs.male.internal_external_std);
    male_hamstring_force = normrnd(inputs.male.hamstring_force_mean, inputs.male.hamstring_force_std);
    male_gastrocnemius_force = normrnd(inputs.male.gastrocnemius_force_mean, inputs.male.gastrocnemius_force_std);
    
    male_fk_ap = male_posterior_grf * cosd(male_shank_angle);
    male_facl_s = male_fk_ap / cosd(male_knee_angle);
    male_facl_vv = male_varus_valgus * (4.179 - 0.1108 * male_knee_angle + 9.43e-4 * male_knee_angle^2);
    male_facl_ie = male_internal_external * (11.815 - 0.7035 * male_knee_angle + 9.47e-3 * male_knee_angle^2);
    male_facl_total = male_facl_s + male_facl_vv + male_facl_ie;
    
    if male_facl_total > threshold_male
        male_injury_count = male_injury_count + 1;
    end

    female_knee_angle = normrnd(inputs.female.knee_angle_mean, inputs.female.knee_angle_std);
    female_shank_angle = normrnd(inputs.female.shank_angle_mean, inputs.female.shank_angle_std);
    female_posterior_grf = normrnd(inputs.female.posterior_grf_mean, inputs.female.posterior_grf_std);
    female_varus_valgus = normrnd(inputs.female.varus_valgus_mean, inputs.female.varus_valgus_std);
    female_internal_external = normrnd(inputs.female.internal_external_mean, inputs.female.internal_external_std);
    female_hamstring_force = normrnd(inputs.female.hamstring_force_mean, inputs.female.hamstring_force_std);
    female_gastrocnemius_force = normrnd(inputs.female.gastrocnemius_force_mean, inputs.female.gastrocnemius_force_std);
    
    female_fk_ap = female_posterior_grf * cosd(female_shank_angle);
    female_facl_s = female_fk_ap / cosd(female_knee_angle);
    female_facl_vv = female_varus_valgus * (4.179 - 0.1108 * female_knee_angle + 9.43e-4 * female_knee_angle^2);
    female_facl_ie = female_internal_external * (11.815 - 0.7035 * female_knee_angle + 9.47e-3 * female_knee_angle^2);
    female_facl_total = female_facl_s + female_facl_vv + female_facl_ie;
    
    if female_facl_total > threshold_female
        female_injury_count = female_injury_count + 1;
    end
end

male_injury_prob = (male_injury_count / num_simulations) * 100;
female_injury_prob = (female_injury_count / num_simulations) * 100;

fprintf('男性ACL损伤概率：%.4f%%\n', male_injury_prob);
fprintf('女性ACL损伤概率：%.4f%%\n', female_injury_prob);
```