c4_qamdemo.m
上传用户:ay_070428
上传日期:2014-12-04
资源大小:11427k
文件大小:1k
源码类别:

语音合成与识别

开发平台:

Matlab

  1. levelx=input('Number of D Levels> ');
  2. levely=input('Number of Q Levels> ');
  3. m=input('Number of symbols> ');
  4. n=input('Number of samples per symbol> ');
  5. bw=input('Filter bandwidth,0<bw<1 > ');%
  6. %
  7. [xd,xq]=qam(levelx,levely,m,n);
  8. %
  9. [b,a]=butter(6,bw);%determine filter coefficients
  10. yd=filter(b,a,xd);%filter direct coefficient
  11. yq=filter(b,a,xq);%filter quadrature coefficient
  12. %
  13. subplot(2,2,1)%first pane
  14. plot(xd,xq,'o')%unfiltered scatterplot
  15. a=1.4;
  16. maxd=max(xd);maxq=max(xq);
  17. mind=min(xd);minq=min(xq);
  18. axis([a*mind a*maxd a*minq a*maxq])
  19. axis equal
  20. xlabel('xd');ylabel('xq');
  21. %
  22. subplot(2,2,2)%second pane
  23. plot(yd,yq)%filtered scatterplot
  24. axis equal
  25. xlabel('xd');ylabel('xq');
  26. %
  27. sym=30;%number of symbols in time plot
  28. nsym=(0:sym*n)/n;%x axis vector for time plots
  29. subplot(2,2,3)%third pane
  30. plot(nsym(1:sym*n),yd(1:sym*n))%filterd direct component
  31. xlabel('symbol index');
  32. ylabel('xd');
  33. %
  34. subplot(2,2,4)%fourth pane
  35. plot(nsym(1:sym*n),yd(1:sym*n))%filterd quadrature component
  36. xlabel('symbol index');
  37. ylabel('xq');
  38. %end of script file