电子科技大学 DSP 实验一多种离散时间信号的产生Word文档下载推荐.docx
- 文档编号:21142040
- 上传时间:2023-01-27
- 格式:DOCX
- 页数:16
- 大小:395.11KB
电子科技大学 DSP 实验一多种离散时间信号的产生Word文档下载推荐.docx
《电子科技大学 DSP 实验一多种离散时间信号的产生Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《电子科技大学 DSP 实验一多种离散时间信号的产生Word文档下载推荐.docx(16页珍藏版)》请在冰豆网上搜索。
在二进制数字调制中,若正弦载波的频率随二进制基带信号在f1和f2两个频率点间变化,则产生二进制移频键控信号(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(DoubleToneMulitiFrequency,双音多频)作为实现电话号码快速可靠传输的一种技术,它具有很强的抗干扰能力和较高的传输速度,因此,可广泛用于电话通信系统中。
但绝大部分是用作电话的音频拨号。
另外,它也可以在数据通信系统中广泛地用来实现各种数据流和语音等信息的远程传输。
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双音多频的组合功能
高频群Hz
功能
低频群Hz
1209
1336
1477
1633
697
1
2
3
A
770
4
5
6
B
852
7
8
9
C
941
*
0
#
D
因此,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、(拓展要求)设计产生数字二进制序列:
101010的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)产生
序列的程序
(2)产生2ASK、2FSK、2PSK调制信号的程序(拓展要求)
(3)产生DTMF信号的程序(拓展要求)
(4)高斯白噪声序列的产生程序(扩展要求)
(4)正弦信号频谱分析的程序(扩展要求)
clearall;
clc;
n=101;
%单位采样序列
x1=zeros(1,n);
x1
(1)=1;
x1=[zeros(1,100),x1];
%单位阶跃序列
x2=ones(1,n);
x2=[zeros(1,100),x2];
%
n1=0:
n-1;
yn1=n1-5;
yn2=n1+8;
%100;
Fs=1000;
n2=0:
100;
%正弦序列
x3=sin(2*pi*n2/32);
x4=sin(n2/2);
x5=sin(pi*n2/12)+cos(3*pi*n2/8);
%指数序列
n3=0:
30;
x61=0.2*(0.8.^n3);
%实指数序列
x62=0.2*(1.2.^n3);
x63=0.2*((-0.8).^n3);
x7=exp((-1/12+1i*pi/6)*n3);
%复指数序列
%画出图形
figure
(1)
subplot(2,2,1),stem(n1,x1),title('
x1'
),axis([-20,20,0,1]);
subplot(2,2,2),stem(n1,x2),title('
x2'
subplot(2,2,3),stem(yn1,x1),title('
y1'
subplot(2,2,4),stem(yn2,x2),title('
y2'
figure
(2)
subplot(3,1,1),stem(n2,x3),title('
x3'
),axis([0,100,-1,1]);
subplot(3,1,2),stem(n2,x4),title('
x4'
subplot(3,1,3),stem(n2,x5),title('
x5'
),axis([0,100,min(x5),max(x5)]);
figure(3)
subplot(3,1,1),stem(n3,x61),title('
x6a=0.8'
),axis([0,30,min(x61),max(x61)]);
subplot(3,1,2),stem(n3,x62),title('
x6a=1.2'
),axis([0,30,min(x62),max(x62)]);
subplot(3,1,3),stem(n3,x63),title('
x6a=-0.8'
),axis([0,30,min(x63),max(x63)]);
figure(4)
subplot(4,1,1),stem(n3,abs(x7)),title('
x7幅值'
),axis([0,30,min(abs(x7)),max(abs(x7))]);
subplot(4,1,2),stem(n3,angle(x7)),title('
x7相角'
),axis([0,30,min(angle(x7)),max(angle(x7))]);
subplot(4,1,3),stem(n3,imag(x7)),title('
x7虚部'
),axis([0,30,min(imag(x7)),max(imag(x7))]);
subplot(4,1,4),stem(n3,real(x7)),title('
x7实部'
),axis([0,30,min(real(x7)),max(real(x7))]);
%调制
x_base=[1,0,1,0,1,0];
Fd=10000;
t=linspace(0,0.6,6*Fd);
if(x_base
(1)==1)
m=ones(1,Fd);
elseif(x_base
(1)==0)
m=zeros(1,Fd);
end
fori=2:
6
if(x_base(i)==1)
m=[m,ones(1,Fd)];
elseif(x_base(i)==0)
m=[m,zeros(1,Fd)];
end
%2ASK
Fc_A=5;
S_ask=m.*cos(2*pi*Fc_A*t);
figure(5)
subplot(4,1,1),stem(0:
0.1:
0.5,x_base),title('
序列An'
),axis([0,0.6,0,1]),xlabel('
s'
);
subplot(4,1,2),plot(t,m),title('
mt'
),axis([0,0.6,0,1.5]),title('
),xlabel('
subplot(4,1,3),plot(t,cos(2*pi*Fc_A*t)),title('
),axis([0,0.6,-1.2,1.2]),title('
载波信号'
subplot(4,1,4),plot(t,S_ask),title('
2ASK调制信号'
%2FSK
F1=5;
F2=1;
s1=m.*cos(2*pi*F1*t);
s2=(1-m).*cos(2*pi*F2*t);
S_fsk=s1+s2;
figure(6)
subplot(4,1,1),plot(t,m),axis([0,0.6,0,1.5]),title('
subplot(4,1,2),plot(t,s1),axis([0,0.6,-1.2,1.2]),title('
载波信号1F=5Hz)'
subplot(4,1,3),plot(t,s2),axis([0,0.6,-1.2,1.2]),title('
载波信号2F=1Hz'
subplot(4,1,4),plot(t,S_fsk),axis([0,0.6,-1.2,1.2]),title('
2FSK调制信号'
%2PSK
Fc_P=1;
S_psk=(2*m-1).*cos(2*pi*Fc_P*t+pi);
figure(7)
subplot(4,1,1),plot(t,2*m-1),axis([0,0.6,-1.5,1.5]),title('
subplot(4,1,2),plot(t,cos(2*pi*Fc_P*t+pi)),axis([0,0.6,-1.2,1.2]),title('
正相载波信号'
subplot(4,1,3),plot(t,-cos(2*pi*Fc_P*t+pi)),axis([0,0.6,-1.2,1.2]),title('
反相载波信号'
subplot(4,1,4),plot(t,S_psk),axis([0,0.6,-1.2,1.2]),title('
2PSK调制信号'
%DTFM
t_dt=linspace(0,0.02,10000);
x_dtfm=cos(2*pi*941*t_dt)+cos(2*pi*1366*t_dt);
plot(t_dt,x_dtfm);
%rand
N=201;
x8=sqrt(0.1)*randn(1,N)+0.15;
x_rnd=x2+x8;
figure(8)
subplot(3,1,1),stem(n1,x2),title('
X2'
subplot(3,1,2),stem(n1,x8),title('
高斯信号'
subplot(3,1,3),stem(n1,x_rnd),title('
加噪声后X2'
%FFT
N_smp=length(n2);
fre=linspace(-1,1,N_smp)*Fs/2;
y3=abs(fftshift(fft(x3)));
y4=abs(fftshift(fft(x4)));
y5=abs(fftshift(fft(x5)));
figure(9)
subplot(3,1,1),plot(fre,y3),xlabel('
Hz'
),title('
X3频谱'
频率Hz'
),axis([-100,100,1.2*min(y3),1.2*max(y3)]);
subplot(3,1,2),plot(fre,y4),xlabel('
X4频谱'
),axis([-200,200,1.2*min(y4),1.2*max(y4)]);
subplot(3,1,3),plot(fre,y5),xlabel('
X5频谱'
),axis([-300,300,1.2*min(y5),1.2*max(y5)]);
结果:
(1)
的时域波形
(2)信号的时移:
答:
相当于
向左平移5个单位,
相当于将
向右平移8个单位
(2)正弦序列
周期的判断
(3)指数序列底数与序列变化的关系总结
当指数序列的底数小于1大于0时,该序列在正数范围能单调递减;
当指数序列的底数大于1时,该序列在正数范围能单调递增;
当指数序列的底数小于0且绝对值小于1时,该序列在正数范围内出现正值与负值的交替,且绝对值单调递减。
(4)2ASK、2FSK、2PSK调制信号时域波形(拓展要求)
(5)数字“0”DTMF时域波形(拓展要求)
(6)高斯白噪声序列的时域波形;
正弦序列加上高斯白噪声后的时域波形(扩展要求)
(6)正弦序列的频谱图(扩展要求)
九、实验结论:
离散信号的频域分布是连续的,且仿真结论与理论值完全吻合
十、总结及心得体会:
通过本次试验,对离散信号的产生有了一定的了解,对各种信号的调制也更加的熟悉
与此同时,对matlab这一仿真工具的使用理解更深。
十一、对本实验过程及方法、手段的改进建议:
无
报告评分:
指导教师签字:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电子科技大学 DSP 实验一多种离散时间信号的产生 实验 多种 离散 时间 信号 产生