实验名称（基于DFT的信号识别系统）
一、实验目的
1)通过实验巩固离散傅立叶变换DFT的认识和理解。
2)熟练掌握应用DFT进行频谱分析的方法。
3)理解DFT离散频谱分析的应用价值。
二、实验内容
在语音识别、雷达信号处理、生物医学信号检测与识别等应用领域广泛使用基于离散傅立叶变换的谱分析技术。一个典型的信号识别系统如图1所示：
图1 信号识别系统框图
设系统的输入信号x(n)是具有单一频谱峰值的正弦信号，短时矩形窗将信号截短为有限长，经过DFT变换得到频谱，频率检测器检测频谱最大峰值的位置，即对应的频率，然后由分类器识别信号的类别。分类器的分类判决规则为：
第一类：最大峰值频率分布范围（Hz）为0≤f≤200。
第二类：最大峰值频率分布范围（Hz）为200≤f≤500。
第三类：最大峰值频率分布范围（Hz）为500≤f≤1000。
第四类：最大峰值频率分布范围（Hz）为f≥1000。
设采样频率fs=10000Hz，短时矩形窗宽度为N=1000，短时加窗信号经过DFT可以得到连续频谱在0≤ω＜2π范围内的1000个取样点。
1)编程实现该系统。
2)输入信号，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察系统
的实际识别结果，分析其正确性。
3)输入信号，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察系统
的实际识别结果，分析其正确性。
4)输入信号，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察系统
的实际识别结果，分析其正确性。
5)输入信号，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察系统
的实际识别结果，分析其正确性。
6)输入信号，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察系统
的实际识别结果，分析其正确性。
三、思考题
1)当矩形窗长度比1000小，例如32，以上实验内容（6）可能出现什么情况？
2)当输入信号时，系统能够得到正确的识别结果吗？为什么？
3)如果输入信号x(n)中含有叠加性宽带噪声e(n)会影响识别结果么？为什么？
4)如果系统中的DFT需要更新为FFT，并且短时窗不变，则FFT计算时应该做哪儿些
考虑？对识别结果会产生什么影响？
四、实验结果与分析
1. 编程实现该系统
本实验使用matlab 软件编程实现。
2. 输入信号
100
200
300
400
500600700
800
900
1000
signal
n
X (n )
frequency spectrum fenquency
M a g n i t u d e
图2 的波形与频谱（fs=10000Hz ，N=1000）
程序输出：最大峰频率：400Hz
最大幅度：1.20 类型：第二类
理论分析：函数中只有一个频率分量，f=fs*w /2π=400Hz ，幅度A=1.2，属于第二类。 程序结果与理论分析相符，程序正确。 3. 输入信号
100
200
300
400
500600700
800
900
1000
signal
n
X (n )
frequency spectrum fenquency
M a g n i t u d e
图3 的波形与频谱（fs=10000Hz ，N=1000）
程序输出：最大峰频率：2500Hz
最大幅度：3.00 类型：第四类
理论分析：函数中有俩个分量，f=fs*w /2π=2500Hz ，幅度A=3，直流分量C=1.5。因为A>C ，
所以属于第四类
程序结果与理论分析相符，程序正确。 4. 输入信号
100
200
300
400
500600700
800
900
1000
signal
n
X (n )
frequency spectrum fenquency
M a g n i t u d e
图4 的波形与频谱（fs=10000Hz ，N=1000）
程序输出：最大峰频率：700Hz
最大幅度：0.70 类型：第三类
理论分析：函数中只有一个频率分量，f=fs*w /2π=700Hz ，幅度A=0.7，属于第三类。 程序结果与理论分析相符，程序正确。 5. 输入信号
signal
n
X (n )
frequency spectrum fenquency
M a g n i t u d e
图5 的波形与频谱（fs=10000Hz ，N=1000）
程序输出：最大峰频率：100Hz
最大幅度：9.50 类型：第一类
理论分析：函数中有俩个分量，f1=fs*w 1/2π=2500Hz ，幅度A1=1.2，f1=fs*w 1/2π=100Hz ，幅度A=9.5。因为A2>A1，所以属于第一类。 程序结果与理论分析相符，程序正确。 6. 输入信号
signal
n
X (n )
frequency spectrum fenquency
M a g n i t u d e
图6 的波形与频谱（fs=10000Hz ，N=1000）
程序输出：最大峰频率：510Hz
最大幅度：1.00 类型：第三类
理论分析：函数中只有一个频率分量，，幅度A=1，属于第三类。 程序结果与理论分析相符，程序正确。
7. 当矩形窗长度比1000小，例如32，以上实验内容（6）可能出现什么情况？
5
10
15
2025
30
35
signal
n
X (n )
frequency spectrum fenquency
M a g n i t u d e
图7 的波形与频谱（fs=10000Hz ，N=32）
因为N 减小，导致?f=fs/N=312.5Hz ，因为中的频率分量为510Hz ，不能够落在采样点上，而
625Hz 距510Hz 最近，所以最大幅度点落在其上面，幅度小于1，由于短时窗过小，导致采样点过少，频谱泄露大，所以在其他频谱上出现了较大的幅值，已经不能够反应信号的准确信息。如果要求准确的反应信号的信息，即要求频率分辨率要小（?f 可取1,2,3,5..），使510Hz 初的k 可以取得整数。
8. 当输入信号时，系统能够得到正确的识别结果吗？为什么？
100
200
300
400
500600700
800
900
1000
signal
n
X (n )
frequency spectrum fenquency
M a g n i t u d e
图8 的波形与频谱（fs=10000Hz ，N=1000）
因为频率分辨率?f=fs/N=10Hz ，信号中的俩个频率分量为f1=995Hz ，幅度A1=1，f2=1020Hz ，幅度A2=0.9，相差25Hz ，从图8中可以区分俩个峰。从理论上讲，最大峰频率应当出现在995Hz 出，但是最小频率分辨率为10Hz ，导致采样点取不到995Hz 这个点，由于频率泄露的原因，导致1020Hz 这点的幅度变大，变成最大峰点。所以可以得出程序输出的结果是错误的，但是可以通过程序绘制的图判断出这个信号用过俩个频率分量。如果希望得到准确的识别结果，可以适当的增加短时窗的长度，降低最小频率分辨率?f 的数值(当最小频率分辨率?f 为5，1，0.5..时，可以准确的进行识别)。
200
400
600
800
100012001400
1600
1800
2000
-2-101
2signal
n
X (n )
frequency spectrum fenquency
M a g n i t u d e
图9 的波形与频谱（fs=10000Hz ，N=2000）
9. 如果输入信号x(n)中含有叠加性宽带噪声e(n)会影响识别结果么？为什么？
a) 噪声信号能量：99.7560信号能量：720.0000实际加入噪声后产生的信噪比值（db ）：
8.5839
100
200
300
400
500600700
800
900
1000
-202n
X (n )
signal
100
200
300
400
500600
700
800
900
1000
-202n
M a g n i t u d e
signal+noise
fenquency
M a g n i t u d e
图10 的波形、加噪声后的波形与频谱（fs=10000Hz ，N=1000）
b) 噪声信号能量：10.1144信号能量：720.0000实际加入噪声后产生的信噪比值：18.5239
100
200
300
400
500600700
800
900
1000
-20
2n
X (n )
signal
100
200
300
400
500600
700
800
900
1000
-20
2n
M a g n i t u d e
signal+noise
fenquency
M a g n i t u d e
图11 的波形、加噪声后的波形与频谱（fs=10000Hz ，N=1000）
c) 噪声信号能量：0.0990信号能量：720.0000实际加入噪声后产生的信噪比值：38.6148
100
200
300
400
500600700
800
900
1000
n
X (n )
signal
n
M a g n i t u d e
signal+noise
fenquency
M a g n i t u d e
图12 的波形、加噪声后的波形与频谱（fs=10000Hz ，N=1000）
噪声信号能量：792.7385信号能量：720.0000实际加入噪声后产生的信噪比值：
-0.4180
n
X
(
n
)
signal
01002003004005006007008009001000
n
M
a
g
n
i
t
u
d
e
signal+noise
fenquency
M
a
g
n
i
t
u
d
e
图13 的波形、加噪声后的波形与频谱（fs=10000Hz，N=1000）
因为信噪比的定义为：SNR=信号功率/噪声功率，信号功率=信号能量/N，噪声功率=噪声能量/N，由公式可得俩者能量，并计算信噪比。
由图10-13可知，加噪声之后主要是对信号的最大幅度产生影响，在信噪比较大的时候，对电路的判读的结果影响不大，在信噪比较小的时候，对最大幅度影响比较大，基本不影响最大峰频率，但如果噪声过大，则对最大峰频率和最大幅度均会产生影响。综上所述，若信噪大，对程序的判断结果影响不大；若信噪比小，对程序的判断结果影响大。
10. 如果系统中的DFT需要更新为FFT，并且短时窗不变，则FFT计算时应该做哪儿些考虑？对识别结果会产生什么影响？
当DFT更新为FFT时，要短时窗不变，在计算的时候应该考虑以下几个方面：
由于FFT要求信号的长度N是2的整数幂，可以采用补零法将N增加为2的整数幂，补零法保持时间分辨率不变，提高了频率分辨率。
①首先应该将信号序列分为偶序列X0(K)和奇序列X1(K)两个长度相等的序列。
②奇序列X1(K)和偶序列X0(K)的短时窗长度都将变为N/2=500点，而不再是N=1000。
③奇序列X1(K)和偶序列X0(K)的K值取值范围为K=0~N/2-1；，当N/2≤K＜N时，应该利用
X0(K)、X1(K)的周期性特征，即X0(K)= X0(K-N/2)，X1(K)=X1(K-N/2)。
对识别结果不会产生什么影响，但是会减少计算量，提高计算速率。
五、实验总结
通过这次实验，我学到了以下几点：
1.巩固了对傅里叶离散变化的认识和价值。
2.了解matlab这款数学软件功能的强大，基本熟悉了使用这款软件进行工程设
计和编程，将数据图形化，使结果更加直观和美观。
3.在实验中，学会合理使用工具书或其他途径查询资料。
4.熟练的使用DFT和FFT进行傅里叶变换，进行频谱分析。
5.懂得了进行信号系统设计时需要考虑很多方面，这样才能是设计的程序或系统
得出正确的结果和结论。
参考文献：
数字信号处理—理论与应用.俞一彪、孙兵.东南大学出版社.2005
《MATLAB函数速查手册》.邓薇.人民邮电出版社.2011
《MATLAB教程及实训》.曹弋.机械工业出版社.2008
附件1 实验源程序
l9.m
function y=l9(EN,N)
%y=l9(使能,短时窗长度)EN=1，则选择内置函数;EN=0,则选择外部输入函数，短时窗长度%返回判别的信号的类型
fs=10000;
n=0:N-1;
if EN==1
disp('1、x(n)=1.2*sin(0.08*pi*n)');
disp('2、x(n)=1.5+3*cos(0.5*pi*n)');
disp('3、x(n)=0.7*sin(0.14*pi*n)');
disp('4、x(n)=1.2*cos(0.5*pi*n)+9.5*sin(0.02*pi*n)');
disp('5、x(n)=cos(0.102*pi*n)');
disp('6、x(n)=cos(0.199*pi*n)+0.9*sin(0.204*pi*n)');
xz=input('选择哪一个函数进行试验\n');
switch(xz)
case 1
x=1.2*sin(0.08*pi*n);
case 2
x=1.5+3*cos(0.5*pi*n);
case 3
x=0.7*sin(0.14*pi*n);
case 4
x=1.2*cos(0.5*pi*n)+9.5*sin(0.02*pi*n);
case 5
x=cos(0.102*pi*n);
case 6
x=cos(0.199*pi*n)+0.9*sin(0.204*pi*n);
otherwise
disp('超过了函数选择范围');
end
else
A1=input('cos前的系数：');
W1=input('W1:');
A2=input('sin前的系数：');
W2=input('W1:');
C=input('直流量：');
x=A1*cos(W1*pi*n)+A2*sin(W2*pi*n)+C;
end
zs=input('是否要加入噪声！Y/N\n','s');
if zs=='y',
snr=input('加入噪声后，产生多少信噪比(db)：');
addx=addnoise(x,snr,N);
end
%x=A1*cos(w1*pi*n)+A2*sin(w2*pi*n)+C;
if zs=='y',
%X=fft(addx,W);%快速傅里叶变换
for i=0:N-1,%DFT系统
X(i+1)=sum(addx.*exp(-j*2*(pi/N)*i*n));
end
else
%X=fft(x,W);
for i=0:N-1,%DFT系统
X(i+1)=sum(x.*exp(-j*2*(pi/N)*i*n));
end
end
%for i=0:W-1,%DFT系统
% X1(i+1)=sum(x.*exp(-j*2*pi/W*i*n));
%end
%而第一个点就是直流分量，它的模值就是直流分量的N倍。
f=fs/N*n;%产生每个采样点的频率
Xm=abs(X)/(N/2);%假设原始信号的峰值为A，那么FFT的结果的每个点（除了第一个点直流分量之外）的模值就是A的N/2倍。
Xm(1)=Xm(1)/2;%第一个点就是直流分量，它的模值就是直流分量的N倍。
%max1=max(Xm);
[max1,wz]=max(Xm(1:N/2));%max1表示返回Xm中的最大值，如果有多个相同最大值时，返回到wz中的将是第一个的索引。
wz=wz-1;%最大值的真实位置要把索引减去1
%wz=find(Xm==max1)-1;
%for i=1:length(wz)
% if wz(i)<(W/2)
% f1=(fs/W)*wz(i);
% end
%end
f1=(fs/N)*wz;%最大值的点所对应的频率点
if f1>=0&&f1<200,
y='第一类';
elseif f1>=200&&f1<500,
y='第二类';
elseif f1>=500&&f1<1000,
y='第三类';
else
y='第四类';
end
if zs=='y',%绘制带噪声的图
subplot(3,1,1);
plot(n,x);
xlabel('n');ylabel('X(n)');
title('signal');
subplot(3,1,2);
plot(n,addx);
xlabel('n');ylabel('Magnitude');
title('signal+noise');
subplot(3,1,3);
%plot(f,Xm);
stem(f,Xm);
xlabel('fenquency');ylabel('Magnitude');
axis([0,fs,0,1.2*max1]);
str=sprintf('最大峰频率：%dHz\n最大幅度：%.2f\n类型：%s',f1,max1,y);
text(f1,max1,str);
else %绘制不添加噪声的图
subplot(2,1,1);
plot(n,x);
title('signal');
xlabel('n');ylabel('X(n)');
subplot(2,1,2);
stem(f,Xm);
title('frequency spectrum');
xlabel('fenquency');ylabel('Magnitude');
axis([0,fs,0,1.2*max1]);
str=sprintf('最大峰频率：%dHz\n最大幅度：%.2f\n类型：%s',f1,max1,y);
text(f1,max1,str); %最大位置处的注解
end
disp(str);%显示最大峰频率，最大幅度，判别类型
end
addnoise.m
function y=addnoise(x,snr,N)
%给x信号中添加高斯白噪声
%y=addnoise(输入的信号，加入噪声后信号拥有的理想信噪比,矩形窗宽度) y=awgn(x,snr);
%y=awgn(x,snr,'measured');
Ezs=sum((y-x).^2);%加入信号的噪声的能量
Exh=sum(x.^2);%信号自身的能量
snr1=10*log10((Exh./N)/(Ezs./N));%加入噪声后，实际的噪声比。
disp('噪声信号能量：');
disp(Ezs);%噪声信号能量
disp('信号能量：');
disp(Exh);%原信号能量
disp('实际加入噪声后产生的信噪比值：'); disp(snr1);
end
《数字信号处理》实验指导书 光电工程学院二○○九年十月
实验一离散时间信号分析 一、实验目的 1．掌握各种常用的序列，理解其数学表达式和波形表示。 2．掌握在计算机中生成及绘制数字信号波形的方法。 3．掌握序列的相加、相乘、移位、反转等基本运算及计算机实现与作用。 4．掌握线性卷积软件实现的方法。 5．掌握计算机的使用方法和常用系统软件及应用软件的使用。 6．通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。 二、实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列来表示，其中代表序列的第n个数字，n代表时间的序列，n的取值范围为的整数，n取其它值没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号进行等间隔采样，采样间隔为T，得到一个有序的数字序列就是离散时间信号，简称序列。 2．常用序列 常用序列有：单位脉冲序列（单位抽样）、单位阶跃序列、矩形序列、实指数序列、复指数序列、正弦型序列等。 3．序列的基本运算 序列的运算包括移位、反转、和、积、标乘、累加、差分运算等。 4．序列的卷积运算 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。 （1）反褶：先将和的变量换成，变成和，再将以纵轴为对称轴反褶成。 （2）移位：将移位，得。当为正数时，右移位；当为负数时，左
移位。 （3）相乘：将和的对应点值相乘。 （4）求和：将以上所有对应点的乘积累加起来，即得。 三、主要实验仪器及材料 微型计算机、Matlab软件6.5或更高版本。 四、实验内容 1．知识准备 认真复习以上基础理论，理解本实验所用到的实验原理。 2．离散时间信号（序列）的产生 利用MATLAB或C语言编程产生和绘制下列有限长序列： （1）单位脉冲序列 （2）单位阶跃序列 （3）矩形序列 （4）正弦型序列 （5）任意序列 3．序列的运算 利用MATLAB编程完成上述两序列的移位、反转、加法、乘法等运算，并绘制运算后序列的波形。 4．卷积运算 利用MATLAB编制一个计算两个序列线性卷积的通用程序，计算上述两序列，并绘制卷积后序列的波形。 5．上机调试并打印或记录实验结果。 6．完成实验报告。 五、实验报告要求 1. 简述实验原理及目的。 2. 给出上述序列的实验结果。 3. 列出计算卷积的公式，画出程序框图，并列出实验程序清单 （可略）（包括必要的程序说明）。 4. 记录调试运行情况及所遇问题的解决方法。 5. 给出实验结果，并对结果做出分析。 6. 简要回答思考题。 1 如何产生方波信号序列和锯齿波信号序列? 2 实验中所产生的正弦序列的频率是多少?是否是周期序列?
（1） 观察高斯序列的时域和幅频特性，固定信号)(n x a 中参数p=8,改变q 的 值，使q 分别等于2、4、8，观察他们的时域和幅频特性，了解当q 取不同值时，对信号序列的时域和幅频特性的影响；固定q=8，改变p,使p 分别等于8、13、14，观察参数p 变化对信号序列的时域和幅频特性的影响，注意p 等于多少时会发生明显的泄漏现象，混叠是否也随之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。 ()() ?????≤≤=-其他0150,2n e n x q p n a 解：程序见附录程序一： P=8,q 变化时： t/T x a (n ) k X a (k ) t/T x a (n ) p=8 q=4 k X a (k ) p=8 q=4 t/T x a (n ) p=8 q=8 k X a (k ) p=8 q=8 幅频特性 时域特性
t/T x a (n ) p=8 q=8 k X a (k ) p=8 q=8 t/T x a (n ) 5 10 15 k X a (k ) p=13 q=8 t/T x a (n ) p=14 q=8 5 10 15 k X a (k ) p=14 q=8 时域特性幅频特性 分析： 由高斯序列表达式知n=p 为期对称轴； 当p 取固定值时，时域图都关于n=8对称截取长度为周期的整数倍，没有发生明显的泄漏现象；但存在混叠，当q 由2增加至8过程中，时域图形变化越来越平缓，中间包络越来越大，可能函数周期开始增加，频率降低，渐渐小于fs/2,混叠减弱； 当q 值固定不变，p 变化时，时域对称中轴右移，截取的时域长度渐渐地不再是周期的整数倍，开始无法代表一个周期，泄漏现象也来越明显，因而图形越来越偏离真实值， p=14时的泄漏现象最为明显，混叠可能也随之出现；
实验一 离散时间信号分析 班级 信息131班 学号 201312030103 姓名 陈娇 日期 一、实验目的 掌握两个序列的相加、相乘、移位、反褶、卷积等基本运算。 二、实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列)}({n x 来表示，其中)(n x 代表序列的第n 个数字，n 代表时间的序列，n 的取值范围为+∞<<∞-n 的整数，n 取其它值)(n x 没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号)(t x a 进行等间隔采样，采样间隔为T ，得到)}({nT x a 一个有序的数字序列就是离散时间信号，简称序列。 2．常用序列 常用序列有：单位脉冲序列（单位抽样）） （n δ、单位阶跃序列)(n u 、矩形序列)(n R N 、实指数序列、复指数序列、正弦型序列等。 3．序列的基本运算 序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。 4．序列的卷积运算 ∑∞ -∞==-= m n h n x m n h m x n y )(*)()()()( 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。 （1）反褶：先将)(n x 和)(n h 的变量n 换成m ，变成)(m x 和)(m h ，再将)(m h 以纵轴为对称轴反褶成)(m h -。
（2）移位：将)(m h -移位n ，得)(m n h -。当n 为正数时，右移n 位；当n 为负数时，左移n 位。 （3）相乘：将)(m n h -和)(m x 的对应点值相乘。 （4）求和：将以上所有对应点的乘积累加起来，即得)(n y 。 三、主要实验仪器及材料 微型计算机、Matlab6.5 教学版、TC 编程环境。 四、实验内容 （1）用Matlab 或C 语言编制两个序列的相加、相乘、移位、反褶、卷积等的程序； （2）画出两个序列运算以后的图形； （3）对结果进行分析； （4）完成实验报告。 五、实验结果 六、实验总结
一、实验目的 1. 通过本次实验回忆并熟悉MATLAB这个软件。 2. 通过本次实验学会如何利用MATLAB进行序列的简单运算。 3. 通过本次实验深刻理解理论课上的数字信号处理的一个常见方法——对时刻n的样本附近的一些样本求平均，产生所需的输出信号。 3. 通过振幅调制信号的产生来理解载波信号与调制信号之间的关系。 二、实验内容 1. 编写程序在MATLAB中实现从被加性噪声污染的信号中移除噪声的算法，本次试验采用三点滑动平均算法，可直接输入程序P1.5。 2. 通过运行程序得出的结果回答习题Q1.31-Q1.33的问题，加深对算法思想的理解。 3. 编写程序在MATLAB中实现振幅调制信号产生的算法，可直接输入程序P1.6。 4. 通过运行程序得出的结果回答习题Q1.34-Q1.35的问题，加深对算法思想的理解。 三、主要算法与程序 1. 三点滑动平均算法的核心程序： %程序P1.5 %通过平均的信号平滑 clf; R=51; d=0.8*(rand(R,1)-0.5);%产生随噪声 m=0:R-1; s=2*m.*(0.9.^m);%产生为污染的信号 x=s+d';%产生被噪音污染的信号 subplot(2,1,1); plot(m,d','r-',m,s,'g--',m,x,'b-.');
xlabel('时间序号n');ylabel('振幅'); legend('d[n]','s[n]','x[n]'); x1=[0 0 x];x2=[0 x 0];x3=[x 0 0]; y=(x1+x2+x3)/3; subplot(2,1,2); plot(m,y(2:R+1),'r-',m,s,'g--'); legend('y[n]','s[n]'); xlabel('时间序号n');ylabel('振幅'); 2. 振幅调制信号的产生核心程序：（由于要几个结果，因此利用subplot函数画图） %程序P1.6 %振幅调制信号的产生 n=0:100; m=0.1;fH=0.1;fL=0.01; m1=0.3;fH1=0.3;fL1=0.03; xH=sin(2*pi*fH*n); xL=sin(2*pi*fL*n); y=(1+m*xL).*xH; xH1=sin(2*pi*fH1*n); xL1=sin(2*pi*fL1*n); y1=(1+m1*xL).*xH; y2=(1+m*xL).*xH1; y3=(1+m*xL1).*xH; subplot(2,2,1); stem(n,y); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.1;fH=0.1;fL=0.01;'); subplot(2,2,2); stem(n,y1); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.3;fH=0.1;fL=0.01;'); subplot(2,2,3); stem(n,y2); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.3;fH=0.3;fL=0.01;'); subplot(2,2,4); stem(n,y3); grid;
实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验容 1、帮助命令 使用 help 命令，查找 sqrt （开方）函数的使用方法； 2、MATLAB 命令窗口 （1）在MATLAB 命令窗口直接输入命令行计算3 1)5.0sin(21+=πy 的值； （2）求多项式 p(x) = x3 + 2x+ 4的根； 3、矩阵运算 （1）矩阵的乘法 已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B
（2）矩阵的行列式 已知A=[1 2 3;4 5 6;7 8 9]，求A （3）矩阵的转置及共轭转置 已知A=[1 2 3;4 5 6;7 8 9]，求A' 已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B' （4）特征值、特征向量、特征多项式 已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征值、特征向量、特征多项式；
（5）使用冒号选出指定元素 已知：A=[1 2 3;4 5 6;7 8 9]；求A 中第3 列前2 个元素；A 中所有列第2，3 行的元素； 4、Matlab 基本编程方法 （1）编写命令文件：计算1+2+…+n<2000 时的最大n 值；
（2）编写函数文件：分别用for 和while 循环结构编写程序，求 2 的0 到15 次幂的和。
5、MATLAB基本绘图命令 （1）绘制余弦曲线 y=cos(t)，t∈[0，2π]
（2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π] （3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求： （a）线形为点划线、颜色为红色、数据点标记为加号； （b）坐标轴控制：显示围、刻度线、比例、网络线 （c）标注控制：坐标轴名称、标题、相应文本； >> clear;
实验一 离散傅里叶变换（DFT ）对确定信号进行谱分析 一．实验目的 1．加深对DFT 算法原理和基本性质的理解。 2．熟悉DFT 算法和原理的编程方法。 3．学习用DFT 对信号进行谱分析的方法，了解可能出现的误差及其原因，以便在实际中正确利用。 二．实验原理 一个连续信号)(t x a 的频谱可以用其傅里叶变换表示，即 dt e t x j X t j a a Ω-∞ ∞ -? = Ω)()( 若对)(t x a 进行理想采样可得采样序列 )(|)()(nT x t x n x a nT t a === 对)(n x 进行DTFT ，可得其频谱为： ∑∞ -∞ =-= n n j j e n x e X ωω )()( 其中数字频率ω与模拟频率Ω的关系为： s f T Ω = Ω=ω )(n x 的DFT 为∑∞ -∞ =-= n nk N j e n x k X π 2)()( 若)(t x a 是限带信号，且在满足采样定理的条件下，)(ω j e X 是)(Ωj X a 的周期延拓， )(k X 是)(ωj e X 在单位圆上的等间隔采样值，即k N j e X k X πωω2| )()(= =。 为在计算机上分析计算方便，常用)(k X 来近似)(ω j e X ，这样对于长度为N 的有限 长序列（无限长序列也可用有限长序列来逼近），便可通过DFT 求其离散频谱。 三．实验内容 1．用DFT 对下列序列进行谱分析。 （1）)()04.0sin(3)(100n R n n x π=
1 （2）]0,0,0,0,0,0,0,0,1,1,1,1[)(=n x 2．为了说明高密度频谱和高分辨率频谱之间的区别，考察序列 )52.0cos()48.0cos()(n n n x ππ+= （1）当0≤n ≤10时，确定并画出x(n)的离散傅里叶变换。 （2）当0≤n ≤100时，确定并画出x(n)的离散傅里叶变换。 四．实验结果 1. （1） （2）
数字信号处理实验答案 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】
实验一熟悉Matlab环境 一、实验目的 1.熟悉MATLAB的主要操作命令。 2.学会简单的矩阵输入和数据读写。 3.掌握简单的绘图命令。 4.用MATLAB编程并学会创建函数。 5.观察离散系统的频率响应。 二、实验内容 认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。在熟悉了MATLAB基本命令的基础上，完成以下实验。 上机实验内容： （1）数组的加、减、乘、除和乘方运算。输入A=[1 2 3 4]，B=[3 4 5 6]，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B并用stem语句画出A、B、C、D、E、F、G。 clear all; a=[1 2 3 4]; b=[3 4 5 6]; c=a+b; d=a-b; e=a.*b; f=a./b; g=a.^b; n=1:4; subplot(4,2,1);stem(n,a); xlabel('n');xlim([0 5]);ylabel('A'); subplot(4,2,2);stem(n,b); xlabel('n');xlim([0 5]);ylabel('B'); subplot(4,2,3);stem(n,c); xlabel('n');xlim([0 5]);ylabel('C'); subplot(4,2,4);stem(n,d); xlabel('n');xlim([0 5]);ylabel('D'); subplot(4,2,5);stem(n,e); xlabel('n');xlim([0 5]);ylabel('E'); subplot(4,2,6);stem(n,f); xlabel('n');xlim([0 5]);ylabel('F'); subplot(4,2,7);stem(n,g); xlabel('n');xlim([0 5]);ylabel('G'); （2）用MATLAB实现下列序列： a) x(n)= 0≤n≤15 b) x(n)=e+3j)n 0≤n≤15 c) x(n)=3cosπn+π)+2sinπn+π) 0≤n≤15 d) 将c)中的x(n)扩展为以16为周期的函数x(n)=x(n+16),绘出四个周期。
实验四 有限长单位脉冲响应滤波器设计 朱方方 0806020433 通信四班 (1) 设计一个线性相位FIR 高通滤波器，通带边界频率为0.6π，阻带边界频率为0.4π，阻 带衰减不小于40dB 。要求给出h(n)的解析式，并用MATLAB 绘出时域波形和幅频特性。 解： (1) 求数字边界频率: 0.6 , .c r ωπωπ== (2) 求理想滤波器的边界频率: 0.5n ωπ= (3) 求理想单位脉冲响应: []d s i n ()s i n [()] () ()1n n n n n n h n n παωαα παωα π?-- -≠??-=? ? -=?? (4) 选择窗函数。阻带最小衰减为-40dB ，因此选择海明窗(其阻带最小衰减为-44dB)；滤 波器的过渡带宽为0.6π-0.4π=0.2π，因此 6.21 0.231 , 152 N N N ππα-=?=== (5) 求FIR 滤波器的单位脉冲响应h(n): []31d sin (15)sin[0.5(15)] 1cos ()15()()()15(15)1 15 n n n R n n h n w n h n n n ππππ?---????-? ?≠? ???==-???? ? ?=? 程序： clear; N=31; n=0:N-1; hd=(sin(pi*(n-15))-sin(0.5*pi*(n-15)))./(pi *(n-15)); hd(16)=0.5; win=hanning(N); h=win'.*hd; figure; stem(n,h); xlabel('n'); ylabel('h(n)'); grid; title('FIR 高通滤波单位脉冲响应h(n)'); [H,w]=freqz(h,1); H=20*log10(abs(H)); figure;3 plot(w/pi,H); axis([0 1 -100 10]); xlabel('\omega/\pi'); ylabel('幅度/dB'); grid; title('FIR 高通滤波器，hanning 窗，N=31');
实验三：离散LSI 系统的频域分析 一、实验内容 2、求以下各序列的z 变换： 12030() ()sin() ()sin()n an x n na x n n x n e n ωω-=== 程序清单如下： syms w0 n z a; x1=n*a^n;X1=ztrans(x1) x2=sin(w0*n);X2=ztrans(x2) x3= exp(-a*n)*sin(w0*n);X3=ztrans(x3) 程序运行结果如下： X1 =z/(a*(z/a - 1)^2) X2 =(z*sin(w0))/(z^2 - 2*cos(w0)*z + 1) X3 =(z*exp(a)*sin(w0))/(exp(2*a)*z^2 - 2*exp(a)*cos(w0)*z + 1) 3、求下列函数的逆z 变换 0 312342 1 1() () () ()() 1j z z z z X z X z X z X z z a z a z e z ω---= = = = ---- 程序清单如下： syms w0 n z a; X1=z/(z-a);x1=iztrans(X1) X2= z/(a-z)^2;x2=iztrans(X2) X3=z/ z-exp(j*w0);x3=iztrans(X3) X4=(1-z^-3)/(1-z^-1);x4=iztrans(X4) 程序运行结果如下： x1 =a^n x2 =n*a^n/a 课程名称 数字信号 实验成绩 指导教师 实 验 报 告 院系 信息工程学院 班级 学号 姓名 日期
x3 =charfcn[0](n)-iztrans(exp(i*w0),w0,n) x4 =charfcn[2](n)+charfcn[1](n)+charfcn[0](n) 4、求一下系统函数所描述的离散系统的零极点分布图，并判断系统的稳定性 （1） (0.3)()(1)(1) z z H z z j z j -= +-++ z1=[0,0.3]';p1=[-1+j,-1-j]';k=1; [b1,a1]=zp2tf(z1,p1,k); subplot(1,2,1);zplane(z1,p1); title('极点在单位圆外); subplot(1,2,2);impz(b1,a1,20); 由图可见:当极点位于单位圆内，系统的单位序列响应随着频率的增大而收敛；当极点位于单位圆上，系统的单位序列响应为等幅振荡；当极点位于单位圆外，系统的单位序列响应随着频率的增大而发散。由此可知系统为不稳定系统。 -1 -0.5 00.51 -2 -1.5-1-0.500.511.5 2Real Part I m a g i n a r y P a r t 极点在单位圆外 n (samples) A m p l i t u d e Impulse Response
实验一 基于Matlab 的数字信号处理基本操作 一、 实验目的：学会运用MA TLAB 表示的常用离散时间信号；学会运用MA TLAB 实现离 散时间信号的基本运算。 二、 实验仪器：电脑一台，MATLAB6.5或更高级版本软件一套。 三、 实验内容： （一） 离散时间信号在MATLAB 中的表示 离散时间信号是指在离散时刻才有定义的信号，简称离散信号，或者序列。离散序列通常用)(n x 来表示，自变量必须是整数。 离散时间信号的波形绘制在MATLAB 中一般用stem 函数。stem 函数的基本用法和plot 函数一样，它绘制的波形图的每个样本点上有一个小圆圈，默认是空心的。如果要实心，需使用参数“fill ”、“filled ”，或者参数“.”。由于MATLAB 中矩阵元素的个数有限，所以MA TLAB 只能表示一定时间范围内有限长度的序列；而对于无限序列，也只能在一定时间范围内表示出来。类似于连续时间信号，离散时间信号也有一些典型的离散时间信号。 1. 单位取样序列 单位取样序列)(n δ，也称为单位冲激序列，定义为 ) 0() 0(0 1)(≠=?? ?=n n n δ 要注意，单位冲激序列不是单位冲激函数的简单离散抽样，它在n =0处是取确定的值1。在MATLAB 中，冲激序列可以通过编写以下的impDT .m 文件来实现，即 function y=impDT(n) y=(n==0); %当参数为0时冲激为1,否则为0 调用该函数时n 必须为整数或整数向量。 【实例1-1】 利用MATLAB 的impDT 函数绘出单位冲激序列的波形图。 解：MATLAB 源程序为 >>n=-3:3; >>x=impDT(n); >>stem(n,x,'fill'),xlabel('n'),grid on >>title('单位冲激序列') >>axis([-3 3 -0.1 1.1]) 程序运行结果如图1-1所示。 图1-1 单位冲激序列
1213225 王聪 微弱信号检测技术概述 在自然现象和规律的科学研究和工程实践中, 经常会遇到需要检测毫微伏量级信号的问题, 比如测定地震的波形和波速、材料分析时测量荧光光强、卫星信号的接收、红外探测以及电信号测量等, 这些问题都归结为噪声中微弱信号的检测。在物理、化学、生物医学、遥感和材料学等领域有广泛应用。微弱信号检测技术是采用电子学、信息论、计算机和物理学的方法, 分析噪声产生的原因和规律, 研究被测信号的特点和相关性, 检测被噪声淹没的微弱有用信号。微弱信号检测的宗旨是研究如何从强噪声中提取有用信号, 任务是研究微弱信号检测的理论、探索新方法和新技术, 从而将其应用于各个学科领域当中。微弱信号检测的不同方法 ( 1) 生物芯片扫描微弱信号检测方法 微弱信号检测是生物芯片扫描仪的重要组成部分, 也是生物芯片技术前进过程中面临的主要困难之一, 特别是在高精度快速扫描中, 其检测灵敏度及响应速度对整个扫描仪的性能将产生重大影响。 随着生物芯片制造技术的蓬勃发展, 与之相应的信号检测方法也迅速发展起来。根据生物芯片相对激光器及探测器是否移动来对生物芯片进行扫读, 有扫描检测和固定检测之分。扫描检测法是将激光器及共聚焦显微镜固定, 生物芯片置于承片台上并随着承片台在X 方向正反线扫描和r 方向步进向前运动, 通过光电倍增管检测激发荧光并收集数据对芯片进行分析。激光共聚焦生物芯片扫描仪就是这种检测方法的典型应用, 这种检测方法灵敏度高, 缺点是扫描时间较长。 固定检测法是将激光器及探测器固定, 激光束从生物芯片侧向照射, 以此解决固定检测系统的荧光激发问题, 激发所有电泳荧光染料通道, 由CCD捕获荧光信号并成像, 从而完成对生物芯片的扫读。CCD 生物芯片扫描仪即由此原理制成。这种方法制成的扫描仪由于其可移动, 部件少, 可大大减少仪器生产中的失误, 使仪器坚固耐用; 但缺点是分辨率及灵敏度较低。根据生物芯片所使用的标记物不同, 相应的信号检测方法有放射性同位素标记法、生物素标记法、荧光染料标记法等。其中放射性同位素由于会损害研究者身体, 所以这种方法基本已被淘汰; 生物素标记样品分子则多用在尼龙膜作载体的生物芯片上, 因为在尼龙膜上荧光标记信号的信噪比较低, 用生物素标记可提高杂交信号的信噪比。目前使用最多的是荧光标记物, 相应的检测方法也最多、最成熟, 主要有激光共聚焦显微镜、CCD 相机、激光扫描荧光显微镜及光纤传感器等。 ( 2) 锁相放大器微弱信号检测 常规的微弱信号检测方法根据信号本身的特点不同, 一般有三条途径: 一是降低传感器与放大器的固有噪声, 尽量提高其信噪比; 二是研制适合微弱检测原理并能满足特殊需要的器件( 如锁相放大器) ;三是利用微弱信号检测技术, 通过各种手段提取信号, 锁相放大器由于具有中心频率稳定, 通频带窄,品质因数高等优点得到广泛应用。常用的模拟锁相放大器虽然速度快, 但是参数稳定性和灵活性差, 而且在与微处理器通信时需要转换电路; 传统数字锁相放大器一般使用高速APDC 对信号进行高速采样, 然后使用比较复杂的算法进行锁相运算, 这对微处理器的速度要求很高。现在提出的新型锁相检测电路是模拟和数字处理方法的有机结合, 这种电路将待测信号和参考信号相乘的结果通过高精度型APDC 采样,
数字信处理上机实验答 案全 Document number【SA80SAB-SAA9SYT-SAATC-SA6UT-SA18】
第十章 上机实验 数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。 实验一系统响应及系统稳定性。 实验二时域采样与频域采样。 实验三用FFT对信号作频谱分析。 实验四 IIR数字滤波器设计及软件实现。 实验五 FIR数字滤波器设计与软件实现 实验六应用实验——数字信号处理在双音多频拨号系统中的应用 任课教师根据教学进度，安排学生上机进行实验。建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在。学习完第六章进行；实验五在学习完第七章后进行。实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。 实验一: 系统响应及系统稳定性 1.实验目的 （1）掌握求系统响应的方法。 （2）掌握时域离散系统的时域特性。 （3）分析、观察及检验系统的稳定性。 2.实验原理与方法 在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。在计算机上适合用递推法求差分方程的解，最简单的方法是采用MATLAB语言的工具箱函数filter函数。也可以用MATLAB语言的工具箱函数conv函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。 系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。 系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。或者系统的单位脉冲响应满足绝对可和的条件。系统的稳定性由其差分方程的系数决定。 实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否都是有界输出，或者检查系统的单位脉冲响应满足绝对可和的条件。可行的方法是在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳定的[19]。系统的稳态输出是指当∞ n时，系统的输出。如果系统稳定，信号加入 → 系统后，系统输出的开始一段称为暂态效应，随n的加大，幅度趋于稳定，达到稳态输出。 注意在以下实验中均假设系统的初始状态为零。 3．实验内容及步骤
数字信号处理实验四 第一题结果： （1）没有增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 %H(3,13) = 0.75;H(5,11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线
（2）增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 H(3) = 0.75;H(13) = 0.75;H(5) = 0.25;H(11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线 第二题结果：
本科生实验报告 实验课程数字信号处理基础 学院名称地球物理学院 专业名称地球物理学 学生姓名 学生学号 指导教师王山山 实验地点5417 实验成绩 二〇一四年十一月二〇一四年十二月
填写说明 1、适用于本科生所有的实验报告（印制实验报告册除外）； 2、专业填写为专业全称，有专业方向的用小括号标明； 3、格式要求： ①用A4纸双面打印（封面双面打印）或在A4大小纸上用蓝黑色水笔书写。 ②打印排版：正文用宋体小四号，1.5倍行距，页边距采取默认形式（上下2.54cm， 左右2.54cm，页眉1.5cm，页脚1.75cm）。字符间距为默认值（缩放100%，间距：标准）；页码用小五号字底端居中。 ③具体要求： 题目（二号黑体居中）； 摘要（“摘要”二字用小二号黑体居中，隔行书写摘要的文字部分，小4号宋体）； 关键词（隔行顶格书写“关键词”三字，提炼3-5个关键词，用分号隔开，小4号黑体)； 正文部分采用三级标题； 第1章××(小二号黑体居中，段前0.5行) 1.1 ×××××小三号黑体×××××（段前、段后0.5行） 1.1.1小四号黑体（段前、段后0.5行） 参考文献（黑体小二号居中，段前0.5行），参考文献用五号宋体，参照《参考文献著录规则（GB/T 7714－2005）》。
实验一生成离散信号并计算其振幅谱 并将信号进行奇偶分解 一、实验原理 单位脉冲响应h(t)=exp(-a*t*t)*sin(2*3.14*f*t)进行离散抽样，分别得到t=0.002s，0.009s，0.011s采样的结果。用Excel软件绘图显示计算结果。并将信号进行奇偶分解，分别得到奇对称信号h(n)-h(-n)与偶对称信号h(n)+h(-n)。用Excel 软件绘图显示计算结果。 二、实验程序代码 （1）离散抽样 double a,t; a=2*f*f*log(m); int i; for(i=0;i 成都理工大学 《信号处理基础》实验 开设时间：2013—2014学年第2学期 题目1：信号的产生和显示 一、实验目的： 认识基本信号 通过使用MATLAB 设计简单程序, 掌握对MATLAB 的基本使用方法 二、实验原理： 找出下列表达式的信号与：正弦信号、最小相位信号、最大相位信号、零相位信号的对应关系。 1、sin60t 2、e-60t sin60t 3、(1- e-60t)sin60t 4、e60t sin60t 三、实验内容： 产生上述信号的信号并显示 （1）t=[-pi/30:0.001:pi/30]; f=sin(60*t); plot(t,f) 产生图形如下： （2）t=[0:0.001:pi/30]; f=exp(-60*t).*sin(60*t); plot(t,f) 产生图形如下: （3）t=[-5*pi/30:0.001:5*pi/30]; f=(1-exp(-60*t)).*sin(60*t); plot(t,f) 产生图形如下： (4) t=[-pi/30:0.001:pi/30]; f=exp(6*t).*sin(60*t); plot(t,f) 产生如下波形: 四、实验结果与讨论： 讨论上述信号的特点 从第一个波形图可以看出，它的波形与正弦函数sin（t）的相像，只是相位上有改变，是一个正弦信号。最大相位信号的能量集中在后面，最小相位能量集中在前面，所以第二个是一个最小相位，第四个是一个最大相位信号。第三个由于波形在t>0时没有，所以是一个零相位信号。 题目2：频谱分析与显示 一、实验目的 初步认识频谱分析 § Z 变换 ? Z 变换的定义及收敛域 【习题】 1. 假如)(n x 的z 变换代数表示式是下式，问)(z X 可能有多少不同的收敛域。 )83451)(411(411)(2122----+++- =z z z z z X 【分析】 ）要单独讨论，（环状、圆外、圆内：有三种收敛域：双边序列的收敛域为：特殊情况有：左边序列的收敛域为：因果序列的收敛域为：右边序列的收敛域为：特殊情况有：有限长序列的收敛域为 0 0 , , 0 0 , , 0 , 0 0 , 0 , 0 22 11 212 1∞==<<≤≤<≤<<≥≥∞≤<≥∞<<≤∞<≤≥∞≤<≤≤∞<<+ -++--z z R z R n n R z n n R z n n z R n n z R n z n z n n n z x x x x x x 解：对X (Z )的分子和分母进行因式分解得 )43 1 )(21 1)(211(2111111----+-+- =Z jZ jZ Z X (Z )的零点为：1/2，极点为：j/2，-j/2，-3/4 ∴ X (Z )的收敛域为： (1) 1/2 < | Z | < 3/4，为双边序列，见图一 (2) | Z | < 1/2，为左边序列，见图二 (3) | Z | > 3/4，为右边序列，见图三 图一 图二 图三 )431)(211)(411()211)(211()(11211-----++++- =Z Z Z Z Z Z X ? Z 反变换 【习题】 2. 有一右边序列 )(n x ，其 z 变换为)1)(211(1 )(11----=z z z X (a) 将上式作部分分式展开(用 1-z 表示)，由展开式求 )(n x 。 (b) 将上式表示成 z 的多项式之比，再作部分分式展开，由展开式求 )(n x ,并说明所得到的序列 与(a)所得的是一样的。 【注意】不管哪种表示法最后求出 x (n ) 应该是相同的。 解：(a) 因为11122 111)(---+--=z z z X 且x(n)是右边序列 所以 )()212()(n u n x n ?? ? ??-= (b) 122 1211 )1)(2 1(21231 )1)(2 1()(2 -+--+=---+=--=z z z z z z z z z X )()212( )1(2)1(21)()( n u n u n u n n x n n ??? ??-=-+-?? ? ??-=δ则 数字信号处理上机实验 答案 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】 第十章上机实验 数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。 实验一系统响应及系统稳定性。 实验二时域采样与频域采样。 实验三用FFT对信号作频谱分析。 实验四 IIR数字滤波器设计及软件实现。 实验五 FIR数字滤波器设计与软件实现 实验六应用实验——数字信号处理在双音多频拨号系统中的应用 任课教师根据教学进度，安排学生上机进行实验。建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在。学习完第六章进行；实验五在学习完第七章后进行。实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。 functiontstem(xn,yn) %时域序列绘图函数 %xn:信号数据序列，yn:绘图信号的纵坐标名称（字符串） n=0:length(xn)-1; stem(n,xn,'.');boxon xlabel('n');ylabel(yn); axis([0,n(end),min(xn),*max(xn)]) 实验一: 系统响应及系统稳定性 1.实验目的 （1）掌握求系统响应的方法。 （2）掌握时域离散系统的时域特性。 （3）分析、观察及检验系统的稳定性。 2.实验原理与方法 在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。在计算机上适合用递推法求差分方程的解，最简单的方法是采用MATLAB语言的工具箱函数filter函数。也可 第十章 上机实验 数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。 实验一 系统响应及系统稳定性。 实验二 时域采样与频域采样。 实验三 用FFT 对信号作频谱分析。 实验四 IIR 数字滤波器设计及软件实现。 实验五 FIR 数字滤波器设计与软件实现 实验六 应用实验——数字信号处理在双音多频拨号系统中的应用 任课教师根据教学进度，安排学生上机进行实验。建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR 数字滤波器设计及软件实现在。学习完第六章进行；实验五在学习完第七章后进行。实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。 10.1 实验一: 系统响应及系统稳定性 1.实验目的 （1）掌握 求系统响应的方法。 （2）掌握时域离散系统的时域特性。 （3）分析、观察及检验系统的稳定性。 2.实验原理与方法 在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。在计算机上适合用递推法求差分方程的解，最简单的方法是采用MA TLAB 语言的工具箱函数filter 函数。也可以用MATLAB 语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。 系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。 系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。或者系统的单位脉冲响应满足绝对可和的条件。系统的稳定性由其差分方程的系数决定。 实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否都是有界输出，或者检查系统的单位脉冲响应满足绝对可和的条件。可行的方法是在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳定的[19]。系统的稳态输出是指当∞→n 时，系统的输出。如果系统稳定，信号加入系统后，系统输出的开始一段称为暂态效应，随n 的加大，幅度趋于稳定，达到稳态输出。 注意在以下实验中均假设系统的初始状态为零。 3．实验内容及步骤 （1）编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用filter 函数或conv 函数求解系统输出响应的主程序。程序中要有绘制信号波形的功能。 （2）给定一个低通滤波器的差分方程为数字信号处理基础实验报告 (2)
习题集-02 数字信号处理习题答案
数字信号处理上机实验答案完整版
数字信号处理上机实验答案(全)1
相关主题
文本预览