| 查看: 580 | 回复: 4 | ||
[求助]
matlab
|
||
| 一幅图像经过小波分解得到它的细节部分，对其细节部分积分怎么用matlab实现 |
回复此楼» 猜你喜欢
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
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
投票:
20