資源簡介
LMD分解方法,采用b樣條包絡(luò),其中根據(jù)原始程序進(jìn)行簡單改變,希望對大家有幫助
代碼片段和文件信息
function[PFASI]=lmd2(m)
%最后一個PF是殘余分量
%A是瞬時賦值
%SI是純調(diào)頻函數(shù),求它的瞬時頻率就是需要的頻率
c=m;
k=0;
wucha1=0.001;
n_l=nengliang(m);
temp=c;
while?1
????k=k+1;
????a=1;
%?????h=c;
%?????[pfasi]=zhaochun1(ahwucha1);
????[pfasi]=zhaochun1(acwucha1);????
????if?length(pos(pf))>length(pos(temp))
????????if?k==1
????????????PF(k:)=temp;
????????else
????????????PF(k:)=c;
????????????break;
????????end
????end
????c=c-pf;
????PF(k:)=pf;
????temp=pf;
????A(k:)=a;
????SI(k:)=si;
????c_pos=pos(c);
????n_c=nengliang(c);
????n_pf=nengliang(pf);
????%停止調(diào)節(jié)
????%1.emd用的是三次樣條求包絡(luò),要求至少3個極值點(diǎn),所以這里c的極值點(diǎn)個數(shù)也應(yīng)該至少為3
????%2.如果上一個PF的極值點(diǎn)數(shù)比下一個PF的極值點(diǎn)數(shù)少,說明結(jié)果也不正確(這個也可以作為停止條件考慮進(jìn)去)
????%上面一句是否可以等價于當(dāng)前PF的極值點(diǎn)個數(shù)一定要大于等于殘量(c)的極值點(diǎn)個數(shù)(目前是用這個作為停止條件的一個參考寫入程序)
????%3.當(dāng)前PF分量的能量應(yīng)該大于殘量c的能量(這個有待商榷)
????%4.
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????1185??2012-11-12?09:17??lmd2.m
?????文件?????????96??2010-12-15?23:40??nengliang.m
?????文件????????185??2010-12-15?23:44??pos.m
?????文件???????1216??2010-12-15?23:43??position.m
?????文件????????921??2012-11-12?09:18??text_lmd.m
?????文件???????6978??2012-11-08?16:15??zhaochun1.m
-----------?---------??----------?-----??----
????????????????10581????????????????????6
評論
共有 條評論