matlab仿真CMMB源代码http:/www.codesoso.net/Search?q=cmmb&l=matlabfunction fredata=change(fre_sum) %转换成频域信号% f=1:4096; for wf=1:4096 if wf=1 fredata(wf)=0; elseif wf=2&wf=2559&wf=1&n=160&n=1&n=1280&n=2047 Xb(n+1)=1-2*PN(n-511); else Xb(n+1)=0; end end Sb_part=ifft(Xb,2048); Sb_all=Sb_part Sb_part; %主循环,53个OFDM信号% %for ofdm_counter=0:52 %每个时隙有53个OFDM符号 %for ofdm_counter=0:1 ofdm_counter=0; %产生频域的3076个符号% if(mod(ofdm_counter,2)=0) Y_before=ofdm_singal(0,ofdmdata(ofdm_counter*2610+1:ofdm_counter*2610+2610); else Y_before=ofdm_singal(1,ofdmdata(ofdm_counter*2610+1:ofdm_counter*2610+2610); end %产生扰码后的3076个符号% Y_after=Y_before.*Pc(ofdm_counter*3076+1:ofdm_counter*3076+3076); %得到频域的4096*53个符号% Y_change=change(Y_after); ofdm_result=ofdm_result,Y_change; %得到时域的4096*53个符号% timedata=ifft(Y_change,4096); %在时域上添加保护间隔% timedata_addcp_part=timedata_addcp(0,timedata); timedata_addcp_part %timedata_addcp_result=timedata_addcp_result timedata_addcp_part; %timedata_result=timedata_result,timedata; %plot(t,abs(timedata_addcp_part); %end %UNTITLED1 Summary of this function goes here %Detailed explanation goes here %CMMB 的频率交织 %clear; function fre_sum=ofdm_singal(qiou,ofdmdata_part) nocarrier = 4096; %共4096个载波 fredata =; %送给IFFT的数据4096个 conpilot_conter =1; data_counter = 1; conpilot_bit = 1; conpilot_direc = 1; scatpilot_conter0 =0; scatpilot_conter1 =192; %数组下标 condata =randint(1,16); %16bit的连续载频的数据信号 %ofdmdata =1/sqrt(2)*(-1+2*round(rand(1,NUM)+i*(-1+2*round(rand(1,NUM); %QPSK星座映射的数据信号，2610个复数据信号 a=0 22 78 92 168 174 244 274 278 344 382 424 426 496 500 564 608; %连续导频的地址 b=650 688 712 740 772 846 848 932 942 950 980 1012 1066 1126 1158; c=1214 1244 1276 1280 1326 1378 1408 1508 1537 1538 1566 1666 1736; d=1748 1794 1798 1830 1860 1916 1948 2008 2062 2094 2124 2132 2142; e=2226 2228 2302 2334 2362 2386 2424 2466 2510 2574 2578 2648 2650; g=2692 2730 2796 2800 2830 2900 2906 2982 2996 3052 3075; aa=22 78 92 168 174 2