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

Đề tài: Báo cáo thí nghiệm xử lý số tín hiệu


BÀI 1: MÔ PHỎNG HỆ THỐNG VÀ TÍN HIỆU RỜI RẠC BẰNG MATLAB

A. Tín hiệu và hệ thống rời rạc ở miền n

1.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 độ 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 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 <= n2
%----------------------------------------------
%[x,n] = r_expseq(a,n1,n2)
n = [n1:n2]; x = randn(size(n));


1.2 Viết các chương trình tạo hàm năng lượng của một dãy. 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ú pháp :
Ex = energy (x,n);

function [Ex] = energy(x,n)
%Tinh Ex = nang luong,
%----------------------------------------------
%[Ex] = energy(x,n)
Ex = sum(abs(x).^2);


1.3 Cho dóy x(n) = {1,2,3,4,5,6,7,6,5,4,3,2,1} . Viết chương trỡnh con thể trờn đồ thị cỏc dóy :
a. x1(n) = 2x(n – 5) – 3x(n+4)
b. x2(n) = x(3 – n) – x(n)x(n – 2)

n = [-2:10]; x = [1:7,6:-1:1];
%
[x11,n11] = sigshift(x,n,5); [x12,n12] = sigshift(x,n,-4);
[x1,n1] = sigadd(2*x11,n11,-3*x12,n12);
subplot(2,1,1); stem(n1,x1);
title('Day so theo dau bai 1.9a');
xlabel('n'); ylabel('x1(n)');
%
[x21,n21] = sigfold(x,n); [x21,n21] = sigshift(x21,n21,3);
[x22,n22] = sigshift(x,n,2); [x22,n22] = sigmult(x,n,x22,n22);
[x2,n2] = sigadd(x21,n21,x22,n22);
subplot(2,1,2); stem(n2,x2);
title('Day so theo dau bai 1.9b');
xlabel('n'); ylabel('x2(n)');