Chương 1: “Cơ sở phép biến đổi Wavelet rời rạc” trình bày cơ sở phép biến đổi Fourier, các phép biến đổi Wavelet và ứng dụng của các phép biến đổi đó trong việc phân tích tín hiệu số. Luận văn trình bày cụ thể cơ sở toán học của phép biến đổi Wavele rời rạc, phép biến đổi Wavelet Packet, Harmonic Wavelet Packet đƣợc sử dụng để xử lý tín hiệu, kết quả của quá trình phân tích này là đầu vào cho mạng nơron. Chương 2: “Mạng Nơron” trình bày cơ sở về mạng nơron, cấu tạo, cấu trúc mạng nơron, các phƣơng pháp học của mạng nơron. Mạng nơron đƣợc sử dụng phổ biến là mạng nơron truyền thẳng đa lớp MLP, với thuật toán học là thuật toán lan truyền ngƣợc.
Bên cạnh đó trình bày mạng Wavelet neural networks, mạng này là sự kết hợp giữa mạng nơron truyền thẳng MLP với đầu vào là giá trị của quá trình phân tích tín hiệu bằng Wavelet packet, mạng đƣợc ứng dụng trong chẩn đoán vết nứt, hƣ hỏng bánh răng. Chương 3: “Phát hiện vết nứt bánh răng nhờ Wavelet neural networks” trình bày các dạng hƣ hỏng chủ yếu của bánh răng, mô hình thí nghiệm đo tín hiệu dao động của hộp số bánh răng với các dạng hỏng đƣợc định trƣớc. Chƣơng này cũng trình bày kết quả phân tích tín hiệu trong miền thời gian, miền tần số, cho thấy đƣợc những hạn chế của các phép phân tích tín hiệu trên trong nhận dạng hƣ hỏng bánh răng, từ đó đƣa ra chẩn đoán tình trạng bánh răng bằng Wavelet neural networks. Với việc trình bày trình tự các bƣớc xử lý, phân tích tín hiệu, nhận dạng bằng Wavelet neural networks, kết quả thu đƣợc là tin cậy, đúng với các sai hỏng định trƣớc có thể thấy đƣợc ƣu điểm và khả năng ứng dụng của phƣơng pháp trong thực tiễn sản xuất.
Luận văn đƣợc thực hiện tại Bộ môn Cơ học ứng dụng- Viện Cơ Khí- Đại học Bách Khoa Hà Nội dƣới sự hƣớng dẫn tận tình của PGS.TS Nguyễn Phong Điền. 12 download by : skknchat@gmail.com Tác giả xin cảm ơn thầy đã hƣớng dẫn, cung cấp những kiến thức quan trọng, quý giá giúp tác giả có thể hoàn thành luận văn này. Hà Nội, ngày tháng 11 năm 2014 Học viên thực hiện: Hà Trung Kiên 13 download by : skknchat@gmail.com CHƢƠNG I CƠ SỞ PHÉP BIẾN ĐỔI WAVELET RỜI RẠC Các phép biến đổi tín hiệu với mục đích tái hiện tín hiệu theo một dạng khác mà không làm thay đổi các thông tin tín hiệu đó ẩn chứa. Cùng một tín hiệu, mỗi cách biến đổi khác nhau sẽ đƣa ra những cách nhìn khác nhau về tín hiệu đó.
Qua đó, những thông tin ẩn chứa trong tín hiệu sẽ đƣợc khai thác và phản ánh hiện trạng của đối tƣợng sinh ra tín hiệu. Có nhiều phƣơng pháp biểu diễn tín hiệu: biểu diễn trong miền thời gian, trong miền tần số và trong miền thời gian- tần số. Ƣu điểm của việc mô tả tín hiệu trong miền thời gian là tính toán tƣơng đối đơn giản, có thể xác định đƣợc các thời điểm xảy ra dao động. Tuy nhiên việc mô tả này lại có nhƣợc điểm là khó đoán biết tần số và khó chẩn đoán.
Việc mô tả tín hiệu trong miền tần số (phân tích phổ) cho phép nhận dạng tần số của tín hiệu nhƣng lại làm mất thông tin về thời gian (ví dụ nhƣ thời điểm xảy ra một hiện tƣợng, khoảng thời gian rung động…). Nhƣ vậy, việc mô tả tín hiệu riêng rẽ trong miền thời gian và trong miền tần số đều có những hạn chế nhất định. Để khắc phục những hạn chế trên, việc mô tả tín hiệu trong miền thời gian - tần số (Time - Frequency Analysis) đã đƣợc đề ra. Cách mô tả tín hiệu này thỏa mãn các yêu cầu của ngành chẩn đoán kĩ thuật là phải thể hiện đƣợc những thông tin về tần số, thời điểm và biên độ của các thành phần tín hiệu.
Wavelet packet là một trong những công cụ để phân tích tín hiệu trong miền thời gian- tần số.1 Phép biến đổi Fourier Phép biến đổi Fourier[1] đƣợc xây dựng trên nền tảng chuỗi Fourier, một hàm tuần hoàn x(t) có chu kỳ T bất kỳ có thể đƣợc khai triển thành tổng của các hàm điều hòa dƣới dạng: x(t ) xk (t) a0 ( ak cos( kt ) bk sin( k t )) (1.1) k 0 k 1 14 download by : skknchat@gmail.com Trong đó, các hệ số 𝑎 0 , 𝑎 𝑘, 𝑏𝑘 đƣợc gọi là các hệ số Fourier và đƣợc xác định: T 1 a0 x(t ) dt (1.4) T0 Ta sử dụng những ký hiệu sau: 1 a a 0 A0 ; Ak ak2 bk2 ; k arctan k ; fk k (1.5) 2 bk Khi đó biểu thức (1.1) đƣợc viết lại: 1 x(t ) xk (t) Ao Ak .6) còn đƣợc gọi là chuỗi Fourier thực một phía, do các tần số f k có giá trị dƣơng và các biên độ Ak nhận các giá trị là số thực. Tần số f1 1/ T (ứng với k=1) đƣợc gọi là tần số cơ bản, các tần số khác có trị số f k k. f1 với k là số nguyên dƣơng. Thành phần A0 / 2 ứng với tần số f 0 0 là hằng số.
Các hàm điều hòa có biên độ Ak và tần số f k đƣợc gọi là điều hòa bậc k. Trong trƣờng hợp tổng quát, hàm x(t) có giá trị phức, ta có dạng biểu diễn khác của chuỗi Fourier: 1 1 x(t)= A0 Ak .t k) ) 2 k 1 2 15 download by : skknchat@gmail.7) k2 k trong đó j là đơn vị ảo.7) đƣợc gọi là chuỗi Fourier phức hai phía, do các tần số f k nhận j các giá trị trên toàn trục số (-∞.+∞) và các biên độ X k Ak e k có giá trị phức. Đây là công thức đƣợc sử dụng trong tính toán chuỗi Fourier, đặc biệt trong các phần mềm tính toán hiện nay. Tuy nhiên, chỉ các tần số dƣơng mới đƣợc sử dụng để phân tích dao động trong các ứng dụng thực tế.
Trong trƣờng hợp tổng quát, với tín hiệu x(t) có thể tuần hoàn hoặc không tuần hoàn, phép biến đổi Fourier (FT) là một công cụ toán học quan trọng để biến đổi tín hiệu x(t) trong miền thời gian t sang miền tần số f : FT x(t ) X ( f ) x(t )e j 2 ft dt (1.8) Việc chuyển đổi tín hiệu X(f) trong miền tần số sang miền thời gian đƣợc thực hiện bởi phép biến đổi Fourier ngƣợc (iFT) iFT X ( f ) x (t ) X ( f ).9) Phép biến đổi Fourier có một số tính chất quan trọng sau: Phép biến đổi Fourier là phép biến đổi tuyến tính, tức : FT a1 x1 (t ) a2 x2 (t ) a1FT x1 (t ) a2 FT x2 (t ) (1.11) 16 download by : skknchat@gmail.com Nếu tín hiệu g(t)= x(t-t0 ) với t0 là hằng số, ta có quan hệ trong miền tần số: G ( f ) e j.13) Phổ tần số X(f) của tín hiệu x(t) cho phép ta xác định đƣợc các thành phần tần số chứa trong tín hiệu.8) ta có thể biểu diễn X(f) dƣới dạng phổ hai phía, gồm các tần số âm và các tần số dƣơng. Trong thực tế, ta chỉ giữ lại các thành phần tần số dƣơng để biểu diễn đồ thị phổ tần số. 1: Phổ tần số- biên độ của một tín hiệu điều biến biên độ có A1 20, A2 25 f1 10, f 2 50( Hz ) 1.2 Biến đổi Wavelet liên tục Phân tích tín hiệu trong miền thời gian (hoặc tần số) chỉ cung cấp cho ta các thông tin của tín hiệu chỉ duy nhất trong miền thời gian (hoặc tần số). Phép biến đổi Wavelet (Wavelet Transform - WT), trong đó có phép biến đổi Wavelet liên tục (Continuous Wavelet Transform- CWT) [2] [6] là phép biến đổi trong miền thời gian- tần số.
Phép biến đổi này đƣợc phát triển để khắc phục những thiếu sót của 17 download by : skknchat@gmail.com phép biến đổi Fourier, từ đó cho ta thấy đƣợc rõ hơn các thông tin chứa trong tín hiệu đồng thời trong miền thời gian và tần số.1 Phép biến đổi Wavelet liên tục Xét phép biến đổi Fourier của một hàm tín hiệu x(t) đƣợc phân tích thành tổng của các hàm điều hòa phức: 1 FT x (t ) X ( f ) x (t )e j 2 ftdt (1. Một cách tƣơng tự, phép biến đổi Wavelet sử dụng các hàm trong đó j Wavelet cơ sở. Một hàm theo biến thời gian (t ) đƣợc gọi là một hàm Wavelet cơ sở nếu thỏa mãn điều kiện sau: 1 2 2 (t ) dt 1 , = (t )dt 0 (1.15) là biến đổi Fourier của (t). Các hàm Wavelet đƣợc tạo ra từ Wavelet cơ sở với các biến tỉ lệ s, và hệ số dịch chuyển τ theo hệ thức sau: -1 t - ,s (t ) s ( 2 ) s , R (1.16) s Biến đổi Wavelet của một tín hiệu dao động x(t) đƣợc định nghĩa dƣới dạng các hàm Wavelet: 1 CWTx ( s,) x (t ), ,s (t ) s x(t ) ( t 2 * (1.17), tín hiệu x(t) đƣợc biến đổi thành hàm của hai biến s và τ.
Bằng cách thay đổi hệ số tỷ lệ s và dịch chuyển hàm 18 download by : skknchat@gmail.com Wavelet dọc theo thời gian một lƣợng τ, hệ số Wavelet CWTx (s ,) có thể tạo ra một đồ thị trong không gian ba chiều biểu diễn sự thay đổi biên độ của tín hiệu đồng thời theo trục thời gian (tỷ lệ với τ) và theo trục tần số (tỷ lệ với s), còn gọi là phân bố thời gian- tần số của tín hiệu. Một Wavelet cơ sở thƣờng hay đƣợc sử dụng đó là hàm Morlet- Wavelet có dạng: 1 t 2 (t ) .18) với 0 là một tham số hằng và có ảnh hƣởng quyết định đến độ phân giải tần số của đồ thị biên độ tín hiệu trong miền thời gian- tần số. Thông thƣờng, giá trị của 0 nằm trong khoảng 5 o 20 khi phân tích tín hiệu dao động. Hàm Morlet- Wavelet tƣơng ứng với hai giá trị khác nhau của 𝜔 0 đƣợc biểu diễn nhƣ hình 1.
2: Hàm Morlet-Wavelet (phần thực) Một đặc điểm quan trọng của phép biến đổi Wavelet với hàm Morlet- Wavelet là hệ số tỷ lệ s có quan hệ với tần số f theo hệ thức: 0 2 02 f (1.19) 4s t Biến đổi Fourier của hàm g (t ) ( ) có dạng: s 19 download by : skknchat@gmail.20) Áp dụng công thức (1.20), biến đổi Wavelet của x(t) theo biểu thức (1.