define_filters.m

matlab中常用电法各种装置的电磁正演计算

function define_filters

% Function to generate the filter.mat file containing the filter 
% coefficients needed by cconvol.m to perform the convolution 
% calculations.
% 
% The filter values are based on:
% Christensen, N. B., 1990. Optimized Fast Hankel Transform filters.
%   Geophysical Prospecting 38(5), 545--568.
%
% Written by:
% Thomas Ingeman-Nielsen
% The Arctic Technology Center, BYG
% Technical University of Denmark
% Email: [e-mail]

J01 = [...
        2.30251282873E-21, 2.89869191107E-21, 3.64923690780E-21,...
        4.59411707689E-21, 5.78365073286E-21, 7.28118488054E-21,...
        9.16646867408E-21, 1.15399003502E-20, 1.45278738005E-20,...
        1.82895095067E-20, 2.30251282873E-20, 2.89869191107E-20,...
        3.64923690780E-20, 4.59411707689E-20, 5.78365073286E-20,...
        7.28118488054E-20, 9.16646867408E-20, 1.15399003502E-19,...
        1.45278738005E-19, 1.82895095067E-19, 2.30251282873E-19,...
        2.89869191107E-19, 3.64923690780E-19, 4.59411707689E-19,...
        5.78365073286E-19, 7.28118488054E-19, 9.16646867408E-19,...
        1.15399003502E-18, 1.45278738005E-18, 1.82895095067E-18,...
        2.30251282873E-18, 2.89869191107E-18, 3.64923690780E-18,...
        4.59411707689E-18, 5.78365073286E-18, 7.28118488054E-18,...
        9.16646867408E-18, 1.15399003502E-17, 1.45278738005E-17,...
        1.82895095067E-17, 2.30251282873E-17, 2.89869191107E-17,...
        3.64923690780E-17, 4.59411707689E-17, 5.78365073286E-17,...
        7.28118488054E-17, 9.16646867408E-17, 1.15399003502E-16,...
        1.45278738005E-16, 1.82895095067E-16, 2.30251282873E-16,...
        2.89869191107E-16, 3.64923690780E-16, 4.59411707689E-16,...
        5.78365073286E-16, 7.28118488054E-16, 9.16646867408E-16,...
        1.15399003502E-15, 1.45278738005E-15, 1.82895095067E-15,... %DATA J02 /						    
        2.30251282873E-15, 2.89869191107E-15, 3.64923690780E-15,...
        4.59411707689E-15, 5.78365073286E-15, 7.28118488054E-15,...
        9.16646867408E-15, 1.15399003502E-14, 1.45278738005E-14,...
        1.82895095067E-14, 2.30251282873E-14, 2.89869191107E-14,...
        3.64923690780E-14, 4.59411707689E-14, 5.78365073286E-14,...
        7.28118488054E-14, 9.16646867408E-14, 1.15399003502E-13,...
        1.45278738005E-13, 1.82895095067E-13, 2.30251282873E-13,...
        2.89869191107E-13, 3.64923690780E-13, 4.59411707689E-13,...
        5.78365073286E-13, 7.28118488054E-13, 9.16646867408E-13,...
        1.15399003502E-12, 1.45278738005E-12, 1.82895095067E-12,...
        2.30251282873E-12, 2.89869191107E-12, 3.64923690780E-12,...
        4.59411707689E-12, 5.78365073286E-12, 7.28118488054E-12,...
        9.16646867408E-12, 1.15399003502E-11, 1.45278738005E-11,...
        1.82895095067E-11, 2.30251282873E-11, 2.89869191107E-11,...
        3.64923690780E-11, 4.59411707689E-11, 5.78365073286E-11,...
        7.28118488054E-11, 9.16646867408E-11, 1.15399003502E-10,...
        1.45278738005E-10, 1.82895095067E-10, 2.30251282873E-10,...
        2.89869191107E-10, 3.64923690780E-10, 4.59411707689E-10,...
        5.78365073286E-10, 7.28118488054E-10, 9.16646867408E-10,... %     DATA J03 /						      
        1.15399003502E-09, 1.45278738005E-09, 1.82895095067E-09,...
        2.30251282873E-09, 2.89869191107E-09, 3.64923690780E-09,...
        4.59411707689E-09, 5.78365073286E-09, 7.28118488054E-09,...
        9.16646867408E-09, 1.15399003502E-08, 1.45278738005E-08,...
        1.82895095067E-08, 2.30251282873E-08, 2.89869191107E-08,...
        3.64923690780E-08, 4.59411707689E-08, 5.78365073286E-08,...
        7.28118488054E-08, 9.16646867408E-08, 1.15399003502E-07,...
        1.45278738005E-07, 1.82895095067E-07, 2.30251282873E-07,...
        2.89869191107E-07, 3.64923690780E-07, 4.59411707689E-07,...
        5.78365073285E-07, 7.28118488052E-07, 9.16646867404E-07,...
        1.15399003501E-06, 1.45278738003E-06, 1.82895095064E-06,...
        2.30251282867E-06, 2.89869191095E-06, 3.64923690758E-06,...
        4.59411707644E-06, 5.78365073195E-06, 7.28118487872E-06,...
        9.16646867045E-06, 1.15399003430E-05, 1.45278737860E-05,...
        1.82895094779E-05, 2.30251282297E-05, 2.89869189958E-05,...
        3.64923688489E-05, 4.59411703117E-05, 5.78365064162E-05,...
        7.28118469850E-05, 9.16646831086E-05, 1.15398996255E-04,...
        1.45278723545E-04, 1.82895066216E-04, 2.30251225306E-04,...
        2.89869076247E-04, 3.64923461604E-04, 4.59411250423E-04,... %DATA J04 /
        5.78364160920E-04, 7.28116667643E-04, 9.16643235217E-04,...
        1.15398278784E-03, 1.45277292006E-03, 1.82892209919E-03,...
        2.30245526274E-03, 2.89857705186E-03, 3.64900773583E-03,...
        4.59365982051E-03, 5.78273840567E-03, 7.27936457548E-03,...
        9.16283685839E-03, 1.15326542506E-02, 1.45134174958E-02,...
        1.82606689882E-02, 2.29675986572E-02, 2.88721706030E-02,...
        3.62635577835E-02, 4.54850275983E-02, 5.69277656746E-02,...
        7.10027008536E-02, 8.80685210410E-02, 1.08264277045E-01,...
        1.31175863070E-01, 1.55158221880E-01, 1.76192710015E-01,...
        1.85876995278E-01, 1.69365111703E-01, 1.03959768545E-01,...
       -3.10930553773E-02,-2.26627317132E-01,-3.63038011380E-01,...
       -2.05495704279E-01, 3.42146510984E-01, 3.05772953116E-01,...
       -5.00859313870E-01, 3.05625273560E-01,-1.35532781438E-01,...
        5.91884875621E-02,-2.64075211274E-02, 1.10378829555E-02,...
       -4.42794704514E-03, 1.76464110345E-03,-7.02632756477E-04,...
        2.79730481886E-04,-1.11363173935E-04, 4.43345073089E-05,...
       -1.76498871074E-05, 7.02654673076E-06,-2.79731864600E-06,...
        1.11363261176E-06,-4.43345128132E-07, 1.76498874541E-07,...
       -7.02654675496E-08, 2.79731864738E-08,-1.11363261185E-08,... %DATA J05 /
        4.43345128138E-09,-1.76498874542E-09, 7.02654675497E-10,...
       -2.79731864738E-10, 1.11363261185E-10,-4.43345128138E-11,...
        1.76498874542E-11,-7.02654675497E-12, 2.79731864738E-12,...
       -1.11363261185E-12, 4.43345128138E-13,-1.76498874542E-13,...
        7.02654675497E-14,-2.79731864738E-14, 1.11363261185E-14,...
       -4.43345128138E-15, 1.76498874542E-15,-7.02654675497E-16,...
        2.79731864738E-16,-1.11363261185E-16, 4.43345128138E-17,...
       -1.76498874542E-17, 7.02654675497E-18,-2.79731864738E-18,...
        1.11363261185E-18,-4.43345128138E-19, 1.76498874542E-19,...
       -7.02654675497E-20, 2.79731864738E-20,-1.11363261185E-20,...
        4.43345128138E-21,-1.76498874542E-21, 7.02654675497E-22,...
       -2.79731864738E-22, 1.11363261185E-22,-4.43345128138E-23,...
        1.76498874542E-23,-7.02654675497E-24, 2.79731864738E-24,...
       -1.11363261185E-24, 4.43345128138E-25,-1.76498874542E-25,...
        7.02654675497E-26,-2.79731864738E-26, 1.11363261185E-26,...
       -4.43345128138E-27, 1.76498874542E-27,-7.02654675497E-28,...
        2.79731864738E-28,-1.11363261185E-28, 4.43345128138E-29,...
       -1.76498874542E-29, 7.02654675497E-30,-2.79731864738E-30,...
        1.11363261185E-30,-4.43345128138E-31, 1.76498874542E-31,... %DATA J06 /
       -7.02654675497E-32, 2.79731864738E-32,-1.11363261185E-32,...
        4.43345128138E-33,-1.76498874542E-33, 7.02654675497E-34,...
       -2.79731864738E-34, 1.11363261185E-34,-4.43345128138E-35,...
        1.76498874542E-35,-7.02654675497E-36, 2.79731864738E-36,...
       -1.11363261185E-36...
    ];

J11 = [ ...
        1.15129254536E-41, 1.82467571768E-41, 2.89191612340E-41,...
        4.58337817714E-41, 7.26416487142E-41, 1.15129254536E-40,...
        1.82467571768E-40, 2.89191612340E-40, 4.58337817714E-40,...
        7.26416487142E-40, 1.15129254536E-39, 1.82467571768E-39,...
        2.89191612340E-39, 4.58337817714E-39, 7.26416487142E-39,...
        1.15129254536E-38, 1.82467571768E-38, 2.89191612340E-38,...
        4.58337817714E-38, 7.26416487142E-38, 1.15129254536E-37,...
        1.82467571768E-37, 2.89191612340E-37, 4.58337817714E-37,...
        7.26416487142E-37, 1.15129254536E-36, 1.82467571768E-36,...
        2.89191612340E-36, 4.58337817714E-36, 7.26416487142E-36,...
        1.15129254536E-35, 1.82467571768E-35, 2.89191612340E-35,...
        4.58337817714E-35, 7.26416487142E-35, 1.15129254536E-34,...
        1.82467571768E-34, 2.89191612340E-34, 4.58337817714E-34,...
        7.26416487142E-34, 1.15129254536E-33, 1.82467571768E-33,...
        2.89191612340E-33, 4.58337817714E-33, 7.26416487142E-33,...
        1.15129254536E-32, 1.82467571768E-32, 2.89191612340E-32,...
        4.58337817714E-32, 7.26416487142E-32, 1.15129254536E-31,...
        1.82467571768E-31, 2.89191612340E-31, 4.58337817714E-31,...
        7.26416487142E-31, 1.15129254536E-30, 1.82467571768E-30,...
        2.89191612340E-30, 4.58337817714E-30, 7.26416487142E-30,... %DATA J12 /
        1.15129254536E-29, 1.82467571768E-29, 2.89191612340E-29,...
        4.58337817714E-29, 7.26416487142E-29, 1.15129254536E-28,...
        1.82467571768E-28, 2.89191612340E-28, 4.58337817714E-28,...
        7.26416487142E-28, 1.15129254536E-27, 1.82467571768E-27,...
        2.89191612340E-27, 4.58337817714E-27, 7.26416487142E-27,...
        1.15129254536E-26, 1.82467571768E-26, 2.89191612340E-26,...
        4.58337817714E-26, 7.26416487142E-26, 1.15129254536E-25,...
        1.82467571768E-25, 2.89191612340E-25, 4.58337817714E-25,...
        7.26416487142E-25, 1.15129254536E-24, 1.82467571768E-24,...
        2.89191612340E-24, 4.58337817714E-24, 7.26416487142E-24,...
        1.15129254536E-23, 1.82467571768E-23, 2.89191612340E-23,...
        4.58337817714E-23, 7.26416487142E-23, 1.15129254536E-22,...
        1.82467571768E-22, 2.89191612340E-22, 4.58337817714E-22,...
        7.26416487142E-22, 1.15129254536E-21, 1.82467571768E-21,...
        2.89191612340E-21, 4.58337817714E-21, 7.26416487142E-21,...
        1.15129254536E-20, 1.82467571768E-20, 2.89191612340E-20,...
        4.58337817714E-20, 7.26416487142E-20, 1.15129254536E-19,...
        1.82467571768E-19, 2.89191612340E-19, 4.58337817714E-19,...
        7.26416487142E-19, 1.15129254536E-18, 1.82467571768E-18,... %DATA J13 /
        2.89191612340E-18, 4.58337817714E-18, 7.26416487142E-18,...
        1.15129254536E-17, 1.82467571768E-17, 2.89191612340E-17,...
        4.58337817714E-17, 7.26416487142E-17, 1.15129254536E-16,...
        1.82467571768E-16, 2.89191612340E-16, 4.58337817714E-16,...
        7.26416487142E-16, 1.15129254536E-15, 1.82467571768E-15,...
        2.89191612340E-15, 4.58337817714E-15, 7.26416487142E-15,...
        1.15129254536E-14, 1.82467571768E-14, 2.89191612340E-14,...
        4.58337817714E-14, 7.26416487142E-14, 1.15129254536E-13,...
        1.82467571768E-13, 2.89191612339E-13, 4.58337817714E-13,...
        7.26416487142E-13, 1.15129254536E-12, 1.82467571767E-12,...
        2.89191612339E-12, 4.58337817712E-12, 7.26416487136E-12,...
        1.15129254535E-11, 1.82467571764E-11, 2.89191612330E-11,...
        4.58337817691E-11, 7.26416487085E-11, 1.15129254522E-10,...
        1.82467571732E-10, 2.89191612249E-10, 4.58337817486E-10,...
        7.26416486570E-10, 1.15129254392E-09, 1.82467571407E-09,...
        2.89191611431E-09, 4.58337815435E-09, 7.26416481407E-09,...
        1.15129253098E-08, 1.82467568149E-08, 2.89191603267E-08,...
        4.58337794883E-08, 7.26416429900E-08, 1.15129240131E-07,...
        1.82467535651E-07, 2.89191521450E-07, 4.58337589830E-07,... %DATA J14 /
        7.26415913665E-07, 1.15129110751E-06, 1.82467209929E-06,...
        2.89190705118E-06, 4.58335534665E-06, 7.26410762971E-06,...
        1.15127814034E-05, 1.82463960074E-05, 2.89182523440E-05,...
        4.58315029698E-05, 7.26359140967E-05, 1.15114876631E-04,...
        1.82431390214E-04, 2.89100899610E-04, 4.58109550314E-04,...
        7.25844219152E-04, 1.14985263625E-03, 1.82106638692E-03,...
        2.88283662699E-03, 4.56062758539E-03, 7.20696761303E-03,...
        1.13697387883E-02, 1.78872970803E-02, 2.80213928499E-02,...
        4.35882514308E-02, 6.70661191541E-02, 1.01311462934E-01,...
        1.48660299844E-01, 2.07323169235E-01, 2.65301367358E-01,...
        2.85244616633E-01, 2.02544074354E-01,-5.73751544738E-02,...
       -3.71353843665E-01,-3.19637336841E-01, 4.69551216345E-01,...
       -2.76990986055E-02,-2.03919684145E-01, 1.81743287906E-01,...
       -1.01182304791E-01, 4.48568788435E-02,-1.80886005471E-02,...
        7.19959869393E-03,-2.86614354212E-03, 1.14102877377E-03,...
       -4.54251518698E-04, 1.80840773115E-04,-7.19940076407E-05,...
        2.86613306248E-05,-1.14102812376E-05, 4.54251477777E-06,...
       -1.80840770536E-06, 7.19940074787E-07,-2.86613306142E-07,...
        1.14102812351E-07,-4.54251477773E-08, 1.80840770536E-08,... %DATA J15 /