ằ b = : N; % bins ằ Ts = 1/128; theo giây % Khoảng lấy mẫu ằ fs = 1/Ts; % Tần số lấy mẫu theo Hz ằ ts = Ts x (b - 1) % Khoảng lấy mẫu ằ a1 = ; f1 = 16; ằ x1 = a1 * sin (2 * pi * f1 * ts) ; % tín hiệu ằ a2 = 3; f2 = 48 ằ x2 = a2 * sin (2 * pi * f2 * ts) ; % tín hiệu thứ hai ằ x = x1 + x2; Nếu bạn nhìn thấy kết tín hiệu (hình 1.13) đa lệnh sau vào ằ plot (ts, x) ằ xlabel (Time, s) , y label (x) Chúng ta xây dựng chấm điểm phổ công suất ằ X = fft(x); % DFT x ằ pwr = x * cosj (X) / N % Công suất tín hiệu ằ frs = (b = 1), N * fs ằ Plot (frs, pwr) % Các tần số % chấm điểm phổ công suất Kết chấm điểm hình 1.14 Tín hiệu x đạt công suất điểm 65, với tần số f = 16H2 (65 = + 512 x 16 / 129 128) bin 449, phần chậm liên hợp số tín hiệu 65 = + 65 đợc cất bin 449 = + 512 - 64 HìnhI.15 Phổ lợng tín hiệu x=x1+x2 Phổ công suất đợc hình 1.14 Ta kiểm tra pwr (65) = (a1/2)2.N Tơng tự nh tín hiệu x2 Công suất bin 193 321 pwr (193) = pwr (321) = (a2/2)2 N Ví dụ 1.5: Nhận dạng tần số thành phần công suất Trong thí dụ phân tích tín hiệu tam giác chu kỳ S = giây điểm nhảy biên độ vào thành phần tần số chúng sử dụng 512 điểm lấy mẫu Chúng ta quan tâm đến việc tìm phần trăm công suất tổng thành phần tín hiệu đợc nhận từ gốc, cách tách thành phần từ thành phần Chúng ta muốn biết cách làm xấp xỉ tín hiệu với tín hiệu chuẩn Đầu tiên, xây dựng phơng án rời rạc x tín hiệu lấy mẫu 512 điểm ằ T = S; ằ N = 512; ằ t = linspace (0,T, N + 1) ; t = (1 : N); ằ x1 = * t/T - 1/2 ; x2 = 2*(T - t) / T - / 2; ằ x = (x1, x2); % tín hiệu tam giác xây dựng phổ công suất chúng: ằb=1:N % Khoảng lấy mẫu tần số ằ X = fft (x); ằ Ts = T / N ; fs = N/T % (b - 1) / N * fs ằ prw = X * conj (X) / N; Để kiểm tra kết chúng ta, dùng đẳng thức Parseval Những số sau phải ằ[sum (pow) norm (x)^ 2] ans = 42.6680 42.6680 Dễ dàng nhận thấy tần số gồm thành phần lớn công suất, sử dụng hàm sort với quay trở lại phần tử pow cách tăng điểm: ằ [spow, spos] = sort (pow); Chúng ta tìm chữ số tần số thành phần công suất lớn nhất: 130 ằ m = 4; spos (N: -1 : (N - m + 1) Chúng ta thấy tần số cấu thành 512, 2, 510 Bây xây dựng tín hiệu xấp xỉ ằ X4 = zesos (X); % Vùng đổ xấp xỉ X ằ h = [512 510 4]; ằ X4 (h) = X (h); suất cao % chép binh cầu thành công Phần trăm công suất tạo thành thành phần đạo đợc đa ằ pere = 100 * (norm (X4) / norm (X))^2 Kết luận, 99,7698 % công suất đợc tạo thành nhóm, tơng ứng với tần số 0.2 Hz; liên quan đến bin số 2, tần số truyền đạt liên quan đến bin 512, giao động thứ 2, 0.6 Hz liên quan đến bin số 4, hệ số truyền nó, liên quan đến bin 510 Chúng ta sử dụng kết qủa ví dụ 1.8 Những dòng sau làm tiến gần đến tín hiệu tam giác góc đợc xấp xỉ, xem hình 1.15 ằ x4 = ift (X4); ằ plot (f, [x; x4]) , grid ằ xlabel ( t ), ylabel ( Tín hiệu tam giác xấp xỉ ý) Hình I.16 Sự xấp xỉ tín hiệu hình tam giác 131