资源简介
集合卡尔曼 以及相关论文。包含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 flowchemical data using an ensembleba
文件 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个文件信息
- 上一篇:美团云运维体系建设与实践
- 下一篇:USB转NRF24L01 UI上位机
评论
共有 条评论