建模与仿真VVA

2020/12/6飞行器工程系单家元博士15建模与仿真VV&AVV&A基本概念模型校核模型验证仿真结果的统计分析实例2020/12/6飞行器工程系单家元博士25.1VV&A基本概念VV&A基本概念校核Verification：证实模型从一种形式转换成另一种形式的过程具有足够的精确度。验证Validation：从预期的角度来确定模型表达实际系统的准确程度，考查模型在其作用域内是否准确地代表了实际系统。包含两层意义：模型是否正确描述系统，模型输出是否接近实际系统行为。确认Accreditation：模型相对特定的研究目的来说是可接受的，相信并接受这一模型的权威性。2020/12/6飞行器工程系单家元博士35.1VV&A基本概念仿真精度与置信度仿真精度：仿真系统实现其规定或期望的动、静态性能技术指标的误差或允许误差。分静态、动态。总精度是分系统精度的综合。仿真置信度：表示仿真系统在特定的建模和仿真意义上，在总体结构和行为上能够复现原型系统的可信性程度。2020/12/6飞行器工程系单家元博士45.1VV&A基本概念影响仿真系统置信度的因素建模的原理和方法建模过程忽略了一些次要因素初始和原始数据的选取上的遗漏或随机性计算机求解过程的误差2020/12/6飞行器工程系单家元博士55.1VV&A基本概念仿真置信度等级1级：结果、概括和结论，鉴定标准、仿真目标、系统组成和功能。2级：子模型、理论模型，比较数据和比较方法3级：硬件分系统，数据收集及其鉴定方法的描述、复杂模型简图4级：硬件在回路仿真5级：实际系统置信度评估理论模型有效性确认、仿真软件校核与验证、运行有效性确认、数据有效性确认、内部安全性确认2020/12/6飞行器工程系单家元博士65.1VV&A基本概念VVA工作模式建模对象目标和要求概念模型（假设、算法、结构等）校验与验证计划技术说明软件开发实现问题求解特殊应用（虚拟、设备、人、仿真器）数据输入特特特特特特特特特特特特特特特特特特特特特特特特2020/12/6飞行器工程系单家元博士75.1VV&A基本概念VVA工作过程导弹制导系统概念模型仿真大纲半实物仿真系统正式模型数学仿真模型数学仿真结果物理效应模型物理仿真结果半实物仿真结果仿真大纲校验数学模型校验物理效应模型校验仿真模型校验数学仿真模型验证物理仿真模型验证半实物仿真系统验证仿真确认实际飞行结果半实物仿真系统校验概念模型校验特特特特特特特特特特特特特特特特特特特特特特特特特特特特特特特特特特2020/12/6飞行器工程系单家元博士85.1VV&A基本概念模型校核与验证的难点模型验证工作是一个过程从理论到实践，再到理论的多次重复的过程。模型验证工作具有模糊性相似程度具有一定的模糊或不确定性。影响模型验证工作的因素的多样性内因：结构：模型的共性；参数：模型的个性外因：激励输入模型验证工作的复杂性大量、复杂的计算实际系统的不可试验性2020/12/6飞行器工程系单家元博士95.2模型校核5.2.1系统仿真误差源校核仿真模型校核：数学模型建立时作了多种假设和简化，忽略了一些次要的、不可观察的因素，而仿真模型则是在此基础上的二次简化。仿真算法校核：编制仿真程序时由所采用的算法本身引起的截断误差和舌入误差。仿真硬件误差：数学仿真：计算机字长。半实物：物理设备、对象实物、接口、时间匹配。2020/12/6飞行器工程系单家元博士105.2模型校核5.2.2程序校验：软件过程中的软件校核方法模型程序校核：零、部件模型校核弹体动力学模型校核控制回路模型校核制导回路模型校核时序控制、设备驱动程序校核程序工作时序是否与武器系统一致？各仿真设备是否同步？驱动程序的坐标转换、比例因子是否准确？2020/12/6飞行器工程系单家元博士115.2模型校核5.2.2程序校验：软件过程中的软件校核方法仿真程序校核方法非正规方法：概念执行、代码会审、设计审查静态分析：词法、语义、结构、数据流动态测试：黑/白箱测试、自上至下/自下至上、跟踪调试符号分析：符号运行、路径分析、因果图示约束分析：边界分析、断言检查理论证明：归纳、演绎、推理2020/12/6飞行器工程系单家元博士125.3模型验证5.3.1概述模型验证：通过模型运行结果与实际系统运行结果的偏差来表示模型的准确度。模型可信性验证：考察演绎过程的可信性前提、假设是否真实：建模是在一定的应用目的和试验条件下。推理过程是否符合思维规律、规则。模型一致性验证方法：考察归纳过程的可信性。主观确认法：专家经验评估法数理统计法：一定的样本数量，样本序列的独立性时－频分析法：比较模型与实际系统的输出时间序列、输出频谱动态关联法：根据先验知识，提出一关联性能指标，利用该指标对模型与实际系统输出进行比较，作出定性判定。定性和定量分析：2020/12/6飞行器工程系单家元博士135.3模型验证主观确认法动态关联分析法数理统计法时－频分析法参数估计法参数假设检验非参数假设检验直观有效评价法TIC不等式系数法点估计t-检验符号检验时间序列分析事件有效评价法灰色关联分析区间估计F-检验秩和检验频谱分析预测有效性确认法相关系数分析最小二乘估计x2检验游程检验小波分析图灵测试法相似系数分析极大似然估计Bayes方法序贯检验回归分析Bayes估计2020/12/6飞行器工程系单家元博士145.3.2置信区间法原理：对于一个可观测的多元响应仿真模型，它在一组给定的试验框架下的精度可以用模型和实际系统输出变量的均值差的置信区间表示。置信区间计算方法一元法：应用一元统计方法和Bonferrioni不等式构造。用t检验法检验均值。多元法：应用多元统计方法直接计算。用Roy-Bose法检验均值。数据要求：满足独立同分布、正态分布、同方差假设。的均值输出变量为矢量jjsjmjsmsm,,2020/12/6飞行器工程系单家元博士155.3.2置信区间法基本步骤确定ARA数据获取自驱动模型：模型输入数据要从与实际系统同分布的随机函数中独立抽样。模型与实际系统输出独立且同分布。协驱动模型：直接采用实际系统输入数据。模型与实际系统输出同分布但不独立。TkjsjmjjUUUUTLkLLLkjUL],...,,1[,],...,,2,[,...,2,1,212020/12/6飞行器工程系单家元博士165.3.2置信区间法选择统计方法：根据观测样本是否满足独立同分布、正态分布、同方差假设条件，选择一元或多元统计方法。统计计算：根据所选统计方法，计算模型和实际系统输出变量的置信区间确认模型的有效性：如果每个变量置信度为(1-r)的置信区间[l,u]在集合[L,U]中，则认为模型有效。2020/12/6飞行器工程系单家元博士175.3.2置信区间法样本假设k,...,2,1j),1n/())yy((S)S...SS(S)S,SNsS,S:nm(k,...,2,1j,n/)y(y)y...yy(Y)y,y,Y:Nsy,y,Y:nm(2n1ijji2njT2k222122Nj2s2nj2mn1ijijTk21sjNsjsmjnmjm：下标，样本数量为实际系统的带下标，样本数量为模型的带样本方差向量：上标，样本数量为实际系统的带上标，样本数量为模型的带样本均值向量：2020/12/6飞行器工程系单家元博士185.3.2置信区间法一元法假设模型和实际系统的输出数据是独立分布的。如果不知道模型和系统的输出方差，但可以认为是相等的，则置信度为(1-r)的同时置信区间为：式8.13如果不知道模型和系统的输出方差，但可以认为是不相等的，则置信度为(1-r)的同时置信区间为：式8.14如果模型和实际系统的输出结果是同分布的（样本为N），但相互不独立。则置信度为(1-r)的同时置信区间为：式8.16如果模型和实际系统的输出结果不是同分布的（样本为N），但相互不独立。则置信度为(1-r)的同时置信区间为：式8.172020/12/6飞行器工程系单家元博士195.3.2置信区间法多元法：如果模型和系统的输出均服从多元正态分布当模型和实际系统的输出数据是相互独立的时，置信度为(1-r)的同时置信区间为：式8.18当模型和实际系统的输出数据相互不独立时，则可利用模型和系统变量观测值的配对来构造同时置信区间，为：式8.21如果已知模型和实际系统变量的置信区间，在模型和系统输出数据相互独立和不独立的情况下，最低置信度水平为(1-r)，则同时置信区间为：式8.26的百分位。分布上的的是自由度为2/)2(,...,2,1,11)(2,2/2,2/jNnrjNnrsjmjrtNntkjNnStyyjj2020/12/6飞行器工程系单家元博士205.3.2置信区间法一元法与多元法的比较：一元法要求模型（系统）观测数量的各分量之间相互独立，即忽略变量之间可能存在地相关性。这种假设与实际情况一般不相符。多元法同时对多个变量的观测数据进行分析，通过样本相关系数和同时置信区间反映输出变量之间的关系。多元法应用条件比较严格，要求影响变量服从多元正态分布。一元法相对宽松。一元法允许模型（系统）的各个输出变量具有不同的置信度水平，样本量也可以不一样。多元法不允许，需系统的置信度水平和样本量。一元法不能得到精确的置信度，也可能得到无意义的结果。2020/12/6飞行器工程系单家元博士215.3.2置信区间法统计分析方法的选择统计分析方法的假设条件各变量的精度即区间的置信度水平是否相同模型和系统的样本量是否相同模型精度即置信度区间的宽度是否合适置信度是否合适2020/12/6飞行器工程系单家元博士225.3.3假设检验法假设检验法包括：t检验法，F检验法，x2检验法，K-S检验，符号检验，游程检验，秩检验等。假设检验法适用：大样本、静态数据。如：脱靶量。原理：使用统计推理，根据仿真数据和实际系统的数据中的信息来研究试验的典型统计量，决定是接受还是拒绝仿真数据与实际输出数据的接近程度的假设。2020/12/6飞行器工程系单家元博士235.3.3假设检验法使用步骤建立初始假设模型确认：H0表示在一定的试验条件和目的下，模型在可接受的精度范围内有效；H1表示在一定的试验条件和目的下，模型在可接受的精度范围内无效。概率Pa1.0α*β*λ*λ2020/12/6飞行器工程系单家元博士245.3.3假设检验法模型确认时会犯两类错误：I：当H0为真时，而样本观测值落入拒绝域，因此拒绝H0，为模型建立者风险；用表示其概率。II：当H1为真时，而样本观测值未落入拒绝域，因此接受H0，为模型用户风险;用表示其概率。显然：＋＝1模型有效性参数：模型精度的有效范围决定了的范围,0*。样本长度：样本量越大，和均越小；小样本情况下相对较小。2020/12/6飞行器工程系单家元博士255.3.3假设检验法使用步骤产生检验信号合理选择施加的测试信号反映仿真误差、降低初始条件对测试信号的影响、尽可能简单选择条件量和假设交换提取合适充分反映建模误差的统计检验参数：均值、方差等将H1H0的假设条件进行转换：转换为检验参数的判断假设检验设计用一种统计算法来估计统计参量（随机量），依据检验规则作出接受还是拒绝假设的判断。2020/12/6飞行器工程系单家元博士265.3.4动态关联分析法判断比较法原理：通过专家对模型和数据系统的输出数据进行比较，定性判断两者是否接近。特性：定性判断。参数数据可以是静态的，也可以是动态的。2020/12/6飞行器工程系单家元博士275.3.4动态关联分析法标量指标法分类：THEIL不等式系数法、灰色关联法原理：通过模型和实际系统输出两个序列的误差计算，可以给出某一性能指标来度量这两个时间序列一