資源簡介
在lmd中加入白噪聲對(duì)信號(hào)進(jìn)行分析,得出仿真圖
代碼片段和文件信息
%???主要調(diào)整參數(shù)為篩選次數(shù)?在程序的?93?行附近,篩選次數(shù)對(duì)信號(hào)結(jié)果影響很大??。。。。?br/>%
function[PFASI]=lmd(m)
[ij]=size(m);
if?i>=2
????m=m‘;
end
xsize=length(m);
c=m;?????%?將原始信號(hào)賦給?c
iter=5;????%自己選定
for?i=1:1:iter
????
k=0;
wucha1=0.001;???%?誤差?0.001
n_l=nengliang(m);???%?計(jì)算信號(hào)的總能量
????
while?1
????k=k+1;????%?調(diào)整存取向量位置
????a=1;??????%?每次循環(huán)開始都為?1?
????h=c;??????%?用來傳遞殘差分量
???
????[pfasi]=zhaochun(ahwucha1);?????
????c=c-pf;???%?減去第一次得到的?pf?分量
????ppf(k:)=pf;???%?存放?pf?分量
????aa(k:)=a;?????%?幅值分量
????ssi(k:)=si;???%?頻率分量
????c_pos=position(c);???%?殘余項(xiàng)的極值點(diǎn)數(shù)
????n_c=nengliang(c);????%?殘余項(xiàng)能量
????n_pf=nengliang(pf);??%?pf?能量
????if?length(c_pos)<3?||?n_c ????????ppf(k+1:)=c;
????????pfz=k+1;?????%分解出來的pf分量k+1為殘差
????????b
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件???????7074??2012-09-21?10:52??lmd.m
?????文件???????1002??2012-09-02?20:45??text_lmd.m
-----------?---------??----------?-----??----
?????????????????8076????????????????????2
評(píng)論
共有 條評(píng)論