# GJB 6289-2008 代入型CEP误差计算方法代码

• JOHNCLOCK
了解作者
• matlab
开发工具
• 625B
文件大小
• zip
文件格式
• 0
收藏次数
• 3 积分
下载积分
• 0
下载次数
• 2022-10-07 20:22
上传日期
GJB 6289-2008 地地弹道式导弹命中精度评定方法中的代入型CEP误差计算方法代码，可以用于命中精度的评定
CEP_calc.m.zip
• CEP_calc.m
921B

clear; xi = [55.85, 136.95, 92.41, 59.15, 64.08, 88.09, 60.57, -16.19, 21.31]; zi = [-97.67, 57.60, 132.20, 2.63, -31.03, 207.75, -205.99, 3.04, -13.22]; N = numel(xi); mu_x_hat = 1/N*sum(xi); mu_z_hat = 1/N*sum(zi); sig_x_hat = sqrt(1/(N-1)*sum((xi-mu_x_hat).^2)); sig_z_hat = sqrt(1/(N-1)*sum((zi-mu_z_hat).^2)); p_hat = sum((xi-mu_x_hat).*(zi-mu_z_hat))/sqrt(sum((xi-mu_x_hat).^2).*sum((zi-mu_z_hat).^2)); mu1 = mu_x_hat; mu2 = mu_z_hat; sig1 = sig_x_hat; sig2 = sig_z_hat; func = @(x,y) exp(-1/2 * ((x-mu1).^2./sig1.^2 + (y-mu2).^2./sig2^2)); polarfunc = @(theta, r) func(cos(theta).*r, sin(theta).*r).*r; % rmax = @(theta) 1./(sin(theta) + cos(theta)); Rhat = 116.97; ry = @(x) sqrt(Rhat.^2-x.^2); rym = @(x) -sqrt(Rhat.^2-x.^2); tic CEP_a = 1/(2*pi*sig1*sig2)*integral2(func, -Rhat, Rhat, rym, ry); toc tic CEP_b = 1/(2*pi*sig1*sig2)*integral2(polarfunc, 0, 2*pi, 0, Rhat); toc

相关推荐