EM算法--基于贝叶斯

EM算法，基于贝叶斯的，matlab程序

%% Adaptive Basis Reconstruction %%
%   A.Kucukelbir 05-Jan-2011
%   Last Revision: 3-Mar-2012

%% Initialize MATLAB and Set Parameters
warning（‘off‘‘Images:imshow:magnificationMustBeFitForDockedFigure‘）;
fig_count = 1;

% Check MATLABPOOL state and initiate multi-core processing
% Note: configuration ‘local‘ should be properly configured for YOUR machine
if（matlabpool（‘size‘）==0）
     matlabpool local
end

% Set path to utilities
path（‘utils‘path）;

% Nesterov Parameters
nest_plot_flag    = 0;          % Plot intermediate slices of iterations of Nesterov loop
basis             = ‘s‘;        % ‘s‘ for swt16 frame ‘w‘ for regular wavelet basis
thr_type          = 0;          % 0 Soft Threshold 1 Hard Threshold
W                 = ‘coif3‘;    % Wavelet family
W_lev             = 2;          % Wavelet tree depth. Note: if basis = ‘s‘ then W_lev must = 2.
nesterov_iter_lim = 30;         % Max no. of iterations for Nesterov
nesterov_stop_lim = 0.02;       % Stopping threshold for Nesterov

%% Load Example Inputs
% proj:      the simulated class means
% data_axes: the projection directions stored as vectors in 3D
% ctfs:      the CTF associated with each class mean （with some error）
% maskR:     integer value of spherical mask used during 
load（‘exampleInputs3QD7simulated.mat‘）;
slice = ceil（size（proj1）/2）;

%% Determine STEP SIZE for Nesterov
% Note: this will depend on the size of your dataset
mu = 1e-3;

%% Run Adaptive Basis Reconstruction using Nesterov‘s Method
tic
x_est = reconstruct_by_nesterov_w_ctf（...
 		proj data_axes ctfs maskR...
 		basis ...      
 		W W_levthr_type...
 		mu nesterov_iter_lim nesterov_stop_lim fig_count slice nest_plot_flag）;
toc

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-04-05 15:11  adaptiveBasis\
     文件     2914650  2012-04-04 20:16  adaptiveBasis\3QD7reconstruction.mat
     文件        1749  2012-04-04 20:14  adaptiveBasis\adaptiveBasisReconstructionExample.m
     文件     5102052  2012-04-04 20:05  adaptiveBasis\exampleInputs3QD7simulated.mat
     文件         627  2012-04-05 15:09  adaptiveBasis\README.txt
     目录           0  2012-04-05 15:07  adaptiveBasis\utils\
     文件        1147  2012-04-02 15:27  adaptiveBasis\utils\get_stat_grad_with_ctf.m
     文件        1143  2012-04-02 15:26  adaptiveBasis\utils\get_wavelet_grad_with_ctf.m
     文件        3519  2012-02-14 20:56  adaptiveBasis\utils\iswt3.m
     文件        3815  2012-04-05 15:13  adaptiveBasis\utils\mex_back_project.c
     文件        3676  2012-04-05 15:13  adaptiveBasis\utils\mex_forward_project.c
     文件         448  2012-04-02 15:26  adaptiveBasis\utils\prox_regularizer.m
     文件        4930  2012-04-05 15:12  adaptiveBasis\utils\reconstruct_by_nesterov_w_ctf.m
     文件        2575  2012-02-14 20:56  adaptiveBasis\utils\swt16dec3.m
     文件        4765  2012-02-14 20:56  adaptiveBasis\utils\swt16rec3.m
     文件        6053  2012-02-14 20:56  adaptiveBasis\utils\swt3.m
     文件        1335  2012-04-05 15:16  license.txt