第二篇：电子科大通信学院DSP数字信号处理实验报告

电子科技大学

实   验   报   告

学生姓名：马露      学号：2010013020012      指导教师：潘晔

一、实验室名称：数字信号处理实验室

二、实验项目名称：多种离散时间信号的产生

三、实验原理：

1、基本离散时间信号

利用MATLAB强大的数值处理工具来实现信号的分析和处理，首先就是要学会应用MATLAB函数来构成信号。常见的基本信号可以简要归纳如下：

（1）．单位采样序列

    

在MATLAB中可以利用zeros()函数实现。

如果在时间轴上延迟了k个单位，得到即：

    

（2）．单位阶跃序列

     

在MATLAB中可以利用ones()函数实现。

 

（3）．正弦序列

采用MATLAB的实现方法，如：

（4）．实指数序列

其中，A、a为实数。采用MATLAB的实现方法，如：

（5）．复指数序列

          

采用MATLAB的实现方法，如：

    为了画出复数信号x[n]，必须要分别画出实部和虚部，或者幅值和相角。MATLAB函数real、imag、abs和angle可以逐次计算出一个复数向量的这些函数。

2、基本数字调制信号

（1）．二进制振幅键控（2ASK）

最简单的数字调制技术是振幅键控（ASK），即二进制信息信号直接调制模拟载波的振幅。二进制幅度键控信号的时域表达式：

   

其中，a_n为要调制的二进制信号，g(t)是单极性脉冲信号的时间波形，Ts 表示调制的信号间隔。典型波形如下：

             图 1 – 1   二进制振幅键控信号时间波形

（2）．二进制频移键控（2FSK）

    在二进制数字调制中，若正弦载波的频率随二进制基带信号在f₁和f₂两个频率点间变化，则产生二进制移频键控信号（2FSK信号）。二进制频域键控已调信号的时域表达式为:

这里，，。典型波形如下：

 

图 1 – 2   二进制频移键控信号时间波形

（3）．二进制相移键控（2PSK或BPSK）

    在二进制数字调制中，当正弦载波的相位随二进制数字基带信号离散变化时，则产生二进制移相键控(2PSK)信号。通常用已调信号载波的0°和 180°分别表示二进制数字基带信号的 1 和 0。二进制移相键控信号的时域表达式为：

 

典型波形如下：

 

图 1 – 3   二进制相移键控信号时间波形

数字调制信号可以根据调制原理采用基本的MATLAB命令实现，也可以用现成的命令：Y = DMOD(X, Fc, Fd, Fs, METHOD...)，其中，X为要调制的二进制信号；Fc为载波频率(Hz)；Fd为符号频率(Hz)；Fs为采样频率。要求：Fs > Fc，同时Fs/Fd 是一个正整数。参数METHOD是可以选择的调制方式：ask 、psk、qask、fsk、msk等。   

3、双音多频DTMF信号

DTMF（Double Tone MulitiFrequency,双音多频）作为实现电话号码快速可靠传输的一种技术，它具有很强的抗干扰能力和较高的传输速度，因此，可广泛用于电话通信系统中。但绝大部分是用作电话的音频拨号。另外，它也可以在数据通信系统中广泛地用来实现各种数据流和语音等信息的远程传输。

DTMF是用两个特定的单音频组合信号来代表数字信号以实现其功能的一种编码技术。两个单音频的频率不同，代表的数字或实现的功能也不同。这种电话机中通常有16个按键，其中有10个数字键0～9和6个功能键*、#、A、B、C、D。由于按照组合原理，一般应有8种不同的单音频信号。因此可采用的频率也有8种，故称之为多频，又因它采用从8种频率中任意抽出2种进行组合来进行编码，所以又称之为“8中取2”的编码技术。

根据CCITT的建议，国际上采用的多种频率为697Hz、770Hz、852Hz、941Hz、1209Hz、1336Hz、1477Hz和1633Hz等8种。用这8种频率可形成16种不同的组合，从而代表16种不同的数字或功能键，具体组合见表1-1。

表1-1 双音多频的组合功能

    因此，DTMF信号可以看作两个有限长度的正弦序列相加，正弦信号的频率由按键数字或字母符号对应的频率决定。如，数字“8”由行频852Hz和列频1336Hz决定。

四、实验目的：

1、掌握几种基本的离散时间信号（包括单位采样序列，单位阶跃序列，单频正弦序列，单频复指数序列，实指数序列等）。

2、能够熟练利用MATLAB产生这些基本的离散时间信号。

3、理解双音多频DTMF信号、ASK、FSK、BPSK等信号的产生原理。

4、学习并运用MATLAB产生各种通信中的调制信号及双音多频信号。

五、实验内容：

1、对几种基本离散时间信号（包括单位采样序列，单位阶跃序列，正弦序列，复指数序列，实指数序列等）在MATLAB中编程产生。

2、（拓展要求）利用MATLAB编程产生2ASK，2FSK，2PSK等数字调制信号。

3、（拓展要求）利用MATLAB编程产生理解双音多频DTFM信号。

4、（拓展要求）利用MATLAB编程产生高斯白噪声序列。

5、（拓展要求）利用MATLAB中的谱分析函数对正弦信号的频谱进行分析。

6、通过硬件（DSP）实验箱演示上述信号的时域（示波器）波形与频域波形（计算结果）。

六、实验器材（设备、元器件）：

安装MATLAB软件的PC机一台，DSP实验演示系统一套。

七、实验步骤：

1、在内，画出单位下列信号：

(a).单位采样序列和单位阶跃序列的时域波形图。

(b).、的波形。说明与、与之间的关系。

2、画出下列信号在内的波形。

     

      观察是否周期信号。如果是周期信号，信号的基波周期是什么？如果不是周期信号，说明原因。

3、在内，画出下列信号：

     

对于复数序列，要求分别画出实部和虚部；幅值和相角。若把中的底数0.8分别改为1.2、-0.8，讨论产生的时域波形有何变化。总结指数序列的底数对序列变化的影响。

4、（拓展要求）设计产生数字二进制序列：1 0 1 0 1 0 的2ASK、2FSK、2PSK调制信号。已知符号速率Fd＝10Hz（即时间间隔Ts为0.1），输出信号的采样频率为20Hz。

(a).2ASK信号的载波频率Fc＝5Hz，

(b).2FSK信号载波1频率F1＝5Hz，载波2频率F2＝1Hz。

(c).2PSK载波频率Fc＝1Hz。

分别画出以上信号调制前后的时域波形图。

5、（拓展要求）利用MATLAB产生DTMF双音多频信号。画出数字“0”的时域波形图。

6、（拓展要求）MATLAB函数randn(1,N)可以产生均值为0，方差为1的高斯随机序列，也就是白噪声序列。试利用randn函数产生均值为0.15，方差为0.1的高斯白噪声序列，要求序列时域范围为。画出时域波形图。同时将实验步骤2中产生的信号与相加，将得到的波形与的波形做比较。

7、（拓展要求）利用MATLAB中的谱分析函数画出的频谱。与理论上根据傅立叶变换的定义计算出的的频谱进行比较。

8、通过硬件（DSP）实验箱演示上述信号的时域（示波器）波形与频域波形（计算结果）。

八、实验数据及结果分析：

程序：

（1）产生序列的程序

1产生x1[n]、x2[n]、y1[n]、y2[n]序列的程序

k=1;

for n1=-20:20

    if n1==0;

        x1(k)=1;

    else

        x1(k)=0;

    end

    k=k+1;

end

subplot(2,2,1),stem(-20:20,x1);title('x1(n)');

k=1;

for n1=-20:20

    if n1>=0;

        x2(k)=1;

    else

        x2(k)=0;

    end

    k=k+1;

end

subplot(2,2,2),stem(-20:20,x2);title('x2(n)');

for k=1:41;

    if k<=36

    y1(k)=x1(k+5);

    else

        y1(k)=0;

    end

end

subplot(2,2,3),stem(-20:20,y1);title('y1(n)');

for k=1:41;

    if k>=28

        y2(k)=x2(k-8);

    else

        y2(k)=0;

    end;

end;

subplot(2,2,4),stem(-20:20,y2);title('y2(n)');

2产生x3[n]、x4[n]、x5[n]、x6[n]、x7[n]序列的程序

N=0:100

x3=sin(pi.*N/16);

x4=sin(N/2);

x5=cos(pi.*N/12)+cos(pi.*N*3/8);

figure(1),

subplot(2,2,1),stem(N,x3);title('x3(n)')

subplot(2,2,2),stem(N,x4);title('x4(n)')

subplot(2,2,3),stem(N,x5);title('x5(n)')

N=0:30;

x6=0.2*(0.8).^N;

x7=exp((-1/12+j*pi/6).*N);

subplot(2,2,4),stem(N,x6);title('x6(n)')

figure(2),

subplot(2,2,1),stem(N,real(x7));title('real{x7(n)}');

subplot(2,2,2),stem(N,imag(x7));title('imag{x7(n)}');

subplot(2,2,3),stem(N,abs(x7));title('Magnitude {x7(n)}');

subplot(2,2,4),stem(N,angle(x7));title('angle {x7(n)}');

（2）产生2ASK、2FSK、2PSK调制信号的程序（拓展要求）

（3）产生DTMF信号的程序（拓展要求）

（4）高斯白噪声序列的产生程序（扩展要求）

（4）正弦信号频谱分析的程序（扩展要求）

结果：

（1）的时域波形

1.x1[n]、x2[n]、x3[n]、x4[n]的时域波形。

X3[n]、x[4]、x[5]、x[6]的时域波形。

复数序列x7[n]的波形：

（2）信号的时移：与、与之间的关系。

         通过比较x1[n]与y1[n]的波形以及x2[n]与y2[n]的波形可知，y1[n]相当于将x1[n]向左搬移5个单位，y1[n]相当于将x1[n]向右搬移8个单位。

    （3）正弦序列周期的判断

         通过观察x3[n]、x4[n]、x5[n]的波形可知，x3[n]是周期信号、x4[n]不      是周期信号、x5[n]是周期信号。

（4）指数序列底数与序列变化的关系总结

对于x6[n]=0.2*(0.8)^N，通过改变底数为1.2、-0.8并画出它们的图像

clear all

N=0:30;

x61=0.2*(0.8).^N;

x62=0.2*(1.2).^N;

x63=0.2*(-0.8).^N;

subplot(3,1,1),stem(N,x61);title('x61(n)'); subplot(3,1,2),stem(N,x62);title('x62(n)');

subplot(3,1,3),stem(N,x63);title('x63(n)'); 

 

通过比较它们的图形，当指数序列的底数小于1大于0时，该序列在正数范围能单调递减；当指数序列的底数大于1时，该序列在正数范围能单调递增；当指数序列的底数小于0且绝对值小于1时，该序列在正数范围内出现正值与负值的交替，且绝对值单调递减。

（4）2ASK、2FSK、2PSK调制信号时域波形（拓展要求）

（5）数字“0” DTMF时域波形（拓展要求）

（6）高斯白噪声序列的时域波形；正弦序列加上高斯白噪声后的时域波形（扩展要求）

（7）正弦序列的频谱图（扩展要求）

九、实验结论：

对多种离散时间信号，进行程序的编写得到时域波形，通过对波形的观察可知，波形符合理论情况，结果正确。

十、总结及心得体会：

了解了多种离散时间信号的产生方法,加深了matlab对该课程的重要作用

十一、对本实验过程及方法、手段的改进建议：

 无

                                                     报告评分：

                                      指导教师签字：