Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#1 |
初级会员
注册日期: 2011-12-04
帖子: 1
声望力: 0 ![]() |
![]()
%得到 db4小波的提升方案
Isdb4=liftwave('db4'); %添加ELS到提升方案中 els={'p',[-0.125,0.125],0}; Isnew=addlift(Isdb4,els); %读取声音文件 FILE='minzu.wav'; y=wavread(FILE); %绘制出原始声音信号图 figure(1);plot(y); title('原始音频'); %用db4小波对读入的声音文件进行2级提升小波分解 xDec=lwt2(y,Isnew,2); %提取2级小波分解的低频系数和高频系数 ca=lwtcoef2('ca',xDec,Isnew,2,2); ch1=lwtcoef2('ch',xDec,Isnew,2,1); cv1=lwtcoef2('cv',xDec,Isnew,2,1); cd1=lwtcoef2('cd',xDec,Isnew,2,1); ch2=lwtcoef2('ch',xDec,Isnew,2,2); cv2=lwtcoef2('cv',xDec,Isnew,2,2); cd2=lwtcoef2('cd',xDec,Isnew,2,2); x=ca; lx=length(x); alpha1=0.02;alpha2=0.1; Arn_Wm1=double(imread('Arn_Wm1.jpg','jpg')); %%%%对离散小波变换的低频系数嵌入水印,嵌入强度根据系数大小定为alpha1或alpha2 I=size(x,lx); sum=0; for j=1:I sum=sum+abs(x); end aver=sum/I;%求离散小波变换的低频系数绝对值的平均值 for k=1:I if (abs(x)<aver)%系数绝对值小于平均值,嵌入强度定为alpha1,使用加性嵌入方法嵌入水印 Aw(k)=double(x+alpha1*Arn_Wm1(k)); else Aw(k)=double(x+alpha2*Arn_Wm1(k));%系数绝对值大于平均值。嵌入强度定为alpha2 ? end end 这个程序在运行的时候总是显示倒数第三行:In an assignment A(I) = B, the number of elements in B and I must be the same. Error in ==> test1 at 47 Aw(k)=double(x+alpha1*Arn_Wm1(k));%系数绝对值大于平均值。嵌入强度定为alpha2 请各位高手帮忙看下,是哪里的错误,好吗?谢谢了哦 |
![]() |
![]() |