ip_07_07.m
上传用户:loeagle
上传日期:2013-03-02
资源大小:1236k
文件大小:1k
源码类别:

通讯编程文档

开发平台:

Matlab

  1. % MATLAB script for Illustrative Problem 7.7. 
  2. echo on
  3. Tb=1;
  4. f1=1000/Tb;
  5. f2=f1+1/Tb;
  6. phi=pi/4;
  7. N=5000;    % number of samples
  8. t=0:Tb/(N-1):Tb;
  9. u1=cos(2*pi*f1*t);
  10. u2=cos(2*pi*f2*t);
  11. % Assuming that u1 is transmitted, the received signal r is
  12. sgma=1; % noise variance
  13. for i=1:N,
  14.   r(i)=cos(2*pi*f1*t(i)+phi)+gngauss(sgma);
  15.   echo off;
  16. end;
  17. echo on ;
  18. % The correlator outputs are computed next.
  19. v1=sin(2*pi*f1*t);
  20. v2=sin(2*pi*f2*t);
  21. r1c(1)=r(1)*u1(1);
  22. r1s(1)=r(1)*v1(1);
  23. r2c(1)=r(1)*u2(1);
  24. r2s(1)=r(1)*v2(1);
  25. for k=2:N,
  26.   r1c(k)=r1c(k-1)+r(k)*u1(k);
  27.   r1s(k)=r1s(k-1)+r(k)*v1(k);
  28.   r2c(k)=r2c(k-1)+r(k)*u2(k);
  29.   r2s(k)=r2s(k-1)+r(k)*v2(k);
  30.   echo off;
  31. end;
  32. echo on;
  33. % decision variables
  34. r1=r1c(5000)^2+r1s(5000)^2;
  35. r2=r2c(5000)^2+r2s(5000)^2;
  36. % Plotting commands follow.