%==========================================================================
%==========================================================================
%Cviceni 10: Radiovy kanal
%==========================================================================
%==========================================================================
close all
clear all
%==========================================================================
%ukol 1
%--------------------------------------------------------------------------
%Prostudujte si základní parametry fyzické vrstvy systému WiMAX(Fixed WiMAX
%OFDM-PHY).Seznamte se s modelem kanálu SUI-3 (Stanford University Interim)
%uvedeném v dokumentu IEEE 802.16.3c-01/29r4 na str. 25.
%==========================================================================
%----------------------------------
%nastaveni vstupnich parametru
%----------------------------------
%fs/b=7/6; pomer vzorkovaciho kmitoctu ku sirce pasma; zname z tab. 9.1.
%b=3.5*10^6;
fs=7*3.5*10e6/6
ts=1/fs
fd=0.4; %doppleruv posuv
k=1; %pomer vykonu mezi jednotlivyma cestama
tau=[0 0.5 1]*1e-6;% zpozdeni
pdb=[0 -5 -10]; %pomer vykonu v dB
%==========================================================================
%ukol 2
%--------------------------------------------------------------------------
% Seznamte se s možností modelování rádiových únikových kanálù
% v programu MATLAB
%==========================================================================
X=ones(1,256); % vysleme vektor samych jednicek
x=ifft(X); %
mm=256/4; %cyklicky prefix
x1=[x(:,end-mm+1:end) x];
chan=ricianchan(ts,fd,k,tau,pdb);% rayleighchan %vytvorim objekt typu kanal
y=filter(chan,x1); %filtruju signal nastavenym kanalem
%==========================================================================
%ukol3
%--------------------------------------------------------------------------
% V programu MATLAB modelujte vliv rádiového kanálu SUI-3 na OFDM systém
% Fixed WiMAX OFDM-PHY. Porovnejte teoretickou hodnotu koherentní šíøky
% pásma kanálu s hodnotou odhadnutou na základì simulace.
%==========================================================================
y1=y(:,mm+1:end);%odstraneni cyklickeho prefixu
Y=fft(y1);
H=Y./X;

%==========================================================================
%ukol 4
%--------------------------------------------------------------------------
% Zmìòte parametry kanálu (maximální èasové rozšíøení, zisk jednotlivých
% cest) a sledujte vliv na koherentní šíøku pásma a únik na jednotlivých
% nosných.
%==========================================================================
%------------------------------------------------
%menime parametry 1
%------------------------------------------------
tau=[0 0.5 2]*1e-6;
pdb=[0 -1 -10]; %pomer vykonu v dB
X=ones(1,256); % vysleme vektor samych jednicek
x=ifft(X); %
mm=256/4; %cyklicky prefix
x1=[x(:,end-mm+1:end) x];
chan=ricianchan(ts,fd,k,tau,pdb);% rayleighchan %vytvorim objekt typu kanal
y=filter(chan,x1);
y1=y(65:end);%odstraneni cyklickeho prefixu
Y=fft(y1);
H=Y./X;
subplot(2,1,1)
plot(abs(H))
axis([0 260 0 2])
%------------------------------------------------
%menime parametry 2
%------------------------------------------------
tau=[0 0.5 2]*1e-6;
pdb=[0 -1 -1]; %pomer vykonu v dB
X=ones(1,256); % vysleme vektor samych jednicek
x=ifft(X); %
mm=256/4; %cyklicky prefix
x1=[x(:,end-mm+1:end) x];
chan=ricianchan(ts,fd,k,tau,pdb);% rayleighchan %vytvorim objekt typu kanal
y=filter(chan,x1);
y1=y(65:end);%odstraneni cyklickeho prefixu
Y=fft(y1);
H=Y./X;