可靠性分析算法：一次二阶矩法的Matlab源代码

function MVFOSM（NE）
%均值一次二阶矩法
%适用范围：随机变量为正态分布，变量间独立
%特点：线性功能函数获得失效概率精确解，非线性程度不高的功能函数获得失效概率近似解
%缺点：对于物理意义相同而数学表达是不同的非线性问题可能得到完全不同的失效概率
%输入参数：NE - 极限状态函数索引，整数
    global Prob
    tic
    if nargin<1 NE = 0; end                                %如果未定义方程索引，默认使用索引0
    
    %1.定义符号变量
    syms x1  x2  x3  x4  x5  x6 x7  x8  x9  xa  xb  xc  xd  xe  xf  xg;
    X=[ x1  x2  x3  x4  x5  x6 x7  x8  x9  xa  xb  xc  xd  xe  xf  xg ];
    XX=[‘x1‘;  ‘x2‘;  ‘x3‘  ;‘x4‘ ; ‘x5‘;  ‘x6‘ ;‘x7‘;  ‘x8‘;  ‘x9‘;  ‘xa‘;  ‘xb‘;  ‘xc‘;  ‘xd‘;  ‘xe‘;  ‘xf‘; ‘xg‘];

    %2.获得方程相关公式、变量数目、变量分布、变量参数
    Prob = FunDist（ NE ）;
    
    %3.判断方法适用性
    for i=1:Prob.Nx
        if strcmp（Prob.Dist{i} ‘norm‘）==0 && strcmp（Prob.Dist{i} ‘Normal‘）==0
        	error（‘This method can only be used for normal dis