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

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

2、 N 63#define 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

3、);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* GetN

4、o();private: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_SIGTRANMI

5、T_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"#incl

6、ude "iostream"#ifdef _DEBUG#undef THIS_FILEstatic 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 SigTra

7、nmit : Basesignal() /基带信号的编码与采样 int d44; int e16; int m=0; int b4; b0=MyID/1000; b1=(MyID%1000)/100; b2=(MyID%100)/10; b3=MyID-b0*1000-b1*100-b2*10; for(int i=0;i<4;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; b

8、reak; 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=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;

9、for(k=0;k<4;k+) for(int j=0;j<4;j+) S0m=dkj; m+; for( i=0;i<16;i+)ei=S0i;if(ei=0)ei=-1;for(int j=0;j<16;j+)for(k=0;k<M;k+)Sa1k+M*j=ej;for(int n=0;n<16*M;n+)San=Sa1n; void SigTranmit :Dft(float x,int m) /DFTfor(int k=0;k<m;k+) ARk=0.0;AIk=0.0; for(int n=0;n<m;n+) ARk=ARk+xn*co

10、s(-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; float wc=(2*PI*400)/fs; float a=(N-1)/2; for(i=0;i<N;i+)hdi=0.0;for(i=0;i<N;i+)if(i=a) hdi=wc/(double)PI;elsehdi=sin(wc*(i-a)/(double)(PI*(i-a); switch(c) case 0:for( i

11、=0;i<N;i+) Wni=0.5*(1.0-cos(2.0*PI*i/(N-1); ;break; case 1:for(i=0;i<N;i+) Wni=0.54-0.46*cos(2.0*PI*i/(N-1); ;break; case 2:for( i=0;i<N;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;i<N;i+) Wni=1; ;break; for(i=0;i<N;i+) H1i=hdi*Wni; void SigTra

12、nmit :Modulate() /调制for(int i=0;i<16*M;i+)Sai=Sai*cos(2*PI*fc*i/fs);void SigTranmit :Channel() /信道for(int i=0;i<16*M;i+)Sai=Sai+GNoise();void SigTranmit :Bpf() /带通滤波 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;ks

13、2=(int)1100*N/fs;for(i=0;i<N;i+)Hbi=0.0;for(i=0;i<N;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-700)/(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/(f

14、loat)N;Hbi=(1300-f)/(float)200;HbN-i=Hbi;for(i=0;i<N;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; float a=(N-1)/2; for(i=0;i<N;i+)hd1i=0.0; hd2i=0.0;for(i=0;i<N;i+)if(i=a) hd1i=wl/(double)PI;elsehd1i=si

15、n(wl*(i-a)/(double)(PI*(i-a); for(i=0;i<N;i+)if(i=a) hd2i=wu/(double)PI;elsehd2i=sin(wu*(i-a)/(double)(PI*(i-a);for(i=0;i<N;i+)hdi=hd2i-hd1i; switch(c) case 0:for( i=0;i<N;i+) Wni=0.5*(1.0-cos(2.0*PI*i/(N-1); ;break; case 1:for(i=0;i<N;i+) Wni=0.54-0.46*cos(2.0*PI*i/(N-1); ;break; case 2

16、:for( i=0;i<N;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;i<N;i+) Wni=1; ;break; for( i=0;i<N;i+) Pi=hdi*Wni; Dft( P,N) ; for( i=0;i<N;i+) Hbi=Ai; break; void SigTranmit :Demodulate() /相干解调for(int i=0;i<16*M;i+)Sai=Sai*cos(2*PI*fc*i/fs); void S

17、igTranmit :Adjust() /恢复判决for(int i=0;i<16;i+)S1i=(float)Sa20+40*i; if(S1i>0)S1i=1;elseS1i=-1;for(int j=0;j<16;j+)for(int k=0;k<M;k+)Sak+M*j=S1j; void SigTranmit : Idft(float h) /IDFT int k; switch(s) case 0: for( k=0;k<N;k+) HbRk=0.0;HbIk=0.0; for(int n=0;n<N;n+) HbRk=(float)HbRk+h

18、n*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;i<I;i+)bi=0;for(int j=0;j<(16*M);j+)if(i-j)<0|(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 :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+)Noj=0;for(int i=0;i<16*M;i+)Noi+=GNoise() ; void SigTr