Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Савёлова Методы решения некорректных задач 2012

.pdf
Скачиваний:
16
Добавлен:
12.11.2022
Размер:
6.26 Mб
Скачать

При вычислении ПФ задаются параметры полуширины гаус-

совского пика, соответствующего каждому зерну, и размер класте-

ра, т.е. участка разбиения ПФ. На рис. 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 2 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 (tl2

(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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]