Matlab小波包分解后如何求各频带信号的能量值？

异步赠书：Kotlin领衔10本好书      SDCC 2017之区块链技术实战线上峰会      程序员9月书讯      每周荐书：Java Web、Python极客编程（评论送书）

window.quickReplyflag = true; var isBole = false; var fasrc="http://my.csdn.net/my/favorite/miniadd?t=Matlab%e5%b0%8f%e6%b3%a2%e5%8c%85%e5%88%86%e8%a7%a3%e5%90%8e%e5%a6%82%e4%bd%95%e6%b1%82%e5%90%84%e9%a2%91%e5%b8%a6%e4%bf%a1%e5%8f%b7%e7%9a%84%e8%83%bd%e9%87%8f%e5%80%bc%ef%bc%9f&u=http://blog.csdn.net/heifan2014/article/details/72811273"

Matlab小波包分解后如何求各频带信号的能量值？

    <div class="article_manage clearfix">
    <div class="article_r">
        <span class="link_postdate">2017-05-30 19:26</span>
        <span class="link_view" title="阅读次数">600人阅读</span>
        <span class="link_comments" title="评论次数"> <a href="#comments" onclick="_gaq.push(['_trackEvent','function', 'onclick', 'blog_articles_pinglun'])">评论</a>(0)</span>
        <span class="link_collect tracking-ad" data-mod="popu_171"> <a href="javascript:void(0);" onclick="javascript:collectArticle('Matlab%e5%b0%8f%e6%b3%a2%e5%8c%85%e5%88%86%e8%a7%a3%e5%90%8e%e5%a6%82%e4%bd%95%e6%b1%82%e5%90%84%e9%a2%91%e5%b8%a6%e4%bf%a1%e5%8f%b7%e7%9a%84%e8%83%bd%e9%87%8f%e5%80%bc%ef%bc%9f','72811273');return false;" title="收藏" target="_blank">收藏</a></span>
         <span class="link_report"> <a href="#report" onclick="javascript:report(72811273,2);return false;" title="举报">举报</a></span>

    </div>
</div>    <style type="text/css">        
        .embody{
            padding:10px 10px 10px;
            margin:0 -20px;
            border-bottom:solid 1px #ededed;                
        }
        .embody_b{
            margin:0 ;
            padding:10px 0;
        }
        .embody .embody_t,.embody .embody_c{
            display: inline-block;
            margin-right:10px;
        }
        .embody_t{
            font-size: 12px;
            color:#999;
        }
        .embody_c{
            font-size: 12px;
        }
        .embody_c img,.embody_c em{
            display: inline-block;
            vertical-align: middle;               
        }
         .embody_c img{               
            width:30px;
            height:30px;
        }
        .embody_c em{
            margin: 0 20px 0 10px;
            color:#333;
            font-style: normal;
        }
</style>
<script type="text/javascript">
    $(function () {
        try
        {
            var lib = eval("("+$("#lib").attr("value")+")");
            var html = "";
            if (lib.err == 0) {
                $.each(lib.data, function (i) {
                    var obj = lib.data[i];
                    //html += '<img src="' + obj.logo + '"/>' + obj.name + "&nbsp;&nbsp;";
                    html += ' <a href="' + obj.url + '" target="_blank">';
                    html += ' <img src="' + obj.logo + '">';
                    html += ' <em><b>' + obj.name + '</b></em>';
                    html += ' </a>';
                });
                if (html != "") {
                    setTimeout(function () {
                        $("#lib").html(html);                      
                        $("#embody").show();
                    }, 100);
                }
            }      
        } catch (err)
        { }

    });
</script>
  <div class="category clearfix">
    <div class="category_l">
       <img src="http://static.blog.csdn.net/images/category_icon.jpg">
        <span>分类：</span>
    </div>
    <div class="category_r">
                <label onclick="GetCategoryArticles('6397818','heifan2014','top','72811273');">
                    <span onclick="_gaq.push(['_trackEvent','function', 'onclick', 'blog_articles_fenlei']);">wavelet<em>（7）</em></span>
                  <img class="arrow-down" src="http://static.blog.csdn.net/images/arrow_triangle _down.jpg" style="display:inline;">
                  <img class="arrow-up" src="http://static.blog.csdn.net/images/arrow_triangle_up.jpg" style="display:none;">
                    <div class="subItem">
                        <div class="subItem_t"><a href="http://blog.csdn.net/heifan2014/article/category/6397818" target="_blank">作者同类文章</a><i class="J_close">X</i></div>
                        <ul class="subItem_l" id="top_6397818">                            
                        </ul>
                    </div>
                </label>                    
    </div>
</div>

clear; clc; n=3; wpname=’db3’; % [b,a]=butter(8,[5/100 99/100]); % load a1_1-1; load k301_1-4; % Data=filter(b,a,Data); Data=Data-mean(Data); Data=detrend(Data); wpt1=wpdec(Data,n,wpname); %对数据进行小波包分解 for i=1:2^n %wpcoef(wpt1,[n,i-1])是求第n层第i个节点的系数 E(i)=norm(wpcoef(wpt1,[n,i-1]),2);%求第i个节点的范数平方，其实也就是平方和 end % disp(‘每个节点的能量E(i)’); % E % disp(‘小波包分解总能量E_total’); E_total=sum(E); %求总能量 for i=1:2^n pfir(i)= E(i)/E_total;%求每个节点的概率 end pfir % load a1_1-2; load k302_1-4; % Data=filter(b,a,Data); Data=Data-mean(Data); Data=detrend(Data); wpt2=wpdec(Data,n,wpname); %对数据进行小波包分解 for i=1:2^n %wpcoef(wpt1,[n,i-1])是求第n层第i个节点的系数 E(i)=norm(wpcoef(wpt2,[n,i-1]),2);%求第i个节点的范数平方，其实也就是平方和 end E_total=sum(E); %求总能量 for i=1:2^n psec(i)= E(i)/E_total;%求每个节点的概率 end psec % load a1_1-4; load k303_1-4; % Data=filter(b,a,Data); Data=Data-mean(Data); Data=detrend(Data); wpt3=wpdec(Data,n,wpname); %对数据进行小波包分解 for i=1:2^n %wpcoef(wpt1,[n,i-1])是求第n层第i个节点的系数 E(i)=norm(wpcoef(wpt3,[n,i-1]),2);%求第i个节点的范数平方，其实也就是平方和 end E_total=sum(E); %求总能量 for i=1:2^n pthi(i)= E(i)/E_total;%求每个节点的概率 end pthi

格式：

y = detrend(x) % 消除时间序列中的线性趋势项
y = detrend(x,’constant’) % 消除时间序列中的均值
y = detrend(x,’linear’,bp) % 分段消除时间序列中的线性趋势项,bp为分段点向量

如

detrend的一个实例

 sig = [0 1 -2 1 0 1 -2 1 0]; % 无线性趋势的信号
trend = [0 1 2 3 4 3 2 1 0]; % 有两段线性的趋势
x = sig+trend; % 将上面趋势叠加到信号上
y = detrend(x,’linear’,5) % 根据指定的分段点去除两段线性趋势


y =


   -0.0000    1.0000   -2.0000    1.0000         0    1.0000   -2.0000    1.0000    0.0000

document.getElementById("bdshell_js").src = "http://bdimg.share.baidu.com/static/js/shell_v2.js?cdnversion=" + Math.ceil(new Date()/3600000)

    <div id="digg" articleid="72811273">
        <dl id="btnDigg" class="digg digg_enable" onclick="btndigga();">

             <dt>顶</dt>
            <dd>0</dd>
        </dl>


        <dl id="btnBury" class="digg digg_enable" onclick="btnburya();">

              <dt>踩</dt>
            <dd>0</dd>               
        </dl>

    </div>
 <div class="tracking-ad" data-mod="popu_222"><a href="javascript:void(0);" target="_blank">&nbsp;</a>   </div>
<div class="tracking-ad" data-mod="popu_223"> <a href="javascript:void(0);" target="_blank">&nbsp;</a></div>
<script type="text/javascript">
    function btndigga() {
        $(".tracking-ad[data-mod='popu_222'] a").click();
    }
    function btnburya() {
        $(".tracking-ad[data-mod='popu_223'] a").click();
    }
        </script>

•