| 输入:
,t,
,
,
,
,数组,n(数组大小) |
| 输出:μ我,
,财政年度(y),F习y |(y),dF习,Wj,W,Rj,圣j |
| (1) |
初始化Nt=2
λmax + 1 |
| (2) |
函数计算gPCE系数(μ我) |
| (3) |
计算gPCE系数μ我与方程(3)和方差与方程(7) |
| (4) |
返回μ我, |
| (5) |
结束函数 |
| (6) |
函数不确定性分析(UA) |
| (7) |
计算无条件的累积分布函数财政年度(y) |
| (8) |
为我= 1到10做 |
| (9) |
集X我作为一个固定的值(平均值) |
| (10) |
计算累积分布函数F习y |(y) |
| (11) |
计算d的距离F习之间的Fy(y),F习y |(y)和方程(11) |
| (12) |
结束了 |
| (13) |
返回财政年度(y),F习y |(y),维F习 |
| (14) |
结束函数 |
| (15) |
函数GSA |
| (16) |
对所有
∈(0,2π),
,t做 |
| (17) |
转换参数x我(t)和方程(9),变换y=f(x我)y=f(t)和方程(10与方程(),进行傅里叶级数展开11) |
| (18) |
结束了 |
| (19) |
为每一个j∈Z,z0=Z−[2],
k=1、2、…Nt做 |
| (20) |
计算造成的方差模型的输入参数的变化W我与方程(15) |
| (21) |
计算的总方差模型W方程(17)和一阶灵敏度指数R我与方程(18) |
| (22) |
结束了 |
| (23) |
为
2
·马克斯{
}做 |
| (24) |
计算的总灵敏度输入项我与方程(18) |
| (25) |
结束了 |
| (26) |
返回Wj,W,Rj,圣我 |
| (27) |
结束函数 |
|
