sigma = 120; % 水平定位的标准差 sigma_z = 40; % 深度定位的标准差 depth_0 = 150; % 潜艇深度定位值 min_depth = 120; % 潜艇深度最小值 kill_radius = 20; % 深弹的杀伤半径 sub_length = 100; % 潜艇长度 sub_width = 20; % 潜艇宽度 % 潜艇中心位置的三维分布 [x, y, z] = ndgrid(-kill_radius:1:kill_radius, -kill_radius:1:kill_radius, min_depth:1:depth_0); % 计算命中概率 hit_prob = normcdf(kill_radius, x, sigma) .* normcdf(kill_radius, y, sigma) .* ... normpdf(z, depth_0, sigma_z) ./ (1 - normcdf(min_depth - depth_0, 0, sigma_z)); % 集成命中概率 total_hit_prob = sum(hit_prob(:)) / numel(hit_prob); % 找到最大命中概率对应的引爆深度 [~, idx] = max(total_hit_prob); optimal_depth = z(idx); % 假设z(idx)是正确的索引 fprintf('最大命中概率为：%f\n', total_hit_prob); fprintf('最佳引爆深度为：%d\n', optimal_depth);模型的建立和求解

1周前