exampl_6.m

仿真模拟多径传播的信道特性

%例六 多径传播

%本程序中时间单位是微秒
%频率单位为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