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

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

数字通信计算机课设  程序代码 
// SigTranmit.cpp: implementation of the SigTranmit class.
 // 
//////////////////////////////////////////////////////////////////////  
#include "stdafx.h" 
#include "ASK.h" 
#include "SigTranmit.h"
 #include   
#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;jmit::Lpf()//低通 {  Conv(Sa,GetH1());  Dft(Sa,16*M);  } 
void SigTranmit::Dft(float x[],int m)//DFT变换 {  for(int k=0;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;ii/(N-1)));};   break;  case 1:   {    for(i=0;iGetHb());  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;j0)    {     Sa[k]=1;     k++;    }    else    {     Sa[k]=-1;     k++;