数字通信计算机仿真课程设计

（1）程序代码之头文件

// SigTranmit.h: interface for the SigTranmit class.

//

//////////////////////////////////////////////////////////////////////

#if !defined(AFX_SIGTRANMIT_H__1E1D11BF_1E65_4834_9496_4E5CAA470F68__INCLUDED_) #define AFX_SIGTRANMIT_H__1E1D11BF_1E65_4834_9496_4E5CAA470F68__INCLUDED_

#if _MSC_VER > 1000

#pragma once

#endif // _MSC_VER > 1000

#define PI 3.1415926

#define N 63

#define T 0.005

#define M 20

#define freq 4000

#define fh 200

#define fc 1000

class SigTranmit

{

public:

SigTranmit();

virtual ~SigTranmit();

void SetDigtal(int ID, int Intensity);

void Basesignal();

void Lpf();

void modulate();

void Channel();

void Bpf();

void Demodulate();

void Adjust();

void Dft(float x[],int m);

void Idft(float h[]);

void Hnc();

void Hmc();

void Blkmc();

void Conv(float array1[],float array2[]);

void Noise(int inten);

void SetHl(float H[]);

void SetHb(float H[]);

void Setc(int Myc);

int GetmyIntensity();

int Getc();

float* GetSa();

float* GetA();

float* GetHl();

float* GetH2();

float* GetHb();

float* GetNo();

private:

int myID;

int myIntensity;

int c;//用来存放创函数选择

int S0[16];

int S1[16];

float Mod[16*M];

float Sa[16*M];//用来存放信号

float Sal[16*M];

float AR[16*M];

float AI[16*M];

float A[16*M];//用来存放频域幅值

float Hl[16*M];

float H2[16*M];

float No[16*M];//用来存放生成的噪声

float Hb[N];//用来存放带通频域数据

float P[N];

float Wn[N];

};

#endif

// !defined(AFX_SIGTRANMIT_H__1E1D11BF_1E65_4834_9496_4E5CAA470F68__INCLUDED_)

（2）程序代码之源文件

// SigTranmit.cpp: implementation of the SigTranmit class.

//

//////////////////////////////////////////////////////////////////////

#include "stdafx.h"

#include "计算机仿真.h"

#include "SigTranmit.h"

#include "math.h"

#include "stdio.h"

#include "stdlib.h"

#include "iostream"

#ifdef _DEBUG

#undef THIS_FILE

static char THIS_FILE[]=__FILE__;

#define new DEBUG_NEW

#endif

//////////////////////////////////////////////////////////////////////

// Construction/Destruction

//////////////////////////////////////////////////////////////////////

SigTranmit::SigTranmit()

{

c=0;

s=1;

}

SigTranmit::~SigTranmit()

{

}

void SigTranmit ::SetDigital(int nMyID,int nIntensity) //初始化{

MyID= nMyID;

myIntensity=nIntensity;

}

void SigTranmit :: Basesignal() //基带信号的编码与采样{

int d[4][4];

int e[16];

int m=0;

int b[4];

b[0]=MyID/1000;

b[1]=(MyID%1000)/100;

b[2]=(MyID%100)/10;

b[3]=MyID-b[0]*1000-b[1]*100-b[2]*10;

for(int i=0;i<4;i++)

{

switch(b[i])

{

case 0:d[i][0]=0;

d[i][1]=0;

d[i][2]=0;

d[i][3]=0;

break;

case 1:d[i][0]=0;

d[i][1]=0;

d[i][2]=0;

d[i][3]=1;

break;

case 2:d[i][0]=0;

d[i][1]=0;

d[i][2]=1;

d[i][3]=0;

break;

case 3:d[i][0]=0;

d[i][1]=0;

d[i][2]=1;

d[i][3]=1;

break;

case 4:d[i][0]=0;

d[i][1]=1;

d[i][2]=0;

d[i][3]=0;

break;

case 5:d[i][0]=0;

d[i][1]=1;

d[i][2]=0;

d[i][3]=1;

break;

case 6:d[i][0]=0;

d[i][1]=1;

d[i][2]=1;

d[i][3]=0;

break;

case 7:d[i][0]=0;

d[i][1]=1;

d[i][2]=1;

d[i][3]=1;

break;

case 8:d[i][0]=1;

d[i][1]=0;

d[i][2]=0;

d[i][3]=0;

break;

case 9:d[i][0]=1;

d[i][1]=0;

d[i][2]=0;

d[i][3]=1;

break;

}

}

int k;

for(k=0;k<4;k++)

{

for(int j=0;j<4;j++)

{

S0[m]=d[k][j];

m++;

}

}

for( i=0;i<16;i++)

{

e[i]=S0[i];

if(e[i]==0)

{

e[i]=-1;

}

}

for(int j=0;j<16;j++)

{

for(k=0;k

{

Sa1[k+M*j]=e[j];

}

}

for(int n=0;n<16*M;n++)

{

Sa[n]=Sa1[n];

}

}

void SigTranmit :: Dft(float x[],int m) //DFT {

for(int k=0;k

{

AR[k]=0.0;AI[k]=0.0;

for(int n=0;n

{

AR[k]=AR[k]+x[n]*cos(-2.0*PI*k*n/m);

AI[k]=AI[k]+x[n]*sin(-2.0*PI*k*n/m);

}

A[k]=sqrt(pow(AR[k],2)+pow(AI[k],2));

}

}

void SigTranmit ::Lpf() //低通滤波{

float hd[N];

int i;

float wc=(2*PI*400)/fs;

float a=(N-1)/2;

for(i=0;i

{

hd[i]=0.0;

}

for(i=0;i

{

if(i==a)

hd[i]=wc/(double)PI;

else

hd[i]=sin(wc*(i-a))/(double)(PI*(i-a));

}

switch(c)

{

case 0:for( i=0;i

{

Wn[i]=0.5*(1.0-cos(2.0*PI*i/(N-1)));

};break;

case 1:for(i=0;i

{

Wn[i]=0.54-0.46*cos(2.0*PI*i/(N-1));

};break;

case 2:for( i=0;i

{

Wn[i]=0.42-0.5*cos(2.0*PI*i/(N-1))+0.08*cos(4.0*PI*i/(N-1)); };break;

case 3:for( i=0;i

{

Wn[i]=1;

};break;

}

for(i=0;i

H1[i]=hd[i]*Wn[i];

}

void SigTranmit ::Modulate() //调制

{

for(int i=0;i<16*M;i++)

Sa[i]=Sa[i]*cos(2*PI*fc*i/fs);

}

void SigTranmit ::Channel() //信道

{

for(int i=0;i<16*M;i++)

Sa[i]=Sa[i]+GNoise();

}

void SigTranmit ::Bpf() //带通滤波

{

switch(s)

{

case 0: int i,ks1,kp1,ks2,kp2; //s=0;选择频率采样法。

float f;

float Z[N] ;

kp1=(int)700*N/fs;

ks1=(int)900*N/fs;

kp2=(int)1300*N/fs;

ks2=(int)1100*N/fs;

for(i=0;i

{

Hb[i]=0.0;

}

for(i=0;i

{

Z[i]=-((N-1)*PI*i)/N; //线性相位}

for(i=0;i<=kp1;i++) //N=127为奇数{

Hb[i]=0;

Hb[N-i]=Hb[i];

}

for(i=kp1+1;i<=ks1;i++)

{

f=i*fs/(float)N;

Hb[i]=(f-700)/(float)200;

Hb[N-i]=Hb[i];

}

for(i=ks1+1;i<=ks2;i++)

{

f=i*fs/(float)N;

Hb[i]=1.0;

Hb[N-i]=Hb[i];

}

for(i=ks2+1;i<=kp2;i++)

{

f=i*fs/(float)N;

Hb[i]=(1300-f)/(float)200;

Hb[N-i]=Hb[i];

}

for(i=0;i

{

HbR[i]=Hb[i]*cos(Z[i]);

HbI[i]=Hb[i]*sin(Z[i]);

};

break;

case 1:

float hd[N], hd1[N],hd2[N]; //s=1,用窗函数法 float wu=(2*PI*1400)/fs;

float wl=(2*PI*600)/fs;

float a=(N-1)/2;

for(i=0;i

{

hd1[i]=0.0;

hd2[i]=0.0;

}

for(i=0;i

{

if(i==a)

hd1[i]=wl/(double)PI;

else

hd1[i]=sin(wl*(i-a))/(double)(PI*(i-a));

}

for(i=0;i

{

if(i==a)

hd2[i]=wu/(double)PI;

else

hd2[i]=sin(wu*(i-a))/(double)(PI*(i-a));

}

for(i=0;i

{

hd[i]=hd2[i]-hd1[i];

}

switch(c)

{

case 0:for( i=0;i

{

Wn[i]=0.5*(1.0-cos(2.0*PI*i/(N-1)));

};break;

case 1:for(i=0;i

{

Wn[i]=0.54-0.46*cos(2.0*PI*i/(N-1));

};break;

case 2:for( i=0;i

{

Wn[i]=0.42-0.5*cos(2.0*PI*i/(N-1))+0.08*cos(4.0*PI*i/(N-1)); };break;

case 3:for( i=0;i

{

Wn[i]=1;

};break;

}

for( i=0;i

P[i]=hd[i]*Wn[i];

Dft( P,N) ;

for( i=0;i

Hb[i]=A[i];

break;

}

}

void SigTranmit ::Demodulate() //相干解调

{

for(int i=0;i<16*M;i++)

Sa[i]=Sa[i]*cos(2*PI*fc*i/fs);

}

void SigTranmit ::Adjust() //恢复判决{

for(int i=0;i<16;i++)

{

S1[i]=(float)Sa[20+40*i];

if(S1[i]>0)

S1[i]=1;

else

S1[i]=-1;

}

for(int j=0;j<16;j++)

{

for(int k=0;k

{

Sa[k+M*j]=S1[j];

}

}

}

void SigTranmit :: Idft(float h[]) //IDFT

{

int k;

switch(s)

{ case 0:

for( k=0;k

{

HbR[k]=0.0;HbI[k]=0.0;

for(int n=0;n

{

HbR[k]=(float)HbR[k]+h[n]*cos((2.0*k-126)*PI*n/N); //?

P[k]=HbR[k]/127;

}

};break;

case 1:break;

}

}

void SigTranmit::Conv(float array1[16*M],float array2[N]) //卷积

{

int I;

float b[16*M+N];

I=16*M+N;

for(int i=0;i

{

b[i]=0;

for(int j=0;j<(16*M);j++)

{

if((i-j)<0||(i-j)>N)

b[i]=b[i]+0;

else

b[i]=b[i]+array1[j]*array2[i-j];

}

}

for(int k=0;k<(16*M);k++)

{

H2[k]=b[k+(N-1)/2];

}

}

float SigTranmit ::GNoise() //高斯白噪声

{

float b,q,randi=0.0;

q=RAND_MAX;

for(int i=0;i<12;i++)

randi+=rand()/q;

b=myIntensity*(randi-6);

b/=10;

return b;

}

void SigTranmit ::Noise()

{

for(int j=0;j<16*M;j++)

{No[j]=0;}

for(int i=0;i<16*M;i++)

{

No[i]+=GNoise() ;

}

}

void SigTranmit ::SetH1(float H[])

{

for(int i=0;i

H1[i]=H[i];

}

void SigTranmit ::SetHb(float H[])

for(int i=0;i

Hb[i]=H[i];

}

void SigTranmit ::SetSa(float H[])

{

for(int i=0;i<16*M;i++)

Sa[i]=H[i];

}

void SigTranmit ::Setc(int Myc) {

c=Myc;

}

void SigTranmit ::Sets(int Mys)

{

s=Mys;

}

int SigTranmit ::GetmyIntensity() {

return myIntensity;

}

int SigTranmit :: Getc()

{

return c;

}

float *SigTranmit ::GetSa()

{

return Sa;

}

float *SigTranmit ::GetA()

{

return A;

}

float *SigTranmit ::GetH1()

{

return H1;

}

float *SigTranmit ::GetH2()

{

return H2;

}

float *SigTranmit ::GetHb()

{

return Hb;

float* SigTranmit ::GetNo() {

return No;

}

float* SigTranmit ::GetP()

{

return P;

}

float * SigTranmit ::GetS1() {

return S1;

}