| 查看: 206 | 回复: 0 | ||
[求助]
反变换出了问题,帮我改改了
|
|
load leleccum; indx=1:2048; x=leleccum(indx);% 将信号中的第1到第2048个采样点赋给x init=2055615866; randn('seed',init);% 产生含噪信号 sn=x+18*randn(size(x)); len=length(sn); L=4; lsbior=liftwave('db4'); els={'p',[-0.125 0.125],0}; lsnew=addlift(lsbior,els); xDec=lwt(sn,lsnew,L); ca=zeros(L,len); cd=zeros(L,len); for i=1:L ca(i, =lwtcoef('a',xDec,lsnew,L,i); %提取各层近似系数与细节系数cd(i, =lwtcoef('d',xDec,lsnew,L,i);end for j=1:L %硬阈值处理细节系数 len=length(cd(j, );sigma=median(abs(cd(j)))/0.6745; r=sigma*sqrt(2*log(len)); for k=1:len if abs(cd(j,k))>r cd(j,k)=sign(cd(j,k))*(abs(cd(j,k))-r); else cd(j,k)=0; end end end ca3=ilwt(ca(4, ,cd(4, ,lsnew,4); %后面是有问题的??ca2=ilwt(ca3,cd(3, ,lsnew,3);ca1=ilwt(ca2,cd(2, ,lsnew,2);S=ilwt(ca1,cd(1, ,lsnew,1);??? In an assignment A(I) = B, the number of elements in B and I must be the same. Error in ==> ilwt at 57 x(firstIdxAPP:2:end) = a; |
» 猜你喜欢
散金币祈福
已经有187人回复
求助,如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
物理学I论文润色/翻译怎么收费?
已经有271人回复
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
深圳大学应用技术学院招聘凝聚态物理博士后
已经有116人回复
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
求助中医院临床中医中标的青年标书学习,谢谢
已经有4人回复
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复










=lwtcoef('a',xDec,lsnew,L,i); %提取各层近似系数与细节系数
回复此楼
点击这里搜索更多相关资源