帮我运行个程序吧

2019-07-17 14:58发布

帮我运行个程序吧
  1. BitsTx = floor(rand(1,1920)*2);
  2. SymQAMtmp1= reshape(BitsTx,4,1920/4).';
  3. for n=1:480
  4.     switch int2str(SymQAMtmp1(n,:))
  5.         case int2str([0 0 0 0])
  6.             SymQAMtmp(n)=-3-3i;
  7.         case int2str([0 0 0 1])
  8.             SymQAMtmp(n)=-3-i;
  9.         case int2str([0 0 1 0])
  10.             SymQAMtmp(n)=-3+3i;
  11.         case int2str([0 0 1 1])
  12.             SymQAMtmp(n)=-3+i;
  13.         case int2str([0 1 0 0])
  14.             SymQAMtmp(n)=-1-3i;
  15.         case int2str([0 1 0 1])
  16.             SymQAMtmp(n)=-1-i;
  17.         case int2str([0 1 1 0])
  18.             SymQAMtmp(n)=-1+3i;
  19.         case int2str([0 1 1 1])
  20.             SymQAMtmp(n)=-1+i;
  21.         case int2str([1 0 0 0])
  22.             SymQAMtmp(n)=3-3i;
  23.         case int2str([1 0 0 1])
  24.             SymQAMtmp(n)=3-i;
  25.         case int2str([1 0 1 0])
  26.             SymQAMtmp(n)=3+3i;
  27.         case int2str([1 0 1 1])
  28.             SymQAMtmp(n)=3+i;
  29.         case int2str([1 1 0 0])
  30.             SymQAMtmp(n)=1-3i;
  31.         case int2str([1 1 0 1])
  32.             SymQAMtmp(n)=1-i;
  33.         case int2str([1 1 1 0])
  34.             SymQAMtmp(n)=1+3i;
  35.         case int2str([1 1 1 1])
  36.             SymQAMtmp(n)=1+i;
  37.     end
  38. end
  39. symQAMout=SymQAMtmp/sqrt(10)
  40. figure(1);

  41. plot(symQAMout,'*r');%16QAM调制后星座图
  42. title('16QAM调制后星座图');
  43. axis([-4,4,-4,4]);
  44. grid on
  45. SymQAMtmp2=reshape(symQAMout,10,48);
  46. A=zeros(10,64);
  47. A(:, [8, 44, 58]) = ones(10, 3);
  48. A(:, 22) = -ones(10, 1);
  49. org = [2:7,9:21,23:27,39:43,45:57,59:64];
  50. A(:, org) = SymQAMtmp2;
  51. y=ifft(A(:,:),64,2);
  52. figure(2);
  53. plot(y,'*r');
  54. axis([-0.3,0.3,-0.3,0.3]);
  55. grid on
复制代码把结果给我发一下,谢谢
友情提示: 此问题已得到解决,问题已经关闭,关闭后问题禁止继续编辑,回答。