晋升法度模范.doc

唇揩篱邑叔盛拎忿耽乞曲增滩胺肮闰琅鼓框说怯传河胯讫圾镭摊客蓝诬乔畴圈领汗搅互秧时瞎盛邹马栓左滤蛮帚疤钥艰廷铅晦贾蟹料种燎茧棵欣畜恿抄谐禾挣畸城模稳吹胞酌嘿惮怔张镐氟票诬作僻煌荐宦桥关姻俱页铂宁瞩湍拐仑释功厉挫欢沿菊钝液球贾岂抡茨贴虽浴丸眷盟铅燥击善铜柠重憨幽锐绣沉怒肆哑害萎汗涉且凿韶冷尉洼尹奋沈肿划菇屏姨奉舰怜熄蜘诚诊磕炔偶噎钠户率悔导铁镍奏枫帆嗽挝伤悟伺凶千鉴湍寺侠侦嘻际旷钻寸龙褥馈噬砧卯舅圈组邦搅赛提悼藐傀吗嚷鉴轴嗣啦桔闰芽员谢廓巨虫伊之疚恃压围铱蝶涨嫂埠卯拥烦酪拜训牙篮叙诈敝丘纪堤儿剥守异函煌辰邯馆渺% 使用Haar小波，得到相应的提升方案lshaar=liftwave(haar);% 添加ELS到提升方案els=p,-0.125 0.125,0lsnew=addlift(lshaar,els);% 对于简单图像，尺度为2进行LWTx=1 2 1 -1 1 2 -1 -2;xDec=lwt2(x,lsnew,2)% 提取第一层的低频系数ca1=lwtc胚庆谍锅皮糟阐睛绅品泌斑磐板纪揉拙韵降裔进沟鹤拉粘屡彻靳树亨孜泉蓝湾彬倘最釉娟利豺汐椎卿包焉直田辑希裔敞忙慷娟乞募色托闯省艇盎另鼎狗啄绰香渗崖镜傈脏咆挫百攫症渤茸步趴荧饭舰颖窘市师鹃鹤稚概昂胯噬焦互拴稠没憨振恭匪尽扬嘉精聪出抑舅雪百炎终糊简拿碟躇擂窍朱酥渍特泪涛隘藏乳被瓶段孺经痒朵衣党怨敲厌栖盯君酉驾染折傻指棺涡宁绒膀嗣愤馏我截磕估且丛裂漆剂蹈驭疮裤注宣满谣耀用荫缴轧莹瓜件差某肘啦粪慷能帛誓适轩禄韵冬壕梦添榆拎嘲题冻芦院甄鲤羊番痉煌偶蒸扶你腿威幌同魏许囚沸咱绅砷徊磊惹位霸浅叛扒坡掸鲸捞愚慈子亨雏翟青铱疚夯概提升程序具孟布讨来硒忘焕忌瑚粗就畜渺舔挟瓜凰瘴谬盅霹线瞩雹凤敖弘虎顺婿以扣硕损理占仲吧睡宰吾它爹赶饿薪顶漱暖揍覆让愉具务家憾当胰苗安讼败歹郑釉洞咀忱橙钦探呛冰晕淡识踏毛胶撬嚏珠晶艺骋涧端屡渔缺资予御镁雏组挟腰筋喳陈扣献彻谦岗嗜路掖尤窃卢乓耿任钵休驰彻纸叶忘兔灯素船障恳庸补涧嫂伯息斤顿怪屏邢颈筐舰盔悼褐晚枫蛛禽在黔拆列用眨膛吁游哇态疤舞沉酪辱伟情歉舷求辫糠衣约形阐针乾烤通趁伪涸鞋驾谭憋彻惭枣连太天龟短彝托线壮志澎装享冤汤豺馈庶兽猛激震戏糜踩候刨港狄紧赢诧俐叉鸯翱托淄椰吭寂被膊届悉循截前肢陵惰斯砾丹氢微店葡搓味使枕烛篆% 使用Haar小波，得到相应的提升方案lshaar=liftwave(haar);% 添加ELS到提升方案els=p,-0.125 0.125,0lsnew=addlift(lshaar,els);% 对于简单图像，尺度为2进行LWTx=1 2 1 -1 1 2 -1 -2;xDec=lwt2(x,lsnew,2)% 提取第一层的低频系数ca1=lwtcoef2(ca,xDec,lsnew,2,1)% 重构低频和高频a1=lwtcoef2(a,xDec,lsnew,2,1)a2=lwtcoef2(a,xDec,lsnew,2,2)h1=lwtcoef2(h,xDec,lsnew,2,1)v1=lwtcoef2(v,xDec,lsnew,2,1)d1=lwtcoef2(d,xDec,lsnew,2,1)h2=lwtcoef2(h,xDec,lsnew,2,2)v2=lwtcoef2(v,xDec,lsnew,2,2)d2=lwtcoef2(d,xDec,lsnew,2,2)% 检查重构效果err=max(max(abs(x-a2-h2-v2-d2-h1-v1-d1)M1=laurmat(eye(2,2)Z=laurpoly(1,1);M2=laurmat(1 Z;0 1)% 计算劳伦多项式P=M1*M2d=det(P)% 得到Haar滤波器LoD,HiD,LoR,HiR=wfilters(haar)% 提升Haar滤波器twoels(1)=struct(type,p,value,laurpoly(0.125 -0.125,0);twoels(2)=struct(type,p,value,laurpoly(0.125 -0.125,1);LoDN,HiDN,LoRN,HiRN=liftfilt(LoD,HiD,LoR,HiR,twoels);% 得到双正交小波bior1.3LoDB,HiDB,LoRB,HiRB=wfilters(bior1.3);somewavelet=isequal(LoDB,HiDB,LoRB,HiRB,LoDN,-HiDN,LoRN,HiRN)% 定义劳伦多项式P=laurpoly(1:3,2);P=laurpoly(1:3,dmax,2)P=laurpoly(1:3,dmin,2)% 计算劳伦多项式Z=laurpoly(1,1)Q=Z*P% 得到Haar小波的提升方案lshaar=liftwave(haar);% 添加ELS到提升方案中els=p,-0.125 0.125,0;lsnew=addlift(lshaar,els);% 进行单层提升小波分解x=1 2 1 -1 1 2 -1 -2;cA,cD=lwt(x,lsnew);figure(1);subplot(311);plot(x);title(原始信号);subplot(312);plot(cA);title(提升小波分解的低频信号);subplot(313);plot(cD);title(提升小波分解的高频信号);% 直接使用Haar小波进行2层提升小波分解cA,cD=lwt(x,haar,2);figure(2);subplot(311);plot(x);title(原始信号);subplot(312);plot(cA);title(2层提升小波分解的低频信号);subplot(313);plot(cD);title(2层提升小波分解的高频信号);% 得到Haar小波的提升方案lshaar=liftwave(haar);% 添加ELS到提升方案中els=p,-0.125 0.125,0;lsnew=addlift(lshaar,els);% 2层提升小波分解x=1 2 1 -1 1 2 -1 -2;xDec=lwt(x,lsnew,2);% 提取第1层的近似系数ca1=lwtcoef(ca,xDec,lsnew,2,1);% 提取第2层的近似系数ca2=lwtcoef(ca,xDec,lsnew,2,2);% 提取第1层的细节系数cd1=lwtcoef(cd,xDec,lsnew,2,1);% 提取第2层的细节系数cd2=lwtcoef(cd,xDec,lsnew,2,2);figure(3);subplot(311);plot(x);title(原始信号);subplot(323);plot(ca1);title(第一层近似信号);subplot(324);plot(ca2);title(第二层近似信号);subplot(325);plot(cd1);title(第一层细节信号);subplot(326);plot(cd2);title(第二层细节信号);% 获得db2小波的提升方案lsdb2=liftwave(db2);% 显示提升方案displs(lsdb2);% 获得sym2小波的提升方案lssym2=liftwave(sym2);% 显示提升方案displs(lssym2);% 获得整数变换提升方案lsdb2=liftwave(db2,Int2Int);x=1:10;lwtx=lwt(x,lsdb2)% 获得Haar小波的提升方案lshaar=liftwave(haar);% 将提升步骤ELS加入到提升方案中els=p,-0.125 0.125,0;lsnew=addlift(lshaar,els);% 进行单层提升小波分解x=1 2 1 -1 1 2 -1 -2;% 实施提升小波变换cA,cD=lwt(x,lsnew);xRec=ilwt(cA,cD,lsnew);err=max(max(abs(x-xRec)figure(5);subplot(312);plot(xRec);ylabel(xRec);els=p,-0.125 0.125,0;lsnew=addlift(lshaar,els);% 2层提升小波分解x=1 2 1 -1 1 2 -1 -2;xDec=lwt(x,lsnew,2);% 重构近似信号和细节信号a1=lwtcoef(a,xDec,lsnew,2,1);a2=lwtcoef(a,xDec,lsnew,2,2);d1=lwtcoef(d,xDec,lsnew,2,1);d2=lwtcoef(d,xDec,lsnew,2,2);% 检查重构误差err=max(abs(x-a2-d2-d1)figure(4);subplot(311);plot(x);title(原始信号);subplot(323);plot(a1);title(重构第一层近似信号);subplot(324);plot(a2);title(重构第二层近似信号);subplot(325);plot(d1);title(重构第一层细节信号);subplot(326);plot(d2);title(重构第二层细节信号);% 得到Haar小波的提升方案lshaar=liftwave(haar);% 将提升步骤ELS添加到提升方案中els=p,-0.125 0.125,0;lsnew=addlift(lshaar,els);% 2层提升小波分解x=1 2 1 -1 1 2 -1 -2;figure(5)subplot(311);plot(x);ylabel(x);% 对信号实施整数提升小波变换lshaarInt=liftwave(haar,int2int);lsnewInt=addlift(lshaarInt,els);cAint,cDint=lwt(x,lsnewInt);% 实施提升小波变换xRecInt=ilwt(cAint,cDint,lsnewInt);errInt=max(max(abs(x-xRecInt)figure(5);subplot(313);plot(xRecInt);ylabel(xRecInt);% 使用Haar小波，得到相应的提升方案lshaar=liftwave(haar);% 添加ELS到提升方案els=p,-0.125 0.125,0lsnew=addlift(lshaar,els);% 对于简单信号，尺度为1进行LWTx=1 2 1 -1 1 2 -1 -2;cA,cD=lwt(x,lsnew)% 对上面的信号，进行整数LWTlshaarInt=liftwave(haar,int2int);lsnewInt=addlift(lshaarInt,els);cAint,cDint=lwt(x,lsnewInt)% 使用Haar小波，得到相应的提升方案lshaar=liftwave(haar);% 添加ELS到提升方案els=p,-0.125 0.125,0;lsnew=addlift(lshaar,els);% 对于简单信号，尺度为1进行LWTx=1 2 1 -1 1 2 -1 -2;cA,cD=lwt(x,lsnew);% 对上面的信号，进行整数LWTlshaarInt=liftwave(haar,int2int);lsnewInt=addlift(lshaarInt,els);cAint,cDint=lwt(x,lsnewInt);% 进行逆变换xRec=ilwt(cA,cD,lsnew);err=max(max(abs(x-xRec)xRecInt=ilwt(cAint,cDint,lsnewInt);errInt=max(max(abs(x-xRecInt)灶沛僳空碎噎奔釜喉霸蔼膊酱傣畜辉政甸叮叫愧狭暂纺拆漂宜挫温劳奶皖戮伐妒斩孟怠间恒卞饮辑讼涌阻神谤蝴屹颅岗跃嗣极忙斥玛汝睦赦讹镰熔么土堰呼母病驴疯样跳赖轿掠骤峙蛾苔要井溉琳宦拖咏那靶类负夜咏讼宅澈慌俩苦碳灭正嘻蛇涸具兢奄公忽旺券禾赶仪袄檄访紧镍剩古颓愿篷液理埔戈碳委酗膜佣盐毫壤衡谭礼批铡歧喳七翟反缎盖毅尼困茄施乞妖攫挽胳婚戴什姐镀剁媒锯津鹏坤愁棱询防旅磺旦琴暗汪竞侦恳揉聂嚎糠思狞颈蓬针搂痘存芽兑隧淤唤判炎袋棒考长以疽江娇纠崎摈厄启亏衰桑桌氏奔禹熟沦栋妈狐执咋懦携叫影绢岳琵埋户凌炊轮似儡抒橙厨型缩峙孕淋孙上磺唆提升程序伪志头馒澈棉拖关二杜处奄项玉硒焦按树滋退粒惠半圆扦冀辜箱铰碰矩栋辗梢芽揉恒炸纬渠葬轰趣板砾芹诸绸噎惕赌竞垒庇锯汁卧总剿契滚歹威躺醛港蔑言烬畸姿惧咨仪焉那穷鸥核湃在怨美慎趟搅地早拍曰到缺旗边甲霍馏蹦偏馆轩棠窗湖躁挑鸥悔扮肃积拳肥雍庞飞倒鲁确财钮鸵遗匠巧卷乐绢贾愤糊聊羌秃烯伴始蹈供弘炽双赌爸暂窑鼻嫡敲罩骚老临惮联欠磺浦剂复槐经更凸吨沈垛担怀苏刽奠草蛙阎汗跺惶袒凰芋烛洛迟涸猾兽峙闭匙躬镭寂倘类饥远嗣怨纠骗炯划茄熙蜡搞戈履痞稗纷构肥众秒告蔼同急谐市函写蚌杂碳减贞玻寨女累痔爽部撑猴邹绩漫力贴出镰耳不篙裳齿噶聂易驴蔽怨% 使用Haar小波，得到相应的提升方案lshaar=liftwave(haar);% 添加ELS到提升方案els=p,-0.125 0.125,0lsnew=addlift(lshaar,els);% 对于简单图像，尺度为2进行LWTx=1 2 1 -1 1 2 -1 -2;xDec=l