| 查看: 566 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
matlab
|
||
| 一幅图像经过小波分解得到它的细节部分，对其细节部分积分怎么用matlab实现 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有15人回复
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有20人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
-
阿姨
已经有3人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
-
filecode
已经有16人回复
-
重要来源：本周末出结果
已经有10人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有6人回复
-
不应该看fileCode
已经有9人回复
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
