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

1、（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#defi

2、ne T 0.005#define M 20#define freq 4000#define fh 200#define fc 1000class 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 Idf

3、t(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();privat

4、e:int myID;int myIntensity; int c;/用来存放创函数选择int S016;int S116;float Mod16*M; float Sa16*M;/用来存放信号 float Sal16*M; float AR16*M; float AI16*M; float A16*M;/用来存放频域幅值 float Hl16*M; float H216*M; float No16*M;/用来存放生成的噪声 float HbN;/用来存放带通频域数据 float PN; float WnN; ;#endif / !defined(AFX_SIGTRANMIT_H_1E1D11

5、BF_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_FILEstatic char THIS_FILE=_FILE_;#define new

6、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 d44; int e16; int m=0; int b4; b0=MyID/1000; b1=(MyID%1000

7、)/100; b2=(MyID%100)/10; b3=MyID-b0*1000-b1*100-b2*10; for(int i=0;i4;i+) switch(bi) case 0:di0=0; di1=0;di2=0; di3=0; break; case 1:di0=0; di1=0;di2=0; di3=1; break; case 2:di0=0; di1=0;di2=1; di3=0; break; case 3:di0=0; di1=0;di2=1; di3=1; break; case 4:di0=0; di1=1;di2=0; di3=0; break; case 5:di0

8、=0; di1=1; di2=0; di3=1; break; case 6:di0=0; di1=1;di2=1; di3=0; break; case 7:di0=0; di1=1;di2=1; di3=1; break; case 8:di0=1; di1=0;di2=0; di3=0; break; case 9:di0=1; di1=0;di2=0; di3=1;break; int k; for(k=0;k4;k+) for(int j=0;j4;j+) S0m=dkj; m+; for( i=0;i16;i+)ei=S0i;if(ei=0)ei=-1;for(int j=0;j1

9、6;j+)for(k=0;kM;k+)Sa1k+M*j=ej;for(int n=0;n16*M;n+)San=Sa1n; void SigTranmit :Dft(float x,int m) /DFTfor(int k=0;km;k+) ARk=0.0;AIk=0.0; for(int n=0;nm;n+) ARk=ARk+xn*cos(-2.0*PI*k*n/m); AIk=AIk+xn*sin(-2.0*PI*k*n/m); Ak=sqrt(pow(ARk,2)+pow(AIk,2);void SigTranmit :Lpf() /低通滤波 float hdN; int i; floa

10、t wc=(2*PI*400)/fs; float a=(N-1)/2; for(i=0;iN;i+)hdi=0.0;for(i=0;iN;i+)if(i=a) hdi=wc/(double)PI;elsehdi=sin(wc*(i-a)/(double)(PI*(i-a); switch(c) case 0:for( i=0;iN;i+) Wni=0.5*(1.0-cos(2.0*PI*i/(N-1); ;break; case 1:for(i=0;iN;i+) Wni=0.54-0.46*cos(2.0*PI*i/(N-1); ;break; case 2:for( i=0;iN;i+)

11、Wni=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;iN;i+) Wni=1; ;break; for(i=0;iN;i+) H1i=hdi*Wni; void SigTranmit :Modulate() /调制for(int i=0;i16*M;i+)Sai=Sai*cos(2*PI*fc*i/fs);void SigTranmit :Channel() /信道for(int i=0;i16*M;i+)Sai=Sai+GNoise();void SigTranmit :Bpf()

12、/带通滤波 switch(s) case 0:int i,ks1,kp1,ks2,kp2; /s=0;选择频率采样法。float f;float ZN ;kp1=(int)700*N/fs;ks1=(int)900*N/fs;kp2=(int)1300*N/fs;ks2=(int)1100*N/fs;for(i=0;iN;i+)Hbi=0.0;for(i=0;iN;i+)Zi=-(N-1)*PI*i)/N; /线性相位for(i=0;i=kp1;i+) /N=127为奇数Hbi=0;HbN-i=Hbi;for(i=kp1+1;i=ks1;i+)f=i*fs/(float)N;Hbi=(f-70

13、0)/(float)200;HbN-i=Hbi;for(i=ks1+1;i=ks2;i+) f=i*fs/(float)N;Hbi=1.0;HbN-i=Hbi;for(i=ks2+1;i=kp2;i+)f=i*fs/(float)N;Hbi=(1300-f)/(float)200;HbN-i=Hbi;for(i=0;iN;i+)HbRi=Hbi*cos(Zi);HbIi=Hbi*sin(Zi);break; case 1: float hdN, hd1N,hd2N; /s=1,用窗函数法 float wu=(2*PI*1400)/fs; float wl=(2*PI*600)/fs; floa

14、t a=(N-1)/2; for(i=0;iN;i+)hd1i=0.0; hd2i=0.0;for(i=0;iN;i+)if(i=a) hd1i=wl/(double)PI;elsehd1i=sin(wl*(i-a)/(double)(PI*(i-a); for(i=0;iN;i+)if(i=a) hd2i=wu/(double)PI;elsehd2i=sin(wu*(i-a)/(double)(PI*(i-a);for(i=0;iN;i+)hdi=hd2i-hd1i; switch(c) case 0:for( i=0;iN;i+) Wni=0.5*(1.0-cos(2.0*PI*i/(N-

15、1); ;break; case 1:for(i=0;iN;i+) Wni=0.54-0.46*cos(2.0*PI*i/(N-1); ;break; case 2:for( i=0;iN;i+) Wni=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;iN;i+) Wni=1; ;break; for( i=0;iN;i+) Pi=hdi*Wni; Dft( P,N) ; for( i=0;iN;i+) Hbi=Ai; break; void SigTranmit :Demodulate

16、() /相干解调for(int i=0;i16*M;i+)Sai=Sai*cos(2*PI*fc*i/fs); void SigTranmit :Adjust() /恢复判决for(int i=0;i0)S1i=1;elseS1i=-1;for(int j=0;j16;j+)for(int k=0;kM;k+)Sak+M*j=S1j; void SigTranmit : Idft(float h) /IDFT int k; switch(s) case 0: for( k=0;kN;k+) HbRk=0.0;HbIk=0.0; for(int n=0;nN;n+) HbRk=(float)Hb

17、Rk+hn*cos(2.0*k-126)*PI*n/N); /? Pk=HbRk/127; ;break; case 1:break; void SigTranmit:Conv(float array116*M,float array2N) /卷积int I; float b16*M+N;I=16*M+N;for(int i=0;iI;i+)bi=0;for(int j=0;j(16*M);j+)if(i-j)N)bi=bi+0;elsebi=bi+array1j*array2i-j;for(int k=0;k(16*M);k+)H2k=bk+(N-1)/2; float SigTranmit

18、 :GNoise() /高斯白噪声 float b,q,randi=0.0; q=RAND_MAX; for(int i=0;i12;i+)randi+=rand()/q;b=myIntensity*(randi-6);b/=10;return b; void SigTranmit :Noise() for(int j=0;j16*M;j+)Noj=0;for(int i=0;i16*M;i+)Noi+=GNoise() ; void SigTranmit :SetH1(float H) for(int i=0;iN;i+)H1i=Hi; void SigTranmit :SetHb(float H) for(int i=0;iN;i+)Hbi=Hi;void SigTranmit :SetSa(float H) for(int i=0;i16*M;i+)Sai=Hi; void SigTranmit :Setc(int Myc)c=Myc; void SigTranmit :Sets(int Mys)s=Mys; int SigTranmit :GetmyIntensity()return myIntensity; int SigTranmit :Getc()retur