




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、chinese journal of mechanical engineering 88 vol.20, no.4, 2007li ning shi tielinschool of mechanical scienceand engineering,huazhong university of scienceand technology,wuhan 430074, chinaestimation of the number of correlated sources with common frequencies based on power spectral density*abstract
2、: blind source separation and estimation of the number of sources usually demand that the number of sensors should be greater than or equal to that of the sources, which, however, is very difficult to satisfy for the complex systems. a new estimating method based on power spectral density (psd) is p
3、resented. when the relation between the number of sensors and that of sources is unknown, the psd matrix is first obtained by the ratio of psd of the observation signals, and then the bound of the number of correlated sources with common frequencies can be estimated by comparing every column vector
4、of psd matrix. the effectiveness of the proposed method is verified by theoretical analysis and experiments, and the influence of noise on the estimation of number of source is simulated. key words: blind signal number of sources power spectral density0introductionblind source separation (bss), a ne
5、w signal processing method, consists of recovering unobserved sources from several observed mixture signals. though bss can be solved using various algorithms, it must be compensated by considering some special assumptions on source signals or mixing system due to the lack of priori knowledge on the
6、 sources. one of the main assumptions is that the number of sensors must be greater than or equal to that of the sources, which is very difficult to satisfy f:ir the unknown blind sources. for this reason, estimating the /itunber of sources becomes an important reseach topic.at present the estimatin
7、g methods foi (he number of sources are mainly based on principal component analysis (pca) and singular value decomposition (svd)m1. in these methods, the number of sources is expected to be equal to the number of non-zero eigenvalues or that of non-zero singular values. that is to say, the number o
8、f sensors must be greater than or equal to that of the sources, which is as difficult as bss to satisfy.in this paper, a new estimating method based on power spectral density (psd) is proposed without imposing any special requirements on the signals or mixing matrix. when the relation between the nu
9、mber of sensors and that of sources is unknown (either greater than or smaller than or equal to), the bound of the number of correlated sources with common frequency can be estimated by comparing the column vector of the psd matrix consisting of the ratio of psd of the observation signals. the simul
10、ation and an example of a pump assembly at normal and fault conditions are used to illustrate the method.1problem formulationin bss, when the noise is not considered, the classical model is the linear instantaneous combinations of sources, that is0)x(t) = as(t)where x(t) is an m-dimensional vector,
11、s(t) is an n-dimensional unknown original source vector, a is an mx dimensional unknown mixing matrix. when m is greater than or equal to , a is with full-column rank, and when m is smaller than n, a is with full-row rank and its column vectors are not proportional. the element s/,(t) of the source
12、signal s(t) can be given by1 this project is supported by national natural science foundation of china (no. 50675076). received september 9, 20o6; received in revised form april 2,2007; accepted april 16, 2007j (0 = 2a exp( j m denotes a non-commor. jreqviency in sh(t), and non-common frequencies di
13、ffrr in vaiue from each other; aydenotes a common frequency of ah sources, and common frequencies also differ in value frc-. each other; 6*, ckv respectively denote coefficient: of non-common and common frequency in sjt); a denotes the number of non-common frequencies in .*(/); kdopjtts the number o
14、f all common frequencies. so, xfc) can be expressed by*,() = zfl* 2a exp(jo0 + zc*v expooyjla a exp( jaj) + a*c*v exp(jd,0 =*=l 1.1tel 1t.|i,zaa exp(jav) + ix p(jfi,)(3)where aih is an element of ith row and ath column in athus, there is the following modelx(t) = am,5,(t) + am2s2(t) = am&(t) + ds2(t
15、) (5)where s,(t) is an n-dimensional (n = nt + n2 +- + nn) sine function vector corresponding to non-common frequencies; m is an nxn dimensional coefficient matrix; s2(t) is a -dimensional sine function vector corresponding to common frequencies; m2 is an nx.v dimensional coefficient matrix; d=av-mi
16、 is an m%v dimensional coefficient matrix. the column vectors in d are not proportional to each other and they are not proportional to the column vectors of a either, for the convenient estimation of the number of sources.formerly, when the number of sensors is smaller than that of sources in real s
17、ignals, it is extremely difficult to estimate the bound of the number of sources, and this is the problem to be solved in the paper.2 estimating the number of correlated sources with common frequencies2.1 power spectral densityconsidering eq. (3), the cross correlation function /?*() of x,t) and xjt
18、) is given by 1994-2007 chin acade c journal electron c ublishing use. a rights res rtp /c k .nchinese journal of mechanical engineering 89 j?;(r) = x,(f+r )*;(/) =jjm 7j0i lla exp( jat + r) +xexp(jav( + r)ixwi)(r-l s=lp=)fix v|uvllwlvi9 cuv 11v111 u1v mu11v oviuvv vi uul, u1vu 1 ml* ibuv9 cuvlimf |
19、xxsaiayf*mexp(j( + r)-j0)+ not reated t0 * me column vectors in a are not proportional,r_ iwi-ir-ii-ithe ratio of psd in fn h t is not ennnl to that in f/i (w(6)zz2xaa exp( - coj + )wv(t + r) +r1 j 1=1z2fl a*# p( jk(f + r) - jr) w = 22fl*ba x exp( ioit) + /,vr) and the power spectral density p(a) is
20、pjico) - 2itj2 .rj(r)exv(-jar)dr =from eq. (13), it can be found that, whether the non-common frequencies are from the same source or not, their psd ratios arenow suppose that av is a non-common frequency of another source s(0- according to eqs. (9)(12), the ratio between pu(,) (i=lm)and pa) atevist
21、he ratio of psd in eq. (12) is not equal to that in eq. (13).(14)now, for all m observation signals, if there are n non-common frequencies, using eq. (13), the psd matrix pi at the n non-common frequencies can be obtainedpm)pim)*;,()p;m)p;m)pn),)pimjpti(fl)km)p;mk)p;m)p;m)pu(o)p;m)p;m)k)*(. (*.)p;,(
22、)(jk = l,2,-,m) is not relatedsuppose wq is a non-common frequency of source st), then to k, and the m values of afcu(a) is equal. but in fact, there arepx, )-2na a b2(9) alwavs some no*se influence and calculation errors when sam- * * pling and processing signals. as a result, usually the m values
23、offor simplicity, the special case with i,j=, 2 is firstly considered, 4u,i(to) is not absolutely equal and must be replaced by the filterand we define two ratios at to,value xu(m) which could be the mean of 4ui() or the mean1 to)k)_2*2,) 2na2palpb2n alpp =(15)from eqs. (9)(11)(14) can be simplified
24、 as follows1 11v,) .,k) :k)km) ik) kmn)= -2t(12)2.lk) =jya(t)_2 a (,) 2napbm a/)km*),)-,km*.,) km*).-/i = l,2,-,n =0*0 fl.=avfrom eq. (16), it can be seen as follows.(1) in matrix pu every element has only relation to mixing matrix a, not to any others.(2) in matrix p, if the denominator of every el
25、ement is not zero, the column vectors which the non-common frequencies coming from the same source correspond to are equal.(3) because the column vectors of mixing matrix a are not proportional to each other, and suppose the denominator of every element in eq. (16) is not zero, the column vectors wh
26、ich the non-common frequencies coming from different sources correspond to are not equal.similarly, the other psd matrix p2, p3,-, pm at nnon-common frequencies which are similar to p can also be obtained, here are not described.13 psd matrix at the common frequencysuppose ) is/?(.) = 2*djja(17)as i
27、n section 2.2, the ratio between /() k, i-lm) andeq. (18) shows that the ratio between/4*(5(i)andp,(q is djdt, which is related only to the matrix d in eq. (5), not to it or any others.now, if there are v common frequency components in m observation signals, similar to the psd matrix at non-common f
28、requency, the psd matrix ft at common frequency can be expressed as follows1 1l*im) km) kmv),(y,) km) k(&y)km,) km) kmv) i i l d*n4,i i i i_i_l_iiiiial.ujiii20 40 60 80 100 120 140160 180 200 frequency /hz fig. 2 power spectrum of the mixing signalsfig. 1 shows that 60 hz and 90 hz are from the same
29、 source; 40 hz, 80 hz and 120 hz are from the same source; 10 hz, 30 hz and 70 hz are from the same source; 70 hz, 100 hz and 130 hz are from the same source. therefore, 70 hz is a common frequency while the others are non-common frequencies. in fig. 2, there are 10 frequency pulses. according to th
30、e method proposed in the paper, the ratios of psd at every frequency pulse are given in table 1.in table 1, the value greater than 20 is written as a, and the value smaller than 0.01 is written as 0. comparing every column vector, it is obvious that the vectors 10 hz and 30 hz correspond to are equa
31、l; the vectors 40 hz, 80 hz and 120 hz, 60 hz and 90 hz, 100 hz and 130 hz correspond to are equal, too. so, the lower-bound of the number of sources is 4. in table 1, the vector not close to any others is the one 70 hz corresponds to. it may be common frequency or non-common with only one frequency
32、. here, it is regarded as the non-common with only one frequency. now, it is certain that the upper-bound of the number of sources is5. the true value is inside of the bound.table 1 ratio of psd at different frequenciesratiofrequency/7hz1030406070*2,i(u)0.6480.620001.6221.0213,|()00000.864003.463rat
33、iofrequency j7hz80901001201302,l(fi)001.607000000*3,i(a)00000000*u()00.622000kl()000000a2j(b)0.889001.7750.8951.630note: a(,/) is the reciprocal value of a/,) which calls for but repetitive calculation for the estimating result, in table 3 only either of af,/o) and i(w) is considered.in table 3, comparing every column vector, it can be seen that there are 3 groups of column vectors that can be regarded as being close, and the 3 groups are c
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 连锁超市转让协议书
- 车位租赁合同协议书
- 顺丰司机合同协议书
- 金融委托贷款协议书
- 造价咨询股东协议书
- Brand KPIs for second-hand apparel online shops IN LOVE AGAIN in Germany-外文版培训课件(2025.2)
- 长期电力交易协议书
- 餐具货物配送协议书
- 闲置资金托管协议书
- 餐具专版定制协议书
- 《比亚迪品牌历史课件》课件
- 小学心理健康家长会课件
- 2025年4月自考00160审计学答案含评分参考
- 购买木地板合同协议
- 严重开放性肢体创伤早期救治专家共识解读
- 2025年公共安全管理考试试题及答案
- 速卖通开店考试最权威答案
- 输液导管相关静脉血栓形成中国专家共识 课件
- 国企岗位笔试题目及答案
- 2024年泉州实验中学初一新生入学考试数学试卷
- 航模课程-飞翔的梦想
评论
0/150
提交评论