相关辨识系统频率响应的方法技术

1、相关辨识系统频率响应的技术方法及若干问题 罗维斌 甘肃省有色地质调查院 甘肃 兰州报告提纲论文工作的目的及意义方法原理循环相关法频率响应的直接辨识方法野外实验结论论文目的及意义电法、电磁法现在已是应用地球物理学科最为活跃的分支电法、电磁法常用信号源方波优化的组合波2N系列伪随机信号最长线性反馈m序列（PRBS）逆重复M序列 观测方法时间域频率域 方波波形及其频谱2N系列伪随机信号及其频谱最长线性反馈m序列m序列 (PRBS) Vs 逆重复M序列 周期序列长度L=2N-1 有直流分量 线状谱 自相关峰单向 周期序列长度L=2*（2N-1） 没有直流 线状谱，低频频率更低 自相关峰正反交替近年来互

2、相关法系统辨识原理应用于地球物理勘探领域,)()()(dsstxsgty.)()()(dsssgxxxy).()()()(KgdsssgKxy如果输入信号是白噪声./ )()(KFGxyMTEM 多通道瞬变电磁法MTEM观测方式输入信号接收信号解卷积获得冲激响应激发电流测量方法国内研究进展国内研究进展国内研究进展文献存在的问题输入m序列信号的自相关函数不是理想的 函数 ，有直流分量冲激响应反算不适定t（ ）本文提出基于循环互相关的频率响应辨识方法( ) ( ) ( )( )yuseuRhhRtttt=*()()()()yuuesSjSjHjHjwwww骣琪琪桫=方法原理基于循环互相关频率响应辨

3、识原理( )( ,)usCtcxcorr u SS( )(,)ysCtcxcorr yout SS( )(t)CysysPFFT C( )(t)CususPFFT C( )H( )( )CysCusPP( )mod( )( ), )PcysPcus 参考信号系统辨识流程图非循环相关法 Vs 循环相关法非循环相关法 Vs 循环相关法非循环自相关循环自相关循环互相关序列特性时间特性 时间延迟表现为互相关峰的向后移位循环互相关序列特性频率特性 只要互相关峰时间间隔相同，其谱线分布相同仿真实验 逆重复M序列激发系统模型信号产生参数如下：寄存器N=7，幅度A=4,分3次激发，时钟频率fc分别为=2Hz，

4、20Hz和508Hz,采样率fs同为200Hz，序列周期P=8220.2909 + 1.5949s+ 0.3418s0.3193 + 3.5228s+ sH(s) =辨识结果系统的阶跃响应叠加周期对辨识结果影响相关法有极强的抗干扰能力，多周期叠加是提高信噪比的有效方法采样率对辨识结果影响采样率应能随频带可调，低频段采样率200-300Hz为宜，高频段观测要提高采样率野外实验实验地：白银矿区仪器选用陈儒军博士湖南强军科技有限公司研制的频率域激电仪。剖面长2000米，点距20米，用了25个4通道接收机首尾相连，一次布置100个并行测量通道。从测线中间开始供电。AB/2最小10米，最大为2600米，

5、共设34个供电极距。每次供电，100个并行测量通道同时测量，并同时采集发送电流信号。测线布置观测方式 观测系统现场观测到的噪声参考信号获取从实测电流波形提取应用本文算法辨识的频率响应在810点获得的幅度谱（黑）和相位谱(红)剖面结果测量精度频率分辨率高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解释参数视电阻率（

12、0.056Hz）相对相位（0.056Hz，0.125Hz）金属因子（0.013889Hz，0.125Hz）在-70号点，视电阻率为92 Ohm.m， 相对相位531.9 mrad， 金属因子824.3，2N序列循环互相关处理结果分两次供电，第一次，时钟频率 Fc=1/16 采样率fs=4096叠加周期数 p=20第二次供电 时钟频率 Fc=1/8Hz 采样率fs=8192 sps叠加周期数 p=202N序列循环互相关处理结果 增强噪声干扰计算参数：2N 7频波，F0=8Hz，采样率4096Hz，叠加周期数 202N序列循环互相关处理结果计算参数：2N 7频波，F0=8Hz，采样率4096Hz，

13、叠加周期数 60恢复的频率响应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恢复的频率响应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恢复的频率响应69120.3417914462.77E-0570400.3417892293.16E-0572960.3417891283.11E-0574240.3417914342.61E-0574640.3417949191.60E-0674960.34179