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
评论
    相关推荐