FMGen.cpp

#include"stdafx.h"#include"mod_demod.h"#include"FMGen.h"#include"head.h"#include"SourceGenClass.h"

#ifdef_DEBUG#undefTHIS_FILEstatic

char

THIS_FILE[]=__FILE__;#definenewDEBUG_NEW#endif

////////////////////////////////////////////////////////////////////////Construction/Destruction//////////////////////////////////////////////////////////////////////

FMGen::FMGen(CStringpnameModel,

CStringptypeModel,

long

vmodelSampleFreq,

long

vinputSampleFreq,

int

vtimeStep,

CStringpnameGetDataFile,

CStringpnamePutDataFile,

double

vampCarrier,

long

vfreqCarrier,

double

vphaseCarrier,

double

vModPara):ModBaseClass(pnameModel,ptypeModel,

vmodelSampleFreq,

vinputSampleFreq,

vtimeStep,

pnameGetDataFile,

pnamePutDataFile,

vampCarrier,

vfreqCarrier,

vphaseCarrier,

vModPara){

}

FMGen::~FMGen(){

}void

FMGen::run(){

double

dt;

float

*fn;

long

num_of_scan;

long

i;

double

phase_n,phase_n_1,dphase;

double

const_var;

dt=1.0/inputSampleFreq;

num_of_scan=(long)timeStep/1000.0*inputSampleFreq;

//begintomodulate

fn=getData(nameGetDataFile,1,num_of_scan);

if(!fn)

{

cout<<"Erroroccuredwhenread"<<nameGetDataFile<<".\n";

getch();

exit(1);

}

phase_n_1=ModPara*(*fn)*dt+phaseCarrier;

*fn=ampCarrier*cos(phase_n_1);

const_var=2*pi*freqCarrier*dt;

for

(i=1;i<num_of_scan;i++)

{

dphase=const_var+ModPara*(*(fn+i))*dt;

while

(dphase>=10000*2*pi)dphase-=2*pi*10000;

while

(dphase<=-2*pi*10000)dphase+=2*pi*10000;

phase_n=phase_n_1+dphase;

while

(phase_n>=2*pi*10000)phase_n-=2*pi*10000;

while

(phase_n<=-2*pi*10000)phase_n+=2*pi*10000;

*(fn+i)=ampCarrier*cos(phase_n);

phase_n_1=phase_n;

}

if(putData(fn,namePutDataFile,num_of_scan)==FALSE)

{

cout<<"whensavemodulationdata,erroroccured.\n";

free(fn);

getch();

exit(1);

}

free(fn);

cout<<"FMmodulationcomplete.\n";

}二．

//FMGen.h:interfacefortheFMGenclass.////////////////////////////////////////////////////////////////////////

#if!defined(AFX_FMGEN_H__119F55AA_0AA6_4CD7_948F_7010ED803ED4__INCLUDED_)#defineAFX_FMGEN_H__119F55AA_0AA6_4CD7_948F_7010ED803ED4__INCLUDED_

#if_MSC_VER>1000#pragmaonce#endif//_MSC_VER>1000

#include"ModBaseClass.h"

class

FMGen:

public

ModBaseClass

{

public:

FMGen(CString,CString,

long,long,

int,CString,

CString,double,

long,double,double);

virtual

~FMGen();

virtual

void

run();

};

#endif//!defined(AFX_FMGEN_H__119F55AA_0AA6_4CD7_948F_7010ED803ED4__INCLUDED_)附錄1：FM調(diào)制源程序清單

%程序名稱：tiaozhixinhao.m

%程序功能：產(chǎn)生模擬基帶信號(hào)并對(duì)調(diào)制信號(hào)進(jìn)行傅立葉變化得到調(diào)制信號(hào)的頻譜圖

%程序代碼

fc

=

150;

%載波頻率

fs

=

400;

%采樣頻率

t

=

(0

:0.001:0.15);

%時(shí)間區(qū)域

x

=

sin(2*pi*30*t);

%輸入信號(hào)

figure(1)

plot(t,x);

%調(diào)制信號(hào)時(shí)域圖xlabel('t');ylabel('x');

axis([0

0.15

-1

1])

a=fft(x,1024);

%對(duì)x進(jìn)行傅利葉變換

f=(0:length(a)-1)*fs/length(a)

-fs/2;

figure（2）

plot(f,abs(a));

%調(diào)制信號(hào)頻譜圖

xlabel('Frequence(Hz)');

%頻率分量

ylabel('powerSpectrum(x)');

%頻率譜能量

附錄2：FM調(diào)制源程序清單

%程序名稱：quyang400.m

%程序功能：調(diào)用函數(shù)modulate實(shí)現(xiàn)FM的調(diào)制

%程序代碼fc

=

150;

%載波頻率

fs

=

400;

%采樣頻率

t

=

(0

:0.001:0.15);

%時(shí)間區(qū)域x

=

sin(2*pi*30*t);

%輸入信號(hào)

y

=

modulate(x,fc,fs,'FM');

%調(diào)制信號(hào)

figure(3)

plot(t,y)

%對(duì)y進(jìn)行傅利葉變換

f=(0:length(b)-1)*fs/length(b)

-fs/2;

figure（4）

plot(f,abs(b));

%FM信號(hào)頻譜圖

xlabel('Frequence(Hz)');

%頻率分量

ylabel('Power

Spectrum(y)');

%頻率譜能量

附錄3：FM調(diào)制源程序清單

%程序名稱：quyang1600.m

%程序功能：調(diào)用函數(shù)modulate實(shí)現(xiàn)FM的調(diào)制

%程序代碼

fc

=

150;

%載波頻率

fs

=

1600;

%采樣頻率

t

=

(0

:0.001:0.15);

%時(shí)間區(qū)域

x

=

sin(2*pi*30*t);

%輸入信號(hào)

y

=

modulate(x,fc,fs,'FM');

%調(diào)制信號(hào)

plot(t,y);

%FM信號(hào)時(shí)域圖

xlabel('t(s)');ylabel('y');

axis([0

0.15

-1

1]);

b=fft(y,1024);

%對(duì)y進(jìn)行傅利葉變換

f=(0:length(b)-1)*fs/length(b)

-fs/2;

figure

plot(f,abs(b));

%FM信號(hào)頻譜圖

xlabel('Frequence(Hz)');

%頻率分量

ylabel('Power

Spectrum(y)');

%頻率譜能量

附錄4：FM調(diào)制源程序清單%程序名稱：xinzaobi1.m

%程序功能：對(duì)FM調(diào)制后的信號(hào)加入高斯噪聲，信噪比為1

%程序代碼

fc

=

150;

%載波頻率

fs

=

400;

%采樣頻率

t

=

(0

:0.001:0.15);

%時(shí)間區(qū)域

x

=

sin(2*pi*30*t);

%輸入信號(hào)

y

=

modulate(x,fc,fs,'FM');

%調(diào)制信號(hào)

y1

=

y

+

awgn(y,1,0);

%

FM信號(hào)加入噪聲

plot(t,y1);

%FM信號(hào)時(shí)域圖

xlabel('t(s)');ylabel('y');

axis([0

0.15

-1

1]);

b=fft(y1,1024);

%對(duì)y進(jìn)行傅利葉變換

f=(0:length(b)-1)*fs/length(b)

-fs/2;

figure

plot(f,abs(b));

%FM信號(hào)頻譜圖

xlabel('Frequence(Hz)');

%頻率分量

ylabel('Power

Spectrum(y)');

%頻率譜能量

附錄5：FM調(diào)制源程序清單

%程序名稱：xinzaobi2.m

%程序功能：對(duì)FM調(diào)制后的信號(hào)加入高斯噪聲，信噪比為0.1

%程序代碼

fc

=

150;

%載波頻率

fs

=

400;

%采樣頻率

t

=

(0

:0.001:0.15);

%時(shí)間區(qū)域

x

=

sin(2*pi*30*t);

%輸入信號(hào)

y

=

modulate(x,fc,fs,'FM');

%調(diào)制信號(hào)

y1

=

y

+

awgn(y,0.1,0);

%FM信號(hào)加入噪聲

plot(t,y1);

%FM信號(hào)時(shí)域圖

xlabel('t(s)');ylabel('y');

axis([0

0.15

-1

1]);

b=fft(y1,1024);

%對(duì)y進(jìn)行傅利葉變換

f=(0:length(b)-1)*fs/length(b)

-fs/2;

figure

plot(f,abs(b));

%FM信號(hào)信號(hào)頻譜

xlabel('Frequence(Hz)');

%頻率分量

ylabel('Power

Spectrum(y)');

%頻率譜能量

附錄6：FM調(diào)制源程序清單

%程序名稱：FM.m

%程序功能：對(duì)FM調(diào)制后的信號(hào)加入高斯噪聲，信噪比為5

%程序代碼

fc

=

150;

%載波頻率

fs

=

400;

%采樣頻率

t

=

(0

:0.001:0.15);

%時(shí)間區(qū)域

x

=

sin(2*pi*30*t);

%輸入信號(hào)

y

=

modulate(x,fc,fs,'FM');

%調(diào)制信號(hào)

y1

=

y