样本量为2的极小样本相容性检验方法

在航空航天领域由于成本、时间周期等原因进行疲劳寿命及可靠性评估时样本量通常极少（m=1或2）,利用相容性检验方法可对样本量进行扩充。常规的Wilcoxon秩和检验和K-S（Kolmogorov-Smirnov）检验适用于小样本情形,而极小样本相容性检验方面研究较少,且缺乏对方法合理性的详细说明和对不同方法检验功效优劣的比较。航空航天产品疲劳寿命多服从正态分布,因此本文主要以正态分布作为研究对象。利用Monte Carlo仿真发现从某一正态分布N（μ,σ²）中随机抽取两个样本x₁、x₂计算均值μ₁和标准差σ₁后构建新正态分布N（μ₁,σ₁²）,其±σ₁、±2σ₁和±3σ₁范围内的点落在原正态分布N（μ,σ²）±3σ范围内的概率依次为99.80%、98.13%和97.37%。在此基础上针对现场试验数据样本量为2的情况,本文提出利用3σ原则对先验信息数据进行相容性检验从而扩充样本量的方法。将该方法与两种文献方法对比后发现其误差率明显更低并呈现出检验性能随先验数据增加而不断提高的优势。     

A comparison of two adaptive detection schemes

Two schemes for adaptive detection are compared: Kelly's generalized likelihood ratio test (GLRT) and the mean level adaptive detector (MLAD). Detection performance, P_D, is predicted for the two schemes under the assumptions that the input noises are zero-mean complex Gaussian random variables that are temporally independent but spatially correlated; and the amplitude of the desired signal is Rayleigh distributed. P_D is computed as a function of the false alarm probability, the number of input channels, the number of independent samples per channel, and the matched filtered output signal-to-noise (S/N) power ratio. In this analysis the GLRT is shown to have better detection performance than the MLAD. The difference in detection performance increases as one uses fewer input samples. However, the required number of samples necessary to have only a 3 dB detection loss for both detection schemes is approximately the same. This is significant since for the present, the MLAD is considerably less complex to implement than the GLRT     

Spatially distributed target detection in non-Gaussian clutter

Two detection schemes for the detection of a spatially distributed, Doppler-shifted target in non-Gaussian clutter are developed. The non-Gaussian clutter is modeled as a spherically invariant random vector (SIRV) distribution. For the first detector, called the non-scatterer density dependent generalized likelihood ratio test (NSDD-GLRT), the detector takes the form of a sum of logarithms of identical functions of data from each individual range cell. It is shown under the clutter only hypothesis, that the detection statistic has the chi-square distribution so that the detector threshold is easily calculated for a given probability of false alarm P_F. The detection probability P_D is shown to be only a function of the signal-to-clutter power ratio (S/C)_opt of the matched filter, the number of pulses N, the number of target range resolution cells J, the spikiness of the clutter determined by a parameter of an assumed underlying mixing distribution, and P_F. For representative examples, it is shown that as N, J, or the clutter spikiness increases, detection performance improves. A second detector is developed which incorporates a priori knowledge of the spatial scatterer density. This detector is called the scatterer density dependent GLRT (SDD-GLRT) and is shown for a representative case to improve significantly the detection performance of a sparsely distributed target relative to the performance of the NSDD-GLRT and to be robust for a moderate mismatch of the expected number of scatterers. For both the NSDD-GLRT and SDD-GLRT, the detectors have the constant false-alarm rate (CFAR) property that P_F is independent of the underlying mixing distribution of the clutter, the clutter covariance matrix, and the steering vector of the desired signal     

Vector subspace detection in compound-Gaussian clutter. Part I: survey and new results

Deals with the problem of detecting subspace random signals against correlated non-Gaussian clutter exploiting different degrees of knowledge on target and clutter statistical characteristics. The clutter process is modeled by the compound-Gaussian distribution. In the first part of the paper, the optimum Neyman-Pearson (NP) detector, the generalized likelihood ratio test (GLRT), and a constant false-alarm rate (CFAR) detector are sequentially derived both for the Gaussian and the compound-Gaussian scenarios. Different interpretations of the various detectors are provided to highlight the relationships and the differences among them. In particular, we show how the GLRT detector may be recast into an estimator-correlator form and into another form, namely a generalized whitening-matched filter (GWMF), which is the GLRT detector against Gaussian disturbance, compared with a data-dependent threshold. In the second part of this paper, the proposed detectors are tested against both simulated data and measured high resolution sea clutter data to investigate the dependence of their performance on the various clutter and signal parameters.     

A CFAR adaptive matched filter detector

An adaptive algorithm for radar target detection using an antenna array is proposed. The detector is derived in a manner similar to that of the generalized likelihood-ratio test (GLRT) but contains a simplified test statistic that is a limiting case of the GLRT detector. This simplified detector is analyzed for performance to signals on boresight, as well as when the signal direction is misaligned with the look direction     

Polarization diversity detection in compound-Gaussian clutter

We present the problem of polarization diversity detection in compound-Gaussian clutter with unknown covariance matrix. To this end we assume that a set of secondary data, free of signal components and with the same covariance structure of the cell under test, is available. Due to the lack of a uniformly most powerful (UMP) detector we resort to a design procedure based upon the Rao and the Wald tests. Specifically we first derive the Rao and the Wald tests assuming that the covariance matrix is known, and then we substitute into the derived decision rules a suitable estimate of the clutter covariance. Interestingly, the newly proposed detectors share the constant false alarm rate (CFAR) property with respect to the texture statistical characterization. Moreover simulation results have shown that the Wald test based detector ensures a performance level higher than the Rao test. We have also conducted a further performance analysis, in the presence of real clutter data and in comparison with the previously proposed generalized likelihood ratio test (GLRT) based receivers, which highlights that, in general, the Wald test receiver outperforms its counterparts. Finally, since the newly proposed decision rules as well as the previously designed GLRTs do not ensure the CFAR property with respect to the clutter covariance matrix, we have developed a sensitivity analysis on the probability of false alarm (P/sub fa/), based on simulated clutter with covariance matrix estimated from real radar data. The results have shown that (P/sub fa/) is only slightly affected by variations in the clutter correlation properties and hence the CFARness is substantially achieved.     

A noncoherent adaptive detection technique

An adaptive detection technique suitable for both stationary and nonstationary noise environments based upon a generalized likelihood ratio test (GLRT) formulation is presented. The detector, which is statistically equivalent to a special form of the Wilks's lambda test, noncoherently combines the information contained in a pulse train of arbitrary length for decision-making purposes. The probability density function of the test under the noise only hypothesis is shown to be central χ². Under the signal plus noise hypothesis, an exact statistical characterization of the test cannot be obtained, and, therefore, a Chernoff bound is derived. Results in terms of the probability of detection versus signal-to-noise ratio (SNR) obtained from Monte Carlo simulation, the Chernoff bound, and the optimal matched filter case are examined. The performance of the noncoherent detector is shown to be a function of the covariance matrix estimate and the number of data samples     

Multi Objective Robust Active Vibration Control for Flexure Jointed Struts of Stewart Platforms via H∞ and μ Synthesis

Active vibration control is needed for future space telescopes, space laser communication and other precision sensitive payloads which require ultra-quiet environments. A Stewart platform based hybrid isolator with 6 hybrid struts is the effective system for active/passive vibration isolation over 5-250 Hz band. Using an identification transfer matrix of the Stewart platform, the coupling analysis of six channels is provided. A dynamics model is derived, and the rigid mode is removed to keep the signal of pointing control. Multi objective robust H∞ and μ synthesis strategies, based on singular values and structured singular values respectively, are presented, which simultaneously satisfy the low frequency pointing and high frequency disturbance rejection requirements and take account of the model uncertainty, parametric uncertainty and sensor noise. Then, by performing robust stability test, it is shown that the two controllers are robust to the uncertainties, the robust stability margin of H, controller is less than that of μ controller, but the order of μ controller is higher than that of H, controller, so the balanced controller reduction is provided. Additionally, the μ controller is compared with a PI controller. The time domain simulation of the μ controller indicates that the two robust control strategies are effective for keeping the pointing command and isolating the harmonic and stochastic disturbances.     

Adaptive resolution of mainlobe and sidelobe detections using model order determination

This work presents the development and performance evaluation of a methodology for distinguishing between mainlobe and sidelobe detections that arise in adaptive radar systems operating in adverse environments. Various adaptive detection test statistics such as the adaptive matched filter (AMF), the generalized likelihood ratio test (GLRT), and adaptive coherence estimate (ACE), and combinations of these, have been previously analyzed with respect to their sidelobe rejection capabilities. In contrast to these methods which are based on detecting a single target with known direction and Doppler, the present method uses model order determination techniques applied to the AMF or GLRT data observed over the range of unknown angle and Doppler parameters. The determination of model order, i.e., the number of signals present in the data, is made by using least-squares model fit error residuals and applying the Akaike Information Criterion (AIC). Comprehensive computer simulation results are presented which demonstrate substantial improvement in sidelobe rejection performance and detections of multiple sources compared with previous methods.     

Vector subspace detection in compound-Gaussian clutter. Part II: performance analysis

For pt. I see ibid., vol. 38, no. 4, p. 1295 (2002). In this second part we deal with the problem of detecting subspace random signals against correlated non-Gaussian clutter modeled by the compound-Gaussian distribution. In the first part of the paper, we derived the optimum Neyman-Pearson (NP) detector, the generalized likelihood ratio test (GLRT), and a constant false-alarm rate (CFAR) detector; we also provided some interesting interpretations of them. In this second part, these detectors are tested against both simulated data and measured high resolution sea clutter data to investigate the dependence of their performance on the various clutter and signal parameters. Numerical examples concern a space-time adaptive processing (STAP) scenario and a ground-based surveillance radar system scenario.     

Covariance matrix estimation for adaptive CFAR detection incompound-Gaussian clutter

We address the estimation of the structure of the covariance matrix and its application to adaptive radar detection of coherent pulse trains in clutter-dominated disturbance modeled as a compound-Gaussian process. For estimation purposes we resort to range cells in spatial proximity with that under test and assume that these cells, free of signal components, can be clustered into groups of data with one and the same value of the texture. We prove that, plugging the proposed estimator of the structure of the covariance matrix into a previously derived detector, based upon the generalized likelihood ratio test (GLRT), leads to an adaptive detector which ensures the constant false alarm rate (CFAR) property with respect to the clutter covariance matrix as well as the statistics of the texture. Finally, we show that this adaptive receiver has an acceptable loss with respect to its nonadaptive counterpart in cases of relevant interest for radar applications     

Nonparametric Radar Extraction Using a Generalized Sign Test

A nonparametric procedure used in a constant false alarm rate (CFAR) radar extractor for detecting targets in a background of noise with unknown statistical properties is described. The detector is based on a generalization of the well-known two-sample sign test and thus requires a set of reference noise observations in addition to the set of observations being tested for signal presence. The detection performance against Gaussian noise is determined for a finite number of observations and asymptotically, for both nonfluctuating and pulse-to-pulse Rayleigh fluctuating target statistics. It is noted that the performance loss, as compared to the optimum parametric detector, depends critically on the number of reference noise observations available when the number of hits per target is not large. In the same case a much larger loss is also found for a pulse-to-pulse fluctuating target even though the asymptotic loss is the same as for a nonfluctuating target. A comparison is finally made with a detector based on the Mann-Whitney test, which usually is considered to be one of the better nonparametric procedures for the two-sample case.     

基于H∞算法的飞机机翼结冰气动参数辨识

针对国内大型飞机结冰防护需求,开展针对大型结冰研究样机的H_∞算法参数辨识结冰探测研究。首先通过参数调节选取一组合适的H_∞算法参数,利用考虑测量噪声的结冰研究样机飞行仿真数据验证H_∞算法的辨识能力,由结果对比发现辨识算法能够跟踪飞机气动导数随结冰累积过程的变化趋势,辨识精度较高,其最大归一化平方根（RMS）误差仅为真值的11%;分析了H_∞算法对81种不同结冰累积过程的辨识能力,通过结果分析发现结冰累积时间较长且结冰速度较慢的情况辨识效果较差,结冰累积时间在100~300 s之间辨识精度较高;最后利用蒙特卡罗仿真分析了不同测量噪声大小对H_∞算法辨识精度和跟踪延时的影响,给出了3个纵向气动导数在随机误差影响下的辨识误差和跟踪延时的统计结果,发现在给定噪声标准差变化范围内,升力和俯仰力矩关于迎角的导数能够得到较为准确的辨识结果,二者的归一化平方根误差均值仅为各自真值的1.8%和4%,其预报延时均值最大仅为3 s和9.5 s。     

GLRT subspace detection for range and Doppler distributed targets

A generalized likelihood ratio test (GLRT) is derived for adaptive detection of range and Doppler-distributed targets. The clutter is modeled as a spherically invariant random process (SIRP) and its texture component is range dependent (heterogeneous clutter). We suppose here that the speckle component covariance matrix is known or estimated thanks to a secondary data set. Thus, unknown parameters to be estimated are local texture values, the complex amplitudes and Doppler frequencies of all scattering centers. To do so, we use superresolution methods. The proposed detector assumes a priori knowledge on the spatial distribution of the target and has the precious property of having a constant false alarm rate (CFAR) with the assumption of a known speckle covariance matrix or by the use of frequency agility.     

Adaptive Radar Detection in Doubly Nonstationary Autoregressive Doppler Spread Clutter

The problem of adaptive radar detection in clutter which is nonstationary both in slow and fast time is addressed. Nonstationarity within a coherent processing interval (CPI) often precludes target detection because of the masking induced by Doppler spreading of the clutter. Across range bins (i.e., fast time), nonstationarity severely limits the amount of training data available to estimate the noise covariance matrix required for adaptive detection. Such difficult clutter conditions are not uncommon in complex multipath propagation conditions where path lengths can change abruptly in dynamic scenarios. To mitigate nonstationary Doppler spread clutter, an approximation to the generalized likelihood ratio test (GLRT) detector is presented wherein the CPI from the hypothesized target range is used for both clutter estimation and target detection. To overcome the lack of training data, a modified time-varying autoregressive (TVAR) model is assumed for the clutter return. In particular, maximum likelihood (ML) estimates of the TVAR parameters, computed from a single snapshot of data, are used in a GLRT for detecting stationary targets in possibly abruptly nonstationary clutter. The GLRT is compared with three alternative methods including a conceptually simpler ad hoc approach based on extrapolation of quasi-stationary data segments. Detection performance is assessed using simulated targets in both synthetically-generated and real radar clutter. Results suggest the proposed GLRT with TVAR clutter modeling can provide between 5–8 dB improvement in signal-to-clutter plus noise ratio (SCNR) when compared with the conventional methods.     

空间非合作航天器抓捕后姿态抗干扰控制

针对空间非合作航天器抓捕后存在未知不确定惯性参数的柔性组合体姿态稳定控制问题,基于中间状态观测器设计方法提出了一种新的姿态稳定抗干扰控制方法,同时考虑了诸多扰动及控制输入受限问题。研究结果表明,传统的姿态稳定控制方法需要已知柔性航天器惯性参数信息及状态信息,上述信息未知情况下会使姿态难以高精度稳定控制,且容易导致控制输入不满足受限要求。针对该问题,考虑控制输入幅值及变化率受限前提,提出了一种基于中间状态观测器的抗干扰控制方法,通过引入辅助变量构造新型中间状态观测器,同时估计组合体状态信息及综合干扰,设计出了一种新的组合体姿态稳定抗干扰控制器。通过Lyapunov稳定性分析方法证明了所设计的控制器能够保证闭环系统的全局渐近稳定性。相比于已有的混合H₂/H_∞控制器,所提出的抗干扰控制器在应用时不需要柔性组合体的姿态及模态信息,并且也不需要惯性参数的辨识过程。最后,通过给定参数进行仿真对比,进一步验证了所设计控制器的有效性和优越性。     

Theory of Adaptive Radar

This paper reviews the principles of adaptive radar in which both the spatial (antenna pattern) and temporal (Doppler filter) responses of the system are controlled adaptively. An adaptive system senses the angular-Doppler distribution of the external noise field and adjusts a set of radar parameters for maximum signal-to-interference ratio and optimum detection performance. A gradient technique for control of the radar array/filter weights is described and shown to generate weights which asymptotically approach optimum values. Simulation results illustrate the convergence rate of adaptive systems and the performance improvement which can be achieved.     

Joint TNT and RDX detection via quadrupole resonance

Quadrupole resonance (QR) technology has been receiving increasing attention for explosive detection in applications including landmine detection. It can be used as a confirming sensor for buried plastic landmine detection by detecting the explosives within the mine. Since a single mine can contain more than one type of explosives (e.g., trinitrotoluene (TNT) and (RDX) compound), a detector designed to detect only one type of explosive may not provide the best detection. We focus herein on the joint detection of TNT and RDX explosives for the landmine detection via the QR sensor. To mitigate the RF interferences (RFIs) associated with the QR measurements, we apply a doubly constrained robust Capon beamformer (DCRCB) to perform the spatial beamforming by exploiting the spatial correlation of the RFIs. Based on the output of the spatial DCRCB, we devise a generalized likelihood ratio test (GLRT) for the joint TNT/RDX detection, which is referred to as the DCRCB-GLRT detector. We show that our detector has the constant false alarm rate (CFAR) property and that the detection variable obeys a x² -mixture distribution in the mine-free scenario. The effectiveness of the proposed DCRCB-GLRT detector is demonstrated with the experimental data collected by Quantum Magnetics, Inc.     

异常值和未知观测噪声鲁棒的非线性滤波器

针对含异常观测值的非线性系统滤波问题,以Huber损失函数替代推导滤波器最大后验准则中观测误差的l₂范数,构造出了一种新的优化准则函数,从而给出了一种对异常值鲁棒的非线性后验线性化滤波器。分析表明：由于Huber损失函数兼具l₁和l₂范数的性质,从而使得由这个新准则推导出的滤波器,不仅具有l₂范数的低误差拟合性,也具备l₁范数对异常值的鲁棒性。而当观测噪声的分布未知时,通过引入箱线图法检测异常值,并对噪声统计分布的参数进行估计,进一步提出了对异常值和未知观测噪声分布鲁棒的非线性后验线性化滤波器。仿真实验验证了分析结果的有效性,并表明本文算法的性能优于现有文献报道的非线性滤波算法。     

三维C-H网格下跨音压气机转子Euler流场数值模拟

介绍了三维C-H网格的生成和全三维转子流场的数值解。以法国宇航院设计的一个高负荷跨音转子为算例, 并与他们的计算结果作了比较。应用这套数值模拟技术得以计算带进口预旋的某跨音压气机转子流场, 并对进口预旋对其下游的转子流场的影响作了数值模拟, 得出了若干有意义的结论。