高手都来看看

dinahlhj

十万火急,请各位大师帮忙看下我的这段程序,小女子在此谢谢了!
load facets
nbc=size(map,1);
subplot(221),image(wcodemat(X,nbc))
%title('原始图像')；
init=2055615866;randn('seed',init);
x=X/18+randn(size(X));
subplot(222),image(wcodemat(x,nbc))
%title('噪声图像')
%使用db1执行2层小波
t=wpdec2(x,2,'db1');
%画出小波包的四叉树
plot(t)
%分解小波包结点（2，0）
[t,cA,cH,cV,cD]=wpsplt(t,[2,0]);
%分解小波包结点（2，1）
[t,cA,cH,cV,cD]=wpsplt(t,[2,1]);
%分解小波包结点（2，2）
[t,cA,cH,cV,cD]=wpsplt(t,[2,2]);
%分解小波包结点（2，3）
[t,cA,cH,cV,cD]=wpsplt(t,[2,3]);
%画出小波包树tree
plot(t)
%计算最优小波包树
bt=besttree(t);%bt=bestlevt(t);
%画出最优小波包树
plot(bt)
%使用wdencmp进行图像降噪
%寻找默认值
[thr,sorh,keepapp,crit]=ddencmp('den','wp'.bt)
%小波包系数的阈值处理
keepapp=1;
T=wpthcoef(bt,keepapp,'s',thr);

问题就出在后面粗体字这块.请大家帮忙指点下.

也可直接发到我邮箱(dinahlzy2002@yahoo.com.cn)或QQ(24744304)联系.
不胜感激!!