武汉理工大学本科生毕业设计（论文）
out(1,:)=real(ifft(dopplertf.*fft(input(1,:)))); out(2,:)=real(ifft(dopplertf.*fft(input(2,:))));
a(:,ptstart:ptend)=a(:,ptstart:ptend)+out; end clear j
alpha=a(1,1:N)+j.*a(2,1:N);
程序4
N_Trials=100; N_number=100; N_snr=10; Q=16; E_M=[];
N=3; % 3 users
bit_rate=100000; Fd_Tb=0.001;
% fd is the maximum doppler frequency y=v/lamda y_Tb=Fd*Tb
M_path=2;
for trials=1:N_Trials trials
%theta=round(rand(1,N))*pi-(1/2)*pi; %a=exp(-j.*L*pi*sin(theta))./sqrt(M);
noise=randn(1,Q*N_number)+j.*randn(1,Q*N_number); s10=round(rand(N,N_number)); pn01=round(rand(N,Q)); pn=(pn01.*2-1)./sqrt(Q); ss=s10*2-1;
phase=exp(-j.*(rand(1,N*M_path)*2*pi)); ph=diag(phase);%相移
sigf_coe=[];
29
武汉理工大学本科生毕业设计（论文）
for k=1:N
for path=1:M_path
[fad,fd,D_Tc]=flat1(N_number,bit_rate,Fd_Tb); SRaley=ss(k,:).*abs(fad);%加衰落系数 sigf_coe=[sigf_coe;SRaley]; end end
spread=[]; for k=1:N
for path=1:M_path
spath=(k-1)*M_path+path; s=kron(sigf_coe(spath,:),pn(k,:)); spread=[spread;s]; end end
delay=fix(rand(N,M_path)*(Q-1)); delay(1,1)=0;%时延
spreadnew=[]; for k=1:N
for path=1:M_path d_path=delay(k,path); spath=(k-1)*M_path+path; s_spath=spread(spath,:);
kp=[zeros(1,d_path),s_spath(1,1:Q*N_number-d_path)]; spreadnew=[spreadnew;kp]; end end
sgma=1; Error_v=[];
for snr_db=0:1:N_snr
snr_v0=[snr_db,5*ones(1,N-1)]; %N-1 interfers with INR=5dB snr_v=kron(snr_v0,ones(1,M_path)); snr=10.^(snr_v./10)/M_path; %uni-power N0=2*sgma.^2; Eb_v=snr.*N0; G=diag(sqrt(Eb_v));
if N*M_path==1
yy=ph*G*spreadnew+noise;
30
武汉理工大学本科生毕业设计（论文）
else
yy=sum(ph*G*spreadnew)+noise; end Y_V=[];
for path=1:M_path; d=delay(1,path);
y_path=[yy(1,d+1:N_number*Q),zeros(1,d)];% 对齐时延信号相加 Y_M=[]; for k=1:N_number
ym=y_path(1,(k-1)*16+1:k*16); Y_M=[Y_M;ym];%拓展 end
ys=Y_M*pn(1,:).'; y=ys.';
Y_V=[Y_V;y];%解扩 end
% for path=1:M_path
phase0=phase(1,1:M_path)'; ys=diag(phase0)*Y_V; %end
if path==1 y=ys; else
y=sum(ys); end
y_real=real(y); s_e=sign(y_real); s_e10=(s_e+1)./2;
Error_snr=sum(abs(s10(1,:)-s_e10(1,:))); Error_v=[Error_v,Error_snr]; end
E_M=[E_M;Error_v] ; end
E_sum=mean(E_M); BER=E_sum./N_number; BER_T=[];
31
武汉理工大学本科生毕业设计（论文）
for snr_db=0:1:N_snr snr=10.^(snr_db./10);
?R_THEROY=Qfunct(sqrt(2.*snr)); BER_THEROY=(1/2)*[1-sqrt(snr/(snr+1))]; BER_T=[BER_T,BER_THEROY]; end i=0:1:10;
semilogy(i,BER,'-r',i,BER_T ,'*g')
32
武汉理工大学本科生毕业设计（论文）
致 谢
四年的读书生活在这个季节即将划上一个句号，而于我的人生却只是一个逗号，我将面对又一次征程的开始。四年的求学生涯在师长、亲友的大力支持下，走得辛苦却也收获满囊，在论文即将付梓之际，思绪万千，心情久久不能平静。 伟人、名人为我所崇拜，可是我更急切地要把我的敬意和赞美献给一位平凡的人，我的导师李平安。我不是您最出色的学生，而您却是我最尊敬的老师。您治学严谨，学识渊博，思想深邃，视野雄阔，为我营造了一种良好的精神氛围。授人以鱼不如授人以渔，置身其间，耳濡目染，潜移默化，使我不仅接受了全新的思想观念，树立了宏伟的学术目标，领会了基本的思考方式，从论文题目的选定到论文写作的指导,经由您悉心的点拨,再经思考后的领悟,常常让我有―山重水复疑无路,柳暗花明又一村‖。
在论文即将完成之际，我的心情无法平静，从开始进入课题到论文的顺利完成，有多少可敬的师长、同学、朋友给了我无言的帮助，在这里请接受我诚挚谢意！
同时也感谢学院为我提供良好的做毕业设计的环境。
最后再一次感谢所有在毕业设计中曾经帮助过我的良师益友和同学，以及在设计中被我引用或参考的论著的作者。
韩立鸿
2010年4月10日
33