自己写的字典学习代码

%% %=================ART_DL算法仿真==================================%
% %===========2014年10月03日===================乔全邦======================%

% 其中清醒的认识各个部分图像矩阵的样子十分重要。

%% ============================================================== %%
function ART_DL_SAM
clc
clear all
close all
tic;
N=256;  % 图像大小
N2=N^2;
D_C_max=0.001;
R_n_max=0.05;
I=phantom（N）;  % 产生头模型图像
theta=linspace（018061）;
theta=theta（1:60）;  % 投影角度
%% ========产生投影数据======== %%
P_num=365;  % 探测器通道个数     
P=My_Func_ParallelBeamForwardProjection（‘shepp-logan‘theta N P_num）;  % 产生投影数据
%% ======获取投影矩阵====== %%
delta=1;  % 网格大小
[W_numberW_length]=My_Func_ProjectionMatrix（thetaNP_numdelta）;
%% ==========ART-DL算法 ========== %%
[Img_Outalpha_Out]=My_Func_ARTDLProcess（PNW_numberW_lengththetaP_numD_C