【正文】 :l for j=1:nalpha P(i,j+i1)=a(j)。 endendfor k=1:Nsb grd=grd+exp((real(w39。*P*Ibq(:,k)))^2/(2*tao^2*w39。*w))*((real(w39。*P*Ibq(:,k))*w)/(w39。*w)P*Ibq(:,k))。endgrd=grd/(2*Nsb*sqrt(2*pi)*tao*sqrt(w39。*w))。endfunction [x,output] = cg(fun,dfun,x0,varargin)epsi=。k=0。iterm=0。funcN=0。rh0=。l=。u=。%線搜索x=x0。f=feval(fun,x,varargin{:})。funcN=funcN+1。n=length(x0)。pe=zeros(1,5000)。pe(1)=f。%%%g=feval(dfun,x,varargin{:})。%梯度函數(shù)while norm(g)epsi amp。amp。 k=10000 iterm=kfloor(k/(1000*n+1))*(1000*n+1)。%重置位 iterm=iterm+1。 k=k+1。 if iterm==1 d=g。 gd=g39。*d。 else beta=g39。*g/(gold39。*gold)。 d=g+beta*dold。 gd=g39。*d。 if gd= d=g。gd=g39。*d。 end end%%%[alpha,funcNk]=armijo(fun,x,d,gd,f)。funcN=funcN+funcNk。%%% s=alpha*d。 x=x/(sqrt(x39。*x))。%%%歸一化 x=x+s。 f=feval(fun,x,varargin{:})。 if(k5000) pe(k+1)=f。 end funcN=funcN+1。 gold=g。 dold=d。 g=feval(dfun,x,varargin{:})。end=pe。=f。=k。=funcN。=norm(g)。endfunction [ x,output] =bfgs1( fun,dfun,x0,varargin )%初始化epsi=。k=0。funcN=0。rh0=。l=。u=。x=x0。f=feval(fun,x,varargin{:})。funcN=funcN+1。n=length(x0)。H=eye(n)。pe=zeros(1,5000)。pe(1)=f。%%%g=feval(dfun,x,varargin{:})。while norm(g)epsi amp。amp。k=9000 itercon=true。 d=H*g。 %%% gd=g39。*d。 [alpha,funcNk,exitflag]=armijo(fun,x,d,gd,f)。 funcN=funcN+funcNk。 if exitflag==1 itercon=false。restart=true。 H=eye(n)。gold=g。 end if itercon s=alpha*d。 x=x/(sqrt(x39。*x))。%%%歸一化 x=x+s。 f=feval(fun,x,varargin{:})。 funcN=funcN+1。 gold=g。 g=feval(dfun,x,varargin{:})。 %%% y=ggold。 hy=H*y。 sy=s39。*y。 yhy=y39。*hy。 if sy*yhy theta=*yhy/(yhysy)。 s=theta*s+()*hy。 sy=*yhy。 end v=sqrt(yhy)*(s/syhy/yhy)。 H=H+s*s39。/syhy*hy39。/yhy+v*v39。 k=k+1。 if(k5000) pe(k+1)=f。 end endend=pe。=f。=k。=funcN。=norm(g)。end%%牛頓法與共軛法的比較load pe_cg。load pe_bfgs1。load pe_steep。k=0:4999。y1=10*log10()*ones(1,5000)。y2=10*log10()*ones(1,5000)。plot(k,10*log10(pe_steep),39。g39。,k,10*log10(pe_cg),39。r39。,k,10*log10(pe_bfgs1),39。k39。)。axis([0 400 70 0])。xlabel(39。迭代次數(shù)i39。)。ylabel(39。10*log10(BER)39。)。hold on。plot(k,y1,39。:39。,k,y2,39。.39。)。legend(39。最速下降法39。,39。共軛梯度法39。,39。擬牛頓法39。,39。MBER39。,39。MMSE39。)。%%%共軛梯度法的迭代情況clearclcload pe_1。%%load pe_2。%%load pe_5。%%load pe_cg。pe_cg=10*log10(pe_cg)。k=0:4999。y1=10*log10()*ones(1,5000)。y2=10*log10()*ones(1,5000)。%plot(k,pe_1,39。c39。,k,pe_5,39。r39。,k,pe_2,39。k39。)。plot(k,pe_1,39。r39。,k,pe_cg,39。k39。)。axis([0 4000 70 0])。xlabel(39。迭代次數(shù)i39。)。ylabel(39。10*log10(BER)39。)。hold on。plot(k,y1,39。k:39。,k,y2,39。.39。)。%legend(39。39。,39。39。,39。39。,39。MBER39。,39。MMSE39。)。legend(39。固定步長39。,39。不精確線搜索39。,39。MBER39。,39。MMSE39。)。function [alpha,funcNk,exitflag] = armijo( fun,xk,dk,gd,f)beta=。sigma=。m=0。mmax=20。funcNk=0。exitflag=1。mk=0。while(mmmax) fn=feval(fun,xk+beta^m*dk)。 if(fn=(f+sigma*beta^m*gd)) mk=m。exitflag=0。 break。 endm=m+1。endfuncNk=mk+1。alpha=beta^mk。end%%最陡下降法function [x,output] = steep(fun,dfun,x0,varargin)epsi=。k=0。iterm=0。funcN=0。x=x0。f=feval(fun,x,varargin{:})。funcN=funcN+1。n=length(x0)。pe=zeros(1,5000)。pe(1)=f。%%%g=feval(dfun,x,varargin{:})。%梯度函數(shù)while norm(g)epsi amp。amp。 k=10000 iterm=kfloor(k/(1000*n+1))*(1000*n+1)。%重置位 iterm=iterm+1。 k=k+1。 if iterm==1 d=g。 gd=g39。*d。 else d=g。 gd=g39。*d。 end%%%[alpha,funcNk]=armijo(fun,x,d,gd,f)。funcN=funcN+funcNk。%%% s=alpha*d。 x=x/(sqrt(x39。*x))。%%%歸一化 x=x+s。 f=feval(fun,x,varargin{:})。 if(k5000) pe(k+1)=f。 end funcN=funcN+1。 g=feval(dfun,x,varargin{:})。end=pe。=f。=k。=funcN。=norm(g)。endfunction I = s_function(d,m )I=de2bi(0:2^m1)39。b=zeros(1,2^(m1))。g=1。for k=1:2^m if I(d+1,k)==1 b(g)=k。g=g+1。 endendI=I(:,b)。for b=1:2^(m1) for g=1:m if I(g,b)==0 I(g,b)=1。 end endendendfunction [ y ] = q( x )y=erfc(x./sqrt(2))./2。endfunction w = mmse()%global tao。tao=。d=3。l=5。%濾波器長度nalpha=3。%輸入信號濾波器長度,信道脈沖響應(yīng)長度m=l+nalpha1。%用戶數(shù)量，對于單用戶是延時長度d=3。%延時Nb=2^m。%bq所有情況數(shù)目Nsb=Nb/2。Ibq=s_function(d,m)。%生成bq矩陣a=[ ]39。%信道脈沖響應(yīng)P=zeros(l,m)。%生成系統(tǒng)矩陣for i=1:l for j=1:nalpha P(i,j+i1)=a(j)。 endendw=inv(P*P39。+2*tao^2*eye(l))*P(:,d+1)。end