相關(guān)辨識(shí)系統(tǒng)頻率響應(yīng)的方法技術(shù)

1、相關(guān)辨識(shí)系統(tǒng)頻率響應(yīng)的技術(shù)方法及若干問(wèn)題 羅維斌 甘肅省有色地質(zhì)調(diào)查院 甘肅 蘭州報(bào)告提綱論文工作的目的及意義方法原理循環(huán)相關(guān)法頻率響應(yīng)的直接辨識(shí)方法野外實(shí)驗(yàn)結(jié)論論文目的及意義電法、電磁法現(xiàn)在已是應(yīng)用地球物理學(xué)科最為活躍的分支電法、電磁法常用信號(hào)源方波優(yōu)化的組合波2N系列偽隨機(jī)信號(hào)最長(zhǎng)線性反饋m序列（PRBS）逆重復(fù)M序列 觀測(cè)方法時(shí)間域頻率域 方波波形及其頻譜2N系列偽隨機(jī)信號(hào)及其頻譜最長(zhǎng)線性反饋m序列m序列 (PRBS) Vs 逆重復(fù)M序列 周期序列長(zhǎng)度L=2N-1 有直流分量 線狀譜 自相關(guān)峰單向 周期序列長(zhǎng)度L=2*（2N-1） 沒(méi)有直流 線狀譜，低頻頻率更低 自相關(guān)峰正反交替近年來(lái)互

2、相關(guān)法系統(tǒng)辨識(shí)原理應(yīng)用于地球物理勘探領(lǐng)域,)()()(dsstxsgty.)()()(dsssgxxxy).()()()(KgdsssgKxy如果輸入信號(hào)是白噪聲./ )()(KFGxyMTEM 多通道瞬變電磁法MTEM觀測(cè)方式輸入信號(hào)接收信號(hào)解卷積獲得沖激響應(yīng)激發(fā)電流測(cè)量方法國(guó)內(nèi)研究進(jìn)展國(guó)內(nèi)研究進(jìn)展國(guó)內(nèi)研究進(jìn)展文獻(xiàn)存在的問(wèn)題輸入m序列信號(hào)的自相關(guān)函數(shù)不是理想的 函數(shù) ，有直流分量沖激響應(yīng)反算不適定t（ ）本文提出基于循環(huán)互相關(guān)的頻率響應(yīng)辨識(shí)方法( ) ( ) ( )( )yuseuRhhRtttt=*()()()()yuuesSjSjHjHjwwww驏琪琪桫=方法原理基于循環(huán)互相關(guān)頻率響應(yīng)辨

3、識(shí)原理( )( ,)usCtcxcorr u SS( )(,)ysCtcxcorr yout SS( )(t)CysysPFFT C( )(t)CususPFFT C( )H( )( )CysCusPP( )mod( )( ), )PcysPcus 參考信號(hào)系統(tǒng)辨識(shí)流程圖非循環(huán)相關(guān)法 Vs 循環(huán)相關(guān)法非循環(huán)相關(guān)法 Vs 循環(huán)相關(guān)法非循環(huán)自相關(guān)循環(huán)自相關(guān)循環(huán)互相關(guān)序列特性時(shí)間特性 時(shí)間延遲表現(xiàn)為互相關(guān)峰的向后移位循環(huán)互相關(guān)序列特性頻率特性 只要互相關(guān)峰時(shí)間間隔相同，其譜線分布相同仿真實(shí)驗(yàn) 逆重復(fù)M序列激發(fā)系統(tǒng)模型信號(hào)產(chǎn)生參數(shù)如下：寄存器N=7，幅度A=4,分3次激發(fā)，時(shí)鐘頻率fc分別為=2Hz，

4、20Hz和508Hz,采樣率fs同為200Hz，序列周期P=8220.2909 + 1.5949s+ 0.3418s0.3193 + 3.5228s+ sH(s) =辨識(shí)結(jié)果系統(tǒng)的階躍響應(yīng)疊加周期對(duì)辨識(shí)結(jié)果影響相關(guān)法有極強(qiáng)的抗干擾能力，多周期疊加是提高信噪比的有效方法采樣率對(duì)辨識(shí)結(jié)果影響采樣率應(yīng)能隨頻帶可調(diào)，低頻段采樣率200-300Hz為宜，高頻段觀測(cè)要提高采樣率野外實(shí)驗(yàn)實(shí)驗(yàn)地：白銀礦區(qū)儀器選用陳儒軍博士湖南強(qiáng)軍科技有限公司研制的頻率域激電儀。剖面長(zhǎng)2000米，點(diǎn)距20米，用了25個(gè)4通道接收機(jī)首尾相連，一次布置100個(gè)并行測(cè)量通道。從測(cè)線中間開(kāi)始供電。AB/2最小10米，最大為2600米，

5、共設(shè)34個(gè)供電極距。每次供電，100個(gè)并行測(cè)量通道同時(shí)測(cè)量，并同時(shí)采集發(fā)送電流信號(hào)。測(cè)線布置觀測(cè)方式 觀測(cè)系統(tǒng)現(xiàn)場(chǎng)觀測(cè)到的噪聲參考信號(hào)獲取從實(shí)測(cè)電流波形提取應(yīng)用本文算法辨識(shí)的頻率響應(yīng)在810點(diǎn)獲得的幅度譜（黑）和相位譜(紅)剖面結(jié)果測(cè)量精度頻率分辨率高800820頻率幅度相位視電阻率頻率幅度相位視電阻率0.0138890.001091-0.06106861.2890.0138890.001256-0.06122976.08190.0277780.001085-0.13808856.04490.0277780.001248-0.1389970.20660.0416670.001088-0.211

6、72858.29470.0416670.001252-0.21349972.71470.0555560.001089-0.28638859.68690.0555560.001254-0.28847974.41070.0694440.001091-0.36146860.55840.0694440.001256-0.36392975.89890.0833330.001095-0.43674863.84010.0833330.001261-0.43946979.9060.0972220.001099-0.50697867.36960.0972220.001266-0.51024984.12970.1

7、111110.001107-0.58074873.39440.1111110.001275-0.58421991.15930.1250.001115-0.64893880.01550.1250.001285-0.65298998.75490.1388890.001122-0.71991885.23210.1388890.001293-0.72431004.8760.1527780.00113-0.78974891.81520.1527780.001303-0.794331012.5820.1666670.001137-0.85962897.41670.1666670.001311-0.8643

8、41019.220.1805560.001147-0.92575904.98170.1805560.001322-0.9311027.8420.1944440.001157-0.99506912.63250.1944440.001334-1.000341036.5950.2083330.001164-1.06158918.81290.2083330.001343-1.067221044.0720.2222220.001173-1.12874925.72860.2222220.001354-1.134661052.0110.2361110.001183-1.19572933.66190.2361

9、110.001365-1.201721061.1810.250.001191-1.26324939.94940.250.001375-1.269351068.6340.2638890.001197-1.32524944.90770.2638890.001383-1.33181074.6090.2777780.00121-1.39041955.23650.2777780.001397-1.397331086.180.2916670.001219-1.45697962.11740.2916670.001408-1.463811094.3440.3055560.001225-1.5204966.91

10、380.3055560.001416-1.527411100.3820.3194440.001238-1.58524977.16730.3194440.00143-1.592541111.7740.3333330.001248-1.6502984.85050.3333330.001442-1.657461120.6620.3472220.001253-1.71637988.87610.3472220.001448-1.72351125.7690.3611110.001264-1.77758997.14510.3611110.001461-1.785531135.3220.3750.001278

11、-1.842181008.2360.3750.001477-1.849971147.7410.3888890.001284-1.908531013.2570.3888890.001485-1.916361154.3190.4027780.001282-1.963811012.0220.4027780.001484-1.972621153.2340.4861110.001336-2.342791053.9840.4861110.001547-2.352211202.7440.50.001359-2.399541072.5490.50.001573-2.409581222.474解釋參數(shù)視電阻率（

12、0.056Hz）相對(duì)相位（0.056Hz，0.125Hz）金屬因子（0.013889Hz，0.125Hz）在-70號(hào)點(diǎn)，視電阻率為92 Ohm.m， 相對(duì)相位531.9 mrad， 金屬因子824.3，2N序列循環(huán)互相關(guān)處理結(jié)果分兩次供電，第一次，時(shí)鐘頻率 Fc=1/16 采樣率fs=4096疊加周期數(shù) p=20第二次供電 時(shí)鐘頻率 Fc=1/8Hz 采樣率fs=8192 sps疊加周期數(shù) p=202N序列循環(huán)互相關(guān)處理結(jié)果 增強(qiáng)噪聲干擾計(jì)算參數(shù)：2N 7頻波，F(xiàn)0=8Hz，采樣率4096Hz，疊加周期數(shù) 202N序列循環(huán)互相關(guān)處理結(jié)果計(jì)算參數(shù)：2N 7頻波，F(xiàn)0=8Hz，采樣率4096Hz，

13、疊加周期數(shù) 60恢復(fù)的頻率響應(yīng)80.3423651020.02258892160.3419435160.011310199240.3418711940.007259607320.3418216010.005609864800.3418305770.001709758880.3418299680.00233052960.3418327440.0019739891040.3418853320.0017252141120.341805560.0015101061200.3416440450.0012313331280.3417906810.0014260681360.3418043150.00130

14、79271440.3417978170.0012024011520.3418097540.001274951600.3417196430.0010470961680.3418029430.0010733351760.3417981540.000949841840.3418609150.001115421920.3417797310.0009561942000.3417643740.0010128962160.3417783160.0007694012240.3418288730.0008435012320.3417700010.0008299572400.3417995160.00057337

15、32480.3417011230.0004975852560.3417919810.0007119182640.3417988980.0006743282720.3417948910.0006440442800.3418079760.000705992880.3417828680.0006244532960.3417980580.0006071993040.3417961990.0004947343120.3418077130.0006134383200.34168370.0005932353280.3417903420.0005681363360.3417935790.00059731634

16、40.3417909880.0005163073520.3418082880.0005241513600.341788770.0005181223680.3417943750.000466423760.3417889660.0004710483840.3417844450.0004771063920.3417478080.0006062324080.3417866380.0004191824160.3418437970.0004629684240.3417885640.0004430144320.3417941780.0003956194400.3417910210.0004056664480

17、.341824050.0003935824560.3418133540.0003312434720.341798290.0003993444800.3417856430.0003764834960.3417943560.0003828765040.3417981250.000370135120.3417929140.0003555755200.3417979720.000336695280.34179410.0003244065440.3417856250.0003324595520.3417981620.000318715680.3418146910.0003662085760.341823

18、4780.0003032935840.3417910180.0003217995920.3417941640.0003383766000.3417882920.00028472恢復(fù)的頻率響應(yīng)6080.3418429150.0003167676160.3417869020.000315736320.3417532960.0001748866400.3417844260.0002869136480.3417888840.0002980676560.3417942290.0003116286640.3417884370.0002565226720.3418080180.0002729636800.3

19、417908930.0002766086880.3417935980.0002202976960.3417901420.0002515717120.3418064410.0002105137200.3417920950.000385437280.3417976630.0002604057360.3417827510.0002464557440.3418079250.0002062647520.3417945640.0002603587600.341798350.0002488357680.3417915930.0002373617920.3417689530.0001640938000.341

20、8266560.0002326718080.3417789090.0002738778240.3417676980.0001519728320.3417798590.0002224848480.3417934270.0003316988560.3418009070.0002322758720.3418081890.0001672938800.3417953830.0002432238880.3418009240.000219698960.3417894650.000204189280.3418209590.0001989499360.3418227560.0001143259520.34180

21、9760.0002233359600.3417853010.0001919489920.3418107540.00018321810080.3417928168.73E-0510400.3417944830.00027596510560.3418107210.00017164410880.3417852720.00016974210960.3418063450.00011561111120.3418240310.00025686411200.3418202370.00016345511520.3417894430.00015909411600.3417996450.00013471711680

22、.3417924550.00011363311760.3418095840.00019856811920.3417999940.00013247712000.3417930694.54E-0512160.3417799350.00015379612240.3417683120.00022971812400.3417748320.00010141612480.3418271210.00014624812560.3417761240.00022127112800.3417915620.00014249712880.3417975760.00012693812960.3417932660.00011

23、869513040.3418086130.00018075113120.3417828020.00014118113200.3417971180.00012652413360.3418091430.00016484713520.3417906130.00014658213600.3417950790.00018249913680.3417902730.0001230613760.3418077320.00012993413840.3417894370.00014808613920.3417929039.90E-0514000.3417882330.00011631414080.34178435

24、0.0001322614320.3417854290.00010634914480.3417892520.00014317514560.3417929619.66E-0514640.3417905230.00011586714720.3418233660.00010187214800.341809275.24E-0514960.3417985940.0001320815040.3417856050.00012311115200.341795040.00013619115280.341798450.00012707515360.3417928040.00011822715440.34179731

25、30.00010370915520.3417931679.78E-0515680.3417855990.00011836515760.3417975890.00010162915920.3418173730.00015335216000.3418223919.13E-0516080.3417912850.00012298216160.3417950460.00014133416240.3417875289.66E-0516400.341787930.00013520916640.3417843950.00011267716720.341789330.00012601716800.3417951

26、980.00014024616880.3417876939.30E-0516960.3418080750.00010354117040.3417912610.00011703517120.3417922095.94E-0517200.3417896999.58E-0517360.341804275.63E-0517520.3417979670.00011204717600.3417827810.00010783317680.3418070835.89E-0517760.3417953620.00012088617840.3417987620.00010996517920.3417915490.

27、00010187418320.3417817410.00015636118560.3417797480.00010342418800.3418015260.00011328718960.3418070764.91E-0519040.3417971590.00013193719120.3418016440.00010765719200.3417894379.62E-0519760.3418114420.00011791819840.3417852729.52E-0520160.3418108198.44E-0520800.3418107688.13E-0521120.3417852178.99E

28、-0521200.3418038743.03E-0521760.3417894148.52E-0521840.3417987835.90E-0521920.3417905744.13E-0522000.3418107620.00012243322160.3417993865.91E-0522400.3417798978.74E-0523040.3417915447.93E-0523120.3417970776.28E-0523200.3417924325.66E-0523280.3418094270.00011578723360.3417829138.46E-0523440.341796771

29、6.45E-0523600.3418106721.00E-0423760.3417909468.88E-0523840.3417962610.00012411723920.3417898926.69E-0524000.3418077086.93E-0524080.3417901369.38E-0524160.3417919484.22E-0524240.3417877566.19E-0524320.341784277.93E-0524560.341784555.63E-0524720.3417899569.28E-0524800.3417921094.38E-0524880.341790198

30、6.48E-0525200.3417989938.07E-0525280.3417856617.69E-0525440.3417957518.75E-0525520.3417988587.76E-0525600.3417927947.06E-0525680.3417968465.51E-0525760.3417924355.09E-0525920.3417856527.53E-0526000.3417971965.42E-0526320.3417915667.91E-0526400.3417958749.69E-0526480.3417869065.51E-0526640.3417889639

31、.50E-0526880.3417843477.25E-0526960.3417897438.54E-0527040.3417961039.84E-0527120.3417870625.40E-0527200.3418083125.91E-0527280.3417916557.76E-0527360.341791031.82E-0527440.3417893155.65E-0527760.3417983097.21E-0527840.3417828857.38E-0528000.3417961648.29E-0528080.3417992017.09E-0528160.3417915366.5

32、0E-0528800.3417795717.10E-0529040.3418021817.68E-0529280.3417988879.74E-0529360.341802367.16E-0529440.3417894196.39E-0530080.3417852256.55E-0531360.3417851416.33E-0532000.3417893865.91E-0532080.3417979162.99E-0532160.3417886841.51E-0532400.3417988263.06E-0532640.3417798026.44E-0533280.3417915335.51E

33、-0533360.341796583.69E-0533440.34179163.22E-0533600.3417830916.46E-0533680.341796443.92E-0534000.3417912796.62E-0534160.3417895234.49E-0534320.3417908747.26E-0534400.3417909931.85E-0534480.341787274.02E-0534560.3417841545.89E-0534800.3417837093.73E-0534960.3417907077.28E-0535040.3417912572.14E-05351

34、20.3417898664.36E-0535440.3417994175.81E-0535520.3417857565.88E-0535680.3417964856.64E-0535760.3417992775.52E-0535840.3417927885.01E-0535920.3417963683.31E-0536000.3417916933.04E-0536160.3417857445.80E-0536240.3417968113.25E-0536560.3417918516.00E-0536640.341796727.72E-0536720.3417862973.77E-0536880

35、.3417900767.81E-0537120.3417842645.56E-0537200.3417901536.78E-0537280.3417970227.94E-0537360.3417864393.74E-0537520.3417920766.01E-0537680.3417889123.88E-0538000.3417986735.28E-0538080.3417830586.02E-0538240.3417969986.48E-0538320.3417996545.16E-0538400.3417915224.79E-05恢復(fù)的頻率響應(yīng)69120.3417914462.77E-0570400.3417892293.16E-0572960.3417891283.11E-0574240.3417914342.61E-0574640.3417949191.60E-0674960.34179