Question

哪位高手能帮忙解答一下此段matlab程序在选择新阈值n时，怎么也不能正确出图是怎么回事？在线急等。

x=wavread('F:\MALTLAB\work\毕设\clean\sp02.wav');
subplot(3,1,1);
plot(x);
axis([0 22521 -0.5 0.5]);

y=wavread('F:\MALTLAB\work\毕设\SNR5_sp02.wav');

subplot(3,1,2);
plot(y);
axis([0 22521 -0.5 0.5]);

%Bark 小波包分解
wptree=wpdec(y,7,'db5');%一维小波包分解采用db5,7层小波包分解
%把 Bark 小波分解包中的第五层的叶结点合并
for k=47:62
wptree=wpjoin(wptree,k);
end
%把 Bark 小波分解包中的第六层的叶结点合并
for k=79:86
wptree=wpjoin(wptree,k);
end
%把 Bark 小波分解包中的第七层的叶结点合并
for k=91:94
wptree=wpjoin(wptree,k);
end
plot(wptree);

%将各节点按相应阈值进行处理
sorh='n'; %所选阈值函数
%-----------
%处理 7 层节点
%-----------
for i=0:31
c70(:,i+1)=wpcoef(wptree,[7 i]);
q70=median(abs(c70(:,i+1)));
arr70(i+1)=(q70/0.6745).*sqrt(2*log(442));
if(i>0)
c70(:,i+1)=wthresh(c70(:,i+1),sorh,arr70(i+1));
wptree = write(wptree,'data',[7 i],c70(:,i+1)); %修改节点系数值
end
for i=48:55
c71(:,i-47)=wpcoef(wptree,[7 i]);
q71=median(abs(c71(:,i-47)));
arr71(i-47)=(q71/0.6745).*sqrt(2*log(442));
c71(:,i-47)=wthresh(c71(:,i-47),sorh,arr71(i-47));
wptree = write(wptree,'data',[7 i],c71(:,i-47));%修改节点系数值
end
%-----------
%处理 6 层节点
%-----------
for j=16:23
c60(:,j-15)=wpcoef(wptree,[6 j]);
q60=median(abs(c60(:,j-15)));
arr60(j-15)=(q60/0.6745).*sqrt(2*log(876));
c60(:,j-15)=wthresh(c60(:,j-15),sorh,arr60(j-15));
wptree = write(wptree,'data',[6 j],c60(:,j-15));%修改节点系数值
end
for j=28:31
c61(:,j-27)=wpcoef(wptree,[6 j]);
q61=median(abs(c61(:,j-27)));
arr61(j-27)=(q61/0.6745).*sqrt(2*log(876));
c61(:,j-27)=wthresh(c61(:,j-27),sorh,arr61(j-27));
wptree = write(wptree,'data',[6 j],c61(:,j-27));%修改节点系数值
end
%-----------
%处理 5 层节点
%-----------
for k=16:31
c5(:,k-15)=wpcoef(wptree,[5 k]);
q5=median(abs(c5(:,k-15)));
arr5(k-15)=(q5/0.6745).*sqrt(2*log(1743));
c5(:,k-15)=wthresh(c5(:,k-15),sorh,arr5(k-15));
wptree = write(wptree,'data',[5 k],c5(:,k-15)); %修改节点系数值
end
%重建信号
z = wprec(wptree);
subplot(3,1,3);
plot(z);
axis([0 22521 -0.5 0.5]);
wavwrite(z,8000,'F:\MALTLAB\work\毕设\WtDeSpeech.wav');
end
axis([0 22521 -0.5 0.5]);单步运行到这时能出现图

wavwrite(z,8000,'F:\MALTLAB\work\毕设\WtDeSpeech.wav');到这儿整个图就全变成蓝色的了

Answer 1

只是因为频率太高。。。放大细节看看

追问

呃，那么请问该怎么改呢？所谓放大细节是什么意思？matlab不太会用。

追答

出来的图上面的工具栏有个放大镜的图标，放大局部，看看特性。

axis([0 22521 -0.5 0.5]);单步运行到这时能出现图 单步的时候到这其实没执行这一条语句。可能是设置坐标轴的问题

追问

请问一下这样一个警告是什么意思呢？

Warning: Data clipped during write to file:F:\MALTLAB\work\毕设\WtDeSpeech.wav
> In wavwrite>PCM_Quantize at 241
In wavwrite>write_wavedat at 267
In wavwrite at 112
In wwl at 87

追答

翻译过来是数据保存的时候长度进行了截取

追问

我知道这个翻译过来的意思，但是不知道要怎么样去修改，能否指导下？

追答

waring一般不影响你的程序运行，如果对你的结果产生影响，你才需要改动的

追问

对呀，出不来全部的结果我要如何完成论文呐，所以才想请教一下怎么样更改才好，我有试着将wavwrite(z,8000,'F:\MALTLAB\work\毕设\WtDeSpeech.wav');
这条语句不要，可是结果还是一样的，出不来图。不过还是感谢你的回答。

追答

axis([0 22521 -0.5 0.5]);单步运行时能出现图 你把这句话不要看看什么结果

追问

这个不要的话就没有警告提示，但是图还是出不来。

追答

但是图还是出不来 什么叫图出不来？？

追问

就是依旧没有增强的语音信号波形图，只是没有了之前的那个警告而已。

追答

你是说plot(z)没有效果？？还有，整图全变成蓝色 你看figure工具栏有个放大镜，你放大局部看看是什么样的曲线。我觉得大概是你增强参数的问题

追问

好吧，是我弄错了，去掉axis([0 22521 -0.5 0.5]);这条语句的时候依然没图有警告。但是我单步运行的话又有波形的，虽然那波形不像增强的。
整图变成蓝色的时候放大看了下，最终也是有个波形的，只是横坐标和纵坐标全不是我设置的那个了，和去掉axis([0 22521 -0.5 0.5]);这条语句出现的一样的结果。
请问这个增强参数的问题是什么问题呢

追答

...没图你是说figure上面看不到图吧。那可能是因为y轴范围大，所以图不明显。

就是参数搭配不适当。axis([0 22521 -0.5 0.5]);改成axis([0 22521 -5 5]);试试

反正看来就是调整纵坐标范围