数字调制仿真.doc

MATLAB的数字调制仿真实验报告

1：实验要求

实验要求通过输入随机信号的长度 ,得到二进制的随机原始信号 ,同时把得到的原始信号用三种不同的方法调制出来。当分别输入各个控件名称时 ,得到原始信号相应的信号输出。

2：实验过程

2.1 实验条件

1：实验的原始信号由MATLAB的randint(n)函数输出 ,需要确定的只是n,就是原始信号的宽度。

2：三种不同的调制函数

原始信号调制信号函数

振幅调制： 0： 0

1： cos(t+pi/3)

频移调制： 0： cos(t+pi/3)

1： cos(2*t+pi/6)

相位移调制： 0： cos(t)

1： cos(t+pi)

时间t为单个信号存在的时间周期 ,为了将图形表达更加清晰 ,这里选择将其选定为2*pi ,并划分为100个具体的时间点,t=0:2*pi/99:2*pi。

2.2 实验步骤

1：首先我要得到原始信号的长度 ,可以通过对s=rindint(n)函数产生的随机矩阵信号用length(s)求取其长度。

2：我们要得到单个的输入原始信号并对其进行调制 ,并同时将其用矩阵进行收集储存 ,最后输出调制后的信号。可以分别求取不同宽度上的信号 ,并将其赋值到对应输出原始信号的时间周期内 ,收集 ,最后输出。

3：调制得到的信号是在每个单个波长时间 ,不同的时间点t应用不同的调制函数的到的。在进行信号调制时 ,需要对这些调制得到的信号信息进行储存。可以在循环内采用矩阵叠加的方法来储存这些信号。

4: 需要的输入只是唯一的信号长度n,输出为得到的三种调制信号

5：编写实现输出全部调制信号的主函数Modulator和三个输出对应的调制信号的子函数ASK,FSK ,PSK.由子函数控制相应的信号输出.

6：编写程序,调试,写实验报告

3: 实验结果

通过输入不同Modulator(n) ,我们得到了调制的信号和相应的图形输出。

4: 附录实验程序及输出图

4.1: 主程序

function Modulator(n)

%定义函数

global ask

global fsk

global psk

global w

global signal

%定义全局变量

ask=[];

fsk=[];

psk=[];

%定义ASK,FSK,PSK调制信号

a=[];

f=[];

p=[];

%定义ASK,FSK,PSK决定信号

signal=[];

%定义输出原始信号和ASK选择信号

dfp=[];

%定义FSK,PSK选择信号

s=randint(n);

%得到原始信号

w=length(s);

%信号长度

t=0:2*pi/99:2*pi;

%划分单个信号的时间周期

for n=1:w

if s(n)==0;

signal1=zeros(1,100);

dfp1=ones(1,100);

f1=cos(t+pi/3);

p1=cos(t);

%产生并收集信号为0的时输出原始信号和调制信号的相应信息else

s(n)==1;

signal1=ones(1,100);

dfp1=ones(1,100);

f1=cos(2*t+pi/6);

p1=cos(t+pi);

%产生并收集信号为1时的输出原始信号和调制信号的相应信息end

signal=[signal signal1];

%得到输出原始信号信息和ASK调制的信号决定信息

a1=cos(t+pi/3);

a=[a a1];

%得到ASK调制的信号决定信息

f=[f f1];

%得到FSK调制的信号决定信息

p=[p p1];

%得到PSK调制的信号决定信息

dfp=[dfp dfp1];

%得到FSK,PSK调制的选择信息

end

%循环结束

ask=signal.*a;

fsk=dfp.*f;

psk=dfp.*p;

%得到ASK ,FSK ,PSK调制信号

ASK

FSK

PSK

%调用子函数

4.2：子函数

4.2.1：ASK调制程序

function ASK()

global ask

global w

global signal

figure(1)

subplot(2,1,1)

plot(signal,'LineWidth',1.5)

axis([0 100*w -1.5 1.5])

ylabel('调制前信号')

title('ASK信号调制图')

grid on

%画出输出原始信号图

subplot(2,1,2)

plot(ask,'LIneWidth',1.5)

axis([0 100*w -1.5 1.5])

xlabel('时间')

ylabel('2ASK调制后信号')

grid on

%画出输出ASK调制信号图

4.2.2: FSK调制程序

function FSK()

global fsk

global w

global signal

figure(2)

subplot(2,1,1)

plot(signal,'LIneWidth',1.5)

axis([0 100*w -1.5 1.5])

ylabel('调制前信号')

title('FSK信号调制图')

grid on

%画出输出原始信号图

subplot(2,1,2)

plot(fsk,'LIneWidth',1.5)

axis([0 100*w -1.5 1.5])

xlabel('时间')

ylabel('2FSK调制后信号')

grid on

%画出输出SFK调制信号图

4.2.3: PSK调制程序

function PSK()

global psk

global w

global signal

figure(3)

subplot(2,1,1)

plot(signal,'LIneWidth',1.5)

axis([0 100*w -1.5 1.5])

ylabel('调制前信号')

title('PSK信号调制图')

grid on

%画出输出信号原始图

subplot(2,1,2)

plot(psk,'LIneWidth',1.5)

axis([0 100*w -1.5 1.5])

xlabel('时间')

ylabel('2PSK调制后信号')

grid on

%画出输出PSK调制信号图