【正文】 spow=spow/(N*NL)。 %加入隨機高斯白噪聲 rece_siganl=add_noise(sgma,trans_o)。 %離散傅里葉變換 receive_signal_y=fft(cut_cp,N)。%信道自相關函數(shù)L=5。 %Rhh信道自相關函數(shù)%三種估計算法 %MMSE估計修正tic modify_signal=mmse(N,NL,interval,count,pilot,Rhh,sgma,pilot_signal,signal)。 modify_signal=invert(N,NL,modify_signal_b)。 %誤碼數(shù)統(tǒng)計 err_mmse(l)=num_of_err_mmse/KL %誤碼率 time_mmse(l)=toc。 %運算所用時間%LS估計修正 tic modify_signal=ls(N,NL,interval,count,pilot,pilot_signal,signal)。 num_of_err_ls=error_stat(qam_out,modify_signal)。 time2=time2+time_ls(l)。modify_signal_b=de_modulation(modify_signal)。 num_of_err_svd=error_stat(qam_out,modify_signal)。 time3=time3+time_svd(l)。ko39。g*39。rd39。誤碼率曲線圖39。信噪比(db)39。誤碼率39。圓圈標記曲線為MMSE39。星號標記曲線為SVD39。菱形標記曲線為LS39。 %輸入的二進制數(shù)據(jù)序列 for j=1:4*N if input_0(j) input(i,j)=1。 end end endx=input。 %qam的輸入 end qam_out(n,j)=qam16(qam_input)。%插入導頻信號子程序function [x,y]=insert_pilot(N,NL,interval,pilot,matrix_1) pilot_0=pilot。 %導頻信號的16qam值 matrix_1=matrix_1.39。 %q表示轉置后的行數(shù)即NL w=1。 %插入導頻信號的次數(shù) while q=realmax for e=1:N insert_qam_out(w,e)=pilot_1。 count=count+1。 for ic=1:interval %每隔interval個插入導頻信號 for e=1:N insert_qam_out(w,e)=matrix_1(q,e)。 w=w+1。 %非共軛轉置.39。%插入循環(huán)前綴子程序function x=insert_cp(num,matrix) [m,n]=size(matrix)。 end end for i=1:m for j=1:n insert_cp_0((i+num),j)=matrix(i,j)。 %Monte Carlo方法生成正弦波疊加法需要的各參數(shù)值子程序function [f_i,c_i,theta_i]=parameter_classical(N_i,var,fmax)sigma=sqrt(var)。f_i=fmax*sin(pi*un/2)。theta_i=rand(1,N_i)*2*pi。%fmax為最大多普勒頻移val_pow=10^(var_power/10)。gauss_u1=zeros(1,samp_num)。N_i=25。[f_1,c_1,theta_1]=parameter_classical(N_i,var_u,fmax)。for k=1:N_i gauss_u1=gauss_u1+c_1(k)*cos(2*pi*f_1(k).*t+theta_1(k))。endray_chan=sqrt(gauss_u1.^2+gauss_u2.^2)。var_pow各徑平均功率；delay各徑延時，單位s；%v運動速率m/s；fc為載頻hz。 %最大doppler頻移t_shift=floor(delay/t_interval)。theta_shift=2*pi*fc*delay。output_sig=zeros(size(input_sig))。ray_chan=zeros(num,chann_l)。for k=1:num atts=10.^(*(var_pow(k)))。 endfor k=1:l input_sig_serial(((k1)*nl+1):k*nl)=input_sig(:,k).39。for f=1:num if t_shift(f)~=0 delay_sig(f,1:t_shift(f))=zeros(1,t_shift(f))。endoutput_sig_serial=zeros(1,chann_l)。 endfor k=1:l output_sig(:,k)=output_sig_serial(((k1)*nl+1):k*nl).39。RH=zeros(N)。for k=1:N for f=1:N if k==f RH(k,f)=1/L。 end endendaa=DFT(N)。)。 for k=1:nl for m=1:n %輸入隨機噪聲 noise=normrnd(0,sgma)+normrnd(0,sgma)*sqrt(1)。 end endx=tran_1。 for i=1:(mnum) for j=1:n cut_cp_0(i,j)=matrix(num+i,j)。%分離導頻信號與數(shù)據(jù)信號子程序function [x,y]=separate_pilot(interval,count,matrix) [m,n]=size(matrix)。 for num=1:count k=1+(num1)*(interval+1)。 end x=pilot_0。 end i=i+1。 %提取的數(shù)據(jù)信號矩陣 else return end endend%MMSE估計算法子程序function out=mmse(N,NL,interval,count,pilot,Rhh,sgma,pilot_signal,signal)pilot_0=qam16(pilot)。 A=sgma^2。 X_P(j)=pilot_0。 X=diag(X_P)。)*(X39。*X39。*Y。 if k=NL for j=1:N X_S(j,k)=signal(j,k)/MMSE(j)。 else return end endend %LS估計算法子程序function out=ls(N,NL,interval,count,pilot,pilot_signal,signal) pilot_0=qam16(pilot)。 for a=1:count for j=1:N Y_P(j)=pilot_signal(j,a)。 end Y=Y_P.39。 LS=inv(X)*Y。 if k=NL for j=1:N X_S(j,k)=signal(j,k)/LS(j)。 else return end endend %SVD估計算法子程序 function out=svd(N,NL,interval,count,pilot,Rhh,p,sgma,pilot_signal,signal) pilot_0=qam16(pilot)。C=4*(1^2+1^2)*(1/16)+2*4*(1^2+3^2)*(1/16)+4*(3^3+3^3)*(1/16)。 for a=1:count for j=1:N Y_P(j)=pilot_signal(j,a)。 end Y=Y_P.39。 LS=inv(X)*Y。 %信道自相關矩陣的特征值分解 UU=UU_0。 for j=1:N if icUU(j,j) ic=UU(j,j)。 end end chara(i)=UU(t,t)。 for j=1:N XX(j,i)=XX_0(j,t)。 end for i=(p+1):N s(i)=0。 SVD=XX*UU_P*(XX.39。 %估計信道特性SVD for i=1:interval k=interval*(a1)+i。 end out=X_S。 for i=1:m ij=1。 for ic=1:4 out_b(i,ij)=y(ic)。 end endend x=out_b。 elseif y==[0 0 1 0] x=1j。 elseif y==[1 0 1 0] x=1j。 elseif y==[0 0 0 1] x=1+3*j。 elseif y==[0 0 1 1] x=13*j。 elseif y==[1 1 0 0] x=3+j。 elseif y==[1 0 1 1] x=13*j。 elseif y==[1 1 0 1] x=3+3*j。 elseif y==[0 1 1 1] x=33*j。 y1=imag(x)。(y=2) y=1。 elseif (y2) y=3。 end if (y1=0)amp。 elseif (y12) y1=3。 else y1=1。 if x==1+j y=[0 0 0 0]。 elseif x==1+j y=[1 0 0 0]。 elseif x==3+j y=[0 1 0 0]。 elseif x==3j y=[0 1 1 0]。 elseif x==1+3*j y=[1 0 0 1]。 elseif x==3j y=[1 1 1 0]。 elseif x==3+3*j y=[0 1 0 1]。 elseif x==33*j y=[1 1 1 1]。 end襖羋蕆袇螀芇蕿蝕聿芆艿蒃肅芅蒁螈羈芄薃薁袆芃芃螆螂芃蒞蕿肁節(jié)蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈螞螂羂薁袈肀肁芀蟻羆肁莃袆袂肀薅蠆袈聿蚇蒂膇肈莇螇肅肇葿薀罿肆薂螆裊肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羈膂莈蚅袇膁蒀袀螃膀薂蚃肂腿節(jié)衿羈腿莄螞襖羋蕆袇螀芇蕿蝕聿芆艿蒃肅芅蒁螈羈芄薃薁袆芃芃螆螂芃蒞蕿肁節(jié)蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈螞螂羂薁袈肀肁芀蟻羆肁莃袆袂肀薅蠆袈聿蚇蒂膇肈莇螇肅肇葿薀罿肆