【正文】 delay_sig(f,:)。 endfor k=1:l output_sig(:,k)=output_sig_serial(((k1)*nl+1):k*nl).39。end%生成信道自相關(guān)函數(shù)子程序function hh=Rhhh(N,L)。RH=zeros(N)。hh=zeros(N)。for k=1:N for f=1:N if k==f RH(k,f)=1/L。 else RH(k,f)=(1exp(j*2*pi*(kf)/N))/(2*pi*j*L*(kf)/N)。 end endendaa=DFT(N)。hh=inv(aa)*RH*inv(aa39。)。%輸入噪聲子程序function x=add_noise(sgma,trans_s) [n,nl]=size(trans_s)。 for k=1:nl for m=1:n %輸入隨機(jī)噪聲 noise=normrnd(0,sgma)+normrnd(0,sgma)*sqrt(1)。 tran_1(m,k)=trans_s(m,k)+noise。 end endx=tran_1。%除去循環(huán)前綴子程序function x=cut_insert_cp(num,matrix) %num是循環(huán)前綴的長(zhǎng)度 [m,n]=size(matrix)。 for i=1:(mnum) for j=1:n cut_cp_0(i,j)=matrix(num+i,j)。 end endx=cut_cp_0。%分離導(dǎo)頻信號(hào)與數(shù)據(jù)信號(hào)子程序function [x,y]=separate_pilot(interval,count,matrix) [m,n]=size(matrix)。 i=1。 for num=1:count k=1+(num1)*(interval+1)。 for j=1:m pilot_0(j,num)=matrix(j,k)。 end x=pilot_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))。 end i=i+1。 y=re_signal_0。 %提取的數(shù)據(jù)信號(hào)矩陣 else return end endend%MMSE估計(jì)算法子程序function out=mmse(N,NL,interval,count,pilot,Rhh,sgma,pilot_signal,signal)pilot_0=qam16(pilot)。 F=DFT(N)。 A=sgma^2。 for a=1:count for j=1:N Y_P(j)=pilot_signal(j,a)。 X_P(j)=pilot_0。 end Y=Y_P.39。 X=diag(X_P)。MMSE=F*Rhh*inv((inv((F39。)*(X39。)*X*F)*A+Rhh))*inv(F39。*X39。*X*F)*F’*X39。*Y。 for i=1:interval k=interval*(a1)+i。 if k=NL for j=1:N X_S(j,k)=signal(j,k)/MMSE(j)。 end out=X_S。 else return end endend %LS估計(jì)算法子程序function out=ls(N,NL,interval,count,pilot,pilot_signal,signal) pilot_0=qam16(pilot)。 F=DFT(N)。 for a=1:count for j=1:N Y_P(j)=pilot_signal(j,a)。 X_P(j)=pilot_0。 end Y=Y_P.39。 X=diag(X_P)。 LS=inv(X)*Y。 for i=1:interval k=interval*(a1)+i。 if k=NL for j=1:N X_S(j,k)=signal(j,k)/LS(j)。 end out=X_S。 else return end endend %SVD估計(jì)算法子程序 function out=svd(N,NL,interval,count,pilot,Rhh,p,sgma,pilot_signal,signal) pilot_0=qam16(pilot)。 F=DFT(N)。C=4*(1^2+1^2)*(1/16)+2*4*(1^2+3^2)*(1/16)+4*(3^3+3^3)*(1/16)。 B=17/9。 for a=1:count for j=1:N Y_P(j)=pilot_signal(j,a)。 X_P(j)=pilot_0。 end Y=Y_P.39。 X=diag(X_P)。 LS=inv(X)*Y。 [XX_0,UU_0]=eig(Rhh)。 %信道自相關(guān)矩陣的特征值分解 UU=UU_0。 for i=1:N %將信道自相關(guān)矩陣的特征值從大到小排列 ic=realmin。 for j=1:N if icUU(j,j) ic=UU(j,j)。t=j。 end end chara(i)=UU(t,t)。 UU(t,t)=realmin。 for j=1:N XX(j,i)=XX_0(j,t)。 end end for i=1:p %只取前p個(gè)較大的特征值，其余設(shè)為零 s(i)=chara(i)/(chara(i)+B/(C/(sgma^2)))。 end for i=(p+1):N s(i)=0。 end UU_P=diag(s)。 SVD=XX*UU_P*(XX.39。)*LS。 %估計(jì)信道特性SVD for i=1:interval k=interval*(a1)+i。 if k=NL for j=1:N X_S(j,k)=signal(j,k)/SVD(j)。 end out=X_S。 else return end endend %16QAM數(shù)據(jù)解調(diào)子程序function x=de_modulation(matrix) [n,m]=size(matrix)。 for i=1:m ij=1。 for j=1:n y=de_qam16(matrix(j,i))。 for ic=1:4 out_b(i,ij)=y(ic)。 ij=ij+1。 end endend x=out_b。%16QAM調(diào)制子程序function x=qam16(y) if y==[0 0 0 0] x=1+j。 elseif y==[0 0 1 0] x=1j。 elseif y==[1 0 0 0] x=1+j。 elseif y==[1 0 1 0] x=1j。 elseif y==[0 1 0 0] x=3+j。 elseif y==[0 0 0 1] x=1+3*j。 elseif y==[0 1 1 0] x=3j。 elseif y==[0 0 1 1] x=13*j。 elseif y==[1 0 0 1] x=1+3*j。 elseif y==[1 1 0 0] x=3+j。 elseif y==[1 1 1 0] x=3j。 elseif y==[1 0 1 1] x=13*j。 elseif y==[0 1 0 1] x=3+3*j。 elseif y==[1 1 0 1] x=3+3*j。 elseif y==[1 1 1 1] x=33*j。 elseif y==[0 1 1 1] x=33*j。end%16QAM解調(diào)子程序function y=de_qam16(x) y=real(x)。 y1=imag(x)。 if (y=0)amp。(y=2) y=1。 elseif (y2) y=3。 elseif (y2) y=3。 else y=1。 end if (y1=0)amp。(y1=2) y1=1。 elseif (y12) y1=3。 elseif (y12) y1=3。 else y1=1。 end x=plex(y,y1)。 if x==1+j y=[0 0 0 0]。 elseif x==1j y=[0 0 1 0]。 elseif x==1+j y=[1 0 0 0]。 elseif x==1j y=[1 0 1 0]。 elseif x==3+j y=[0 1 0 0]。 elseif x==1+3*j y=[0 0 0 1]。 elseif x==3j y=[0 1 1 0]。 elseif x==13*j y=[0 0 1 1]。 elseif x==1+3*j y=[1 0 0 1]。 elseif x==3+j y=[1 1 0 0]。 elseif x==3j y=[1 1 1 0]。 elseif x==13*j y=[1 0 1 1]。 elseif x==3+3*j y=[0 1 0 1]。 elseif x==3+3*j y=[1 1 0 1]。 elseif x==33*j y=[1 1 1 1]。 elseif x==33*j y=[0 1 1 1]。 end襖羋蕆袇螀芇蕿蝕聿芆艿蒃肅芅蒁螈羈芄薃薁袆芃芃螆螂芃蒞蕿肁節(jié)蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈螞螂羂薁袈肀肁芀蟻羆肁莃袆袂肀薅蠆袈聿蚇蒂膇肈莇螇肅肇葿薀罿肆薂螆裊肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羈膂莈蚅袇膁蒀袀螃膀薂蚃肂腿節(jié)衿羈腿莄螞襖羋蕆袇螀芇蕿蝕聿芆艿蒃肅芅蒁螈羈芄薃薁袆芃芃螆螂芃蒞蕿肁節(jié)蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈螞螂羂薁袈肀肁芀蟻羆肁莃袆袂肀薅蠆袈聿蚇蒂膇肈莇螇肅肇葿薀罿肆薂