Sobol灵敏性分析.rar

Sobol灵敏性分析.rar

function [y] = calY（x）
%CALY 仿真函数
%x是m*n矩阵，m是数据量，n是自变量数量

%   测试函数
%     y=sin（x（:1））+7*（sin（x（:2）））.^2+0.1*x（:3）.^4.*sin（x（:1））;
%     x（1）=dx（2）=Gx（3）=qwx（4）=Tsat
%测试函数2
    rhol=958.59;
    rhov=0.5902;
    Cpl=4.1929;
    Hfg=2258.62;
    sigma=0.0579;
    mu=25.73e-5;
    v=x（:2）/rhol;
    Ja=rhol*Cpl*x（:4）/rhov/Hfg;
%     Re=1e4;
     Re=rhol*v.*x（:1）/mu;
    Bo=x（:3）./x（:2）/Hfg;
    ratio=（rhol/rhov）^1.326;
    y=0.0683*ratio./（Re.^0.314）./（Ja+149.2*ratio./（Bo.^0.487）./（Re.^1.6））.*sqrt（sigma/9.8/（rhol-rhov））;

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        591  2018-07-20 21:32  calY.m

     文件       1366  2018-07-20 21:33  Sensitivity.m

----------- ---------  ---------- -----  ----

                 1957                    2