| 查看: 490 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
matlab
|
||
| 一幅图像经过小波分解得到它的细节部分，对其细节部分积分怎么用matlab实现 |
回复此楼» 猜你喜欢
北京林业大学硕导招生广告
已经有4人回复
-
材料调剂
已经有7人回复
-
化学工程085602 305分求调剂
已经有10人回复
-
289求调剂
已经有15人回复
-
291 求调剂
已经有7人回复
-
274求调剂
已经有14人回复
-
309求调剂
已经有5人回复
-
292求调剂
已经有8人回复
-
求调剂
已经有4人回复
-
一志愿 西北大学 总分282 英语一62 求调剂
已经有3人回复
4楼2012-11-07 18:19:44
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
5楼2012-11-08 09:05:43
