%例六 多徑傳播

%本程序中時間單位是微秒
%頻率單位為MHz
%碼速率單位是Mb/s

global dt t f df N
close all

N=2^13; %采樣點數
L=16;   %每碼元的采樣點數
M=N/L   %碼元數
Rb=2;		%碼速率是2Mb/s
Ts=1/Rb;        %碼元間隔
dt=Ts/L;        %時域采樣間隔
df=1/(N*dt)     %頻域采樣間隔
T=N*dt          %截短時間
Bs=N*df/2       %系統帶寬
Na=6;           %示波器掃描寬度為6個碼元
alpha=1.

t=linspace(-T/2,T/2,N); %時域橫坐標
f=linspace(-Bs,Bs,N)+eps;   %頻域橫坐標

%升余弦
hr1=sin(pi*t/Ts)./(pi*t/Ts);
hr2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2);
hr=hr1.*hr2;		 
HR=t2f(hr);

figure(1)
 set(1,'Position',[300,10,300,200])
                %設定窗口位置及大小
figure(2)
 set(2,'Position',[50,500,400,150])
                %設定窗口位置及大小 
figure(3)
 set(3,'Position',[10,10,250,200])
                %設定窗口位置及大小 

for ii=1:100
  a=sqrt(rand^2+rand^2);
  tao=rand*Ts;
  phi=rand*2*pi;
  [a,tao,phi]
  C=1-a*exp(-j*(2*pi*f*tao+phi));
  aC=10*log10(abs(C));

  figure(1)
  plot(f,aC)
  grid
  axis([-2/Ts,2/Ts,-30,10])
  xlabel('f (MHz)')
  ylabel('信道特性 (dB)')

  code=sign(randn(1,M));
  imp=zeros(1,N);    %產生沖激序列
  imp(L/2:L:N)=code/dt;
  IMP=t2f(imp);

  H=C.*HR;		%總體特性
  h=real(f2t(H));

  figure(3)
  plot(t,h)
  grid
  axis([-Ts,3*Ts,-1,2])
  xlabel('t (us)')
  ylabel('總體時域特性 ')

  S=IMP.*H;  		%接收信號
  s=real(f2t(S));

  figure(2) 
  clf
  hold on
  grid
  xlabel('t (us)')
  ylabel('V')
  tt=[0:dt:Na*L*dt];
  for jj=1:Na*L:N/2
   plot(tt,s(jj:jj+Na*L),'g');
  end
  axis([0,Na*Ts,min(s)-0.3,max(s)+0.3])
  '按任意鍵繼續'
  pause     %暫停等待按鍵
end