user7chp.m

rake接收机的MATLAB仿真程序 希望能对大家有些帮助

                function output=user7chip()

%//input siganl form transmitter side for chip based\\
                  
clear all;
close all;
clc;
us=3;
% nb1=10;
% nb2=10;
% nb3=10;
% nb4=10;
% nb5=10;
% nb6=10;
% nb7=10;
% m1=[1 0 0 1 0 1 0 1 0 0 ];
% m2=[1 0 0 0 1  0 1 1 0 0];
% m3=[1 0 1 0 0 1 0 0 1 1];
% m4=[1 1 0 0 0  1 0 0 0 1];
% m5=[1 0 0 1 0  0 1 1 0 0];
% m6=[1 0 1 0 0  1 0 0 0 1 ];
% m7=[1 0 0 0 1 0 1 0 0 0];
% nb1=20;
% nb2=20;
% nb3=20;
% nb4=20;
% nb5=20;
% nb6=20;
% nb7=20;
% m1=[1 1 0 0 0 0 0 0 0 0 1 1 1 0 1 0 1 1 0 0 ];
% m2=[1 0 1 0 1  0 1 1 0 1 0 0 0 0 1 0 1 0 1 0];
% m3=[1 0 1 0 1 1 0 1 1 1 1 0 1 0 1 0 1 0 1 1 ];
% m4=[1 1 0 1 0  1 0 1 0 1 0 0 0 0 1 1 0 1 0 1];
% m5=[1 0 1 1 0  0 1 1 0 1 1 0 1 1 1 0 0 1 0 0];
% m6=[1 0 1 0 1  1 0 1 0 1 1 1 0 0 1 0 1 0 1 0];
% m7=[1 1 0 1 0 0 1 0 1 0 1 0 1 0 1 1 0 1 1 0];
% m1=[1 1 0 0 0 0 0 0 0 0 1 1 1 0 1 0 1 1 0 0 ];
% m2=[1 0 1 0 1  0 1 1 0 1 0 0 0 0 1 0 1 0 1 0];
% m3=[1 0 1 0 1 1 0 1 0 0 0 0 1 0 1 0 1 0 1 1 ];
% m4=[1 1 0 1 0  1 0 1 0 1 0 0 0 0 1 1 0 1 0 1];
% m5=[1 0 1 1 0  0 1 1 0 1 1 0 1 1 1 0 0 1 0 0];
% m6=[1 0 1 0 1  1 0 1 0 1 1 1 0 0 1 0 1 0 1 0];
% m7=[1 1 0 1 0 0 1 0 1 0 1 0 1 0 1 1 0 1 1 0];

nb1=20;
nb2=20;
nb3=20;
nb4=20;
nb5=20;
nb6=20;
nb7=20;
m1=[1 1 1 1 1 0 1 1 1 0 1 0 1 1 1 1 0 1 1 0];
m2=[1 1 0 1 1 1 1 0 1 1 1 0 1 1 0 1 1 1 1 0 ];
m3=[1 0 1 1 0 1 1 0 1 1 1 0 1 0 1 1 0 1 1 0];
m4=[1 1 1 1 0 1 1 1 0 1 1 0 1 1 1 1 1 1 0 1 ];
m5=[1 1 1 1 0 1 1 0 1 1 1 0 1 1 1 1 0 1 1 1 ];
m6=[1 0 1 1 0 1 0 1 1 1 1 1 1 1 0 1 1 1 1 1 ];
m7=[1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1];

                %//message1 into polar form\\
                    
for i=1:nb1
    if m1(i)==1
        mb1(i)=m1(i);
    else
        mb1(i)=-1;
    end
end 
  % disp(mb1);
   
                %//message2 into polar form\\
                
for i=1:nb2
    if m2(i)==1
        mb2(i)=m2(i);
    else
        mb2(i)=-1;
    end
end 
  %disp(mb2);
  
                %//message3 into polar form\\
                
for i=1:nb3
    if m3(i)==1
        mb3(i)=m3(i);
    else
        mb3(i)=-1;
    end
end 
  %disp(mb3);
  
    %//message4 into polar form\\
                
for i=1:nb4
    if m4(i)==1
        mb4(i)=m4(i);
    else
        mb4(i)=-1;
    end
end 
  %disp(mb4);
  
    %//message5 into polar form\\
                
for i=1:nb5
    if m5(i)==1
        mb5(i)=m5(i);
    else
        mb5(i)=-1;
    end
end 
  %disp(mb5);
  
             %//message6 into polar form\\
                
for i=1:nb6
    if m6(i)==1
        mb6(i)=m6(i);
    else
        mb6(i)=-1;
    end
end 
  %disp(mb6);
  
            %//message7 into polar form\\
                
for i=1:nb7
    if m7(i)==1
        mb7(i)=m7(i);
    else
        mb7(i)=-1;
    end
end 
  %disp(mb7);
                  %//generation of maximal length sequence1\\
                  
  f1=1;f2=0;f3=1;f4=0;f5=1;
m=(2^5)-1;
for j=1:m
    p=xor(f2,f5);
    f5=f4;
    op1(j)=f5;
    f4=f3;
    f3=f2;
    f2=f1;
    f1=p;
    
end
  %disp(op1);
  
                    %//msequence1 into polar form\\
                    
    for i=1:31
        
         if op1(i)==1
             pb1(i)=1;
         else
         pb1(i)=-1;
        end
    end 
    %disp(op1);
    %disp(pb1);
    
                %//generation of maximal length sequence2\\
      
  f1=1;f2=1;f3=1;f4=0;f5=1;
m=(2^5)-1;
for j=1:m
    p=xor(f3,f4);
    f5=f4;
    op2(j)=f5;
    f4=f3;
    f3=f2;
    f2=f1;
    f1=p;
    
end
  %disp(op2);
  
                    %//msequence2 into polar form\\
                    
    for i=1:31
        
         if op2(i)==1
             pb2(i)=1;
         else
         pb2(i)=-1;
        end
    end 
    %disp(op2);
    %disp(pb2);
    
             %//generation of maximal length sequence3\\
      
  f1=1;f2=0;f3=0;f4=1;f5=1;
m=(2^5)-1;
for j=1:m
    p=xor(f3,f5);
    f5=f4;
    op3(j)=f5;
    f4=f3;
    f3=f2;
    f2=f1;
    f1=p;
    
end
  %disp(op3);
  
       %//msequence3 into polar form\\
                    
    for i=1:31
        
         if op3(i)==1
             pb3(i)=1;
         else
         pb3(i)=-1;
        end
    end 
    %disp(op3);
    %disp(pb3);
  
   %//generation of maximal length sequence4\\
                  
  f1=1;f2=1;f3=1;f4=0;f5=1;
m=(2^5)-1;
for j=1:m
    p=xor(f2,f5);
    f5=f4;
    op4(j)=f5;
    f4=f3;
    f3=f2;
    f2=f1;
    f1=p;
    
end
  %disp(op4);
  
                    %//msequence4 into polar form\\
                    
    for i=1:31
        
         if op4(i)==1
             pb4(i)=1;
         else
         pb4(i)=-1;
        end
    end 
    %disp(op4);
    %disp(pb4);
    
             %//generation of maximal length sequence5\\
                  
  f1=0;f2=1;f3=1;f4=0;f5=1;
m=(2^5)-1;
for j=1:m
    p=xor(f2,f5);
    f5=f4;
    op5(j)=f5;
    f4=f3;
    f3=f2;
    f2=f1;
    f1=p;
    
end
  %disp(op5);
  
                    %//msequence5 into polar form\\
                    
    for i=1:31
        
         if op5(i)==1
             pb5(i)=1;
         else
         pb5(i)=-1;
        end
    end 
    %disp(op5);
    %disp(pb5);
               
           % //generation of maximal length sequence6\\
  f1=1;f2=1;f3=1;f4=1;f5=1;
m=(2^5)-1;
for j=1:m
    p=xor(f1,f5);
    f5=f4;
    op6(j)=f5;
    f4=f3;
    f3=f2;
    f2=f1;
    f1=p;
    
end
  %disp(op6);
  
                    %//msequence6 into polar form\\
                    
    for i=1:31
        
         if op6(i)==1
             pb6(i)=1;
         else
         pb6(i)=-1;
        end
    end 
    %disp(op6);
    %disp(pb6);
    
           % //generation of maximal length sequence7\\              
  f1=0;f2=1;f3=1;f4=0;f5=1;
m=(2^5)-1;
for j=1:m
    p=xor(f3,f5);
    f5=f4;
    op7(j)=f5;
    f4=f3;
    f3=f2;
    f2=f1;
    f1=p;
    
end
  %disp(op7);
  
                    %//msequence7 into polar form\\
                    
    for i=1:31
        
         if op7(i)==1
             pb7(i)=1;
         else
         pb7(i)=-1;
        end
    end 
    %disp(op7);
    %disp(pb7);
               
    
               
         %//first transmit bit\\
         
k=1;
for i=1:nb1
    for j=1:31
        tb1(k)=mb1(i)*pb1(j);
        k=k+1;
    end
end
%disp(tb1);

       %//second transmit bit\\
k=1;
for i=1:nb2
    for j=1:31
        tb2(k)=mb2(i)*pb2(j);
        k=k+1;
    end
end
%disp(tb2);
            
         %//third transmit bit\\
         
k=1;
for i=1:nb3
    for j=1:31
        tb3(k)=mb3(i)*pb3(j);
        k=k+1;
    end
end
%disp(tb3);
 
%// fourth transmit bit\\
         
k=1;
for i=1:nb4
    for j=1:31
        tb4(k)=mb4(i)*pb4(j);
        k=k+1;
    end
end
%disp(tb4);

        %//fifth transmit bit\\
         
k=1;
for i=1:nb5
    for j=1:31
        tb5(k)=mb5(i)*pb5(j);
        k=k+1;
    end
end
%disp(tb5);

   %//sixth transmit bit\\
         
k=1;
for i=1:nb6
    for j=1:31
        tb6(k)=mb6(i)*pb6(j);
        k=k+1;
    end
end
%disp(tb6);

   %//seventh transmit bit\\
         
k=1;
for i=1:nb7
    for j=1:31
        tb7(k)=mb7(i)*pb7(j);
        k=k+1;
    end
end
%disp(tb7);

      
     %//spreading of three  signals transmitted in the channel\\
       
n=1;
for i=1:k-1
    tb(n)=tb1(n)+tb2(n)+tb3(n)+tb4(n)+tb5(n)+tb6(n)+tb7(n);
    n=n+1;
end
%disp(tb);
 figure(1);
grid on;
subplot(2,2,1);
plot(tb);
title('Received signal with out noise for CB');
xlabel('Time');
ylabel('Amplitude');
g=awgn(tb,10);
%disp(g);
subplot(2,2,2);
plot(g);
title('Received signal with noise for CB');
xlabel('Time');
ylabel('Amplitude');

                    %//optimum weight for CB configuration\\
                    
%q=input('enter the no of bits');
q1=length(tb);
s=randsrc(1,q1);
c=complex(tb,s);
%disp(c);
d=conj(c);
d1=d*d';
%disp(d1);
m=mean(d1);
lamda=170;
teta=45;
k=((2*pi)/lamda);
K=50;
for i=0:K-1
    v=(exp(j*k*i*d*sin(teta)))';
end             
%disp(v);
v1=conj(v)';
teta=30;
k=((2*pi)/lamda);
K=5;
for i=0:K-1
    eta=exp(sqrt(-1)*k*i*d*sin(teta))';
end  
%u1=[2+2j 3-1j 4+3j 5+2j]'
%m=randsrc(1,124);
%a=randsrc(1,124);
%l=input('enter the no of bits');
l1=length(tb);
u1=complex(tb,l1);
u=eta*u1;
u2=conj(u)';
u3=u*u2;
Ru=mean(u3);

%Ru1=imresize(Ru,[124 124]);
%a1=inv(Ru1);
%a1=Ru1';
a1=Ru';
%disp(a1);
%calculating beta value
%g=1;
dem=v*v1*a1;
beta=dem.^-1;
 %optimum weight
%Wopt=beta*a1*v';
temp=beta*v';
Wopt=temp*a1;
length(Wopt);
%disp('optimum weights for CB')
%disp(Wopt);
subplot(2,2,3);
plot(Wopt);
title('Received signal with weight for CB');
xlabel('Time');
ylabel('weight');

                 %//beamformer output\\
                 
o=conj(Wopt)*c;
%disp('Beamformer output for CB')
%disp(o);
subplot(2,2,4);
plot(o);
title('Beam former output with weight for CB');
xlabel('Time');
ylabel('weight');

   %//%receiver side\\
        %//reconstruction of user2 signal\\
        
s2=0;
t2=1;
for i=0:31:k-32
    for j=1:31
        ob2(t2)=o(i+j)*pb2(j)+s2;
        s2=ob2(t2);
        if ob2(t2)>0
            ob2(t2)=1;
        else ob2(t2)=0;
        end
    end
    t2=t2+1;
    s2=0;
end
%disp('second message is')
%disp(ob2);

      %//reconstruction of received signal from noise\\
            
for i=1:nb2
     if g(i)>0
            g2(i)=1;
        else g2(i)=0;
        end
    end
    
% disp(g2);

                  %//SINR CALCULATION FOR CB CONFIGURATION\\
hop1=conj(pb2)';
hx=conj(c)';
nr=hop1*hx';
 sum=nr*Wopt; 
 %disp('sum');
 %disp(sum);
 cm=sum'*sum;
 signal=mean(cm);
 % disp('cm');
 %disp(cm);
 hop2=conj(pb2)';
hx1=conj(g)';
nr1=hop2*hx1';
 sum1=nr1*Wopt; 
 %disp('sum');
 %disp(sum);
  cm1=sum1'*sum1;
  noise=mean(cm1);
 %disp('cm1');
 %disp(cm1);
 CSINR=signal/noise;
 length(CSINR);
 disp('CSINR');
 disp(CSINR);
  
                 %//GRAPH FOR DOA VERSUS SINR\\
              
     DOA=-80:20:80;
     figure(2);
     grid on;
     plot(DOA,CSINR,'R-*');
     title('simulation for DOA versus CSINR');
     xlabel('DOA in degree');
     ylabel('CSINR in db');
     
                    %// BER CALCULATION FOR CB CONFIGURATION\\
                    
    k=biterr(m2,g2);
    output=k;
    disp('biterror rate for cb configuration');
    disp(k);
    
     %//GRAPH FOR NUMBER USERS VERSUS BER\\
            
     figure(3);
     plot(us,k,'R-*');
     
     title('simulation for No. of users versus BER CHIP BASED ');
     xlabel('NUMBER USERS');
     ylabel('BER');