⭐ 欢迎来到虫虫下载站! | 📦 资源下载 📁 资源专辑 ℹ️ 关于我们
⭐ 虫虫下载站

📄 rereiver.cpp

📁 单用户多载波mccdma的仿真
💻 CPP
📖 第 1 页 / 共 2 页
字号:
  //cout<<min<<" ";
 if(min==dic[0]) 
	 fI[k][i]=1*sqrt(2)/2,fQ[k][i]=1*sqrt(2)/2;
 else if(min==dic[1]) 
	 fI[k][i]=-1*sqrt(2)/2,fQ[k][i]=1*sqrt(2)/2;
 else if(min==dic[2]) 
	 fI[k][i]=1*sqrt(2)/2,fQ[k][i]=-1*sqrt(2)/2;
 else if(min==dic[3]) 
	 fI[k][i]=-1*sqrt(2)/2,fQ[k][i]=-1*sqrt(2)/2;
 //cout<<min<<" ";
}
//cout<<"\n";
}
  /////////////////////////////////////////////
 /*  cout<<"解调"<<"\n";
   cout<<"I支路"<<"\n";
for(i=0;i<factor;i++)
  
{ for(j=0;j<numberofuserbit/2+numberofuserbit%2;j++) 
    cout<<fI[i][j]<<" ";
    cout<<"\n";
} 
 cout<<"Q支路"<<"\n";
for(i=0;i<factor;i++)
{ for(j=0;j<numberofuserbit/2+numberofuserbit%2;j++) 
    cout<<fQ[i][j]<<" ";
    cout<<"\n";
} */
////////////////////合并IQ两路///////////////////////////////

for(k=0;k<factor;k++)
  for(i=1;i<numberofuserbit/2+numberofuserbit%2+1;i++)
  {
	  //cout<<int(fI[k][i-1]*sqrt(2)+2*fQ[k][i-1]*sqrt(2))<<" ";
   
	  if(int(fI[k][i-1]*sqrt(2)+2*fQ[k][i-1]*sqrt(2))==3)  Demodulation[k][i*2-2]=-1,Demodulation[k][i*2-1]=-1;
         
	  else if(int(fI[k][i-1]*sqrt(2)+2*fQ[k][i-1]*sqrt(2))==1) Demodulation[k][i*2-2]=-1,Demodulation[k][i*2-1]=1;
	  else if(int(fI[k][i-1]*sqrt(2)+2*fQ[k][i-1]*sqrt(2))==-1) Demodulation[k][i*2-2]=1,Demodulation[k][i*2-1]=-1;
	  else if(int(fI[k][i-1]*sqrt(2)+2*fQ[k][i-1]*sqrt(2))==-3) Demodulation[k][i*2-2]=1,Demodulation[k][i*2-1]=1;
	  //finde[k][i*2-2]=fI[k][i-1]*sqrt(2); 
      //finde[k][i*2-1]=fQ[k][i-1]*sqrt(2);  
   //cout<<finde[k][i*2-2]<<" ";

 }
  //cout<<"\n";
  /*cout<<"解调合并"<<"\n";
for(i=0;i<Fspreadfactor;i++)
{ for(j=0;j<numberofuserbit;j++) 
    cout<<Demodulation[i][j]<<" ";
 cout<<"\n";
}*/
for(i=0;i<numberofuserbit;i++)
{for(k=1;k<factor;k++) 
{   
//	Demodulation[0][i]= Despreadout[0][i]+Despreadout[k][i];
 Demodulation[0][i]+=Demodulation[k][i];
}
if(Demodulation[0][i]>0 ) Demodulation[0][i]=1;
else Demodulation[0][i]=0;

}
/*for(j=0;j<numberofuserbit;j++) 
    cout<<Demodulation[0][j]<<" ";*/
DeAllocate2DArray(fI,x1);                                   //释放内存
DeAllocate2DArray(fQ,x1);                                   //释放内存
}
/////////////////////////////////////////////////////////
//                                                  /////
//                    产生m序列扩频码             ///////
//////////////////////////////////////////////////////////
	void receiver::M_seq(int n)
    {
		int zita[8200],i,k,N;
		N=int(pow(2,n)-1);
	for(i=0;i<n;i++)                    //对zita[k]赋初值
	{
      zita[i]=1;
	   Mseq[i]=1;
    }
	
	switch(n)
	{
	case 2: 
		    for(k=n+1;k<=2041;k++)
			{
               zita[k]=zita[k-n]^zita[k-1];	
               Mseq[k]=(2*zita[k]-1);
			 // cout<<Mseq[k]<<"  ";
			 
            }
		    break;
case 3: 
		    for(k=n;k<=2041;k++)
			{
               zita[k]=zita[k-n]^zita[k-2];	
               Mseq[k]=(2*zita[k]-1);
			   //cout<<Mseq[k]<<"  ";
			    
            }
		    break;
	case 4:
		    for(k=n;k<=2041;k++)
			{
               zita[k]=zita[k-n]^zita[k-3];	
               Mseq[k]=(2*zita[k]-1);
			  // cout<<Mseq[k]<<"  ";
			   
            }
		    break;
	case 5:
		    for(k=n+1;k<=2041;k++)
			{
               zita[k]=zita[k-n]^zita[k-3];	
               Mseq[k]=(2*zita[k]-1);
			  // cout<<Mseq[k]<<"  ";
			  
            }
		    break;
	case 6:
		     for(k=n+1;k<=2041;k++)
			{
               zita[k]=zita[k-n]^zita[k-5];	
               Mseq[k]=(2*zita[k]-1);
			  // cout<<Mseq[k]<<"  ";
			    
            }
		    break;
	case 7:
		    for(k=n+1;k<=2041;k++)
			{
               zita[k]=zita[k-n]^zita[k-4];	
               Mseq[k]=(2*zita[k]-1);
			  // cout<<Mseq[k]<<"  ";
			    
            }
		    break;
	case 9: 
		     for(k=n+1;k<=2041;k++)
			{
               zita[k]=zita[k-n]^zita[k-5];	
               Mseq[k]=(2*zita[k]-1);
			  // cout<<Mseq[k]<<"  ";
			   
            }
		    break;
	case 10: 
		     for(k=n+1;k<=2041;k++)
			{
               zita[k]=zita[k-n]^zita[k-7];	
               Mseq[k]=(2*zita[k]-1);
			 //  cout<<Mseq[k]<<"  ";
			   
            }
		    break;
	case 11:
		    for(k=n+1;k<=2041;k++)
			{
               zita[k]=zita[k-n]^zita[k-9];	
               Mseq[k]=(2*zita[k]-1);
			 //  cout<<Mseq[k]<<"  ";
			   
            }
		    break;   
	case 8:
		    for(k=n+1;k<=2047;k++)
			{
               zita[k]=zita[k-n]^zita[k-7]^zita[k-6]^zita[k-1];
               Mseq[k]=(2*zita[k]-1);
			 //  cout<<Mseq[k]<<"  "; 
			    
			}
			break;
	case 12:
		    for(k=n+1;k<=N;k++)
			{
               zita[k]=zita[k-n]^zita[k-8]^zita[k-5]^zita[k-1];
               Mseq[k]=(2*zita[k]-1);
			 //  cout<<Mseq[k]<<"  "; 
			   
			}
			break;
	case 13:
		    for(k=n+1;k<=N;k++)
			{
               zita[k]=zita[k-n]^zita[k-8]^zita[k-6]^zita[k-1];
               Mseq[k]=(2*zita[k]-1);
			//   cout<<Mseq[k]<<"  "; 
			    
			}
			break;
	default: 
		    cout<<"请重新输入2到13之间的数:";
			}
	}
	/////////////////////////////////////////////////////////
//                                                  /////
//                    进行频域解扩                  ///////
//////////////////////////////////////////////////////////
void receiver::FdomainDespread(int factor,double **a)
{
	int i,j,k,**sp;
//sp=intAllocate2DArray(x1,1000);                       //动态2维数组分配内存
for(k=0;k<factor;k++) 
for(i=0;i<numberofuserbit/2+numberofuserbit%2;i++)
{ Despreadout[k][i]=a[k][i]*Mseq[k];
  
  //cout<<Despreadout[k][i];
}
 //for(k=0;k<n2-1;k++) 
//for(j=0;j<10;j++)
 //for(i=0;i<r;i++)
   // sp2[k][j]=sp2[k][j]+sp[k][i+r*j];
 //for(i=0;i<10;i++)
//	 if(de2[i]>0) de2[i]=1;
//	 else de2[i]=0;
////////////////////////
/*  cout<<"解扩"<<"\n";
 for(i=0;i<factor;i++)
{ for(j=0;j<numberofuserbit;j++) 
    cout<<sp[i][j]<<" ";
 cout<<"\n";
}*/
 //////////////////
/*for(i=0;i<numberofuserbit;i++)
{Despreadout[i]=0;                              
for(k=0;k<factor;k++) 
{   
	Despreadout[i]= Despreadout[i]+sp[k][i];
//cout<<Despreadout[i];
}
}
intDeAllocate2DArray(sp,x1);                                   //释放内存
for(i=0;i<numberofuserbit;i++)
if( Despreadout[i]>0) Despreadout[i]=1;
   else  Despreadout[i]=0;
*/
}

// --------------------------------------------------------------------------
// 动态分配二维数组
// 使用方法(以int类型8行10列为例): int **dyarr = Allocate2DArray<int>(8,10)
//template<double>
 double ** receiver::  Allocate2DArray(size_t pm_Rows, size_t pm_Cols)
{
    double **arr = new double*[pm_Rows];
    for (int i = 0; i < pm_Rows; i++)
        arr[i] = new double[pm_Cols];
    return arr;
}

// --------------------------------------------------------------------------
// 释放动态分配的二维数组
// 使用方法(以int类型8行10列为例): Allocate2DArray<int>(dyarr,8)
//template<double>
void receiver::DeAllocate2DArray(double **arr, size_t pm_Rows)
{
    for (int i = 0; i < pm_Rows; i++)
        delete[] arr[i];
    delete[] arr;
}

int ** receiver::intAllocate2DArray(size_t pm_Rows, size_t pm_Cols)
{
    int **arr = new int*[pm_Rows];
    for (int i = 0; i < pm_Rows; i++)
        arr[i] = new int[pm_Cols];
    return arr;
}

// --------------------------------------------------------------------------
// 释放动态分配的二维数组
// 使用方法(以int类型8行10列为例): Allocate2DArray<int>(dyarr,8)
//template<double>
void receiver::intDeAllocate2DArray(int **arr, size_t pm_Rows)
{
    for (int i = 0; i < pm_Rows; i++)
        delete[] arr[i];
    delete[] arr;
}

⌨️ 快捷键说明

复制代码 Ctrl + C
搜索代码 Ctrl + F
全屏模式 F11
切换主题 Ctrl + Shift + D
显示快捷键 ?
增大字号 Ctrl + =
减小字号 Ctrl + -