数字信号处理期末试卷及答案 本文关键词：答案，数字信号处理，期末试卷
数字信号处理期末试卷及答案 本文简介：A一、选择题（每题3分，共5题）1、，该序列是。A.非周期序列B.周期C.周期D.周期2、序列，则的收敛域为。A.B.C.D.3、对和分别作20点DFT，得和，，，n在范围内时，是和的线性卷积。A.B.C.D.4、，，用DFT计算二者的线性卷积，为使计算量尽可能的少，应使DFT的长度N满足。A.B.
数字信号处理期末试卷及答案 本文内容：
A
一、
选择题（每题3分，共5题）
1、
，该序列是
。
A.非周期序列B.周期C.周期D.
周期
2、
序列，则的收敛域为
。
A.B.C.D.
3、
对和分别作20点DFT，得和，，，
n在
范围内时，是和的线性卷积。
A.B.C.D.
4、
，，用DFT计算二者的线性卷积，为使计算量尽可能的少，应使DFT的长度N满足
。
A.B.C.D.
5.已知序列Z变换的收敛域为｜z｜2Ωc
B.Ωs>Ωc
C.Ωs|z|>3，则该序列为（
）
A.有限长序列
B.右边序列
C.左边序列D.双边序列
4.实偶序列傅里叶变换是（
）
A.实偶序列
B.实奇序列
C.虚偶序列
D.虚奇序列
5.已知x(n)=δ(n)，其N点的DFT［x(n)］=X(k)，则X(N-1)=（
）
A.N-1
B.1
C.0D.-N+1
6.设两有限长序列的长度分别是M与N，欲通过计算两者的圆周卷积来得到两者的线性卷积，则圆周卷积的点数至少应取（
）
A.M+N
B.M+N-1
C.M+N+1D.2(M+N)
7.下面说法中正确的是（
）
A.连续非周期信号的频谱为周期连续函数
B.连续周期信号的频谱为周期连续函数
C.离散非周期信号的频谱为周期连续函数
D.离散周期信号的频谱为周期连续函数
8.下列各种滤波器的结构中哪种不是IIR滤波器的基本结构？（
）
A.直接型
B.级联型
C.频率抽样型
D.并联型
9.下列关于FIR滤波器的说法中正确的是（
）
A.FIR滤波器容易设计成线性相位特性
B.FIR滤波器的脉冲响应长度是无限的
C.FIR滤波器的脉冲响应长度是确定的
D.对于相同的幅频特性要求，用FIR滤波器实现要比用IIR滤波器实现阶数低
10.下列关于冲激响应不变法的说法中错误的是（
）
A.数字频率与模拟频率之间呈线性关系
B.能将线性相位的模拟滤波器映射为一个线性相位的数字滤波器
C.具有频率混叠效应
D.可以用于设计低通、高通和带阻滤波器
三、填空题(本大题共5小题，每空2分，共20分)。
16.线性移不变系统是因果系统的充分必要条件是________。
17.傅里叶变换的四种形式________，________，________和________。
18.使用DFT分析模拟信号的频谱时，可能出现的问题有________、栅栏效应和________。
19.下图所示信号流图的系统函数为________。
20.对于N点（N＝2L）的按时间抽取的基2FFT算法，共需要作________次复数乘和________次复数加。
四、计算题
23．(10分)考虑一个具有系统函数的稳定系统。
1）求系统的零点和极点，并作出图表示；
2）画出系统的级联型结构图。
24.(10分)有一用于频谱分析的FFT处理器，其抽样点数必须是2的整数次幂，假定没有采用任何特殊的数据处理措施，已知条件为：1）频率分辨率小于10Hz；2）信号最高频率小于4kHz。试确定以下参量：
1）最小记录长度tp；
2）最大抽样间隔T；
3）在一个记录中的最少点数N。
25．(10分)将双线性变换应用于模拟巴特沃兹滤波器，设计一个3dB截止频率的一阶数字滤波器。（注：式中模拟巴特沃兹滤波器的3dB截止频率为
Ωc）
D
一、单项选择题(每小题3分，共24分)
1、在对连续信号均匀采样时，要从离散采样值不失真恢复原信号，则采样周期Ts与信号最高截止频率fh应满足关系
A.Ts>2/fhB.Ts>1/fh
C.TsM），二者线性卷积的长度为
N点循环卷积中混叠的点有
个，循环卷积与线性卷积的关系是
5．全通系统的极零点分布特点是
三、分析计算题：（共
50分）
1．（15分）已知序列，n=0,1…,4
(1)
该序列是否可以作为线性相位FIR滤波器的单位脉冲响应？为什么？
(2)
设序列的傅立叶变换用表示，不用求，分别计算
、、、。
(3)
求与序列的线性卷积及7点循环卷积。
2．（15分）已知一因果系统的系统函数为
试完成下列问题：
(1)
系统是否稳定？为什么？
(2)
求单位脉冲响应
(3)
写出差分方程；
(4)
画出系统的极零图；
(5)
画出系统的所需存储器最少的实现结构。
3．（5分）已知模拟滤波器的传输函数：式中，a、b为常数，设因果稳定，试用脉冲响应不变法将其转换成数字滤波器。
F
一、单项选择题(在每小题的四个备选答案中，选出一个正确答案，并将正确答案的序号填在括号内。
1.若一模拟信号为带限，且对其抽样满足奈奎斯特采样定理，则只要将抽样信号通过(
)即可完全不失真恢复原信号。
A.理想低通滤波器
B.理想高通滤波器
C.理想带通滤波器
D.理想带阻滤波器
2．下列系统（其中y(n)为输出序列，x(n)为输入序列）中哪个属于线性系统？(
)
A.y(n)=x3(n)B.y(n)=x(n)x(n+2)
C.y(n)=x(n)+2D.y(n)=x(n2)
3.．设两有限长序列的长度分别是M与N，欲用圆周卷积计算两者的线性卷积，则圆周卷积的长度至少应取(
)。
A．M+NB.M+N-1
C.M+N+1D.2(M+N)
4.若序列的长度为M，要能够由频域抽样信号X(k)恢复原序列，而不发生时域混叠现象，则频域抽样点数N需满足的条件是(
)。
A.N≥M
B.N≤M
C.N≤2M
D.N≥2M
5.直接计算N点DFT所需的复数乘法次数与(
)成正比。
A.N
B.N2
C.N3
D.Nlog2N
6.下列各种滤波器的结构中哪种不是FIR滤波器的基本结构(
)。
A.直接型
B.级联型
C.并联型
D.频率抽样型
7.第二种类型线性FIR滤波器的幅度响应H(w)特点(
)：
A
关于、、偶对称
B
关于、、奇对称
C
关于、偶对称
关于奇对称
D关于、奇对称
关于偶对称
8.适合带阻滤波器设计的是：
（
）
A
N为偶数
B
N为奇数
C
N为偶数
D
N为奇数
9.以下对双线性变换的描述中不正确的是(
)。
A.双线性变换是一种非线性变换
B.双线性变换可以用来进行数字频率与模拟频率间的变换
C.双线性变换把s平面的左半平面单值映射到z平面的单位圆内
D.以上说法都不对
10.关于窗函数设计法中错误的是：
A窗函数的截取长度增加，则主瓣宽度减小；
B窗函数的旁瓣相对幅度取决于窗函数的形状，与窗函数的截取长度无关；
C为减小旁瓣相对幅度而改变窗函数的形状，通常主瓣的宽度会增加；
D窗函数法不能用于设计高通滤波器；
二、填空题(每空2分，共20分)
1.
用DFT近似分析连续信号频谱时,_________效应是指DFT只能计算一些离散点上的频谱。
2.有限长序列X(z)与X（k）的关系
X（k）与的关系
3.下图所示信号流图的系统函数为：
4.如果通用计算机的速度为平均每次复数乘需要４μs，每次复数加需要1μs，则在此计算机上计算210点的基2FFT需要__________级蝶形运算,总的运算时间是__________μs。
5.单位脉冲响应不变法优点,缺点____________,适合_______________________滤波器设计
6.已知FIR滤波器具有线性相位，则a＝______,冲激响应h（2）＝___,相位___
7.的周期__________________
8.用频率采样法设计数字滤波器，对第二类型相位滤波器H(k)应具有的约束条件：幅值__________，相位_____________
9．两序列h(n)=δ(n)+2δ(n-1)+3δ(n-2),x(n)=δ(n)+δ(n-1)，两者的线性卷积为y(n)，则y(2)_____
________；若两者3点圆周卷积为y1(n)，则y1(0)＝__________________y1(2)＝__________________。
三
计算题
1.
有一个线性移不变的系统，其系统函数为：
1）用直接型结构实现该系统
2)讨论系统稳定性，并求出相应的单位脉冲响应
4．试用冲激响应不变法与双线性变换法将以下模拟滤波器系统函数变换为数字滤波器系统函数：
H(s)=其中抽样周期T=1s。
G
三、有一个线性移不变的因果系统，其系统函数为：
1用直接型结构实现该系统
2)讨论系统稳定性，并求出相应的单位脉冲响应
七、用双线性变换设计一个三阶巴特沃思数字低通虑波器，采样频率为（即采样周期为）,其3dB截止频率为。三阶模拟巴特沃思滤波器为：
答案
一、
选择题（10分，每题1分）
1.A
2.D
3.B
4.A
5.B
6.C
7.C
8.D
9.D
10.D
二、填空题（共25分
3、4、7、9每空2分；其余每空1分）
1.栅栏效应
2.x(z)|z=wN-k
x(k)＝X(ejw)|w＝
3.
4.
8
6144us
5.线性相位
频谱混迭、低通带通
6.
2、5
、－2w
7、14
9.
、
10、5、
4
、
5
三计算题
1.（15分）
解1)
……………………………
2分
当时：
收敛域包括单位圆……………………………6分
系统稳定系统。……………………………….10分
………………………………12分
………………………………….15分
4.（10分）解：
………………1分
……………………3分
……………5分
2)……8分
……………………………
10分
三、（15）
1.解1)
……………………………
2分
2)当时：
收敛域包括单位圆……………………………6分
系统稳定系统。……………………………….10分
………………………………12分
………………………………….15分
七、（12分）解：
………………………………………3分
………………………………………5分
……………………………8分
篇2：数字信号处理习题集大题及答案
数字信号处理习题集大题及答案 本文关键词：习题集，大题，答案，数字信号处理
数字信号处理习题集大题及答案 本文简介：1设序列x(n)={4，3，2，1}，另一序列h(n)={1，1，1，1}，n=0,1,2,3（1）试求线性卷积y(n)=x(n)*h(n)（2）试求6点圆周卷积。（3）试求8点圆周卷积。解：1．y(n)=x(n)*h(n)={4,7,9,10,6,3,1}2．6点圆周卷积={5,7,9,10,6,
数字信号处理习题集大题及答案 本文内容：
1设序列x(n)={4，3，2，1}
，
另一序列h(n)
={1，1，1，1}，n=0,1,2,3
（1）试求线性卷积
y(n)=x(n)*h(n)
（2）试求6点圆周卷积。
（3）试求8点圆周卷积。
解：1．y(n)=x(n)*h(n)={4,7,9,10,6,3,1}
2．6点圆周卷积={5,7,9,10,6,3}
3．8点圆周卷积={4,7,9,10,6,3,1,0}
2二．数字序列
x(n)如图所示.
画出下列每个序列时域序列：
(1)
x(n-2);
(2)x(3-n);(3)x[((n-1))6],(0≤n≤5);(4)x[((-n-1))6],(0≤n≤5);
3．已知一稳定的LTI
系统的H(z)为
试确定该系统H(z)的收敛域和脉冲响应h[n]。
解:
系统有两个极点，其收敛域可能有三种形式，|z|2
因为稳定，收敛域应包含单位圆，则系统收敛域为：0.5(5+3-1),所以y3(n)=
x(n)⑧h(n)＝{-15,4,-3,13,-4,3,2，0}
y3(n)与y(n)非零部分相同。
6一个因果线性时不变离散系统，其输入为x[n]、输出为y[n]，系统的差分方程如下：
y（n）-0.16y(n-2)=
0.25x(n-2)＋x(n)
(1)
求系统的系统函数
H(z)=Y(z)/X(z);
(2)
系统稳定吗?
(3)
画出系统直接型II的信号流图;
(4)
画出系统幅频特性。
解：(1)方程两边同求Z变换：
Y(z)-0.16z-2Y(z)=
0.25z-2X(z)＋X(z)
(2)系统的极点为：0.4和－0.4,在单位圆内，故系统稳定。
(3)
(4)
7．如果需要设计FIR低通数字滤波器，其性能要求如下：
(1)阻带的衰减大于35dB,(2)过渡带宽度小于p/6.
请选择满足上述条件的窗函数，并确定滤波器h(n)最小长度N
解：根据上表，我们应该选择汉宁窗函数，
8两个有限长的复序列x[n]和h[n]，其长度分别为N
和M，设两序列的线性卷积为y[n]=x[n]*h[n]，回答下列问题：.
(1)
序列y[n]的有效长度为多长？
(2)
如果我们直接利用卷积公式计算y[n]
，那么计算全部有效y[n]的需要多少次复数乘法？
(3)
现用FFT
来计算y[n]，说明实现的原理，并给出实现时所需满足的条件，画出实现的方框图，计算该方法实现时所需要的复数乘法计算量。
解：(1)
序列y[n]的有效长度为：N+M-1；
(2)
直接利用卷积公式计算y[n]，
需要MN次复数乘法
(3)
需要次复数乘法。
9用倒序输入顺序输出的基2
DIT-FFT
算法分析一长度为N点的复序列x[n]
的DFT，回答下列问题：
(1)
说明N所需满足的条件，并说明如果N不满足的话，如何处理？
(2)
如果N=8,那么在蝶形流图中，共有几级蝶形？每级有几个蝶形？确定第2级中蝶形的蝶距(dm)和第2级中不同的权系数(WNr
)。
(3)
如果有两个长度为N点的实序列y1[n]和y2
[n]，能否只用一次N点的上述FFT运算来计算出y1[n]和y2
[n]的DFT，如果可以的话，写出实现的原理及步骤，并计算实现时所需的复数乘法次数；如果不行，说明理由。
解(1)N应为2的幂，即N＝2m，（m为整数）；如果N不满足条件，可以补零。
(2)3级，4个，蝶距为2，WN0
，WN2
(3)
y[n]=y1[n]+jy2[n]
10已知系统函数，求其差分方程。
解：
11已知，画系统结构图。
解：
直接型I：
直接型II：
级联型：
并联型：
12若x
(n)=
{3，2，1，2，1，2
}，0≤n≤5，
1)
求序列x(n)的6点DFT，X
(k)=？
2)
若，试确定6点序列g(n)=?
3)
若y(n)
=x(n)⑨x(n)，求y(n)=？
1）
2）
3）
13用DFT对连续信号进行谱分析的误差问题有哪些？
答：混叠失真；截断效应（频谱泄漏）；栅栏效应
14画出模拟信号数字化处理框图，并简要说明框图中每一部分的功能作用。
第1部分：滤除模拟信号高频部分；第2部分：模拟信号经抽样变为离散信号；第3部分：按照预制要求对数字信号处理加工；第4部分：数字信号变为模拟信号；第5部分：滤除高频部分，平滑模拟信号。
15简述用双线性法设计IIR数字低通滤波器设计的步骤。
答：确定数字滤波器的技术指标；将数字滤波器的技术指标转变成模拟滤波器的技术指标；按模拟滤波器的技术指标设计模拟低通滤波器；将模拟低通滤波器转换成数字低通滤波器。
16
8点序列的按时间抽取的（DIT）基-2
FFT如何表示？
17已知，求x(n)。（6分）
解：由题部分分式展开
求系数得
A=1/3
，
B=2/3
所以
（3分）
收敛域÷z÷>2，故上式第一项为因果序列象函数，第二项为反因果序列象函数，
则
（3分）
18写出差分方程表示系统的直接型和级联型结构。（8分）
解：(8分)
19计算下面序列的N点DFT。
（1）（4分）
（2）
（4分）
解：（1）
（4分）
（2）
（4分）
20设序列x(n)={1，3，2，1；n=0,1,2,3
}，另一序列h(n)
={1，2，1，2；n=0,1,2,3}，
（1）求两序列的线性卷积
yL(n)；
（4分）
（2）求两序列的6点循环卷积yC(n)。
（4分）
（3）说明循环卷积能代替线性卷积的条件。（2分）
解：（1）
yL(n)={1，5，9，10，10，5，2；n=0,1,2…6}
（4分）
（2）
yC(n)=
{3，5，9，10，10，5；n=0,1,2,4,5}
（4分）
（3）c≥L1+L2-1
（2分）
21设系统由下面差分方程描述：
（1）求系统函数H（z）；（2分）
（2）限定系统稳定，写出H（z）的收敛域，并求出其单位脉冲响应h(n)。（6分）
解：（1）
（2分）
（2）
（2分）；
（4分）
23求
，
的反变换。
25有一频谱分析用的FFT处理器，其抽样点数必须是2的整数幂，假定没有任何特殊的数据处理措施，已给条件为频率分辨力≤10Hz,信号最高频率≤4kHz.
试确定一下参量：（1）最小记录长度T0；
（2）抽样点间的最大时间间隔（即最小抽样频率）；
（3）在一个记录中最小点数N0
解：(1)T0≥=0.1(s)
(2)fs＞2fh
T＜==0.125×(s)
（3）N＞==800
N===1024＞800
26有一个线性移不变的系统，其系统函数为：
1）用直接型结构实现该系统
2)讨论系统稳定性，并求出相应的单位脉冲响应
解1)
……………………………
2分
当时：
收敛域包括单位圆……………………………6分
系统稳定系统。……………………………….10分
………………………………
………………………………….15分
27试用冲激响应不变法与双线性变换法将以下模拟滤波器系统函数变换为数字滤波器系统函数：
H(s)=其中抽样周期T=1s。
解：
……………………3分
……………5分
2)……8分
……………………………
10分
28用长除法、留数定理法、部分分式法分别求以下X(Z)的Z反变换：
(1)
;
(2)
;
(3)
a.
长除法
b．留数法
c．部分分式法
30设系统差分方程
y(n)=ay(n-1)+x(n)
其中x(n)为输入，y(n)为输出。当边界条件选为y(-1)=0时，是判断系统是否线性的、移不变的。
①
令
则
同样可求得
所以
②令
则
同样可求得
所以
因为与为移1位关系，而且与也是移1位关系，所以在y(-1)=0条件下，系统是移不变系统。
③令
n=0时，
综上，可得
所以系统是线性系统。
31用级联型结构实现以下系统函数，试问一共能构成几种级联型网络，并画出结构图。
32已知，求x(n)。
解：由题部分分式展开
求系数得
A=1/3
，
B=2/3
所以
（3分）
收敛域÷z÷>2，故上式第一项为因果序列象函数，第二项为反因果序列象函数，
则
33．写出差分方程表示系统的直接型和级联型结构。（8分）
34计算下面序列的N点DFT。（1）（2）
1
2
35设系统由下面差分方程描述：（1）求系统函数H（z）；（2分）
（2）限定系统稳定，写出H（z）的收敛域，并求出其单位脉冲响应h(n)。（6分）
（1）
（2）
篇3：维纳滤波器设计随机信号处理
维纳滤波器设计随机信号处理 本文关键词：滤波器，信号处理，随机，维纳，设计
维纳滤波器设计随机信号处理 本文简介：随机信号处理实验报告维纳滤波器设计随机信号处理学号：姓名：实验一维纳滤波器设计1实验内容设计一个维纳滤波器：（1）产生三组观测数据，首先根据产生信号，将其加噪，（信噪比分别为），得到观测数据。（2）估计，的AR模型参数。假设信号长度为L，AR模型阶数为N，分析实验结果，并讨论改变L，N对实验结果的影
维纳滤波器设计随机信号处理 本文内容：
随机信号处理实验报告
维纳滤波器设计
随机信号处理
学号：
姓名：
实验一
维纳滤波器设计
1实验内容
设计一个维纳滤波器：
（1）
产生三组观测数据，首先根据产生信号，将其加噪，（信噪比分别为），得到观测数据。
（2）
估计，的AR模型参数。假设信号长度为L，AR模型阶数为N，分析实验结果，并讨论改变L，N对实验结果的影响。
2实验原理
滤波目的是从被噪声污染的信号中分离出有用的信号来，最大限度地抑制噪声。对信号进行滤波的实质就是对信号进行估计。滤波问题就是设计一个线性滤波器，使得滤波器的输出信号是期望响应的一个估计值。下图就是观测信号的组成和信号滤波的一般模型。
观测信号
信号滤波的一般模型
维纳滤波解决从噪声中提取信号的滤波问题，并以估计的结果与真值之间的误差均方值最小作为最佳准则。它根据估计信号的当前值，它的解以系统的系统函数或单位脉冲形式给出，这种系统常称为最佳线性滤波器。
维纳滤波器设计的任务就是选择，使其输出信号与期望信号误差的均方值最小。
假设滤波系统是一个线性时不变系统，它的和输入信号都是复函数，设
考虑系统的因果性，可得到滤波器的输出
设期望信号，误差信号及其均方误差分别为
要使均方误差为最小，需满足：
整理得，等价于
上式说明，均方误差达到最小值的充要条件使误差信号与任一进入估计的输入信号正交，这就是正交性原理。
将展开，得
整理得
等价于
此式称为维纳-霍夫（Wiener-Holf）方程。解此方程可得到最优权系数，此式是Wiener滤波器的一般方程。
定义
则维纳-霍夫方程可写成矩阵形式
求逆，得
此式表明，已知期望信号与观测数据的互相关函数及观测数据的自相关函数时，可以通过矩阵求逆运算，得到维纳滤波器的最佳解。
3实验结果及分析
（1）当L=200,N=6
信噪比为20dB的滤波效果
信噪比为10dB的滤波效果
信噪比为6dB的滤波效果
（2）当L=200,N=60
信噪比为20dB的滤波效果
信噪比为10dB的滤波效果
信噪比为6dB的滤波效果
（3）当L=600,N=6
信噪比为20dB的滤波效果
信噪比为10dB的滤波效果
信噪比为6dB的滤波效果
实验分析：别取信号长度为200、600，滤波器长度为6、60，加噪信噪比为20dB、10dB、6dB，组合进行实验。每组实验得到的最小均方误差统计如下表。由此表可以看出，信号长度越长，最小均方误差（绝对值）越大，精度越差；在信噪比较大（误差影响较小）的滤波过程中，滤波器长度约长，最小均方误差（绝对值）越小，精度越好。
表1
最小均方误差统计表
加噪
信号
信噪比
L=200
N=6
L=200
N=60
L=600
N=6
20dB
-1785.3971
-1574.1107
-6888.4103
10dB
-845.5714
-1727.5260
-10219.3397
6dB
-1213.6367
-2135.4596
-7743.9358
对于相同信号和滤波器（这里取L=200,N=6），信噪比越大，最小距离误差约小；而当信噪比较小时，信号与噪声值接近，导致滤波效果受到影响，最小距离误差变大。
4源代码
clear;
clc;
%初始化变量
L=200;
%信号长度
N=6;
%滤波器的阶次
a=0.96;
wn=randn(L,1);
%wn为用于生成信号的噪声信号，随机生成一个L*1矩阵，生成矩阵的元素值在%区间(0.0,1.0)之间
sn=zeros(L,1);
%sn为信号，生成一个L*1的零矩阵
hn=zeros(N,1);
%hn为系统单位脉冲响应
生成一个N*1的零矩阵
rxx=zeros(N,1);
%rxx为自相关函数，生成一个N*1的零矩阵
rxd=zeros(N,1);
%rxd为互相关函数，生成一个N*1的零矩阵
yn=zeros(L,1);
%yn为输出信号，生成一个L*1的零矩阵
xt=zeros(L+N,1);
%生成一个(L+N)*1的零矩阵
gn=zeros(L,1);
%gn为yn与sn最小距离误差信号，生成一个L*1的零矩阵
%根据给定公式s(n)=as(n-1)+w(n),生成理想信号
for
i=2:L
sn(i,1)=a*sn(i-1,1)+wn(i,1);
end
sn(1,1)=wn(1,1);
subplot(2,2,1);
plot(sn,r
),axis([0
200
-10
10]),xlabel(
时间
),ylabel(
幅度
),title(
sn
)
%生成期望信号方差cd
cd=(var(wn))/(1-a^2);
%对信号加噪
x1=awgn(sn,20);
x2=awgn(sn,10);
x3=awgn(sn,6);
subplot(2,2,2)
plot(x3,g
),axis([0
200
-10
10]),xlabel(
时间
),ylabel(
幅度
),title(
x3
);
%生成输入信号与理想信号的互相关函数,此处x1为输入信号，sn为期望信号
for
i=1:N,for
m=i:1:L,rxd(i,1)=rxd(i,1)+x3(m,1)*sn(m-i+1,1);
end
end
%生成输入信号的自相关函数
for
i=1:N,for
m=i:1:L,rxx(i,1)=rxx(i,1)+x3(m,1)*x3(m-i+1,1);
end
end
%将自相关函数生成托普勒斯矩阵
rxx1=toeplitz(rxx);
%生成逆矩阵
irxx=inv(rxx1);
%生成滤波器系数h(n)
hn=irxx*rxd;
for
i=1:L
xt(i+N,1)=x3(i,1);
end
%实现滤波
for
i=1:L,for
m=1:N,yn(i,1)=yn(i,1)+xt(i+N+1-m,1)*hn(m,1);
end
end
%计算最小均方误差信号en
en=0;
en=cd-(rxd
)*hn;
%生成最小距离误差信号gn
gn=yn-sn;
%画出滤波后的信号时域图
subplot(2,2,3);
plot(yn),axis([0
200
-10
10]),xlabel(
时间
),ylabel(
幅度
),title(
yn
);
%画出理想信号与输出信号对比图
subplot(2,2,4);
plot(sn,r
),axis([0
200
-10
10]),xlabel(
时间
),ylabel(
幅度
),title(
sn与yn对比
);
hold
on;
plot(yn,b
),axis([0
200
-10
10]);
hold
off;
%画出最小距离误差信号图
figure;
plot(gn),axis([0
200
-2
2]),xlabel(
时间
),ylabel(
幅度
),title(
gn
);
10