• 大小: 21.1MB
    文件类型: .rar
    金币: 2
    下载: 1 次
    发布日期: 2023-06-14
  • 语言: 其他
  • 标签: ENKF  

资源简介

集合卡尔曼 以及相关论文。包含2中简单模型florenz63以及florenz96

资源截图

代码片段和文件信息

function Output = Assimilate( paramXaFaXfFfXRHOlYRsig )


    CL1 = Taper( param.NxOlparam.lc ‘PH‘ ) ; 
    CL2 = Taper( param.NxOlparam.lc‘HPH‘ ) ; 
    
    Fam = mean( Fa2 ) ;
    Xam = mean( Xa2 ) ;

    t = 0 ; 
    s = 0 ;
    
    while( t < param.Nt )

        t = t + param.do ;
        s = s + 1 ;

        % Propagate ensemble members to next observation time:
        for e = 1:param.Ne
            Ff( :e )= Fa( :e ) ;
            Xf( :e )= modelrun( paramXa( :e )‘Fa( :e )param.do )‘ + sqrt(param.me)*randn(param.Nx1) ;
        end      
        Xfm= mean( Xf2 ) ;
        Ffm= mean( Ff2 ) ;

        % Perturb Observations for each ensemble member:
        Yp= zeros( param.Noparam.Ne ) ;
        for e = 1:param.Ne
            Yp( :e )= Y( :t ) + sig.*randn( param.No1 ) ;
        end

        % Inflate forecast members around the mean
        for e = 1:param.Ne
            Xf( :e )= param.ga(1) * ( Xf( :e ) - Xfm ) + Xfm ; 
            Ff( :e )= param.ga(2) * ( Ff( :e ) - Ffm ) + Ffm ; 
        end
        Xfp= Xf - repmat( Xfm1param.Ne ) ;
        Ffp= Ff - repmat( Ffm1param.Ne ) ;

        % Forecast statistics:
        PXf= 1/( param.Ne-1 ) * ( Xfp * Xfp‘ ) ;
        PFf= 1/( param.Ne-1 ) * ( Ffp * Xfp‘ ) ; 

        if ( strcmp(param.aa‘joint‘) == 1 )
            
            Fa = Update_param( paramXfFfPXfPFfCL2HYpR );        %Update param ensmeble
            [ GainXXa ] = Update_state( paramXfPXfCL1CL2HYpR ); %Update state ensemble
            

        elseif ( strcmp(param.aa‘dual‘) == 1 )

            Fa = Update_param( paramXfFfPXfPFfCL2HYpR );
            
            for e = 1:param.Ne
                Ff( :e )= Fa( :e ) ;
                Xf( :e )= modelrun( paramXa( :e )‘Fa( :e )param.do )‘ + sqrt(param.me)*randn(param.Nx1) ;
            end      
            Xfm= mean( Xf2 ) ;

            for e = 1:param.Ne
                Xf( :e )= param.ga(1) * ( Xf( :e ) - Xfm ) + Xfm ; 
            end
            Xfp= Xf - repmat( Xfm1param.Ne ) ;
            PXf= 1/( param.Ne-1 ) * ( Xfp * Xfp‘ ) ;

            [ GainXXa ] = Update_state( paramXfPXfCL1CL2HYpR );


        elseif ( strcmp(param.aa‘joint-osa‘) == 1 )

            % Most recent analysis ensemble:
            Xap= Xa - repmat( Xam1param.Ne ) ;
            Fap= Fa - repmat( Fam1param.Ne ) ;
            
            PFa= 1/( param.Ne-1 ) * ( Fap * Xfp‘ ) ;
            PXa= 1/( param.Ne-1 ) * ( Xap * Xfp‘ ) ;

            Fs = Update_param( paramXfFaPXfPFaCL2HYpR );
            [ GainXXs ] = Smooth_state( paramXfXaPXfPXaCL1CL2HYpR );

            for e = 1:param.Ne
                Ff( :e )= Fs( :e ) ;
                Xf( :e )= modelrun( paramXs( :e )‘Fs( :e )param.do )‘ + sqrt(param.me)*randn(param.Nx1) ;
            end  
            Ffm= mean( Ff2 );
            Xfm= mean( Xf2 );
            
            for e = 1:param.Ne
                Xf( :e )= para

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

     文件    4856557  2017-10-01 20:28  Code\A Bayesian Consistent Dual Ensemble Kalman Filter for State-Parameter Estimation in Subsurface Hydrology.pdf

     文件    3488184  2017-10-01 20:47  Code\An iterative ensemble Kalman filter with one-step-ahead smoothing for state-parameters estimation of contaminant transport models.pdf

     文件       6148  2016-08-03 12:26  Code\Code\.DS_Store

     文件       7194  2016-06-16 14:40  Code\Code\Assimilate.m

     文件       2273  2017-10-01 15:32  Code\Code\enkf.m

     文件        206  2015-11-18 18:23  Code\Code\florenz63.m

     文件        341  2015-11-18 18:22  Code\Code\florenz96.m

     文件        467  2017-10-23 17:08  Code\Code\GetIE.m

     文件        466  2015-11-19 17:14  Code\Code\GetIE.m.bak

     文件        598  2015-11-13 16:02  Code\Code\Hmatrix.m

     文件        340  2015-11-18 18:27  Code\Code\lorenz.m

     文件        354  2015-11-18 18:19  Code\Code\modelrun.m

     文件       7167  2016-08-03 11:51  Code\Code\MSPEL.m

     文件        590  2015-11-13 16:02  Code\Code\obs_operator.m

     文件        369  2015-11-19 22:11  Code\Code\rankHist.m

     文件       1657  2016-08-03 14:06  Code\Code\README.txt

     文件       1495  2016-05-30 15:16  Code\Code\Smooth_state.m

     文件       2793  2016-03-23 10:10  Code\Code\Taper.m

     文件       1381  2016-05-27 16:28  Code\Code\Update_param.m

     文件       1481  2016-05-27 16:28  Code\Code\Update_state.m

     文件       9353  2017-10-27 10:05  Code\Code\VisualOutput.m

     文件   24514830  2017-10-01 20:48  Code\Constraining a compositional flow model with flowchemical data using an ensemblebased Kalman filter.pdf

     文件      24099  2017-10-01 15:06  Code\sample_figure_errors.pdf

     文件      28359  2017-10-01 15:07  Code\sample_figure_influence.pdf

     文件        239  2016-08-03 14:07  Code\__MACOSX\._Code

     文件        120  2016-08-03 12:26  Code\__MACOSX\Code\._.DS_Store

     文件        239  2016-06-16 14:40  Code\__MACOSX\Code\._Assimilate.m

     文件        239  2016-05-27 16:09  Code\__MACOSX\Code\._enkf.m

     文件        239  2015-11-18 18:23  Code\__MACOSX\Code\._florenz63.m

     文件        239  2015-11-18 18:22  Code\__MACOSX\Code\._florenz96.m

............此处省略20个文件信息

评论

共有 条评论

相关资源