数字信号处理课设报告.xWord格式.docx
《数字信号处理课设报告.xWord格式.docx》由会员分享,可在线阅读,更多相关《数字信号处理课设报告.xWord格式.docx(24页珍藏版)》请在冰点文库上搜索。
3.1设计题目和题目设计要求................................................................................
3.2设计思想和系统功能分析................................................................................
3.3测试数据、测试输出结果,及理论分析和比较............................................
三、设计心得与体会......................................................................................................
四、参考资料..................................................................................................................
五、程序源代码..............................................................................................................
二、程序分析与设计
1.1设计题目和题目设计要求
一、
1)生成信号发生器:
能产生频率(或基频)为10Hz的周期性正弦波、三角波和方波信号。
绘出它们的时域波形
2)为避免频谱混叠,试确定各信号的采样频率。
说明选择理由。
3)对周期信号进行离散傅立叶变换,为了克服频谱泄露现象,试确定截取数据的长度,即信号长度。
分析说明选择理由。
4)绘出各信号频域的幅频特性和相频特性
5)以正弦周期信号为例,观察讨论基本概念(频谱混叠、频谱泄漏、整周期截取等)。
1.2设计思想和系统功能分析
(1)首先调用简单函数生成三种信号绘出时域波形,
(2)然后对连续信号进行采样,将其离散化得到离散信号,由题意可知信号最高频率分量为10HZ,由采样定理得知采样频率应该大于2fs=20Hz,一般取3到5倍,所以采样频率可取100Hz,采样点数可取128。
采样直接将t变换为1/fs即可,由简单的matlab语句即可实现,在采样结束后对三个序列进行离散傅里叶变换。
通过改变各信号的采样频率来观察频谱混叠现象;
(3)加窗将原信号截短形成有限长序列,通过改变截取数据的长度来观察分析频谱泄露现象;
(4)绘出各信号频域的幅频特性和相频特性;
(5)通过整周期截取与非整周期截取信号后波形及频谱来分析整周期截取的原因等。
1.3测试数据、测试输出结果,及理论分析和比较
(1)基频为10Hz的周期性正弦波、三角波和方波信号的时域波形正弦波信号:
x=sin(2*pi*f*t-0.02*pi);
三角波信号:
y=sawtooth(2*pi*f*t,0.5);
方波信号:
z=square(2*pi*f*t,50);
□三三三三三三三
1
-1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
2
-2
□三三三三三三三三
(2) 观察混叠现象
fs=100Hz时,混叠图像为
1 60
0 0.5 1
40
20
0 20 40 60
60
40
100
50
fs=50Hz时,混叠图像为
1 30
10
0 10 20 30
30
0 20
0
采样频率Fs只有在满足采样定理(Fs大于两倍的最大频率)时,采样信号通过一个理想低通滤波器恢复出原连续信号,所以信号的采样频率
Fs=100Hz。
如图中当Fs越来越小时,在f=Fs/2附近产生了频谱混叠现象,这时用DFT分析的结果必然在f=Fs/2产生较大的误差。
对于Fs确定的情况下,我们一般在采样前进行预滤波,滤除高于折叠频率Fs/2的频率成分,以免发生混叠现象。
(3)观察截断效应
Timedomain
Frequencydomain
0.5
-0.5
-10
-20
40 60 80
Samples
100 120
0.2
0.4
0.6
0.8
NormalizedFrequency(´
prad/sample)
Amplitude
Magnitude(dB)
下图为当N=40时,正弦波的截断信号与频谱图:
1 40
下图为当N=100时,正弦波的截断信号与频谱图:
-30
下图为当N=40时,三角波的截断信号与频谱图:
下图为当N=100时,三角波的截断信号与频谱图:
下图为当N=40时,方波的截断信号与频谱图:
下图为当N=100时,方波的截断信号与频谱图:
由以上图对比可以看出,增大N使信号幅度谱的主瓣变窄,减小泄露,提高了频率分辨率,但旁瓣的相对幅度并没有减小。
所以确定信号截取数据的长
度为N=100,能够更好地克服频谱泄露现象。
(4)绘出各信号频域的幅频特性和相频特性
□三三三三三
0 10 20 30 40 50 60 70 80 90 100
4
-4
□三三三三三三
80
(5)任以一数字正弦信号为例,观察讨论整周期截取
xn1
0 20 40 60 80 100 120 140
n1
□三三三三128三dft三三
□三
xn2
0 1 2 3 4 5 6 7
0 10 20 30 40 50 60 70
n2
□三三三三三128三dft三三
整周期截取实现不必加窗采到整周期的信号,上图N=128由于截取了整个周期得到单一谱线的频谱,下图N=60由于时域中对信号的非整周期截断产生频谱泄露。
2.1设计题目和题目设计要求
3
二、已知三个信号aip(n),经调制产生信号s(n)=å
aip(n)cos(ipn/4),其中ai为常数,p(n)为具有窄带特性的Hanning信号。
将此已调信号通过信道传输,描述该信道的差分方程为
y(n)-1.1172y(n-1)+0.9841y(n-2)-0.4022y(n-3)+0.2247y(n-4)
=0.2247x(n)-0.4022x(n-1)+0.9841x(n-2)-1.1172x(n-3)+x(n-4)
得到接收信号y(n)=s(n)*h(n)
1)分析Hanning信号p(n)的时域与频域特性
2)分析已调信号s(n)的时域与频域特性
3)分析系统的单位脉冲响应h(n)
4)分析接收信号y(n)的频谱
5)设计带通滤波器从接收信号y(n)中还原出三个已调信号。
2.2设计思想和系统功能分析
(1)首先调用简单的汉宁窗函数wn=hanning(N)或利用公式pn=0.5*[1-cos((2*pi*n)/(N-1))].*Rn绘出Hanning信号p(n)的时域波形与频谱图来分析
p(n)的时域与频域特性;
(2)与上一小问思路类似,由题中给出的信号s(n)的公式,给定ai值后绘出其时域波形与频谱图然后分析;
(3)构造单位脉冲函数xn,利用滤波器函数filter(b,a,xn),由题中差分方程确定参数a、b后即可得到单位脉冲响应h(n);
(4)与上一小问思路相同,可以将信道当作类似于滤波器的作用,经过信道传输后的信号y(n)也可以利用滤波器函数filter(b,a,yn)得到;
(5)最后利用带通滤波器从接收信号y(n)中分别还原出三个已调信号即可,但是试验后发现要滤出靠左和靠右的两个信号时,分别使用低通滤波器和高通滤波器滤出的结果比使用带通滤波器时得出的信号可靠一些,所以做出了两组结果图以供比较使用。
2.3测试数据、测试输出结果,及理论分析和比较
(1)Hanning信号p(n)的时域波形与频谱图
□三三
0 5 10 15 20 25 30 35 40
-50
-100
-150
4 30
-40
5 10 15
25 30 35 40
(2)已调信号s(n)的时域波形与频谱图
(3)系统的单位脉冲响应h(n)的时域波形与频谱图
-0.2
-0.4
0 5 10 15 20 25 30 35
0.02
0.01
-0.01
-0.02
-0.03
(4)接收信号y(n)的时域波形与频谱图
3 30
20 25
30 35 40
(5)设计带通滤波器从接收信号y(n)中还原出三个已调信号
5).1.1 设计带通滤波器从y(n)中还原出左边的信号
0.8 15
5
-5
-0.6
-15
5).1.2 设计低通滤波器从y(n)中还原出左边的信号
-0.8
-25
5 10
15 20
5).2 设计带通滤波器从y(n)中还原出中间的已调信号
5).3.1 设计带通滤波器从y(n)中还原出右边的已调信号
2 20
1.5
15
-1.5
2 25
5).3.2 设计高通滤波器从y(n)中还原出右边的已调信号
由以上两图可以看出,带通滤波器滤出的信号相比原信号主峰的高度下降很多,高通滤波器滤出的结果相对可靠一些。
3.1设计题目和题目设计要求
三、图像信号相关处理
1)读入一幅彩色图像
2)将