Савёлова Методы решения некорректных задач 2012
.pdfПри вычислении ПФ задаются параметры полуширины гаус-
совского пика, соответствующего каждому зерну, и размер класте-
ра, т.е. участка разбиения ПФ. На рис. 7.10 и 7.11 приведены ПФ с параметрами 8º – полуширина гауссовского пика, 8º – размер кла-
стера (рис. 7.10), 4 и 8º (рис. 7.11).
Рис. 7.10. ПФ с полушириной гауссовского пика 8º
Рис. 7.11. ПФ с полушириной гауссовского пика 4º
161
На рис. 7.10 ПФ получились более заглаженными, чем на рис. 7.11, где более детально прорисованы подробности текстуры материала.
При вычислении ФРО задаются параметры полуширины гаус-
совского пика, соответствующего каждому зерну, и длина ряда разложения. На рис. 7.12, 7.13 представлены сечения ФРО (посто-
янным был взят угол φ1 ) с параметрами Lmax = 28, полуширина – 8˚ (рис. 7.12) и Lmax= 34, полуширина равна 4˚ (рис. 7.13).
Рис. 7.12. ФРО с полушириной пика 8º
162
Рис. 7.13. ФРО с полушириной гауссовского пика 4º
На рис. 7.12, 7.13 можно видеть, что с уменьшением параметра полуширины гауссовского пика, насаживаемого на каждое зерно,
происходит большая детализация, ФРО становится более изрезан-
ной.
В данной главе рассмотрены вопросы точности проведения
EBSD-измерений ориентировок зерен поликристаллических мате-
риалов. На точность определения величин, характеризующих тек-
стуру материалов, влияют несколько факторов, а именно: выбор шага измерения ориентировок, количество измеряемых зерен и связанный с ним факт зависимости ориентировок, вопрос о выборе угла разориентации, разграничивающий соседние зерна, определе-
163
ние границ и формы зерен, вопрос усреднения неиндексируемых
областей.
Контрольные вопросы
1.Получение ориентаций отдельных зерен поликристаллов мето-
дами электронной микроскопии.
2.Перечислите основные погрешности EBSD-измерений.
3.Как влияют параметры EBSD-измерений на точность вычисле-
ния ФРО и ПФ?
4.Следствием чего является зависимость ориентаций при EBSD-
измерениях?
164
Глава 8. Робастные методы восстановления ФРО
по набору отдельных ориентировок
Функция распределения зерен по ориентациям широко исполь-
зуется в текстурном анализе для исследования физико-
механических свойств поликристаллических материалов (напри-
мер, [30, 58]). В последние годы автоматизация ЕВSD (Electron Backscattering Diffraction) систем позволяет получать большое (по-
рядка 107 ) количество измерений отдельных кристаллографиче-
ских ориентаций для изучения глобальной текстуры.
В данной главе рассмотрены аналитические и численные оценки точности восстановления плотности распределения ориентаций на
SO(3):
в различных метрических пространствах L2 , L1 , С;
зависимость точности ядерных методов от объема выборки,
ядра сглаживания, параметра ядра сглаживания;
точность проекционных методов приближения плотности рас-
пределения ориентаций по базису из обобщенных шаровых функ-
ций, в частности характеров представлений группы SO(3) [12], за-
висимость точности проекционных методов от параметров – объе-
ма выборки, количества членов ряда Фурье;
зависимость приближения от свойств искомой плотности рас-
пределения, в частности от аналога дисперсии распределения (в
текстурном анализе принят термин «острота текстуры»);
165
оценки точности для зависимых элементов выборки, влияние размеров зерен на точность вычислений.
Численные эксперименты проводятся на модельных примерах,
полученных путем статистического моделирования нормальных распределений (НР), удовлетворяющих центральной предельной теореме на группе SO(3) [5, 7].
Для оценки точности сходимости ядерных и проекционных ап-
проксимаций плотности распределения на SO(3) по выборке из от-
дельных ориентировок в L2 , L1 используют результаты, получен-
ные в евклидовом пространстве [4, 14, 18, 45], а для сглаживания ориентировок зерен в ядерных методах – ядра, указанные в работах
[14, 104, 105].
8.1. Вычисление статистических характеристик ядерных
оценок ФРО
Пусть плотность распределения ориентаций |
f (g) , g SO(3) , |
представляется в виде ряда по обобщенным шаровым функциям |
|
l |
|
f (g) ClmnTl mn g , |
(8.1) |
l 0 m, n l
где
166
|
|
|
|
Clmn 2l 1 |
|
|
|
f g |
|
|
|
|
|
dg , |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
Tl mn g |
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
SO(3) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
f (g) |
|
2 dg |
|
|
|
|
|
|
|
|
Clmn |
2 . |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
2l |
|
|
|
|
|||||||||||||||||||||||||||
|
|
SO(3) |
l 0 |
|
1 m,n l |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
По выборке |
|
из независимых ориентировок |
g1, , gN , отве- |
|||||||||||||||||||||||||||||||||
чающих |
плотности распределения |
|
|
|
|
f (g) , |
|
рассмотрим ядерную |
||||||||||||||||||||||||||||
оценку |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
π |
|
|
|
|
|
1 |
|
|
|
t |
|
|
|||||
fNk (g) |
|
q(gi 1g) , |
q(t) 0 , |
|
|
q(t) |
sin |
2 |
dt 1, |
(8.2) |
||||||||||||||||||||||||||
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
N i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
π |
|
|
|
|
|
π |
|
2 |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
sin l |
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
где q(t) Ql χl (t) , χl t |
|
|
|
|
|
|
|
2 |
|
|
|
|
, есть ядро сглаживания. |
|||||||||||||||||||||||
|
|
|
|
|
t |
|
|
|
|
|||||||||||||||||||||||||||
|
l 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
sin |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Оценка (8.2) может быть преобразована к виду |
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
k |
|
|
Ql |
|
|
|
|
l |
|
|
ˆ mn |
mn |
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
fN (g) |
|
2l |
|
|
|
|
|
|
|
|
Cl Tl |
|
(g) , |
(8.3) |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
l 0 |
1 m,n l |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
ˆ mn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
mn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где Cl |
– оценки коэффициентов Cl |
: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
ˆ mn |
|
|
2l 1 |
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
mn |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
Cl |
|
|
|
|
|
|
|
|
|
|
|
|
|
Tl |
|
(gi ). |
|
|
|
(8.4) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N i 1
Очевидно, получаем
167
|
Ql |
l |
|
fN (g) MfNk (g) |
ClmnTl mn (g) . (8.5) |
||
|
|||
l 0 |
2l 1 m,n l |
Из выражений (8.1), (8.5) находим смещение оценки (8.2):
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Q |
|
|
l |
|
|
|
|
|
sN (g) f (g) fN (g) 1 |
|
|
|
l |
|
|
ClmnTl mn (g) . (8.6) |
|||||||||||||||
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
l 0 |
|
|
2l 1 m,n l |
|
|
|
|
||||||||
Дисперсия оценки (8.2) или (8.3), (8.4) имеет вид |
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
Ql |
|
2 |
l |
|
|
|
|
|
|
2 |
|
|
||||
|
k |
|
|
|
|
|
|
ˆ l |
l |
|
|
|
||||||||||
DfN (g) |
|
|
|
|
|
|
DCmn |
Tmn (g) |
|
|
||||||||||||
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||
|
|
|
l 0 |
|
2l 1 |
|
m,n l |
|
|
|
|
|
|
|
|
|
||||||
|
Ql1 |
|
Ql2 |
|
|
l1 |
|
|
|
|
l2 |
|
|
|
|
ˆ l1 |
ˆ l2 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
2l 1 2l |
1 |
|
|
|
|
|
|
|
cov Cm1n1 |
, Cm2n2 |
|
||||||||||
l1 ,l2 0 |
m1 |
,n1 l1 |
m2 ,n2 l2 |
|
|
|
|
|
|
|
|
|||||||||||
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|||||||||||
l1 l2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T l1 |
|
(g)T l2 |
|
(g), |
|
|
|
|
(8.7) |
|||||||
|
|
|
|
|
|
|
m n |
|
|
|
m n |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
1 |
1 |
|
|
|
2 |
2 |
|
|
|
|
|
|
|
|
где
ˆ mn |
|
1 |
|
2l |
|
|
|
|
|
|
|
2 |
C(l, l, l , m, m, 2m)C(l, l, l , n, n, 2n) |
||||||||
DCl |
|
|
(2l 1) |
|
||||||
|
|
|||||||||
|
|
N |
|
l 0 |
|
|
|
|
|
|
|
|
|
|
|
2m2n 1 |
mn 2 |
|
|
||
|
|
|
Cl |
|
|
(Cl ) |
|
, |
||
|
|
|
2l 1 |
|||||||
|
|
|
|
|
|
|
|
|
168
ˆ m1n1 |
ˆ m2n2 |
|
1 |
|
|
l1 l2 |
1 |
|
|
||
|
|
|
|
|
|
|
|
||||
cov Cl1 |
, Cl2 |
|
N |
(2l1 1) (2l2 1) |
|
2l 1 |
|
||||
|
|
|
|
l |
l1 l2 |
|
|
||||
C(l1, l2 , l , m1, m2 , m1 m2 ) C(l1,l2 ,l , n1 |
, n2 , n1 n2 ) |
||||||||||
|
Clm1 m2 ,n1 n2 Clm1n1 |
Clm2n2 , |
|
|
|
||||||
|
|
|
|
1 |
2 |
|
|
|
|
|
|
C(l1,l2 ,l , m, n, m n) – коэффициенты Клебша–Гордана (приложе-
ние 4) [12].
Из формул (8.6), (8.7) получаем оценку точности ядерной оценки (8.2) в L2 :
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
Ql |
|
2 |
1 |
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
||||||||
M |
f (g) fN (g) |
|
|
1 |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
L2 |
l 0 |
|
|
2l 1 |
|
2l 1 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
l |
|
mn |
|
2 |
|
|
1 |
|
|
|
Ql |
|
2 |
l |
(2l 1) |
2 |
|
||||
|
|
|
|
|
|
||||||||||||||||
|
|
Cl |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
m,n l |
|
|
|
|
|
|
N l 0 |
|
2l 1 |
|
m,n l |
|
|
|
|
|
2l
C(l,l,l , m, m, 2m)C(l,l,l , n, n, 2n)Cl2m,2nl
l 0
|
1 |
2 |
1 |
|
|
|
|
|
Clmn |
|
|
. |
(8.8) |
2l 1 |
2l 1 |
В силу компактности группы SO(3) и неравенства Коши– Буняковского находим, что оценки в пространстве L2 справедливы
в пространстве L1 сверху:
f (g) L1 22π f (g) L2 .
169
8.2. Статистические характеристики для центральных
функций
Для центральных функций на SO(3) имеем разложение по ха-
рактерам представлений:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f t Cl χl t , |
|
|
|
(8.9) |
||||||
|
|
|
|
|
|
|
|
|
|
l 0 |
|
|
|
|
|
|
|
|
l |
|
1 |
π |
f t χl (t)sin |
2 |
t |
|
t |
|
|
θ |
|
ψ |
|
||
C |
|
|
|
|
|
|
|
dt , cos |
|
|
cos |
|
cos |
|
. |
||
|
|
|
|
|
|
|
|||||||||||
|
|
|
π |
π |
|
|
2 |
|
2 |
|
|
2 |
|
2 |
|
Ядерная оценка по выборке t1, , tN , отвечающей независимым вращениям g1, , gN , имеет вид
|
|
|
|
|
|
|
|
|
|
|
|
Ql |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
N |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
k |
t |
|
|
|
|
|
|
ˆ l |
χl t |
|
|
ˆ l |
|
|
χl (ti ) . |
|
|
|
|
|
|||||||||||||||||||
|
|
fN |
|
2l |
|
|
C |
, C |
|
|
|
|
|
|
|
(8.10) |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
l 0 |
|
1 |
|
|
|
|
|
|
|
|
|
|
N i 1 |
|
|
|
|
|
|
|
|
|
|
||||||||||||
Смещение в этом случае вычисляется по формуле |
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Q |
|
|
|
|
|
|
|
|
|
|
||||
|
|
sN (t) f (t) MfNk (t) 1 |
|
|
|
|
l |
Cl |
χl (t). |
|
(8.11) |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l 0 |
|
|
|
|
2l 1 |
|
|
|
|
|
|
|
|
|
|||||||
Из (8.10) находим дисперсию оценки |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ql |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
ˆ l |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
DfN (t) |
|
|
|
|
|
|
DC |
|
χl (t) |
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l 0 |
|
|
2l 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
Ql1 |
|
|
|
|
Ql2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ˆ l1 |
ˆ l2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
(2l |
1) (2l |
|
|
|
|
|
|
|
cov(C |
, C |
)χl1 (t)χl2 |
(t), |
|
|
(8.12) |
|||||||||||||||||||||||||
|
|
l1 ,l2 |
0 |
|
2 |
1) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
l1 l2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ˆ l |
|
1 |
|
|
2l |
l |
|
|
|
|
l |
|
|
2 |
|
|
|
|
|
|
ˆ l1 |
|
ˆ l2 |
|
|
|
|
|
|
1 |
|
|
l1 l2 |
l |
l1 |
|
l2 |
|
|||||
DC |
|
|
|
C |
|
|
C |
|
|
|
, |
|
cov(C |
,C |
) |
|
|
|
C |
|
C |
C |
|
. |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
N l 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
170 |
|
|
|
|
|
|
|
|
|
l |
l1 l2 |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|