clc; clear; close all; % 参数定义 L = 100; % 潜艇长度 (m) W = 20; % 潜艇宽度 (m) H = 25; % 潜艇高度 (m) d_s = 150; % 潜艇上表面深度 (m) Rk = 20; % 深弹杀伤半径 (m) sigma = 120; % 水平定位误差标准差 (m) N = 9; % 深弹数量 depth_fuze = 160; % 定深引信引爆深度 % 投弹间隔设定 d_x = 60; % 水平方向间隔 (m) d_y = 60; % 垂直方向间隔 (m) % 投弹位置网格 [x_grid, y_grid] = meshgrid(-1:1, -1:1); x_positions = x_grid * d_x; % 3x3阵列的X坐标 y_positions = y_grid * d_y; % 3x3阵列的Y坐标 % 计算命中概率 P_total = 1; % 初始化总命中概率 for i = 1:N % 计算每一枚深弹到潜艇中心的水平距离 distance_x = x_positions(i); distance_y = y_positions(i); horizontal_distance = sqrt(distance_x^2 + distance_y^2); % 平面距离 % 计算深弹的杀伤范围 vertical_distance = abs(depth_fuze - d_s); % 深弹与潜艇上表面的垂直距离 total_distance = sqrt(horizontal_distance^2 + vertical_distance^2); % 总距离 % 判断是否在潜艇平面范围内 if abs(distance_x) <= L/2 && abs(distance_y) <= W/2 if depth_fuze >= d_s && depth_fuze <= (d_s + H) % 在潜艇的上下表面范围内 P_hit = 1; % 触发引信直接命中 elseif vertical_distance <= Rk % 在定深引信引爆且潜艇在杀伤半径内 P_hit = 1; else P_hit = exp(-horizontal_distance^2 / (2 * sigma^2)); % 计算命中概率 end elseif total_distance <= Rk % 如果深弹落点在潜艇外，但潜艇在杀伤半径内 P_hit = 1; else P_hit = exp(-horizontal_distance^2 / (2 * sigma^2)); % 正态分布的命中概率 end % 更新未命中的概率 P_total = P_total * (1 - P_hit); end % 计算总命中概率 P_total = 1 - P_total; % 输出结果 fprintf('投弹间隔 (dx = %.2f m, dy = %.2f m), 定深引信深度 %.2f m\n', d_x, d_y, depth_fuze); fprintf('至少一枚深弹命中的概率为: %.4f\n', P_total); % 可视化投弹方案和潜艇 figure; hold on; grid on; axis equal; % 绘制潜艇 rectangle('Position', [-L/2, -W/2, L, W], 'EdgeColor', 'b', 'LineWidth', 2); text(0, 0, 'Submarine', 'HorizontalAlignment', 'center'); % 绘制深弹落点 scatter(x_positions, y_positions, 100 ... , 'r', 'filled'); text(x_positions(:), y_positions(:), 'Depth Charge', 'HorizontalAlignment', 'center', 'VerticalAlignment', 'bottom'); xlabel('X (m)'); ylabel('Y (m)'); title('Depth Charge Drop Array and Submarine Position'); hold off; 修改代码

1周前