单载波频域均衡MMSE的代码

function s=modulation(sym,flag_m) % sym=randi([0 63],Num,1); % sym=0:Num-1; %Num==64,测试星座图用 Num=length(sym); s=zeros(1,Num); for i=1:Num switch sym(i) case 0 s(i)= 1.000000+1i*0.000000; case 1 s(i)= 0.822878+1i*0.568218; case 2 s(i)= 0.821137+1i*0.152996; case 3 s(i)= 0.932897+1i*0.360142; case 4 s(i)= 0.000000-1i*1.000000; case 5 s(i)= 0.822878-1i*0.568218; case 6 s(i)= 0.821137-1i*0.152996; case 7 s(i)= 0.932897-1i*0.360142; case 8 s(i)= 0.568218+1i*0.822878; case 9 s(i)= 0.588429+1i*0.588429; case 10 s(i)= 0.588429+1i*0.117686; case 11 s(i)= 0.588429+1i*0.353057; case 12 s(i)= 0.568218-1i*0.822878; case 13 s(i)= 0.588429-1i*0.588429; case 14 s(i)= 0.588429-1i*0.117686; case 15 s(i)= 0.588429-1i*0.353057; case 16 s(i)= 0.152996+1i*0.821137; case 17 s(i)= 0.117686+1i*0.588429; case 18 s(i)= 0.117686+1i*0.117686; case 19 s(i)= 0.117686+1i*0.353057; case 20 s(i)= 0.152996-1i*0.821137; case 21 s(i)= 0.117686-1i*0.588429; case 22 s(i)= 0.117686-1i*0.117686; case 23 s(i)= 0.117686-1i*0.353057; case 24 s(i)= 0.360142+1i*0.932897; case 25 s(i)= 0.353057+1i*0.588429; case 26 s(i)= 0.353057+1i*0.117686; case 27 s(i)= 0.353057+1i*0.353057; case 28 s(i)= 0.360142-1i*0.932897; case 29 s(i)= 0.353057-1i*0.588429; case 30 s(i)= 0.353057-1i*0.117686; case 31 s(i)= 0.353057-1i*0.353057; case 32 s(i)= 0.000000+1i*1.000000; case 33 s(i)=-0.822878+1i*0.568218; case 34 s(i)=-0.821137+1i*0.152996; case 35 s(i)=-0.932897+1i*0.360142; case 36 s(i)=-1.000000+1i*0.000000; case 37 s(i)=-0.822878-1i*0.568218; case 38 s(i)=-0.821137-1i*0.152996; case 39 s(i)=-0.932897-1i*0.360142; case 40 s(i)=-0.568218+1i*0.822878; case 41 s(i)=-0.588429+1i*0.588429; case 42 s(i)=-0.588429+1i*0.117686; case 43 s(i)=-0.588429+1i*0.353057; case 44 s(i)=-0.568218-1i*0.822878; case 45 s(i)=-0.588429-1i*0.588429; case 46 s(i)=-0.588429-1i*0.117686; case 47 s(i)=-0.588429-1i*0.353057; case 48 s(i)=-0.152996+1i*0.821137; case 49 s(i)=-0.117686+1i*0.588429; case 50 s(i)=-0.117686+1i*0.117686; case 51 s(i)=-0.117686+1i*0.353057; case 52 s(i)=-0.152996-1i*0.821137; case 53 s(i)=-0.117686-1i*0.588429; case 54 s(i)=-0.117686-1i*0.117686; case 55 s(i)=-0.117686-1i*0.353057; case 56 s(i)=-0.360142+1i*0.932897; case 57 s(i)=-0.353057+1i*0.588429; case 58 s(i)=-0.353057+1i*0.117686; case 59 s(i)=-0.353057+1i*0.353057; case 60 s(i)=-0.360142-1i*0.932897; case 61 s(i)=-0.353057-1i*0.588429; case 62 s(i)=-0.353057-1i*0.117686; case 63 s(i)=-0.353057-1i*0.353057; end end % scatterplot(s);%画星座图用 end