数字通信计算机仿真课设.docVIP

数字通信计算机课设 程序代码 // SigTranmit.cpp: implementation of the SigTranmit class. // ////////////////////////////////////////////////////////////////////// #include stdafx.h #include ASK.h #include SigTranmit.h #include math.h #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;j4;j++) { i=0; temp=a[j]; for(i;i4;i++) { S0[k]=temp 0x1; temp = temp1; k--; } } //采样 k=0; for(i=0;i16;i++) { for(j=0;jM;j++) { Sal[k]=S0[i]; k++; } } for(i=0;i16*M;i++) { if(Sal[i] == 0) Sa[i] = -1; else Sa[i] = Sal[i]; } Dft(Sa,16*M); } void SigTranmit::Lpf()//低通 { Conv(Sa,GetH1()); Dft(Sa,16*M); } void SigTranmit::Dft(float x[],int m)//DFT变换 { for(int k=0;km;k++) { AR[k]=0; AI[k]=0; for(int i=0;im;i++) { AR[k] =AR[k] + x[i]*cos(2*PI*k*i/m); AI[k] =AI[k] + x[i]*sin(2*PI*k*i/m); } A[k] = sqrt(AR[k]*AR[k] + AI[k]*AI[k]); } } void SigTranmit::modulate()//调制 { float a=0.0; for(int i=0;i16*M;i++) { a=2.0*PI*i*0.15; Mod[i] = Sa[i]*cos(a); } for(i=0;i16*M;i++) Sa[i] = Mod[i]; Dft(Sa,16*M); } void SigTranmit::Channel()//信道 { Noise(myIntensity); for(int i=0;i16*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; k16*M; k