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?
????????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個文件信息