4.64 Lọc và cộng tín hiệu

Tạo ra tín hiệu sin trong hai chu kỳ có biên độ bằng 1, cộng vào tín hiệu nhiễu trắng có biên độ bằng 25% biên độ tín hiệu. Cho tín hiệu qua bộ lọc trung bình động có kích thước cửa sổ bằng 10. Vẽ tín hiệu trước và sau khi cộng nhiễu và lọc.

Giải: Kích thước cửa sổ lọc là 10.

t = linspace(-2*pi,2*pi,100);
x = sin(t);
n = 0.25*randn(size(t));
s = x + n;
% Filter
windowSize = 10;
b = (1/windowSize)*ones(1,windowSize);
a = 1;
y = filter(b,a,x);
plot(t,x,'-',t,s,'--',t,y,':','LineWidth',2)
legend('Tín hiệu vào','Tín hiệu vào với nhiễu','Tín hiệu sau khi lọc');
set(gcf,'color','w');

Tín hiệu trước và sau khi lọc nhiễu bằng bộ lọc trung bình động. 

4.63 Biểu đồ chòm sao điều chế 16-QAM

Vẽ biểu đồ chòm sao của điều chế 16-QAM tại máy thu của kênh truyền nhiễu trắng với tỷ lệ tín hiệu trên nhiễu là 20 dB trong hai trường hợp:
 a. Không có chuẩn hóa công suất.
 b. Có chuẩn hóa công suất.

 Giải: Trong ví dụ này, chúng ta không dùng hàm modnorm(), mà sử dụng tham số “UnitAveragePower” trong hàm qammod() để thực hiện chức năng chuẩn hóa công suất. Chế độ mặc định của hàm qammod() là không chuẩn hóa công suất. Để vẽ biểu đồ chòm sao, chúng ta dùng hàm scatterplot().

% Tạo dữ liệu vào
M = 16;
data = randi([0 M-1],1000,1);
SNR_dB = 20;

%% Trường hợp 1: Không chuẩn hóa công suất
% Điều chế tín hiệu 16-QAM
txSig = qammod(data,M);
% Cộng nhiễu vào tín hiệu -3dB
rxSig = awgn(txSig,SNR_dB);
scatterplot(rxSig)
set(gcf,'color','w');

%% Trường hợp 2: Chuẩn hóa công suất
txSig_N =qammod(data,16,'UnitAveragePower', true);
rxSig_N = awgn(txSig_N,SNR_dB);
scatterplot(rxSig_N)
set(gcf,'color','w');

Biểu đồ chòm sao của 16-QAM khi không chuẩn hóa công suất.

 

 Biểu đồ chòm sao của 16-QAM khi chuẩn hóa công suất.