【正文】 elseif x==33*j y=[0 1 1 1]。 elseif x==3+3*j y=[1 1 0 1]。 elseif x==13*j y=[1 0 1 1]。 elseif x==3+j y=[1 1 0 0]。 elseif x==13*j y=[0 0 1 1]。 elseif x==1+3*j y=[0 0 0 1]。 elseif x==1j y=[1 0 1 0]。 elseif x==1j y=[0 0 1 0]。 end x=plex(y,y1)。 elseif (y12) y1=3。(y1=2) y1=1。 else y=1。 elseif (y2) y=3。 if (y=0)amp。end%16QAM解調(diào)子程序function y=de_qam16(x) y=real(x)。 elseif y==[1 1 1 1] x=33*j。 elseif y==[0 1 0 1] x=3+3*j。 elseif y==[1 1 1 0] x=3j。 elseif y==[1 0 0 1] x=1+3*j。 elseif y==[0 1 1 0] x=3j。 elseif y==[0 1 0 0] x=3+j。 elseif y==[1 0 0 0] x=1+j。%16QAM調(diào)制子程序function x=qam16(y) if y==[0 0 0 0] x=1+j。 ij=ij+1。 for j=1:n y=de_qam16(matrix(j,i))。 else return end endend %16QAM數(shù)據(jù)解調(diào)子程序function x=de_modulation(matrix) [n,m]=size(matrix)。 if k=NL for j=1:N X_S(j,k)=signal(j,k)/SVD(j)。)*LS。 end UU_P=diag(s)。 end end for i=1:p %只取前p個(gè)較大的特征值，其余設(shè)為零 s(i)=chara(i)/(chara(i)+B/(C/(sgma^2)))。 UU(t,t)=realmin。t=j。 for i=1:N %將信道自相關(guān)矩陣的特征值從大到小排列 ic=realmin。 [XX_0,UU_0]=eig(Rhh)。 X=diag(X_P)。 X_P(j)=pilot_0。 B=17/9。 F=DFT(N)。 end out=X_S。 for i=1:interval k=interval*(a1)+i。 X=diag(X_P)。 X_P(j)=pilot_0。 F=DFT(N)。 end out=X_S。 for i=1:interval k=interval*(a1)+i。*X*F)*F’*X39。)*X*F)*A+Rhh))*inv(F39。MMSE=F*Rhh*inv((inv((F39。 end Y=Y_P.39。 for a=1:count for j=1:N Y_P(j)=pilot_signal(j,a)。 F=DFT(N)。 y=re_signal_0。 %提取的導(dǎo)頻信號(hào)矩陣 for ic=1:interval if (k+ic)=n for j=1:m re_signal_0(j,i)=matrix(j,(k+ic))。 for j=1:m pilot_0(j,num)=matrix(j,k)。 i=1。 end endx=cut_cp_0。%除去循環(huán)前綴子程序function x=cut_insert_cp(num,matrix) %num是循環(huán)前綴的長(zhǎng)度 [m,n]=size(matrix)。 tran_1(m,k)=trans_s(m,k)+noise。%輸入噪聲子程序function x=add_noise(sgma,trans_s) [n,nl]=size(trans_s)。hh=inv(aa)*RH*inv(aa39。 else RH(k,f)=(1exp(j*2*pi*(kf)/N))/(2*pi*j*L*(kf)/N)。hh=zeros(N)。end%生成信道自相關(guān)函數(shù)子程序function hh=Rhhh(N,L)。for f=1:num output_sig_serial= output_sig_serial+ray_chan(f,:).*delay_sig(f,:)。 end delay_sig(f,(t_shift(f)+1):chann_l)=input_sig_serial(1:(chann_lt_shift(f)))。enddelay_sig=zeros(num,chann_l)。 ray_chan(k,:)=atts*rayleigh(0,chann_l,t_interval,fmax)/sqrt(total_pow)。total_pow=sum(10.^(var_pow)/10)。chann_l=nl*l。[nl,l]=size(input_sig)。t_delay=delay/t_interval。%t_interval為離散信道抽樣時(shí)間間隔，OFDM符號(hào)長(zhǎng)度/(子信道數(shù)N+cp長(zhǎng)度lp)；%output_sig為經(jīng)過多徑信道的輸出信號(hào)矢量fmax=v*fc/3e+8。%多徑信道子程序function output_sig=multipath_channg(input_sig,num,var_pow,delay, v,fc,t_interval,N,counter)%input_sig,加了cp后的輸入信號(hào)矩陣，NL(N+cp長(zhǎng)度lp)；%num多徑數(shù)。 gauss_u2=gauss_u2+c_2(k)*cos(2*pi*f_2(k).*t+theta_2(k))。[f_2,c_2,theta_2]=parameter_classical(N_i,var_u,fmax)。t=(0:samp_num1)*t_interval。 gauss_u2=zeros(1,samp_num)。var_u=var_power/。 %正弦疊加法生成瑞利過程function ray_chan=rayleigh(var_power,samp_num,t_interval,fmax)%var_power為該徑平均功率單位db；%samp_num為離散信道采樣數(shù)%t_interval為抽樣時(shí)間間隔。 c_i=sigma*sqrt(2/N_i)*ones(size(un))。 un=rand(1,N_i)。 end end x=insert_cp_0。 for i=1:num %循環(huán)前綴CP的長(zhǎng)度 for j=1:n insert_cp_0(i,j)=matrix((mnum+i),j)。 if qNL return end end end y=count。 x=insert_qam_out.39。 end q=q+1。 y=count。 end w=w+1。 count=0。 %非共軛轉(zhuǎn)置 q=1。 %導(dǎo)頻信號(hào)的二進(jìn)制值 pilot_1=qam16(pilot_0)。 %每一列表示一個(gè)OFDM符號(hào) end end inputqam_0=qam_out。%將二進(jìn)制數(shù)據(jù)流進(jìn)行16QAM調(diào)制子程序function inputqam_0=invert(N,NL,input_1)for j=1:NL for n=1:N for ic=1:4 qam_input(ic)=input_1(j,(n1)*4+ic)。 else input(i,j)=0。) grid on%產(chǎn)生二進(jìn)制數(shù)據(jù)流子程序function x=input_b(N,NL) for i=1:NL input_0=rand(1,4*N)。) text(,39。) text(,39。) text(,39。) ylabel(39。) xlabel(39。) axis([0,1e3,1]) title(39。,SNR_dB,err_svd,39。,SNR_dB,err_ls,39。 %運(yùn)算所用時(shí)間end%時(shí)間統(tǒng)計(jì) time_mmse_add=time1 time_ls_add=time2 time_lmmse_add=time3%畫圖 semilogy(SNR_dB,err_mmse,39。 %誤碼數(shù)統(tǒng)計(jì) err_svd(l)=num_of_err_svd/KL %誤碼率 time_svd(l)=toc。modify_signal=invert(N,NL,modify_signal_b)。 %運(yùn)算所用時(shí)間% SVD估計(jì)修正 ticmodify_signal=svd(N,NL,interval,count,pilot,Rhh,p,sgma,pilot_signal,signal)。 %誤碼數(shù)統(tǒng)計(jì) err_ls(l)=num_of_err_ls/KL %誤碼率 time_ls(l)=toc。 modify_signal=invert(N,NL,modify_signal_b)。 time1=time1+time_mmse(l)。num_of_err_mmse=error_stat(qam_out,modify_signal)。 modify_signal_b=de_modulation(modify_signal)。 %最大時(shí)延對(duì)符號(hào)間隔的歸一化Rhh=Rhhh(N,L)。 %分離導(dǎo)頻信號(hào)和有用信號(hào) [pilot_signal,signal]=separate_pilot(interval,count,receive_signal_y)。 %去掉循環(huán)前綴 cut_cp=cut_insert_cp(num,rece_siganl)。 %每個(gè)OFDM符號(hào)的平均能量 sgma=sqrt(spow/(2*snr))。for k=1:nnl for b=1:mml spow=spow+real(trans_s(k,b))^2+imag(trans_s(k,b))^2。[nnl,mml]=size(trans_s)。trans_s=multipath_channg(trans_o,nn,var_po