Báo cáo xử lý số tín hiệu

Viết chương trình con tạo một dãy thực ngẫu nhiên xuất phát từ n1 đến n2 và có giá trị của biên độ theo phân bố Gauss với trung bình bằng 0, phương sai bằng 1. Yêu cầu chương trình con có các tham số đầu vào và đầu ra được nhập theo câu lệnh với cú pháp: [x,n] = randnseq(n1,n2); function [x,n] = randnseq(n1,n2) %Generates x(n) = a^n; n1
Nội dung trích xuất từ tài liệu:
Báo cáo xử lý số tín hiệuBáo cáo xử lý số tín hiệuBáo cáo xử lý số tín hiệu ĐHBKHN BÀI 1. MÔ PHỎNG HỆ THỐNG VÀ TÍN HIỆU RỜI RẠC BẰNG MATLABA. Tín hiệu và hệ thống rời rạc ở miền n1.1. Viết chương trình con tạo một dãy thực ngẫu nhiên xuất phát từ n1 đến n2 và có giá trị của biên độ theophân bố Gauss với trung bình bằng 0, phương sai bằng 1. Yêu cầu chương trình con có các tham s ố đầu vàovà đầu ra được nhập theo câu lệnh với cú pháp:[x,n] = randnseq(n1,n2); function [x,n] = randnseq(n1,n2) %Generates x(n) = a^n; n1 Báo cáo xử lý số tín hiệu ĐHBKHN1.4. Cho hệ thống được mô tả bởi phương trình sai phân tuyến tính hệ số hằng như sau: y(n) – y(n-1) + 0.9y(n-2) = x(n)Sử dụng hàm filter của MATLAB, viết chương trình thực hiện các công việc sau: a. Biểu diễn bằng đồ thị hàm đáp ứng xung đơn vị của hệ thống với -20 ≤ n ≤ 100 b. Biểu diễn bằng đồ thị dãy đáp ứng của hệ thống với -20 ≤ n ≤ 100 khi dãy đầu vàolà dãy nhảy đơn vị. b = [1]; a = [1, -1, 0.9]; % x= impseq(0,-20,120); n = [-20:120]; h= filter(b,a,x); subplot(2,1,1); stem(n,h); title(Impulse Response); xlabel(n); ylabel(h(n)); % x = stepseq(0,-20,120); s = filter(b,a,x); subplot(2,1,2); stem(n,s); title(Step Response); xlabel(n); ylabel(s(n));nang_trong_dem_90@yahoo.com Page 2Báo cáo xử lý số tín hiệu ĐHBKHNB. Tín hiệu và hệ thống rời rạc ở miền Z, miền tần số liên tục ω, và miền tần số rời rạc k1.5. Cho dãy x(n)=0.5nu(n)a. Dựa trên định nghĩa của biến đổi Z, tìm biến đổi Z của dãy trênb. Kiểm chứng lại kết quả câu a bằng hàm ztransc. Từ kết quả trên, tìm biến đổi Fourier của x(n)Dùng MATLAB thể hiện trên đồ thị phổ X(ejω) tại 501 điểm rời rạc trong khoảng [0,π] %Bien doi Z cua day syms n; ztrans(0.5^n) ans = 2*z/(2*z-1) %Do thi pho w = [0:1:500]*pi/500; X = exp(j*w) ./ (exp(j*w)- 0.5*ones(1,501)); magX = abs(X); angX = angle(X); realX = real(X); imagX = imag(X); % subplot(2,2,1); plot(w/pi,magX); grid; title(Magnitude Part); xlabel(frequency in pi units); ylabel(Magnitude); subplot(2,2,3); plot(w/pi,angX); grid; title(Angle Part); xlabel(frequency in pi units); ylabel(Radians); subplot(2,2,2); plot(w/pi,realX); grid; title(Real Part); xlabel(frequency in pi units); ylabel(Real); subplot(2,2,4); plot(w/pi,imagX); grid; title(Imaginary Part); xlabel(frequency in pi units); ylabel(Imaginary);nang_trong_dem_90@yahoo.com Page 3Báo cáo xử lý số tín hiệu ĐHBKHN1.6. Cho dãy: x(n) = {..,0,0,1,2,3,4,5,0,0,....} Đây là một dãy số xác định trong một khoảng hữu hạn từ -1 đến 3. Dựa trên công thức định nghĩa của biến đổi Fourier, viết chương trình tính và thể hiệnphổ của dãy x(n) tại 501 điểm rời rạc trong khoảng [0,π] n = -1:3; x = 1:5; k = 0:500; w = (pi/500)*k; X = x*(exp(-j*pi/500)).^(n*k); magX = abs(X); angX = angle(X); realX = real(X); imagX = imag(X); % subplot(2,2,1); plot(k/500,magX); grid; title(Magnitude Part); xlabel(frequency in pi units); ylabel(Magnitude); subplot(2,2,3); plot(k/500,angX); grid; title(Angle Part); xlabel(frequency in pi units); ylabel(Radians); subplot(2,2,2); plot(k/500,realX); grid; title(Real Part); xlabel(frequency in pi units); ylabel(Real); subplot(2,2,4); plot(k/500,imagX); grid; title(Imaginary Part); xlabel(frequency in pi units); ylabel(Imaginary);nang_trong_dem_90@yahoo.com Page 4Báo cáo xử lý số tín hiệu ĐHBKHN1.7. Một hàm ở miền Z được cho với công thức sau đây: Hàm số X(z) có thể viết dưới dạng tỷ số của hai đa thức theo z-1 như sau:a. Sử dụng lệnh residuez của MATLAB, tính các điểm cực, thặng dư tại các điểm cựcb. Từ kết quả câu trên, viết công thức khai triển X(z) thành tổng các phân thức đơn giản, từ đó tìm biến đổi Zngược của X(z), cho biết x(n) là một dãy nhân quảc. Kiểm chứng lại kết quả câu b bằng hàm iztrans %Tinh cac diem cuc, thang du tai cac diem cuc b = [0 1]; a = [3 -4 1]; [R,p,C] = residuez(b,a) ...