| 查看: 562 | 回复: 4 | ||
[求助]
matlab
|
||
| 一幅图像经过小波分解得到它的细节部分，对其细节部分积分怎么用matlab实现 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
3楼2012-11-07 09:23:57
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
4楼2012-11-07 18:19:44
5楼2012-11-08 09:05:43
100