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

книги / Теория сварочных процессов

..pdf
Скачиваний:
11
Добавлен:
12.11.2023
Размер:
10.02 Mб
Скачать

изображении, кодированные числом, меньшим k, на бинарном изображении станут белыми, в матрице – кодировка 0; все темные участки на реальном изображении, кодированные числом, большим k, на бинарном изображении станут черными, в матрице – кодировка 1;

2) выделить на изучаемом шлифе квадратный участок площадью 300×300 пикселей (в зависимости от увеличения и типа металлографического микроскопа это может соответствовать площади от 10 мкм2 до 1 мм2). Это значение площади в пикселях является оптимальным. Участки меньших размеров не желательны, поскольку могут охватывать только отдельный структурный фрагмент и не дать полного представления о периодичности и фрагментарности структуры в целом. Верхняя граница площади обрабатываемого фрактальным анализом участка шлифа определяется только возможностями процессора компьютера;

3) полученные изображения структур образцов сохранить в формате точечного рисунка bmp или jpg.

После этих подготовительных действий можно приступать собственно к проведению анализа, последовательные этапы которого задаются алгоритмом, разработанным на кафедре СП и ТКМ ПГТУ.

При фрактальном описании структур сварных швов алгоритм, реализованный программными средствами Mathcad 2001 Professional, MATLAB 6.1, заключался в следующем:

1.Дискретная аппроксимация исследуемых структур: разбиение соответствующих бинарных изображений квадратной сеткой, состоящей из одинаковых ячеек, и присвоение ячейкам, приходящимся на область перлита, значений 0, а на однофазную светлую область видманштеттова феррита – значений 1. Таким образом, полученное оцифрованное изображение структуры представляет собой матрицу, состоящую из нулей и единиц.

2.Обработка массивов цифровых значений (матриц), полученных

путем разбиения матриц на более крупные ячейки с размерами Li×Li

(i = 1…32).

3.Построение для каждого разбиения характеристической меры

ввиде вероятности распределения единиц Pi, необходимых для покрытия поля феррита.

4.Аппроксимация зависимостей ln(Pi) от ln(Li) методом наимень-

ших квадратов и определение фрактальных размерностей Df, из соот-

ношения

Pi = const(Li)D.

141

Данный алгоритм был реализован для анализа структурообразования металла швов опытных электродов серии 17 (по нумерации УИС), с покрытием на основе уральского ильменита. В качестве базового состава использовали электроды марки МР-3, с заменой ферромарганца на силикомарганец и добавкой ферротитана.

Для исследования микроструктуры использовали микроскоп МИМ-8, цифровую видеокамеру AverMedia (разрешение 640×420), транслирование изображения на дисплей компьютера осуществлялось системой Amcap.

Фрактальный анализ структур сварных швов проводился по основным зонам сварного соединения:

по центру шва;

зоне термовлияния (ЗТВ);

зоне сплавления;

основному металлу.

Изучали фотоизображения всех зон сварного соединения размером 440×440 пикселей, что соответствовало участку сварного шва 150×150 мкм (при увеличении ×150). Такой размер охватываемой области вполне корректен, так как толщина зон, в которых имеются структурные отличия может достигать 400 мкм (за исключением зоны сплавления). Дополнительно исследовали фотоизображения центра шва при большем увеличении (×450), размер изображения составлял 300×300 пикселей, что соответствовало участку образца размером 40×40 мкм. На рис. 4.18–4.22 представлен листинг программы бинаризации структуры центра шва опытных электродов серии 17.

Приведенный фрагмент матрицы М соответствует реальному изображению структуры, а фрагмент матрицы Mb – бинарному изображению. Полученное бинарное изображение обрабатывалось программой фрактального анализа, приведенной на рис. 4.23, листинг обработки приведен на рис. 4.24–4.27.

В результате реализации программного алгоритма получено значение фрактальной размерности структуры Df .

Результаты фрактального анализа структуры металла шва приведены в табл. 4.3, где представлены структуры различных зон наиболее характерных из сварных швов и соответствующие им фрактальные размерности.

Как видно из сопоставления значений фрактальной размерности Df с топографией структуры (см. табл. 4.3), более грубая структура имеет меньшее значение фрактальной размерности, а более гомогенная структура имеет большее значение фрактальной размерности (исключение составляет только зона сплавления, результаты анализа которой объяснены ниже).

142

Считывание изображения M READBMP ( "1793center3.bmp")

 

Число строк rows (M) = 640 и столбцов cols (M) = 480 матрицы М

 

 

 

 

 

 

 

 

 

 

 

 

Выделение квадратной матрицы из

 

 

 

 

Трансформирование матрицы в бинарный

 

считываемого массива (изображения)

 

 

 

массив (0; 479)

---> (0;1)

 

 

 

 

 

 

 

 

 

M

 

submatrix( M ,i1,i2,j1, j2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

150

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

submatrix ( M ,40,370,170,500)

 

 

 

 

 

 

 

 

c

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mb

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

for

i 0.. rows ( M )

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

for j 0.. cols ( M )

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dj ,i

0 if Mj ,i <k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

if Mj ,i

k

 

 

 

 

WRITEBMP ("1793center3_bin.bmp"

)

Mb .255

 

 

 

 

 

c

c

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

6

7

8

9

10

11

12

13

14

15

 

 

 

.

 

 

5 6

7

8

9

10

11

12

 

13

14

15

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

0

255

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

20

17

18

19

18

22

13

22

23

24

 

 

 

0

0

0

0

0

0

0

0

 

 

0

0

0

 

 

 

 

1

20

20

20

26

24

21

21

19

28

24

46

 

1

 

0

0

0

0

0

0

0

0

 

 

0

0

0

 

 

 

 

2

20

19

19

21

16

13

17

30

41

51

51

 

2

 

0

0

0

0

0

0

0

0

 

 

0

0

0

 

 

 

 

3

20

24

23

12

15

23

28

60

57

53

51

 

3

 

0

0

0

0

0

0

0

0

 

 

0

0

0

 

 

 

 

4

20

19

19

5

17

37

42

79

60

40

33

 

4

 

0

0

0

0

0

0

0

0

 

 

0

0

0

 

 

 

 

5

20

24

18

24

31

53

67

97

82

57

53

 

5

 

0

0

0

0

0

0

0

0

 

 

0

0

0

 

 

M =

6

20

19

14

63

57

70

84

89

83

63

62

Mb =

6

 

0

0

0

0

0

0

0

0

 

 

0

0

0

 

 

 

 

7

20

29

38

138

128

125

123

88

80

61

61

 

7

 

0

0

0

0

0

0

0

0

 

 

0

0

0

 

 

 

 

8

20

44

58

150

137

133

126

91

74

52

51

 

8

 

0

0

0

1

0

0

0

0

 

 

0

0

0

 

 

 

 

9

15

146

161

85

63

128

111

115

109

63

79

 

9

 

0

0

1

0

0

0

0

0

 

 

0

0

0

 

 

 

 

10

67

184

174

58

68

119

101

158

113

81

111

 

10

0

1

1

0

0

0

0

1

 

 

0

0

0

 

 

 

 

11

150

137

116

113

148

79

64

202

142

126

146

 

11

1

0

0

0

0

0

0

1

 

 

0

0

0

 

 

 

 

12

141

107

92

141

182

58

41

192

140

140

137

 

12

0

0

0

0

1

0

0

1

 

 

0

0

0

 

 

 

 

13

114

136

100

126

163

51

38

169

110

104

95

 

13

0

0

0

0

1

0

0

1

 

 

0

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

б

Рис. 4.18. Листинг программы бинаризации структуры центра сварного шва электрода 1793 (вверху); фрагменты матриц изображения структуры (в центре); структуры центра шва (внизу): реальная (а) и бинарная (б). ×150

143

Считывание изображения M READBMP ( "1793ztv.bmp")

Число строк rows (M) = 640 и столбцов cols (M) = 480 матрицы М

 

 

 

 

 

 

Выделение квадратной матрицы из

 

Трансформирование матрицы в бинарный

считываемого массива (изображения)

массив (0; 479)

---> (0;1)

 

 

 

 

M

 

submatrix( M ,i1,i2,j1, j2)

 

 

 

 

 

 

 

k

 

150

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

submatrix ( M ,40,370,170,500)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mb

c

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

for

i 0.. rows ( M )

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

for j 0.. cols ( M )

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dj ,i

0 if Mj ,i <k

 

 

 

 

 

 

 

 

 

1 if Mj ,i

k

WRITEBMP ("1793ztv_bin.bmp"

)

Mb .255

 

c

c

1

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

а

б

Рис. 4.19. Листинг программы бинаризации структуры ЗТВ сварного шва электрода 1793 (вверху); матрицы изображения структуры (в центре); структуры ЗТВ (внизу): реальная (а) и бинарная (б). ×150

144

Считывание изображения M READBMP ( "1793line_spl.bmp")

Число строк rows (M) = 640 и столбцов cols (M) = 480 матрицы М

Выделение квадратной матрицы из

 

 

Трансформирование матрицы в бинарный

 

 

 

считываемого массива (изображения)

 

 

массив (0; 479)

---> (0;1)

 

 

 

 

 

 

 

 

 

M

 

 

 

submatrix( M ,i1,i2,j1,

j2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

150

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

submatrix ( M ,40,370,170,500)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mb

 

 

 

c

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

for i 0.. rows ( M )

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

for

j 0.. cols ( M )

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dj ,i

0 if Mj ,i <k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

if Mj ,i

k

 

 

 

 

WRITEBMP ("1793line_spl_bin.bmp" )

 

Mb .255

 

 

 

 

 

D

 

c

c

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

2

3

4

5

6

7

8

9

10

 

 

 

 

0

1

2

3

4

5

6

 

7

8

9

10

 

 

0

74

66

67

73

66

58

59

54

64

74

63

 

0

0

 

 

0

 

0

0

0

0

0

 

0

 

0

0

0

 

 

 

1

56

49

68

95

73

55

75

71

67

72

58

 

1

0

 

 

0

 

0

0

0

0

0

 

0

 

0

0

0

 

 

 

2

52

42

70

101

64

53

98

92

59

76

67

 

2

0

 

 

0

 

0

0

0

0

0

 

0

 

0

0

0

 

 

 

3

80

60

87

113

66

70

127

99

68

89

86

 

3

0

 

 

0

 

0

1

0

0

1

 

0

 

0

0

0

 

 

 

4

97

66

89

116

71

76

118

76

86

92

78

 

4

0

 

 

0

 

0

1

0

0

1

 

0

 

0

0

0

 

 

 

5

104 57

75

94

57

60

87

64

77

71

56

 

5

0

 

 

0

 

0

0

0

0

0

 

0

 

0

0

0

 

 

M =

6

96

56

79

82

46

53

76

77

68

60

48

Mb =

6

0

 

 

0

 

0

0

0

0

0

 

0

 

0

0

0

 

 

 

7

71

54

91

80

43

54

66

77

87

66

52

 

7

0

 

 

0

 

0

0

0

0

0

 

0

 

0

0

0

 

 

 

8

63

62

80

73

42

48

72

66

88

66

55

 

8

0

 

 

0

 

0

0

0

0

0

 

0

 

0

0

0

 

 

 

9

57

54

79

77

45

57

92

90

89

64

48

 

9

0

 

 

0

 

0

0

0

0

0

 

0

 

0

0

0

 

 

 

10

59

50

71

72

46

64

98

94

89

66

55

 

10

0

 

 

0

 

0

0

0

0

0

 

0

 

0

0

0

 

 

 

11

70

52

53

54

50

71

89

83

72

64

67

 

11

0

 

 

0

 

0

0

0

0

0

 

0

 

0

0

0

 

 

 

12

67

57

43

42

58

73

78

82

62

67

80

 

12

0

 

 

0

 

0

0

0

0

0

 

0

 

0

0

0

 

 

 

13

55

60

49

43

58

62

62

80

69

75

84

 

13

0

 

 

0

 

0

0

0

0

0

 

0

 

0

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЗТВ

К центру шва

 

а

б

Рис. 4.20. Листинг программы бинаризации структуры зоны сплавления сварного шва электрода 1793 (вверху); матрицы изображения структуры (в центре); структуры зоны сплавления (внизу): реальная (а) и бинарная (б). ×150

145

Считывание изображения M READBMP ( "1793osn.bmp" )

Число строк rows (M) = 640 и столбцов

Выделение квадратной матрицы из считываемого массива (изображения)

M submatrix( M ,i1,i2,j1, j2)

M submatrix ( M ,40,370,170,500)

WRITEBMP ("1793osn_bin.bmp" ) Mb

cols (M) = 480 матрицы М

Трансформирование матрицы в бинарный

массив (0; 479)

---> (0;1)

 

 

 

 

k

 

150

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mb

c

0

 

 

 

 

 

 

 

 

 

for

i 0.. rows ( M )

 

1

 

 

 

 

 

 

 

for j 0.. cols ( M )

 

 

1

 

 

 

 

 

 

Dj ,i

0

if Mj ,i <k

 

 

 

 

 

1

if Mj ,i

k

.255

 

D

c

c

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

6

7

8

9

10

11

12

13

14

15

 

 

 

 

0

1

2

3

4

5

6

7

8

9

10

 

 

 

 

0

39

63

38

45

39

39

34

92

133

124

37

 

0

 

0

0

0

0

0

0

0

0

0

0

 

0

 

 

 

 

1

55

44

57

37

56

35

35

129

91

54

70

 

1

 

0

0

0

0

0

0

0

0

0

0

 

0

 

 

 

 

2

47

44

49

46

53

44

39

120

86

78

102

 

2

 

0

0

0

0

0

0

0

0

0

0

 

0

 

 

 

 

3

51

57

45

42

18

106

89

38

54

170

149

 

3

 

1

1

0

0

0

0

0

0

0

0

 

0

 

 

 

 

4

86

69

45

42

0

148

139

24

84

178

170

 

4

 

1

1

0

0

0

0

0

0

0

0

 

1

 

 

 

 

5

60

155

129

32

59

134

122

104

155

129

76

 

5

 

0

0

0

0

0

0

1

1

0

0

 

1

 

 

 

M =

6

23

173

164

20

78

97

81

103

137

86

71

Mb =

6

 

0

0

0

0

0

0

1

1

0

0

 

0

 

 

 

 

7

23

82

66

39

44

42

41

42

52

61

97

 

7

 

0

0

0

0

0

0

0

0

0

0

 

0

 

 

 

 

8

51

44

38

39

40

54

53

28

29

85

96

 

8

 

1

1

0

0

0

0

0

0

0

0

 

0

 

 

 

 

9

48

42

31

54

106

116

79

45

61

128

86

 

9

 

0

0

0

0

0

0

0

0

0

0

 

1

 

 

 

 

10

46

54

36

92

145

109

68

50

73

117

55

 

10

0

0

0

0

0

0

0

0

0

1

 

0

 

 

 

 

11

55

47

34

76

103

54

59

42

51

94

47

 

11

0

0

0

0

0

0

0

0

0

0

 

0

 

 

 

 

12

78

31

53

58

52

44

50

44

45

72

84

 

12

0

0

0

0

0

0

0

0

0

0

 

0

 

 

 

 

13

122

56

45

55

68

62

32

80

119

59

32

 

13

0

0

0

0

0

1

0

0

0

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а б

Рис. 4.21. Листинг программы бинаризации структуры основного металла (вверху); матрицы изображения структуры (в центре); структуры основного металла (внизу):

реальная (а) и бинарная (б). ×150

146

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

2

3

4

5

6

7

8

9

10

 

 

 

 

0

 

1

 

2

3

4

5

 

6

 

7

 

8

 

9

10

 

 

0

 

21

20

20

21

21

22

15

15

29

20

21

 

 

0

0

 

0

 

0

 

0

0

0

0

 

0

 

0

 

0

 

0

 

 

1

 

21

20

20

21

21

21

18

18

23

15

27

 

1

0

0

0

 

0

0

0

0

0

0

0

 

0

 

 

2

 

21

21

21

21

21

21

18

15

25

19

25

 

2

0

0

0

 

0

0

0

0

0

0

0

 

0

 

 

3

 

21

21

21

21

21

21

27

18

16

7

13

 

3

0

0

0

 

0

0

0

0

0

0

0

 

0

 

 

4

 

21

21

21

21

21

21

20

19

4

7

39

 

4

0

0

0

 

0

0

0

0

0

0

0

 

0

 

 

5

 

21

21

21

21

21

21

19

22

29

58

113

 

5

0

0

0

 

0

0

0

0

0

0

0

 

0

 

M =

6

 

21

22

22

21

21

21

21

22

98

119

115

Mb =

6

0

0

0

 

0

0

0

0

0

0

0

 

0

 

 

7

 

21

22

22

21

21

20

37

50

232

219

88

 

7

0

0

0

 

0

0

0

0

0

1

1

 

0

 

 

8

 

23

15

17

2

28

11

40

208

239

133

156

 

8

0

0

0

 

0

0

0

0

1

1

0

 

0

 

 

9

 

20

15

25

20

76

79

133

142

111

142

228

 

9

0

0

0

 

0

0

0

0

0

0

0

 

1

 

 

 

10

22

16

26

38

121

113

101

70

76

168

208

 

 

0

0

0

 

0

0

0

0

0

0

0

 

1

 

 

 

11

27

0

56

133

202

136

111

98

94

120

165

 

 

0

0

0

 

0

1

0

0

0

0

0

 

0

 

 

 

12

35

4

136

228

231

151

154

142

85

56

163

 

 

0

0

0

 

1

1

0

0

0

0

0

 

0

 

 

 

13

34

52

229

241

184

163

168

149

115

123

203

 

 

0

0

1

 

1

1

0

0

0

0

0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

б

Рис. 4.22. Листинг программы бинаризации структуры центра сварного шва электрода 1807 (вверху); матрицы изображения структуры (в центре); структуры центра шва (внизу): реальная (а) и бинарная (б). ×150

147

clear a

 

end;%for

clear N

 

if cnt>0;

clear NL

 

break;

clear EL

 

end;

a = imread('1802center_bin.bmp','bmp');

 

end;%for

B1=uint8(a);

 

if cnt>0

SS = B1 ~= 0;

 

N(L)=N(L)+1;

figure, imshow(SS, 'notruesize');

 

end;%if

%pause;

 

end;%for

p=430;

 

end;%for

S(1:p,1:p)=1;

 

end;

for i=1:p;

 

%

for j=1:p;

 

%

S(i,j)=SS(i,j);

 

for L=1:Lmax;

end;

 

EL(L)=log(L);

end;

 

NL(L)=log(N(L));

Lmax=round(p/16);

 

end;%for

N(1:Lmax)=0.0;

 

pP=polyfit(EL(1:Lmax),NL(1:Lma

NL(1:Lmax)=0.0;

 

x),1);

EL(1:Lmax)=0.0;

 

ff=polyval(pP,EL);

NN=sum(S);

 

D=pP(1)

NMN=sum(NN);

 

axis auto;

N(1)=p^2-sum(NMN)

 

figure,

for L=2:Lmax;

 

plot(EL(1:Lmax),NL(1:Lmax));

B=floor(p/L);

 

title(['D = ',num2str(D),'']);

for i=1:B;

 

xlabel(['D = ',num2str(D),'']);

for j=1:B;

 

l=1

cnt=0.0;

 

DD=N(1)

xb=L.*(i-1)+1;

 

l=2

xe=L.*i;

 

DD=N(2)

yb=L.*(j-1)+1;

 

l=4

ye=L.*j;

 

DD=N(4)

for x=xb:xe;

 

l=8

for y=yb:ye;

 

DD=N(8)

if S(x,y)==0

 

 

cnt=cnt+1;

 

 

end;%if

 

 

if cnt>0;

 

 

 

 

 

Рис. 4.23. Реализация алгоритма программы фрактального анализа структуры металла сварного шва на языке Pascal

148

Рис. 4.24. Листинг программы фрактального анализа бинарного изображения структуры центра шва электрода 1793

149

Рис. 4.25. Листинг программы фрактального анализа бинарного изображения структуры центра шва электрода 1807

150