CN105068931B - Single-particle soft error reliability calculation method for analyzing DSP software system - Google Patents

Single-particle soft error reliability calculation method for analyzing DSP software system Download PDF

Info

Publication number
CN105068931B
CN105068931B CN201510518615.8A CN201510518615A CN105068931B CN 105068931 B CN105068931 B CN 105068931B CN 201510518615 A CN201510518615 A CN 201510518615A CN 105068931 B CN105068931 B CN 105068931B
Authority
CN
China
Prior art keywords
module
software system
path
dsp software
paths
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510518615.8A
Other languages
Chinese (zh)
Other versions
CN105068931A (en
Inventor
王健
郭宝龙
贾亮
闫允一
高翔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN201510518615.8A priority Critical patent/CN105068931B/en
Publication of CN105068931A publication Critical patent/CN105068931A/en
Application granted granted Critical
Publication of CN105068931B publication Critical patent/CN105068931B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Stored Programmes (AREA)

Abstract

A single-particle soft error reliability calculation method for analyzing a DSP software system is characterized by utilizing results of module division and module topological relation of the DSP software system, searching all paths between an initial module and a termination module in the module topological relation by using depth-first search, extracting a plurality of paths by using an importance sampling method and solving the single-particle soft error reliability of the DSP software system. The invention adopts a method of dividing and searching paths based on modules for engineering files, solves the problem of probability calculation of correct results output by a DSP system under the single event effect, selects partial path paths by using an importance sampling method, endows the paths with large integral contribution with larger importance weight, extracts more path samples in the part with larger importance weight, inhibits the influence of small weight paths on the system, increases the calculation efficiency, and ensures that the calculation results are more reliable because the extracted sample mean value is unbiased estimation of the overall mean value.

Description

一种分析DSP软件系统的单粒子软错误可靠性计算方法A Reliability Calculation Method for Single Event Soft Errors for Analyzing DSP Software Systems

技术领域technical field

本发明涉及一种可靠性计算方法,特别是一种分析DSP软件系统的单粒子软错误可靠性计算方法,涉及单粒子效应的DSP系统的可靠性计算,可用于指导提高系统可靠性设计,属于系统可靠性领域。The invention relates to a reliability calculation method, in particular to a single particle soft error reliability calculation method for analyzing a DSP software system, and relates to the reliability calculation of a single particle effect DSP system, which can be used to guide and improve system reliability design, belonging to System reliability field.

背景技术Background technique

航天领域中的电子器件常常会受到空间高能粒子的冲击而发生电路逻辑状态的改变,这种现象被称为单粒子效应。单粒子效应会对航天器的安全造成极大的威胁,因此人们对此十分关注。对此人们希望得到在单粒子效应下电子器件的可靠性。Electronic devices in the aerospace field are often impacted by high-energy particles in space to change the logic state of the circuit. This phenomenon is called the single event effect. The single event effect will pose a great threat to the safety of spacecraft, so people are very concerned about it. In this regard, it is desirable to obtain the reliability of electronic devices under the single event effect.

目前,有许多论文对如何评价由于单粒子效应造成的可靠性进行了研究。例如,论文“A systematic methodology to compute the architectural vulnerabilityfactors for a high-performance microprocessor”研究了发生软错误之后导致系统失效的可能性,为了量化体系结构和微体系结构的掩盖因素,论文中引入了体系结构脆弱性因子来表示某个处理器结构发生软错误后系统失效的可能性。但是这种方法只跟踪指令,而忽略了数据,分析粒度太粗。At present, there are many papers on how to evaluate the reliability due to single event effects. For example, the paper "A systematic methodology to compute the architectural vulnerabilityfactors for a high-performance microprocessor" examines the likelihood of system failure following the occurrence of soft errors. In order to quantify the masking factors of architecture and microarchitecture, the paper introduces architecture Vulnerability factor to express the probability of system failure after a soft error occurs in a certain processor structure. But this method only tracks the instructions and ignores the data, and the analysis granularity is too coarse.

发明内容SUMMARY OF THE INVENTION

本发明的技术解决问题是:克服现有技术的不足,提供了一种分析DSP软件系统的单粒子软错误可靠性计算方法,这种方法利用DSP软件系统的模块划分和模块拓扑关系的结果,使用深度优先搜索寻找模块拓扑关系中起始模块和终止模块之间的所有路径,并使用重要性抽样的方法抽取若干路径,求取DSP软件系统的单粒子软错误可靠性;本发明对工程文件采用基于模块划分并搜索路径的方法,解决了单粒子效应下DSP系统输出正确结果的概率计算问题,利用重要性抽样方法选取部分路径路径,对整体贡献大的路径赋予更大的重要性权值,在更大的重要性权值部分抽取更多的路径样本,抑制了小权值路径对系统的影响,增加了计算效率,并且抽取的样本均值为总体均值的无偏估计,使得计算结果更为可靠。The technical solution of the present invention is to overcome the deficiencies of the prior art and provide a single particle soft error reliability calculation method for analyzing a DSP software system. This method utilizes the results of the module division of the DSP software system and the module topology relationship, The depth-first search is used to find all paths between the start module and the end module in the module topology relationship, and the importance sampling method is used to extract several paths to obtain the single-particle soft error reliability of the DSP software system; The method based on module division and path search is adopted to solve the problem of probability calculation of the correct output of the DSP system under the single event effect. The importance sampling method is used to select some paths, and the paths that contribute more to the whole are given greater importance weights. , extract more path samples in the larger importance weight part, suppress the influence of the small weight path on the system, increase the calculation efficiency, and the sample mean is an unbiased estimate of the overall mean, making the calculation result more accurate to be reliable.

本发明的技术解决方案是:一种分析DSP软件系统的单粒子软错误可靠性计算方法,步骤如下:The technical solution of the present invention is: a single particle soft error reliability calculation method for analyzing a DSP software system, the steps are as follows:

(1)对DSP软件系统进行模块划分;(1) Divide the DSP software system into modules;

(2)建立步骤(1)中划分出的各模块之间的拓扑关系图G(V,E);(2) establishing the topological relationship graph G(V, E) between the modules divided in step (1);

(3)利用深度优先遍历搜索模块拓扑关系图G(V,E)中起始模块到终止模块之间所有的路径;(3) Use depth-first traversal to search all paths between the start module and the end module in the module topology graph G(V, E);

(4)设置并计算模块正确输出结果的相关参数,所述相关参数包括SERF值、当前模块本身受到单粒子翻转的错误率、传播错误率、模块从失效状态转移到正常状态的概率μ、马尔可夫模型的连续时间转移矩阵和各个模块输出正确结果的概率;(4) Set and calculate the relevant parameters of the correct output result of the module, the relevant parameters include the SERF value, the error rate of the current module itself subjected to single event flip, the propagation error rate, the probability of the module transitioning from the failure state to the normal state μ, Mal The continuous-time transition matrix of the Kov model and the probability of each module outputting the correct result;

(5)采用重要性抽样法选取路径并计算DSP系统整体输出正确结果的概率。(5) The importance sampling method is used to select the path and calculate the probability that the DSP system outputs the correct result as a whole.

所述步骤(1)中对DSP软件系统进行模块划分;具体步骤为:In the described step (1), the DSP software system is divided into modules; the concrete steps are:

(1-1)使用TI反汇编工具分析软件系统的可执行.out文件,生成反汇编文件;进而读取反汇编文件,对DSP软件系统的汇编代码进行分析;(1-1) Use the TI disassembly tool to analyze the executable .out file of the software system to generate a disassembly file; then read the disassembly file and analyze the assembly code of the DSP software system;

(1-2)对DSP软件系统的汇编代码进行第一次扫描,根据汇编代码中函数的命名规则识别出DSP软件系统工程下的各个函数,并根据每个函数内部的函数调用指令得到各个函数之间的调用关系;(1-2) Scan the assembly code of the DSP software system for the first time, identify each function under the DSP software system engineering according to the naming rules of the functions in the assembly code, and obtain each function according to the function call instruction inside each function The calling relationship between;

(1-3)对DSP软件系统的汇编代码进行第二次扫描,分析第一遍扫描得到的每个函数中的汇编指令,确定每个函数中各个模块的入口指令,并记录其入口指令对应的地址;(1-3) Scan the assembly code of the DSP software system for the second time, analyze the assembly instructions in each function obtained by the first scan, determine the entry instructions of each module in each function, and record the corresponding entry instructions the address of;

(1-4)对DSP软件系统的汇编代码进行第三次扫描,划分出各个函数中的模块,根据步骤(1-3)中得到的各个模块入口指令地址信息确定出模块的范围,从而完成DSP软件系统模块划分。(1-4) Scan the assembly code of the DSP software system for the third time, divide the modules in each function, and determine the scope of the module according to the entry instruction address information of each module obtained in step (1-3), thereby completing the DSP software system module division.

所述步骤(2)中建立步骤(1)中划分出的各模块之间的拓扑关系图G(V,E);具体为:首先获得每个函数内的模块之间的关系,然后根据步骤(1-2)中获得的各个函数间的调用关系得到整个DSP软件系统工程中模块之间的拓扑关系;用数据结构中的有向图G(V,E)表示模块的拓扑关系,其中点集V={vi|i=1,2,,...,M}表示模块,E表示两个模块之间的连接关系。In the step (2), establish the topological relationship graph G(V, E) between the modules divided in the step (1); specifically: first obtain the relationship between the modules in each function, and then according to the step The calling relationship between the functions obtained in (1-2) obtains the topological relationship between the modules in the entire DSP software system engineering; the directed graph G(V, E) in the data structure is used to represent the topological relationship of the modules, where the point The set V={v i |i=1,2,,...,M} represents a module, and E represents a connection relationship between two modules.

所述步骤(3)中利用深度优先遍历搜索模块拓扑关系图G(V,E)中起始模块到终止模块之间所有的路径;步骤为:In the described step (3), depth-first traversal is used to search all paths between the start module and the end module in the module topology relation graph G(V, E); the steps are:

(3-1)确定模块拓扑关系图G(V,E)的起始模块和终止模块;所述拓扑关系图G(V,E)的起始模块为main函数的第一个模块;所述拓扑关系图G(V,E)的终止模块为没有后继模块的模块;(3-1) Determine the start module and the end module of the module topology relation graph G(V, E); the start module of the topological relation graph G(V, E) is the first module of the main function; the described The termination module of the topological relationship graph G(V, E) is a module without a successor module;

(3-2)使用C语言标准库函数malloc函数在计算机内存中申请一个模块结构体大小的内存空间,用来存放模块的数据和指针信息;(3-2) Use the C language standard library function malloc function to apply for a memory space of the size of a module structure in the computer memory to store the data and pointer information of the module;

(3-3)将起始模块入栈,并设置为已访问;(3-3) Push the start module into the stack and set it as accessed;

(3-4)判断当前栈顶模块是否为终止模块,若是终止模块,则执行步骤(3-5),否则执行步骤(3-6);(3-4) judge whether the current stack top module is a termination module, if it is a termination module, then execute step (3-5), otherwise execute step (3-6);

(3-5)新建一个数组用于保存栈中的各个元素,数组中的所有元素构成一条路径,执行步骤(3-7);(3-5) Create a new array to save each element in the stack, all elements in the array form a path, and execute step (3-7);

(3-6)查看栈顶模块是否有后继且未访问的模块,若有,则随机选择其中一个后继且未访问的模块并将其入栈,并将该入栈模块设置为已访问,返回步骤(3-4),否则,执行步骤(3-7);(3-6) Check whether the module at the top of the stack has a subsequent and unvisited module, if so, randomly select one of the succeeding and unvisited modules and push it into the stack, and set the pushed module as visited, and return Step (3-4), otherwise, execute step (3-7);

(3-7)查看栈中元素是否为空,若是,则结束执行,否则将栈顶模块出栈,并设置为未访问,执行步骤(3-6)。(3-7) Check whether the element in the stack is empty, if so, end the execution, otherwise, pop the top module of the stack from the stack and set it as not accessed, and execute step (3-6).

所述SERF值为软错误鲁棒性因子,表示的是单粒子效应下每个bit翻转引起可见错误的概率,具体由公式:The SERF value is a soft error robustness factor, which represents the probability of visible errors caused by each bit flip under the single event effect, specifically by the formula:

Figure GDA0002353651170000051
Figure GDA0002353651170000051

给出,其中,ACE表示执行结果会对应用程序的输出结果产生影响的位;BH表示硬件结构H的位数;N表示CPU的周期数;分子表示在N个指令周期内硬件结构H中ACE位的总个数;分母表示在N个指令周期内硬件结构H中的总位数;Given, where ACE represents the bit whose execution result will affect the output result of the application; B H represents the number of bits of the hardware structure H; N represents the number of cycles of the CPU; the numerator represents the hardware structure H in N instruction cycles The total number of ACE bits; the denominator represents the total number of bits in the hardware structure H in N instruction cycles;

所述当前模块本身受到单粒子翻转的错误率具体由公式:The error rate of the current module itself subject to single event flipping is specifically defined by the formula:

λseu=SERF×Rp×Bλ seu = SERF × R p × B

给出,其中,Rp为单粒子翻转率,即DSP软件系统受到高能粒子撞击以后其数据发生翻转的概率;B为当前模块程序存储区的代码大小和数据大小之和。is given, where R p is the single-particle turnover rate, that is, the probability that the data of the DSP software system will be turned over after being hit by high-energy particles; B is the sum of the code size and data size of the current module program storage area.

所述传播错误率,具体由公式:The propagation error rate is specified by the formula:

Figure GDA0002353651170000052
Figure GDA0002353651170000052

给出,其中,λj表示进入当前模块的所有模块中第j个模块的错误率,nj为其相应的出度;所述出度为一个模块直接到达的其他模块的总数目。is given, where λj represents the error rate of the jth module in all modules entering the current module, and nj is its corresponding out-degree; the out-degree is the total number of other modules directly reached by a module.

所述模块从失效状态转移到正常状态的概率μ具体由公式:The probability μ of the module transitioning from the failed state to the normal state is specified by the formula:

Figure GDA0002353651170000053
Figure GDA0002353651170000053

给出,其中,Trc表示模块周期修复间隔;τ表示执行一次DSP软件系统重构所需要的时间。is given, where T rc represents the module cycle repair interval; τ represents the time required to perform a DSP software system reconstruction.

所述马尔可夫模型的连续时间转移矩阵具体由公式:The continuous-time transition matrix of the Markov model is specified by the formula:

Figure GDA0002353651170000061
Figure GDA0002353651170000061

给出,其中T表示连续时间转移矩阵,Δt表示时间间隔,其中λ为当前模块的错误率且λ=λseupropis given, where T represents the continuous-time transition matrix, Δt represents the time interval, where λ is the error rate of the current module and λ=λ seuprop .

所述各个模块输出正确结果的概率由公式:The probability of each module outputting the correct result is given by the formula:

Figure GDA0002353651170000062
Figure GDA0002353651170000062

给出。given.

所述步骤(5)中采用重要性抽样法选取路径并计算DSP系统整体输出正确结果的概率,具体为:In the step (5), the importance sampling method is adopted to select the path and calculate the probability of the DSP system outputting the correct result as a whole, specifically:

(5-1)记步骤(3)得到的起始模块和终止模块之间的所有路径的条数为S;(5-1) The number of all paths between the start module and the end module obtained in step (3) is S;

(5-2)对起始模块和终止模块之间的每条路径进行编号,依次为1,2,…,S;(5-2) Number each path between the start module and the end module, 1, 2, ..., S in sequence;

(5-3)计算每一条路径的重要性权值,具体由公式:(5-3) Calculate the importance weight of each path, specifically by the formula:

Figure GDA0002353651170000063
Figure GDA0002353651170000063

给出,其中,wi为第i条路径的重要性权值;li表示第i条路径中包含的模块个数,nj表示第i条路径中第j个模块的出度;is given, where w i is the importance weight of the i-th path; l i represents the number of modules contained in the i-th path, and n j represents the out-degree of the j-th module in the i-th path;

(5-4)根据步骤(5-3)中路径的重要性权值对所有的路径进行排序;(5-4) Sort all the paths according to the importance weights of the paths in step (5-3);

(5-5)将排序后的路径均匀分成k段,则每段的路径条数为:

Figure GDA0002353651170000064
(5-5) Divide the sorted paths into k segments evenly, then the number of paths in each segment is:
Figure GDA0002353651170000064

(5-6)计算每一段路径的采样率,具体由公式:(5-6) Calculate the sampling rate of each path, specifically by the formula:

Figure GDA0002353651170000071
Figure GDA0002353651170000071

给出,其中ηi为第i段路径的采样率,wj为第i段路径中第j个路径的重要性权值;is given, where η i is the sampling rate of the i-th path, and w j is the importance weight of the j-th path in the i-th path;

(5-7)对每一段路径进行随机采样,采样的样本个数由公式:(5-7) Randomly sample each path, and the number of samples to be sampled is determined by the formula:

ni=S×ηi i=1,2,..,k;n i =S×η i i =1,2,..,k;

给出,其中ni为第i段路径的采样样本个数;given, where n i is the number of sampling samples of the i-th path;

(5-8)读取样本中每条路径中每个模块正常工作的概率;(5-8) Read the probability that each module in each path in the sample works normally;

(5-9)计算起始模块到终止模块之间所有路径输出正确结果的概率,具体由公式:(5-9) Calculate the probability that all paths between the start module and the end module output the correct result, specifically by the formula:

Figure GDA0002353651170000072
Figure GDA0002353651170000072

给出,其中nm为第m段采样的样本个数,即路径条数;Pi m为第m段中起始模块到终止模块之间第i条路径输出正确结果的概率;

Figure GDA0002353651170000073
为m段中第i条路径的第j个模块输出正确结果的概率;given, where n m is the number of samples sampled in the mth segment, that is, the number of paths; P i m is the probability that the i-th path between the start module and the termination module in the mth segment outputs the correct result;
Figure GDA0002353651170000073
The probability of outputting the correct result for the jth module of the ith path in the m segment;

(5-10)计算起始模块到终止模块之间抽取的n条路径输出正确结果的加权均值

Figure GDA0002353651170000074
即为DSP软件系统整体输出正确结果的概率,具体由公式:(5-10) Calculate the weighted mean of the correct results of the n paths extracted between the start module and the end module
Figure GDA0002353651170000074
It is the probability that the DSP software system outputs the correct result as a whole, which is specified by the formula:

Figure GDA0002353651170000075
其中
Figure GDA0002353651170000076
Figure GDA0002353651170000075
in
Figure GDA0002353651170000076

给出。given.

本发明与现有技术相比的有益效果是:The beneficial effects of the present invention compared with the prior art are:

(1)与传统的计算方法相比,本发明在DSP拓扑关系图的基础上采用深度优先遍历算法进行路径的搜索,可以充分利用拓扑关系图的信息,使得路径搜索更为准确,解决了单粒子软错误造成的系统功能评估。(1) Compared with the traditional calculation method, the present invention adopts the depth-first traversal algorithm to search the path on the basis of the DSP topology relationship diagram, which can make full use of the information of the topology relationship diagram, making the path search more accurate, and solving the problem of single Evaluation of system function due to particle soft errors.

(2)本发明引入了重要性抽样方法选取部分路径,对整体贡献大的路径赋予更大的重要性权值,在更大的重要性权值部分抽取更多的路径样本,抑制了小权值路径对系统的影响,增加了计算效率,并且抽取的样本均值为总体均值的无偏估计,使得计算结果更为可靠。(2) The present invention introduces the importance sampling method to select partial paths, assigns greater importance weights to the paths with greater overall contribution, and extracts more path samples from the larger importance weights, suppressing small weights The influence of the value path on the system increases the computational efficiency, and the sample mean is an unbiased estimate of the population mean, making the calculation result more reliable.

附图说明Description of drawings

图1是本发明采用重要性抽样的DSP系统单粒子软错误可靠性概率计算方法流程图;Fig. 1 is the flow chart of the method for calculating the reliability probability of single particle soft error reliability of the DSP system that adopts importance sampling in the present invention;

图2是模块结构体示意图;Fig. 2 is a schematic diagram of a module structure;

图3模块单链表存储形式;Figure 3 module singly linked list storage form;

图4是工程的模块关系示意图;Figure 4 is a schematic diagram of the module relationship of the project;

图5搜寻工程的模块关系流程图。Figure 5 is a flow chart of the module relationship of the search project.

具体实施方式Detailed ways

下面结合附图对本发明的具体实施方式进行进一步的详细描述。The specific embodiments of the present invention will be further described in detail below with reference to the accompanying drawings.

本发明的核心思想是对划分后的模块拓扑关系图进行基于深度优先搜索的路径寻找,然后使用重要性采样对所有的路径进行采样,使用样本值来评估DSP系统单粒子软错误可靠性。因此,首先根据现有的模块划分方法进行模块划分和模块拓扑关系建立;然后,使用深度优先遍历搜索模块拓扑关系的所有路径;最后,采用重要性采样方法从所有路径中抽取样本,求解样本路径的加权平均值,得到DSP系统单粒子软错误可靠性。The core idea of the present invention is to find the path based on depth-first search on the divided topological relationship graph of modules, then use importance sampling to sample all paths, and use the sample value to evaluate the single-particle soft error reliability of the DSP system. Therefore, firstly, the module division and module topology relationship are established according to the existing module division method; then, all paths of the module topology relationship are searched using depth-first traversal; finally, the importance sampling method is used to extract samples from all paths to solve the sample path. The weighted average of , obtains the single-event soft error reliability of the DSP system.

图1为本发明的方法流程图,从图1可知,本发明提出的一种分析DSP软件系统的单粒子软错误可靠性计算方法,其实现步骤如下:Fig. 1 is the method flow chart of the present invention, as can be seen from Fig. 1, a kind of single particle soft error reliability calculation method of analyzing DSP software system proposed by the present invention, its realization steps are as follows:

(1)对DSP软件系统进行模块划分;具体为:(1) Divide the DSP software system into modules; specifically:

(1-1)使用TI反汇编工具分析软件系统的可执行.out文件,生成反汇编文件;进而读取反汇编文件,对DSP软件系统的汇编代码进行分析;(1-1) Use the TI disassembly tool to analyze the executable .out file of the software system to generate a disassembly file; then read the disassembly file and analyze the assembly code of the DSP software system;

(1-2)对DSP软件系统的汇编代码进行第一次扫描,根据汇编代码中函数的命名规则识别出DSP软件系统工程下的各个函数,并根据每个函数内部的函数调用指令得到各个函数之间的调用关系,函数之间的调用关系示意图如图4所示,图中各个字母代表不同的函数;(1-2) Scan the assembly code of the DSP software system for the first time, identify each function under the DSP software system engineering according to the naming rules of the functions in the assembly code, and obtain each function according to the function call instruction inside each function The calling relationship between the functions and the calling relationship between the functions are shown in Figure 4, and each letter in the figure represents different functions;

(1-3)对DSP软件系统的汇编代码进行第二次扫描,分析第一遍扫描得到的每个函数中的汇编指令,确定每个函数中各个模块的入口指令,并记录其入口指令对应的地址;(1-3) Scan the assembly code of the DSP software system for the second time, analyze the assembly instructions in each function obtained by the first scan, determine the entry instructions of each module in each function, and record the corresponding entry instructions the address of;

(1-4)对DSP软件系统的汇编代码进行第三次扫描,划分出各个函数中的模块,根据步骤(1-3)中得到的各个模块入口指令地址信息确定出模块的范围,从而完成DSP软件系统模块划分,建立结构体存储模块的信息,结构体示意图如图2所示,数据区表示对应模块的信息,指针区存储邻接模块的地址。(1-4) Scan the assembly code of the DSP software system for the third time, divide the modules in each function, and determine the scope of the module according to the entry instruction address information of each module obtained in step (1-3), thereby completing the The DSP software system is divided into modules, and a structure is established to store the information of the module. The schematic diagram of the structure is shown in Figure 2. The data area represents the information of the corresponding module, and the pointer area stores the address of the adjacent module.

(2)建立步骤(1)中划分出的各模块之间的拓扑关系图G(V,E);具体为:(2) Establish the topological relationship graph G(V, E) between the modules divided in step (1); specifically:

建立步骤(1)中划分出的各模块之间的拓扑关系图G(V,E);具体为:首先获得每个函数内的模块之间的关系,然后根据步骤(1-2)中获得的各个函数间的调用关系得到整个DSP软件系统工程中模块之间的拓扑关系;算法示意图如图5所示,首先获得各个函数的模块拓扑关系,然后从main函数开始分析各个模块,若有模块调用其它函数,则由图4得到的函数调用关系,将被调用函数的模块拓扑关系插入到此模块处,继续分析被调用函数模块的信息,直到将整个软件的函数模块分析完,得到DSP软件系统工程中的模块之间的拓扑关系;用数据结构中的有向图G(V,E)表示模块的拓扑关系。其中点集V={vi|i=1,2,,...,M}表示模块,E表示两个模块之间的连接关系,拓扑图的存储结构示意图如图3所示,以邻接链表的形式存储了每一个模块的数据信息以及该模块的所有邻接模块。Establish a topological relationship graph G(V, E) between the modules divided in step (1); specifically: first obtain the relationship between modules in each function, and then obtain according to step (1-2) The calling relationship between the various functions of the DSP can obtain the topology relationship between modules in the entire DSP software system engineering; the schematic diagram of the algorithm is shown in Figure 5. First, the module topology relationship of each function is obtained, and then each module is analyzed from the main function. If there are modules To call other functions, the function call relationship obtained in Figure 4, insert the module topology relationship of the called function into this module, continue to analyze the information of the called function module, until the function module of the entire software is analyzed, and the DSP software is obtained. Topological relationship between modules in systems engineering; the topological relationship of modules is represented by a directed graph G(V, E) in the data structure. The point set V={v i |i=1,2,,...,M} represents the module, and E represents the connection between the two modules. The schematic diagram of the storage structure of the topology map is shown in Figure 3. The data information of each module and all adjacent modules of the module are stored in the form of a linked list.

(3)利用深度优先遍历搜索模块拓扑关系图G(V,E)中起始模块到终止模块之间所有的路径;具体为:(3) Use depth-first traversal to search for all paths between the start module and the end module in the module topology graph G(V, E); specifically:

(3-1)确定模块拓扑关系图G(V,E)的起始模块和终止模块;所述拓扑关系图G(V,E)的起始模块为main函数的第一个模块;所述拓扑关系图G(V,E)的终止模块为没有后继模块的模块;(3-1) Determine the start module and the end module of the module topology relation graph G(V, E); the start module of the topological relation graph G(V, E) is the first module of the main function; the described The termination module of the topological relationship graph G(V, E) is a module without a successor module;

(3-2)使用C语言标准库函数malloc函数在计算机内存中申请一个模块结构体(模块结构体用于存储模块的各种信息,包括该模块在汇编代码中的起始地址,终止地址以及拓扑图中该模块的后继模块等)大小的内存空间,用来存放模块的数据和指针信息;(3-2) Use the C language standard library function malloc function to apply for a module structure in the computer memory (the module structure is used to store various information of the module, including the start address, end address and The memory space of the size of the successor module of the module in the topology diagram is used to store the data and pointer information of the module;

(3-3)将起始模块入栈,并设置为已访问;(3-3) Push the start module into the stack and set it as accessed;

(3-4)判断当前栈顶模块是否为终止模块,若是,则执行步骤(3-5),否则执行步骤(3-6);(3-4) judge whether the current stack top module is a termination module, if so, execute step (3-5), otherwise execute step (3-6);

(3-5)新建一个数组用于保存栈中的各个元素(模块),数组中的所有元素构成一条路径,执行步骤(3-7);(3-5) Create a new array to save each element (module) in the stack, all elements in the array form a path, and execute step (3-7);

(3-6)查看栈顶模块是否有后继且未访问的模块,若有,则随机选择其中一个后继且未访问的模块并将其入栈,并将该入栈模块设置为已访问,返回步骤(3-4),否则,执行步骤(3-7);(3-6) Check whether the module at the top of the stack has a subsequent and unvisited module, if so, randomly select one of the succeeding and unvisited modules and push it into the stack, and set the pushed module as visited, and return Step (3-4), otherwise, execute step (3-7);

(3-7)查看栈中元素是否为空,若是,则结束执行,否则将栈顶模块出栈,并设置为未访问,执行步骤(3-6)。(3-7) Check whether the element in the stack is empty, if so, end the execution, otherwise, pop the top module of the stack from the stack and set it as not accessed, and execute step (3-6).

(4)设置并计算模块正确输出结果的相关参数,所述相关参数包括SERF值、当前模块本身受到单粒子翻转的错误率、传播错误率、模块从失效状态转移到正常状态的概率μ、马尔可夫模型的连续时间转移矩阵和各个模块输出正确结果的概率;(4) Set and calculate the relevant parameters of the correct output result of the module, the relevant parameters include the SERF value, the error rate of the current module itself subjected to single event flip, the propagation error rate, the probability of the module transitioning from the failure state to the normal state μ, Mal The continuous-time transition matrix of the Kov model and the probability of each module outputting the correct result;

所述SERF(Soft Error Robust Factor)值为软错误鲁棒性因子,表示的是单粒子效应下每个bit翻转引起可见错误(指能够被检测出来的输出错误)的概率,具体由公式:The SERF (Soft Error Robust Factor) value is a soft error robustness factor, which represents the probability that each bit flip causes a visible error (referring to an output error that can be detected) under the single event effect, specifically by the formula:

Figure GDA0002353651170000121
Figure GDA0002353651170000121

给出,其中,ACE表示执行结果会对应用程序的输出结果产生影响的位;BH表示硬件结构H的位数;N表示CPU的周期数;分子表示在N个指令周期内硬件结构H中ACE位的总个数;分母表示在N个指令周期内硬件结构H中的总位数;Given, where ACE represents the bit whose execution result will affect the output result of the application; B H represents the number of bits of the hardware structure H; N represents the number of cycles of the CPU; the numerator represents the hardware structure H in N instruction cycles The total number of ACE bits; the denominator represents the total number of bits in the hardware structure H in N instruction cycles;

所述当前模块本身受到单粒子翻转的错误率具体由公式:The error rate of the current module itself subject to single event flipping is specifically defined by the formula:

λseu=SERF×Rp×Bλ seu = SERF × R p × B

给出,其中,Rp为单粒子翻转率,即DSP软件系统受到高能粒子撞击以后其数据发生翻转的概率;B为当前模块程序存储区的代码大小和数据大小之和。is given, where R p is the single-particle turnover rate, that is, the probability that the data of the DSP software system will be turned over after being hit by high-energy particles; B is the sum of the code size and data size of the current module program storage area.

所述传播错误率,具体由公式:The propagation error rate is specified by the formula:

Figure GDA0002353651170000122
Figure GDA0002353651170000122

给出,其中,λj表示进入当前模块的所有模块中第j个模块的错误率,nj为其相应的出度;所述出度为一个模块直接到达的其他模块的总数目。is given, where λj represents the error rate of the jth module in all modules entering the current module, and nj is its corresponding out-degree; the out-degree is the total number of other modules directly reached by a module.

所述模块从失效状态转移到正常状态的概率μ具体由公式:The probability μ of the module transitioning from the failed state to the normal state is specified by the formula:

Figure GDA0002353651170000123
Figure GDA0002353651170000123

给出,其中,Trc表示模块周期修复间隔;τ表示执行一次DSP软件系统重构所需要的时间,对于型号为C6701的DSP,其值为6.378ms。is given, where T rc represents the module cycle repair interval; τ represents the time required to perform a DSP software system reconstruction, and for the DSP of the model C6701, its value is 6.378ms.

所述马尔可夫模型的连续时间转移矩阵具体由公式:The continuous-time transition matrix of the Markov model is specified by the formula:

Figure GDA0002353651170000131
Figure GDA0002353651170000131

给出,其中T表示连续时间转移矩阵,Δt表示时间间隔,其中λ为当前模块的错误率且λ=λseupropis given, where T represents the continuous-time transition matrix, Δt represents the time interval, where λ is the error rate of the current module and λ=λ seuprop .

基于转移矩阵获取模块在时刻t+Δt处于正常工作状态和失效状态的概率。数学表达式为:Based on the transition matrix, the probability that the module is in the normal working state and the failure state at time t+Δt is obtained. The mathematical expression is:

[p0i(t+Δt),p1i(t+Δt)]=[p0i(t),p1i(t)]T[p 0i (t+Δt),p 1i (t+Δt)]=[p 0i (t),p 1i (t)]T

其中,p0i(t)表示t时刻模块在单粒子效应下仍然正常工作状态的概率,p1i(t)表示t时刻模块在单粒子效应下处于失效状态的概率。Among them, p 0i (t) represents the probability that the module is still working normally under the single event effect at time t, and p 1i (t) represents the probability that the module is in a failure state under the single event effect at time t.

将上述方程经过变形得到概率计算公式如下:The probability calculation formula is obtained by deforming the above equation as follows:

Figure GDA0002353651170000132
Figure GDA0002353651170000132

Figure GDA0002353651170000133
Figure GDA0002353651170000133

p0i(t)随着时间t的增加,这个由单粒子翻转引起的模块输出错误的概率会增大,而正常工作的概率会减小。根据上述的公式可以计算得到各个模块输出正确结果的概率,将其记为p0i(t),其中i代表每个划分的模块。p 0i (t) As time t increases, the probability of this module output error due to single event flipping increases and the probability of normal operation decreases. According to the above formula, the probability of each module outputting the correct result can be calculated, which is recorded as p 0i (t), where i represents each divided module.

于是得到各个模块输出正确结果的概率:So the probability of each module outputting the correct result is obtained:

Figure GDA0002353651170000134
Figure GDA0002353651170000134

给出。given.

(5)采用重要性抽样法选取路径并计算DSP系统整体输出正确结果的概率。具体为:(5) The importance sampling method is used to select the path and calculate the probability that the DSP system outputs the correct result as a whole. Specifically:

(5-1)记步骤(3)得到的起始模块和终止模块之间的所有路径的条数为S;(5-1) The number of all paths between the start module and the end module obtained in step (3) is S;

(5-2)对起始模块和终止模块之间的每条路径进行编号,依次为1,2,…,S;(5-2) Number each path between the start module and the end module, 1, 2, ..., S in sequence;

(5-3)计算每一条路径的重要性权值,具体由公式:(5-3) Calculate the importance weight of each path, specifically by the formula:

Figure GDA0002353651170000141
Figure GDA0002353651170000141

给出,其中,wi为第i条路径的重要性权值;li表示第i条路径中包含的模块个数,nj表示第i条路径中第j个模块的出度;is given, where w i is the importance weight of the i-th path; l i represents the number of modules contained in the i-th path, and n j represents the out-degree of the j-th module in the i-th path;

(5-4)根据步骤(5-3)中路径的重要性权值对所有的路径进行排序;(5-4) Sort all the paths according to the importance weights of the paths in step (5-3);

(5-5)将排序后的路径均匀分成k段,则每段的路径条数为:

Figure GDA0002353651170000142
(5-5) Divide the sorted paths into k segments evenly, then the number of paths in each segment is:
Figure GDA0002353651170000142

(5-6)计算每一段路径的采样率,具体由公式:(5-6) Calculate the sampling rate of each path, specifically by the formula:

Figure GDA0002353651170000143
Figure GDA0002353651170000143

给出,其中ηi为第i段路径的采样率,wj为第i段路径中第j个路径的重要性权值;is given, where η i is the sampling rate of the i-th path, and w j is the importance weight of the j-th path in the i-th path;

(5-7)对每一段路径进行随机采样,采样的样本个数由公式:(5-7) Randomly sample each path, and the number of samples to be sampled is determined by the formula:

ni=S×ηi i=1,2,..,k;n i =S×η i i =1,2,..,k;

给出,其中ni为第i段路径的采样样本个数;given, where n i is the number of sampling samples of the i-th path;

(5-8)读取样本中每条路径中每个模块正常工作的概率;(5-8) Read the probability that each module in each path in the sample works normally;

(5-9)计算起始模块到终止模块之间所有路径输出正确结果的概率,具体由公式:(5-9) Calculate the probability that all paths between the start module and the end module output the correct result, specifically by the formula:

Figure GDA0002353651170000151
Figure GDA0002353651170000151

给出,其中nm为第m段采样的样本个数,即路径条数;Pi m为第m段中起始模块到终止模块之间第i条路径输出正确结果的概率;

Figure GDA0002353651170000152
为m段中第i条路径的第j个模块输出正确结果的概率;is given, where n m is the number of samples sampled in the mth segment, that is, the number of paths; P i m is the probability that the i-th path between the start module and the termination module in the mth segment outputs the correct result;
Figure GDA0002353651170000152
The probability of outputting the correct result for the jth module of the ith path in the m segment;

(5-10)计算起始模块到终止模块之间抽取的n条路径输出正确结果的加权均值

Figure GDA0002353651170000153
即为DSP软件系统整体输出正确结果的概率,具体由公式:(5-10) Calculate the weighted mean of the correct results of the n paths extracted between the start module and the end module
Figure GDA0002353651170000153
It is the probability that the DSP software system outputs the correct result as a whole, which is specified by the formula:

Figure GDA0002353651170000154
其中
Figure GDA0002353651170000155
Figure GDA0002353651170000154
in
Figure GDA0002353651170000155

给出。given.

本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。Contents that are not described in detail in the specification of the present invention belong to the well-known technology of those skilled in the art.

Claims (10)

1.一种分析DSP软件系统的单粒子软错误可靠性计算方法,其特征在于步骤如下:1. a single particle soft error reliability calculation method of analyzing DSP software system, is characterized in that step is as follows: (1)对DSP软件系统进行模块划分;(1) Divide the DSP software system into modules; (2)建立步骤(1)中划分出的各模块之间的拓扑关系图G(V,E);(2) establishing the topological relationship graph G(V, E) between the modules divided in step (1); (3)利用深度优先遍历搜索模块拓扑关系图G(V,E)中起始模块到终止模块之间所有的路径;(3) Use depth-first traversal to search all paths between the start module and the end module in the module topology graph G(V, E); (4)设置并计算模块正确输出结果的相关参数,所述相关参数包括SERF值、当前模块本身受到单粒子翻转的错误率、传播错误率、模块从失效状态转移到正常状态的概率μ、马尔可夫模型的连续时间转移矩阵和各个模块输出正确结果的概率;(4) Set and calculate the relevant parameters of the correct output result of the module, the relevant parameters include the SERF value, the error rate of the current module itself subjected to single event flip, the propagation error rate, the probability of the module transitioning from the failure state to the normal state μ, Mal The continuous-time transition matrix of the Kov model and the probability of each module outputting the correct result; (5)采用重要性抽样法选取路径并计算DSP系统整体输出正确结果的概率。(5) The importance sampling method is used to select the path and calculate the probability that the DSP system outputs the correct result as a whole. 2.根据权利要求1所述的一种分析DSP软件系统的单粒子软错误可靠性计算方法,其特征在于:所述步骤(1)中对DSP软件系统进行模块划分;具体步骤为:2. a kind of single particle soft error reliability calculation method of analyzing DSP software system according to claim 1, is characterized in that: in described step (1), DSP software system is carried out module division; Concrete steps are: (1-1)使用TI反汇编工具分析软件系统的可执行.out文件,生成反汇编文件;进而读取反汇编文件,对DSP软件系统的汇编代码进行分析;(1-1) Use the TI disassembly tool to analyze the executable .out file of the software system to generate a disassembly file; then read the disassembly file and analyze the assembly code of the DSP software system; (1-2)对DSP软件系统的汇编代码进行第一次扫描,根据汇编代码中函数的命名规则识别出DSP软件系统工程下的各个函数,并根据每个函数内部的函数调用指令得到各个函数之间的调用关系;(1-2) Scan the assembly code of the DSP software system for the first time, identify each function under the DSP software system engineering according to the naming rules of the functions in the assembly code, and obtain each function according to the function call instruction inside each function The calling relationship between; (1-3)对DSP软件系统的汇编代码进行第二次扫描,分析第一遍扫描得到的每个函数中的汇编指令,确定每个函数中各个模块的入口指令,并记录其入口指令对应的地址;(1-3) Scan the assembly code of the DSP software system for the second time, analyze the assembly instructions in each function obtained by the first scan, determine the entry instructions of each module in each function, and record the corresponding entry instructions the address of; (1-4)对DSP软件系统的汇编代码进行第三次扫描,划分出各个函数中的模块,根据步骤(1-3)中得到的各个模块入口指令地址信息确定出模块的范围,从而完成DSP软件系统模块划分。(1-4) Scan the assembly code of the DSP software system for the third time, divide the modules in each function, and determine the scope of the module according to the entry instruction address information of each module obtained in step (1-3), thereby completing the DSP software system module division. 3.根据权利要求2所述的一种分析DSP软件系统的单粒子软错误可靠性计算方法,其特征在于:所述步骤(2)中建立步骤(1)中划分出的各模块之间的拓扑关系图G(V,E);具体为:首先获得每个函数内的模块之间的关系,然后根据步骤(1-2)中获得的各个函数间的调用关系得到整个DSP软件系统工程中模块之间的拓扑关系;用数据结构中的有向图G(V,E)表示模块的拓扑关系,其中点集V={vi|i=1,2,,...,M}表示模块,E表示两个模块之间的连接关系。3. the single particle soft error reliability calculation method of a kind of analysis DSP software system according to claim 2, is characterized in that: in described step (2), establishes between each module divided in step (1) Topological relationship graph G(V, E); specifically: first obtain the relationship between the modules in each function, and then obtain the entire DSP software system engineering according to the calling relationship between the functions obtained in step (1-2). Topological relationship between modules; the topological relationship of modules is represented by a directed graph G(V, E) in the data structure, where the point set V={v i |i=1,2,,...,M} represents module, E represents the connection relationship between two modules. 4.根据权利要求1或2或3所述的一种分析DSP软件系统的单粒子软错误可靠性计算方法,其特征在于:所述步骤(3)中利用深度优先遍历搜索模块拓扑关系图G(V,E)中起始模块到终止模块之间所有的路径;步骤为:4. the single particle soft error reliability calculation method of a kind of analysis DSP software system according to claim 1 or 2 or 3, is characterized in that: in described step (3), utilize depth-first traversal to search for module topology relation graph G All paths between the start module and the end module in (V,E); the steps are: (3-1)确定模块拓扑关系图G(V,E)的起始模块和终止模块;所述拓扑关系图G(V,E)的起始模块为main函数的第一个模块;所述拓扑关系图G(V,E)的终止模块为没有后继模块的模块;(3-1) Determine the start module and the end module of the module topology relation graph G(V, E); the start module of the topological relation graph G(V, E) is the first module of the main function; the described The termination module of the topological relationship graph G(V, E) is a module without a successor module; (3-2)使用C语言标准库函数malloc函数在计算机内存中申请一个模块结构体大小的内存空间,用来存放模块的数据和指针信息;(3-2) Use the C language standard library function malloc function to apply for a memory space of the size of a module structure in the computer memory to store the data and pointer information of the module; (3-3)将起始模块入栈,并设置为已访问;(3-3) Push the start module into the stack and set it as accessed; (3-4)判断当前栈顶模块是否为终止模块,若是终止模块,则执行步骤(3-5),否则执行步骤(3-6);(3-4) judge whether the current stack top module is a termination module, if it is a termination module, then execute step (3-5), otherwise execute step (3-6); (3-5)新建一个数组用于保存栈中的各个元素,数组中的所有元素构成一条路径,执行步骤(3-7);(3-5) Create a new array to save each element in the stack, all elements in the array form a path, and execute step (3-7); (3-6)查看栈顶模块是否有后继且未访问的模块,若有,则随机选择其中一个后继且未访问的模块并将其入栈,并将该入栈模块设置为已访问,返回步骤(3-4),否则,执行步骤(3-7);(3-6) Check whether the module at the top of the stack has a subsequent and unvisited module, if so, randomly select one of the succeeding and unvisited modules and push it into the stack, and set the pushed module as visited, and return Step (3-4), otherwise, execute step (3-7); (3-7)查看栈中元素是否为空,若是,则结束执行,否则将栈顶模块出栈,并设置为未访问,执行步骤(3-6)。(3-7) Check whether the element in the stack is empty, if so, end the execution, otherwise, pop the top module of the stack from the stack and set it as not accessed, and execute step (3-6). 5.根据权利要求1所述的一种分析DSP软件系统的单粒子软错误可靠性计算方法,其特征在于:所述SERF值为软错误鲁棒性因子,表示的是单粒子效应下每个bit翻转引起可见错误的概率,具体由公式:5. the single particle soft error reliability calculation method of a kind of analysis DSP software system according to claim 1, is characterized in that: described SERF value is soft error robustness factor, what represents is each under single particle effect. The probability of a bit flip causing a visible error is given by the formula:
Figure FDA0002353651160000031
Figure FDA0002353651160000031
给出,其中,ACE表示执行结果会对应用程序的输出结果产生影响的位;BH表示硬件结构H的位数;N表示CPU的周期数;分子表示在N个指令周期内硬件结构H中ACE位的总个数;分母表示在N个指令周期内硬件结构H中的总位数;Given, where ACE represents the bit whose execution result will affect the output result of the application; B H represents the number of bits of the hardware structure H; N represents the number of cycles of the CPU; the numerator represents the hardware structure H in N instruction cycles The total number of ACE bits; the denominator represents the total number of bits in the hardware structure H in N instruction cycles; 所述当前模块本身受到单粒子翻转的错误率具体由公式:The error rate of the current module itself subject to single event flipping is specifically defined by the formula: λseu=SERF×Rp×Bλ seu = SERF × R p × B 给出,其中,Rp为单粒子翻转率,即DSP软件系统受到高能粒子撞击以后其数据发生翻转的概率;B为当前模块程序存储区的代码大小和数据大小之和。is given, where R p is the single-particle turnover rate, that is, the probability that the data of the DSP software system will be turned over after being hit by high-energy particles; B is the sum of the code size and data size of the current module program storage area.
6.根据权利要求5所述的一种分析DSP软件系统的单粒子软错误可靠性计算方法,其特征在于:所述传播错误率,具体由公式:6. the single particle soft error reliability calculation method of a kind of analysis DSP software system according to claim 5, is characterized in that: described propagation error rate, specifically by formula:
Figure FDA0002353651160000041
Figure FDA0002353651160000041
给出,其中,λj表示进入当前模块的所有模块中第j个模块的错误率,nj为其相应的出度;所述出度为一个模块直接到达的其他模块的总数目。is given, where λj represents the error rate of the jth module in all modules entering the current module, and nj is its corresponding out-degree; the out-degree is the total number of other modules directly reached by a module.
7.根据权利要求6所述的一种分析DSP软件系统的单粒子软错误可靠性计算方法,其特征在于:所述模块从失效状态转移到正常状态的概率μ具体由公式:7. a kind of single particle soft error reliability calculation method of analyzing DSP software system according to claim 6, is characterized in that: the probability μ that described module transfers from failure state to normal state is specifically by formula:
Figure FDA0002353651160000042
Figure FDA0002353651160000042
给出,其中,Trc表示模块周期修复间隔;τ表示执行一次DSP软件系统重构所需要的时间。is given, where T rc represents the module cycle repair interval; τ represents the time required to perform a DSP software system reconstruction.
8.根据权利要求7所述的一种分析DSP软件系统的单粒子软错误可靠性计算方法,其特征在于:所述马尔可夫模型的连续时间转移矩阵具体由公式:8. the single particle soft error reliability calculation method of a kind of analysis DSP software system according to claim 7, is characterized in that: the continuous time transition matrix of described Markov model is specifically by formula:
Figure FDA0002353651160000051
Figure FDA0002353651160000051
给出,其中T表示连续时间转移矩阵,Δt表示时间间隔,其中λ为当前模块的错误率且λ=λseupropis given, where T represents the continuous-time transition matrix, Δt represents the time interval, where λ is the error rate of the current module and λ=λ seuprop .
9.根据权利要求8所述的一种分析DSP软件系统的单粒子软错误可靠性计算方法,其特征在于:所述各个模块输出正确结果的概率由公式:9. the single particle soft error reliability calculation method of a kind of analysis DSP software system according to claim 8, is characterized in that: the probability that described each module outputs correct result by formula:
Figure FDA0002353651160000052
Figure FDA0002353651160000052
给出。given.
10.根据权利要求1所述的一种分析DSP软件系统的单粒子软错误可靠性计算方法,其特征在于:所述步骤(5)中采用重要性抽样法选取路径并计算DSP系统整体输出正确结果的概率,具体为:10. a kind of single particle soft error reliability calculation method for analyzing DSP software system according to claim 1, is characterized in that: in described step (5), adopt importance sampling method to select path and calculate DSP system overall output is correct The probability of the outcome, specifically: (5-1)记步骤(3)得到的起始模块和终止模块之间的所有路径的条数为S;(5-1) The number of all paths between the start module and the end module obtained in step (3) is S; (5-2)对起始模块和终止模块之间的每条路径进行编号,依次为1,2,…,S;(5-2) Number each path between the start module and the end module, 1, 2, ..., S in sequence; (5-3)计算每一条路径的重要性权值,具体由公式:(5-3) Calculate the importance weight of each path, specifically by the formula:
Figure FDA0002353651160000053
Figure FDA0002353651160000053
给出,其中,wi为第i条路径的重要性权值;li表示第i条路径中包含的模块个数,nj表示第i条路径中第j个模块的出度;is given, where w i is the importance weight of the i-th path; l i represents the number of modules contained in the i-th path, and n j represents the out-degree of the j-th module in the i-th path; (5-4)根据步骤(5-3)中路径的重要性权值对所有的路径进行排序;(5-4) Sort all the paths according to the importance weights of the paths in step (5-3); (5-5)将排序后的路径均匀分成k段,则每段的路径条数为:
Figure FDA0002353651160000061
(5-5) Divide the sorted paths into k segments evenly, then the number of paths in each segment is:
Figure FDA0002353651160000061
(5-6)计算每一段路径的采样率,具体由公式:(5-6) Calculate the sampling rate of each path, specifically by the formula:
Figure FDA0002353651160000062
Figure FDA0002353651160000062
给出,其中ηi为第i段路径的采样率,wj为第i段路径中第j个路径的重要性权值;is given, where η i is the sampling rate of the i-th path, and w j is the importance weight of the j-th path in the i-th path; (5-7)对每一段路径进行随机采样,采样的样本个数由公式:(5-7) Randomly sample each path, and the number of samples to be sampled is determined by the formula: ni=S×ηii=1,2,..,k;n i =S×η i i =1,2,..,k; 给出,其中ni为第i段路径的采样样本个数;given, where n i is the number of sampling samples of the i-th path; (5-8)读取样本中每条路径中每个模块正常工作的概率;(5-8) Read the probability that each module in each path in the sample works normally; (5-9)计算起始模块到终止模块之间所有路径输出正确结果的概率,具体由公式:(5-9) Calculate the probability that all paths between the start module and the end module output the correct result, specifically by the formula:
Figure FDA0002353651160000063
Figure FDA0002353651160000063
给出,其中nm为第m段采样的样本个数,即路径条数;Pi m为第m段中起始模块到终止模块之间第i条路径输出正确结果的概率;
Figure FDA0002353651160000064
为m段中第i条路径的第j个模块输出正确结果的概率;
is given, where n m is the number of samples sampled in the mth segment, that is, the number of paths; P i m is the probability that the i-th path between the start module and the termination module in the mth segment outputs the correct result;
Figure FDA0002353651160000064
The probability of outputting the correct result for the jth module of the ith path in the m segment;
(5-10)计算起始模块到终止模块之间抽取的n条路径输出正确结果的加权均值
Figure FDA0002353651160000065
即为DSP软件系统整体输出正确结果的概率,具体由公式:
(5-10) Calculate the weighted mean of the correct results of the n paths extracted between the start module and the end module
Figure FDA0002353651160000065
It is the probability that the DSP software system outputs the correct result as a whole, which is specified by the formula:
Figure FDA0002353651160000071
其中
Figure FDA0002353651160000072
Figure FDA0002353651160000071
in
Figure FDA0002353651160000072
给出。given.
CN201510518615.8A 2015-08-21 2015-08-21 Single-particle soft error reliability calculation method for analyzing DSP software system Active CN105068931B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510518615.8A CN105068931B (en) 2015-08-21 2015-08-21 Single-particle soft error reliability calculation method for analyzing DSP software system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510518615.8A CN105068931B (en) 2015-08-21 2015-08-21 Single-particle soft error reliability calculation method for analyzing DSP software system

Publications (2)

Publication Number Publication Date
CN105068931A CN105068931A (en) 2015-11-18
CN105068931B true CN105068931B (en) 2020-07-14

Family

ID=54498308

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510518615.8A Active CN105068931B (en) 2015-08-21 2015-08-21 Single-particle soft error reliability calculation method for analyzing DSP software system

Country Status (1)

Country Link
CN (1) CN105068931B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109815507B (en) * 2017-11-21 2023-04-07 中国商用飞机有限责任公司 Fault sample selection method of flight control system based on sign directed graph
CN108304278B (en) * 2018-01-18 2021-03-26 西安电子科技大学 Protection method for optimal redundancy allocation of engineering code modules based on improved modular factorial algorithm

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103646129A (en) * 2013-11-22 2014-03-19 中国科学院计算技术研究所 Reliability assessment method and device applied to FPGA
CN103678123A (en) * 2013-11-29 2014-03-26 西安空间无线电技术研究所 Method applied to recognizing vulnerabilities of single-event soft errors in processor systems
CN104035834A (en) * 2014-07-02 2014-09-10 东南大学 Buffering reliability analytical method considering safeguard measures
CN104598352A (en) * 2015-01-08 2015-05-06 西安空间无线电技术研究所 Rapid reliability evaluation method for SRAM (Static Random Access Memory) type FPGA (Field Programmable Gate Array)
CN104778324A (en) * 2015-04-16 2015-07-15 合肥工业大学 Integrated circuit selective reinforcement method with aging delaying and soft error tolerance functions

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103646129A (en) * 2013-11-22 2014-03-19 中国科学院计算技术研究所 Reliability assessment method and device applied to FPGA
CN103678123A (en) * 2013-11-29 2014-03-26 西安空间无线电技术研究所 Method applied to recognizing vulnerabilities of single-event soft errors in processor systems
CN104035834A (en) * 2014-07-02 2014-09-10 东南大学 Buffering reliability analytical method considering safeguard measures
CN104598352A (en) * 2015-01-08 2015-05-06 西安空间无线电技术研究所 Rapid reliability evaluation method for SRAM (Static Random Access Memory) type FPGA (Field Programmable Gate Array)
CN104778324A (en) * 2015-04-16 2015-07-15 合肥工业大学 Integrated circuit selective reinforcement method with aging delaying and soft error tolerance functions

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Technology-aware system failure analysis in the presence of soft errors by Mixture Importance Sampling";Veit B. Kleeberger 等;《2013 IEEE International Symposium on Defect and Fault Tolerance in VLSI and Nanotechnology Systems (DFTS)》;20131004;第118-124页 *
"一种针对软错误的程序可靠性定量分析方法";徐建军 等;《电子学报》;20110331;第39卷(第3期);第675-679页 *

Also Published As

Publication number Publication date
CN105068931A (en) 2015-11-18

Similar Documents

Publication Publication Date Title
US11741373B2 (en) Turbulence field update method and apparatus, and related device thereof
US20180082215A1 (en) Information processing apparatus and information processing method
CN101650687B (en) A Realization Method of Performance Prediction for Large-Scale Parallel Programs
CN111930903A (en) System anomaly detection method and system based on deep log sequence analysis
CN107229563A (en) A kind of binary program leak function correlating method across framework
CN107103332A (en) A kind of Method Using Relevance Vector Machine sorting technique towards large-scale dataset
CN114584379B (en) Log anomaly detection method based on optimized feature extraction granularity
CN108921424B (en) Power data anomaly detection method, device, equipment and readable storage medium
CN110686633A (en) Landslide displacement prediction method and device and electronic equipment
CN105630458B (en) The Forecasting Methodology of average throughput under a kind of out-of order processor stable state based on artificial neural network
CN111343127A (en) Method, device, medium and equipment for improving crawler recognition recall rate
CN105068931B (en) Single-particle soft error reliability calculation method for analyzing DSP software system
CN112948193B (en) FPGA comprehensive tool defect detection method based on difference test
CN111770053B (en) A Malicious Program Detection Method Based on Improved Clustering and Self-similarity
CN115185920B (en) Method, device and equipment for detecting log type
CN111767546A (en) A deep learning-based input structure inference method and device
CN108133234B (en) Community detection method, device and device based on sparse subset selection algorithm
CN114580239B (en) A bridge damage identification method considering uncertainty
CN106326123A (en) Method and system for detecting array bound-crossing defect
CN117035012A (en) Fine granularity defect positioning method, system and storage medium based on abstract syntax tree path
CN116881172B (en) Software defect prediction method based on graph convolution network
Mao et al. Extracting the representative failure executions via clustering analysis based on Markov profile model
CN110232023A (en) A kind of software defect positioning method based on convolutional neural networks, device and medium
CN110928757B (en) Performance analysis method for positioning HDFS (Hadoop distributed File System) key low-efficiency function based on Bayesian network
Bhattacharyya et al. Adaptive, online models to detect and estimate gross error in SPNDs

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant