-
Notifications
You must be signed in to change notification settings - Fork 4
/
getHankelFC.m
144 lines (143 loc) · 9.33 KB
/
getHankelFC.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
function [fc, nc, nc0] = getHankelFC(name)
% TODO: Update help text.
%
% [FC,NC,NC0]=HANKELFC(ORDER)
% Filter coefficients for Hankel transform
% FC(NC0) refers to zero argument
% NC number of coefficients
% ORDER=1: NY=+0.5 (SIN)
% ORDER=2: NY=-0.5 (COS)
% ORDER=3: NY=0.0 (J0)
% ORDER=4: NY=1.0 (J1)
% 10 data points per decade
switch name
case 'sin'
fc = [...
2.59526236E-07 3.66544843E-07 5.17830795E-07 7.31340622E-07 ...
1.03322805E-06 1.45918500E-06 2.06161065E-06 2.91137793E-06 ...
4.11357863E-06 5.80876420E-06 8.20798075E-06 1.15895083E-05 ...
1.63778560E-05 2.31228459E-05 3.26800649E-05 4.61329334E-05 ...
6.52101085E-05 9.20390575E-05 1.30122935E-04 1.83620431E-04 ...
2.59656626E-04 3.66311982E-04 5.18141184E-04 7.30717340E-04 ...
1.03392184E-03 1.45742714E-03 2.06292302E-03 2.90599911E-03 ...
4.11471902E-03 5.79042763E-03 8.20004722E-03 1.15192930E-02 ...
1.63039133E-02 2.28257757E-02 3.22249222E-02 4.47864328E-02 ...
6.27329625E-02 8.57059100E-02 1.17418314E-01 1.53632655E-01 ...
1.97717964E-01 2.28849849E-01 2.40311038E-01 1.65409220E-01 ...
2.84701476E-03 -2.88016057E-01 -3.69097406E-01 -2.50107514E-02 ...
5.71811256E-01 -3.92261572E-01 7.63280044E-02 5.16233994E-02 ...
-6.48012082E-02 4.89047141E-02 -3.26936331E-02 2.10539842E-02 ...
-1.33862549E-02 8.47124695E-03 -5.35123972E-03 3.37796651E-03 ...
-2.13174466E-03 1.34513833E-03 -8.48749612E-04 5.35531006E-04 ...
-3.37898780E-04 2.13200109E-04 -1.34520273E-04 8.48765787E-05 ...
-5.35535069E-05 3.37899801E-05 -2.13200365E-05 1.34520337E-05 ...
-8.48765949E-06 5.35535110E-06 -3.37899811E-06 2.13200368E-06 ...
-1.34520338E-06 8.48765951E-07 -5.35535110E-07 3.37899811E-07].';
nc = 80;
nc0 = 40;
case 'cos'
fc = [...
1.63740363E-07 1.83719709E-07 2.06136904E-07 2.31289411E-07 ...
2.59510987E-07 2.91176117E-07 3.26704977E-07 3.66569013E-07 ...
4.11297197E-07 4.61483045E-07 5.17792493E-07 5.80972733E-07 ...
6.51862128E-07 7.31401337E-07 8.20645798E-07 9.20779729E-07 ...
1.03313185E-06 1.15919300E-06 1.30063594E-06 1.45933752E-06 ...
1.63740363E-06 1.83719709E-06 2.06136904E-06 2.31289411E-06 ...
2.59510987E-06 2.91176117E-06 3.26704977E-06 3.66569013E-06 ...
4.11297197E-06 4.61483045E-06 5.17792493E-06 5.80972733E-06 ...
6.51862128E-06 7.31401337E-06 8.20645798E-06 9.20779729E-06 ...
1.03313185E-05 1.15919300E-05 1.30063594E-05 1.45933752E-05 ...
1.63740363E-05 1.83719709E-05 2.06136904E-05 2.31289411E-05 ...
2.59510987E-05 2.91176117E-05 3.26704977E-05 3.66569013E-05 ...
4.11297197E-05 4.61483045E-05 5.17792493E-05 5.80972733E-05 ...
6.51862128E-05 7.31401337E-05 8.20645798E-05 9.20779729E-05 ...
1.03313185E-04 1.15919300E-04 1.30063594E-04 1.45933752E-04 ...
1.63740363E-04 1.83719709E-04 2.06136904E-04 2.31289411E-04 ...
2.59510987E-04 2.91176117E-04 3.26704976E-04 3.66569013E-04 ...
4.11297197E-04 4.61483045E-04 5.17792493E-04 5.80972733E-04 ...
6.51862127E-04 7.31401337E-04 8.20645797E-04 9.20779730E-04 ...
1.03313185E-03 1.15919300E-03 1.30063593E-03 1.45933753E-03 ...
1.63740362E-03 1.83719710E-03 2.06136901E-03 2.31289411E-03 ...
2.59510977E-03 2.91176115E-03 3.26704948E-03 3.66569003E-03 ...
4.11297114E-03 4.61483003E-03 5.17792252E-03 5.80972566E-03 ...
6.51861416E-03 7.31400728E-03 8.20643673E-03 9.20777603E-03 ...
1.03312545E-02 1.15918577E-02 1.30061650E-02 1.45931339E-02 ...
1.63734419E-02 1.83711757E-02 2.06118614E-02 2.31263461E-02 ...
2.59454421E-02 2.91092045E-02 3.26529302E-02 3.66298115E-02 ...
4.10749753E-02 4.60613861E-02 5.16081994E-02 5.78193646E-02 ...
6.46507780E-02 7.22544422E-02 8.03873578E-02 8.92661837E-02 ...
9.80670729E-02 1.07049506E-01 1.13757572E-01 1.18327217E-01 ...
1.13965041E-01 1.00497783E-01 6.12958082E-02 -1.61234222E-04 ...
-1.11788551E-01 -2.27536948E-01 -3.39004453E-01 -2.25128800E-01 ...
8.98279919E-02 5.12510388E-01 -1.31991937E-01 -3.35136479E-01 ...
3.64868100E-01 -2.34039961E-01 1.32085237E-01 -7.56739672E-02 ...
4.52296662E-02 -2.78297002E-02 1.73727753E-02 -1.09136894E-02 ...
6.87397283E-03 -4.33413470E-03 2.73388730E-03 -1.72477355E-03 ...
1.08821012E-03 -6.86602007E-04 4.33213523E-04 -2.73338487E-04 ...
1.72464733E-04 -1.08817842E-04 6.86594042E-05 -4.33211523E-05 ...
2.73337984E-05 -1.72464607E-05 1.08817810E-05 -6.86593962E-06 ...
4.33211503E-06 -2.73337979E-06 1.72464606E-06 -1.08817810E-06 ...
6.86593961E-07 -4.33211503E-07 2.73337979E-07 -1.72464606E-07].';
nc = 164;
nc0 = 122;
case 'j0'
fc = [...
2.89878288E-07 3.64935144E-07 4.59426126E-07 5.78383226E-07 ...
7.28141338E-07 9.16675639E-07 1.15402625E-06 1.45283298E-06 ...
1.82900834E-06 2.30258511E-06 2.89878286E-06 3.64935148E-06 ...
4.59426119E-06 5.78383236E-06 7.28141322E-06 9.16675664E-06 ...
1.15402621E-05 1.45283305E-05 1.82900824E-05 2.30258527E-05 ...
2.89878259E-05 3.64935186E-05 4.59426051E-05 5.78383329E-05 ...
7.28141144E-05 9.16675882E-05 1.15402573E-04 1.45283354E-04 ...
1.82900694E-04 2.30258630E-04 2.89877891E-04 3.64935362E-04 ...
4.59424960E-04 5.78383437E-04 7.28137738E-04 9.16674828E-04 ...
1.15401453E-03 1.45282561E-03 1.82896826E-03 2.30254535E-03 ...
2.89863979E-03 3.64916703E-03 4.59373308E-03 5.78303238E-03 ...
7.27941497E-03 9.16340705E-03 1.15325691E-02 1.45145832E-02 ...
1.82601199E-02 2.29701042E-02 2.88702619E-02 3.62691810E-02 ...
4.54794031E-02 5.69408192E-02 7.09873072E-02 8.80995426E-02 ...
1.08223889E-01 1.31250483E-01 1.55055715E-01 1.76371506E-01 ...
1.85627738E-01 1.69778044E-01 1.03405245E-01 -3.02583233E-02 ...
-2.27574393E-01 -3.62173217E-01 -2.05500446E-01 3.37394873E-01 ...
3.17689897E-01 -5.13762160E-01 3.09130264E-01 -1.26757592E-01 ...
4.61967890E-02 -1.80968674E-02 8.35426050E-03 -4.47368304E-03 ...
2.61974783E-03 -1.60171357E-03 9.97717882E-04 -6.26275815E-04 ...
3.94338818E-04 -2.48606354E-04 1.56808604E-04 -9.89266288E-05 ...
6.24152398E-05 -3.93805393E-05 2.48472358E-05 -1.56774945E-05 ...
9.89181741E-06 -6.24131160E-06 3.93800058E-06 -2.48471018E-06 ...
1.56774609E-06 -9.89180896E-07 6.24130948E-07 -3.93800005E-07 ...
2.48471005E-07 -1.56774605E-07 9.89180888E-08 -6.24130946E-08].';
nc = 100;
nc0 = 60;
case 'j1'
fc = [...
1.84909557E-13 2.85321327E-13 4.64471808E-13 7.16694771E-13 ...
1.16670043E-12 1.80025587E-12 2.93061898E-12 4.52203829E-12 ...
7.36138206E-12 1.13588466E-11 1.84909557E-11 2.85321327E-11 ...
4.64471808E-11 7.16694771E-11 1.16670043E-10 1.80025587E-10 ...
2.93061898E-10 4.52203829E-10 7.36138206E-10 1.13588466E-09 ...
1.84909557E-09 2.85321326E-09 4.64471806E-09 7.16694765E-09 ...
1.16670042E-08 1.80025583E-08 2.93061889E-08 4.52203807E-08 ...
7.36138149E-08 1.13588452E-07 1.84909521E-07 2.85321237E-07 ...
4.64471580E-07 7.16694198E-07 1.16669899E-06 1.80025226E-06 ...
2.93060990E-06 4.52201549E-06 7.36132477E-06 1.13587027E-05 ...
1.84905942E-05 2.85312247E-05 4.64449000E-05 7.16637480E-05 ...
1.16655653E-04 1.79989440E-04 2.92971106E-04 4.51975783E-04 ...
7.35565435E-04 1.13444615E-03 1.84548306E-03 2.84414257E-03 ...
4.62194743E-03 7.10980590E-03 1.15236911E-02 1.76434485E-02 ...
2.84076233E-02 4.29770596E-02 6.80332569E-02 9.97845929E-02 ...
1.51070544E-01 2.03540581E-01 2.71235377E-01 2.76073871E-01 ...
2.16691977E-01 -7.83723737E-02 -3.40675627E-01 -3.60693673E-01 ...
5.13024526E-01 -5.94724729E-02 -1.95117123E-01 1.99235600E-01 ...
-1.38521553E-01 8.79320859E-02 -5.50697146E-02 3.45637848E-02 ...
-2.17527180E-02 1.37100291E-02 -8.64656417E-03 5.45462758E-03 ...
-3.44138864E-03 2.17130686E-03 -1.36998628E-03 8.64398952E-04 ...
-5.45397874E-04 3.44122545E-04 -2.17126585E-04 1.36997597E-04 ...
-8.64396364E-05 5.45397224E-05 -3.44122382E-05 2.17126544E-05 ...
-1.36997587E-05 8.64396338E-06 -5.45397218E-06 3.44122380E-06 ...
-2.17126543E-06 1.36997587E-06 -8.64396337E-07 5.45397218E-07].';
nc = 100;
nc0 = 60;
otherwise
error('Unrecognized name.');
end