燕山大学数字通信计算机仿真课程设计报告模板.doc

. .

数字通信计算机课设  程序代码 
// SigTranmit.cpp: implementation of the SigTranmit class.
 // 
//////////////////////////////////////////////////////////////////////  
#include "stdafx.h" 
#include "ASK.h" 
#include "SigTranmit.h"
 #include <math.h>  
#ifdef _DEBUG
 #undef THIS_FILE 
static char THIS_FILE[]=__FILE__; 
#define new DEBUG_NEW 
#endif  
////////////////////////////////////////////////////////////////////// 
// Construction/Destruction 
//////////////////////////////////////////////////////////////////////  
SigTranmit::SigTranmit() 
{  c=0; 
/* S0[16]={0}; 
 float S1[16]={0};
Mod[16*M]={0}; 
 Sa[16*M]={0}; 
 Sal[16*M]={0};   
 H1[16*M]={0}; 
 H2[16*M]={0}; 
 No[16*M]={0};  
Hb[N]={0}; 
 P[N]={0}; 
 Wn[N]={0};*/
 }  
SigTranmit::~SigTranmit() 
{ 
 } 
void SigTranmit::SetDigital(int ID,int Intensity)//赋值 
{ 
   myID=ID; 
   myIntensity=Intensity;
 } 
void SigTranmit::Basesignal()//基带信号
 { 
 //ID变成int型数组
int a[4] = {0}; 
 int i=0;  
int k=15;
  int temp=0;  
a[0]=myID%10; 
 a[1]=myID/10%10; 
 a[2]=myID/100%10; 
 a[3]=myID/1000;  
for(int j=0;j<4;j++)  
{   i=0; 
  temp=a[j]; 
  for(i;i<4;i++)   {    S0[k]=temp & 0x1;    temp = temp>>1;    k--;   }  }  //采样  k=0;  for(i=0;i<16;i++)  {   for(j=0;j<M;j++)   {    Sal[k]=S0[i];    k++;   }  }  for(i=0;i<16*M;i++)  {   if(Sal[i] == 0)    Sa[i] = -1;   else    Sa[i] = Sal[i];  }  Dft(Sa,16*M);  }  

void SigTranmit::Lpf()//低通 {  Conv(Sa,GetH1());  Dft(Sa,16*M);  } 
void SigTranmit::Dft(float x[],int m)//DFT变换 {  for(int k=0;k<m;k++)  {   AR[k]=0;   AI[k]=0;   for(int i=0;i<m;i++)   {    AR[k] =AR[k] + x[i]*cos(2*PI*k*i/m);     AI[k] =AI[k] + x[i]*sin(2*PI*k*i/m);   }   A[k] = sqrt(AR[k]*AR[k] + AI[k]*AI[k]);  }  } 
void SigTranmit::modulate()//调制 {  float a=0.0;  for(int i=0;i<16*M;i++)  {   a=2.0*PI*i*0.15;   Mod[i] = Sa[i]*cos(a);  }  for(i=0;i<16*M;i++)   Sa[i] = Mod[i];  Dft(Sa,16*M); } 
void SigTranmit::Channel()//信道 {  Noise(myIntensity);  for(int i=0;i<16*M;i++)   Sa[i] = No[i] + Sa[i];  Dft(Sa,16*M); } 
void SigTranmit::Noise(int inten)//噪声 {  int i,k;  float r0=0.0;  float r[12]; 

 for(k=0; k<16*M; k++)  {   for(i=0; i<12; i++)   {    r[i] = (float)(2*rand()-RAND_MAX)/RAND_MAX;    r0 += r[i];    r0 = float(inten*r0/12.0);   }   No[k] = r0;  }  Dft(No,16*M); } 
void SigTranmit::Bpf()//带通 {  int i; 
    float a=(N-1)/2;   //生成理想带通滤波器的单位取样响应  for(i=0;i<N;i++)  {   if(i==a)     Hb[i]=(float)0.4;   else     Hb[i]=sin(0.7*PI*(i-a))/(float)(PI*(i-a))-sin(0.3*PI*(i-a))/(float)(PI*(i-a));  }   switch(c)  {  case 0:   {for(i=0;i<N;i++)    Wn[i]=0.5*(1.0-cos(2.0*PI*i/(N-1)));};   break;  case 1:   {    for(i=0;i<N;i++)     Wn[i]=0.54-0.46*cos(2.0*PI*i/(N-1));};    break;  case 2:  //生成布莱克曼窗   {    for(i=0;i<N;i++)     Wn[i]=0.42-0.5*cos(2.0*PI*i/(N-1))+0.08*cos(4.0*PI*i/(N-1));   };  }  for( i=0;i<N;i++)  {   Hb[i]=Hb[i]*Wn[i]; 

 }  Conv(Sa,GetHb());  Dft(Sa,16*M);  } 
void SigTranmit::Demodulate()//解调 {  float a=0.0;  for(int i=0;i<16*M;i++)  {   a=2.0*PI*i*0.15;   Mod[i] = Sa[i]*cos(a);  }  for(i=0;i<16*M;i++)   Sa[i] = Mod[i];  Dft(Sa,16*M); } 
void SigTranmit::Adjust()//判决 {  int k=0;  int j=0;  for(int i=0;i<16;i++)  {   S1[i]=0;   for(j=0;j<M;j++)   {    S1[i] += Sa[k];    k++;   }  }  k=0;  for(i=0;i<16;i++)  {   for(j=0;j<M;j++)   {    if(S1[i]>0)    {     Sa[k]=1;     k++;    }    else    {     Sa[k]=-1;     k++; 

. .word.zl.