| 查看: 571 | 回复: 4 | |||
[求助]
matlab
|
| 一幅图像经过小波分解得到它的细节部分，对其细节部分积分怎么用matlab实现 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有46人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有16人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有4人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
xiaoshui_cd
木虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 1645.5
- 帖子: 12
- 在线: 13.4小时
- 虫号: 1857770
- 注册: 2012-06-13
- 性别: MM
- 专业: 信号理论与信号处理
【答案】应助回帖
感谢参与，应助指数 +1
|
你指的是低频部分吗？有个函数可以的。 function deccoef = wavelet_dec(signal,scales,wavelet); %---------------------------------------------------------------------- % wavelet_dec.m, 函数子程序，用于多分辨率分解，使用 A'trous 算法； % 在分解的过程中，同时给出各级的细节和概貌，它们 % 和原数据有着同样的长度。 %---------------------------------------------------------------------- coeflength = size(signal,2); [Lo_D,Hi_D,Lo_R,Hi_R] = wfilters(wavelet); Lo_D = Lo_D(find(Lo_D)); Hi_D = Hi_D(find(Hi_D)); Lo_R = Lo_R(find(Lo_R)); Hi_R = Hi_R(find(Hi_R)); deccoef = []; figure(gcf); j = 0; A0 = signal; A0length = length(A0); while j < scales LoDlength = length(Lo_D); HiDlength = length(Hi_D); Dtemp = conv(Hi_D,A0); Atemp = conv(Lo_D,A0); A = wkeep(Atemp,coeflength); D = wkeep(Dtemp,coeflength); deccoef(j+1,  = D;%--------------------------------------------------------------- % Plot the Approximated coefficients and Detail coefficients subplot(scales+1,2,2*j+3); plot(A); label = strcat('a ',num2str(j+1)); Ylabel(label); subplot(scales+1,2,2*j+4); plot(D); label = strcat('d ',num2str(j+1)); Ylabel(label); %-------------------------------------------------------------- Lo = zeros(1,2*LoDlength-1); Hi = zeros(1,2*HiDlength-1); Lo(1:2:2*LoDlength-1)=Lo_D; Hi(1:2:2*HiDlength-1)=Hi_D; Lo_D = Lo; Hi_D = Hi; A0 = A; j=j+1; end deccoef(j+1, = A; |
2楼2012-11-07 08:04:30
3楼2012-11-07 09:23:57
4楼2012-11-07 18:19:44
5楼2012-11-08 09:05:43
