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

数字通信计算机课设  程序代码 

// 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;j 

void SigTranmit::Lpf()//低通 {  Conv(Sa,GetH1());  Dft(Sa,16*M);  } 

void SigTranmit::Dft(float x[],int m)//DFT变换 {  for(int k=0;kvoid 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 

 }  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;j0)    {     Sa[k]=1;     k++;    }    else    {     Sa[k]=-1;     k++;