matlab code in digital communication

3 次查看(过去 30 天)
can you please explain the code for the above questions using the below code
clear all
% solution 6.5-1
% ==============
Given_Pr = [0.36 -1.2 0.21];
length_Pr = length(Given_Pr);
% Zero- clc
Forcing-Equalizer:
Pr_zf = toeplitz([Given_Pr([2:end]) zeros(1,2*1+1-length_Pr+1)], ...
[ Given_Pr([2:-1:1]) zeros(1,2*1+1-length_Pr+1) ]);
Po_zf = zeros(1,2*1+1);
Po_zf(2) = 1;
fprintf('solution 6.5-1 outputs are :\n');
C_val_ZF = [inv(Pr_zf)*Po_zf.'].'
mse_ZF = mse(C_val_ZF)
% solution 6.5-2
% ==============
Given_Pr = [0.36 -1.2 0.21];
length_Pr = length(Given_Pr);
K = 2;
N = 1;
% MMSE-Equalizer:
Pr_AutoCorr = conv(Given_Pr,fliplr(Given_Pr));
Pr_mmse = toeplitz([Pr_AutoCorr([3:end]) zeros(1,2*K+1-length_Pr)], ...
[ Pr_AutoCorr([3:end]) zeros(1,2*K+1-length_Pr) ]);
Pr_mmse = Pr_mmse + 1/2*10^(-1/10)*eye(2*K+1)*N;
Po_mmse = zeros(1,2*K+1);
Po_mmse([-1:1]+K+1) = fliplr(Given_Pr);
fprintf('\nsolution 6.5-2 outputs are :\n');
C_val_MMSE = [inv(Pr_mmse)*Po_mmse.'].'
mse_MMSE = mse(min(C_val_MMSE))

回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Manage Products 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by