[tplee923]

I = imread('lena.png');
X = im2double(I);
x = noise(X, 'gaussian', 0.01);
wave={'harr','db1','db2'};
figure; imshow(I); title('the original');
range=1;
index=1;
for i=1:3 %一共有45个字符串元素
for ii=1:wmaxlev(size(I),(wave(index)))
[C, S] = wavedec2(x, range, wave(index));
thr_lvd = Birge_Massart(C, S); % 计算Bige-Massa策略阈值
alpha = 0.5;

x_soft_lvd = wdenoise(x, 'lvd', 's', thr_lvd, wave(index), range); % 软阈值方法 Birge-Massart策略 计算的阈值 然后重构出的图像
figure; imshow(x_soft_lvd); title(wave(index));
range=range+1;
end
index = index +1;
end


MATLAB提示错误
??? Function 'ne' is not defined for values of class 'cell'.

Error in ==> deblankl at 13
s = s(s~=' ');

Error in ==> wmaxlev at 31
wname = deblankl(wname);