数字信号处理应用的作业

合集下载

数字信号处理作业

数字信号处理作业
(2) 用h1(n)和h2(n)分别构成的低通滤波器是否具有线性 相位?群延时为多少?
题8图
解: (1) 由题8图可以看出h2(n)与h1(n)是循环移位 关系:
h2(n)=h1((n+4))8R8(n) 由DFT的循环移位性质可得
H 2 (k) W8k 4 H1(k) e jk H1(k) (1)k H1(k) | H 2 (k) || W8k 4 H1 (k) || H1 (k ) |
(2) 由题8图可知, h1(n)和h2(n)均满足线性相位条件: h1(n)=h1(N-1-n) h2(n)=h2(N-1-n)
x(t)=cos(2πfct+jc)[1+cos(2πfmt+jm)]
所以, 已调AM信号x(t) 只有3个频率: fc、 fc+fm、 fc-fm。 x(t)的最高频率fmax=1.1 kHz, 频率分辨率F≤100 Hz(对本题 所给单频AM调制信号应满足100/F=整数, 以便能采样到这
三个频率成分)。 故
(1)
Tp m in

1 F
1 100
0.01s
10ms
(2) Fsmin 2 fmax 2.2 kHz
(3)
N m in

Tp Tm a x

Tp

f m in
10 103

2.2 103

22
5-1. 已知系统用下面差分方程描述:
y(n)=3 y(n 1)-1 y(n 2)+x(n) 1 x(n 1)

j
3
2

j 3
j 3
H(z)
3
3
1 e z

数字信号处理实验作业完全版

数字信号处理实验作业完全版

实验1:理想采样信号的序列,幅度谱,相位谱,以及改变参数后的图像。

源程序: clc;n=0:50;A=444.128;a=50*sqrt(2.0*pi;T=0.001;w0=50*sqrt(2.0*pi;x=A*exp(-a*n*T.*sin(w0*n*T;close allsubplot(3,2,1;stem(x,’.’;title('理想采样信号序列';k=-25:25;W=(pi/12.5*k;X=x*(exp(-j*pi/12.5.^(n'*k;magX=abs(X;s ubplot(3,2,2;stem(magX,’.’;title('理想采样信号序列的幅度谱';angX=angle(X;subplot(3,2,3;stem(angX;title('理想采样信号序列的相位谱'n=0:50;A=1;a=0.4,w0=2.0734;T=1; x=A*exp(-a*n*T.*sin(w0*n*T;subplot(3,2,4;stem(x,’.’; title('理想采样信号序列'; k=-25:25; W=(pi/12.5*k;X=x*(exp(-j*pi/12.5.^(n'*k; magX=abs(X; subplot(3,2,5; stem(magX,’.’title('理想采样信号序列的幅度谱';0204060-2000200理想采样信号序列020406005001000理想采样信号序列的幅度谱0204060-505理想采样信号序列的相位谱0204060-11理想采样信号序列020406012理想采样信号序列的幅度谱上机实验答案:分析理想采样信号序列的特性产生在不同采样频率时的理想采样信号序列Xa(n,并记录各自的幅频特性,观察频谱‚混淆‛现象是否明显存在,说明原因。

源程序:A=444.128;a=50*pi*sqrt(2.0;W0=50*pi*sqrt(2.0;n=-50:1:50; T1=1/1000;Xa=A*(exp(a*n*T1.*(sin(W0*n*T1;subplot(3,3,1;plot(n,Xa;title('Xa序列';xlabel('n';ylabel('Xa';k=-25:25;X1=Xa*(exp(-j*pi/12.5.^(n'*k;subplot(3,3,2; stem(k,abs(X1,'.';title('Xa的幅度谱';xlabel('k';ylabel('〃幅度';subplot(3,3,3;stem(k,angle(X1,'.';title('Xa的相位谱';xlabel('k';ylabel('相位';T2=1/300;Xb=A*(exp(a*n*T2.*(sin(W0*n*T2;subplot(3,3,4;plot(n,Xb;title('Xb序列';xlabel('n';ylabel('相位';k=-25:25;X2=Xb*(exp(-j*pi/12.5.^(n'*k;subplot(3,3,5; stem(k,abs(X2,'.'; title('Xb 的幅度谱';xlabel('k';ylabel('〃幅度';subplot(3,3,6;stem(k,angle(X2,'.'; title(' Xb 的相位谱';xlabel('k';ylabel('相位';T3=1/200;Xc=A*(exp(a*n*T3.*(sin(W0*n*T3; subplot(3,3,7;plot(n,Xc;title('Xc 序列'; xlabel('n';ylabel('Xc';k=-25:25;X3=Xc*(exp(-j*pi/12.5.^(n'*k;subplot(3,3,8; stem(k,abs(X3,'.'; title('Xc 的幅度谱'; xlabel('k';ylabel('幅度';subplot(3,3,9;stem(k,angle(X3,'.'; title('Xc 的相位谱';xlabel('k';ylabel('相位';-50050-5057X a 序列n X a-500500128X a 的幅度谱k 幅度-50050-55X a 的相位谱k相位-50050-50518X b 序列n 相位-50050051018X b 的幅度谱k 幅度-50050-55X b 的相位谱k相位-50050-505x 1026X c 序列nX c-500500510x 1026X c 的幅度谱k幅度-50050-505X c 的相位谱k相位由图可以看出:当采样频率为1000Hz时,采样序列在折叠频率附近处,无明显混叠。

数字信号处理作业答案(参考版-第一章)

数字信号处理作业答案(参考版-第一章)

1-2习题1-2图所示为一个理想采样—恢复系统,采样频率Ωs =8π,采样后经过理想低通G jΩ 还原。

解:(1)根据余弦函数傅里叶变换知:)]2()2([)]2[cos(πδπδππ-Ω++Ω=t F ,)]6()6([)]6[cos(πδπδππ-Ω++Ω=t F 。

又根据抽样后频谱公式:∑∞-∞=∧Ω-Ω=Ωk s a a jk j X T j X )(1)(,得到14T= ∑∞-∞=∧--Ω+-+Ω=Ωk a k k j X )]82()82([4)(1ππδππδπ∑∞-∞=∧--Ω+-+Ω=Ωk a k k j X )]86()86([4)(2ππδππδπ所以,)(1t x a ∧频谱如下所示)(2t x a ∧频谱如下所示(2))(1t y a 是由)(1t x a ∧经过理想低通滤波器)(Ωj G 得到,)]2()2([)()()]([11πδπδπ-Ω++Ω=ΩΩ=∧j G j X t y F a a ,故)2cos()(1t t y a π=(4π) (4π) (4π)(4π)(4π) (4π) Ω-6π-10π-2π 2π0 6π10π)(1Ω∧j X a Ω10π-10π -6π-2π 0 2π6π-14π 14π(4π)(4π) (4π)(4π) (4π) (4π)(4π) (4π))(2Ω∧j X a同理,)]2()2([)()()]([22πδπδπ-Ω++Ω=ΩΩ=∧j G j X t y F a a 故)2cos()(2t t y a π=(3)由题(2)可知,无失真,有失真。

原因是根据采样定理,采样频率满足信号)(1t x a 的采样率,而不满足)(2t x a 的,发生了频谱混叠。

1-3判断下列序列是否为周期序列,对周期序列确定其周期。

(1)()5cos 86x n A ππ⎛⎫=+ ⎪⎝⎭(2)()8n j x n eπ⎛⎫- ⎪⎝⎭=(3)()3sin 43x n A ππ⎛⎫=+ ⎪⎝⎭解:(1)85πω=,5162=ωπ为有理数,是周期序列,.16=N (2)πωπω162,81==,为无理数,是非周期序列; (3)382,43==ωππω,为有理数,是周期序列,8=N 。

数字信号处理作业(附答案)1-郑佳春

数字信号处理作业(附答案)1-郑佳春

习题一1.2 在过滤限带的模拟数据时,常采用数字滤波器,如图中T 表示采样周期(假设T 足够小,足以防止混迭效应),把从)()(t y t x 到的整个系统等效为一个模拟滤波器。

(a ) 如果kHz rad n h 101,8)(=π截止于,求整个系统的截止频率。

(b )对于kHz T 201=,重复(a )的计算。

解 (a )因为当0)(=≥ωπωj e H rad 时,在数-模变换中)(1)(1)(Tj X T j X T e Y a a j ωω=Ω= 所以)(n h 得截止频率πω=c对应于模拟信号的角频率c Ω为8π=ΩT c因此 Hz Tf c c 6251612==Ω=π 由于最后一级的低通滤波器的截止频率为Tπ,因此对T 8π没有影响,故整个系统的截止频率由)(ωj eH 决定,是625Hz 。

(b )采用同样的方法求得kHz 201=,整个系统的截止频率为 Hz Tf c 1250161==1.3 一模拟信号x(t)具有如图所示的带通型频谱,若对其进行采样,试确定最佳采样频率,并绘制采样信号的频谱。

解:由已知可得:==35,25H L f kHz f kHz ,10k H L B f f Hz =-=,为使无失真的恢复原始信号,采样频率应满足:2f 21c c s B f Bf m m+-≤≤+且220s f B kHz >=、0/12H m f B ≤≤-=⎡⎤⎣⎦ 当m=1时,2501c s f Bf kHz -==,满足: 3550s kHz f kHz ≤≤ 当m=2时,2252c s f Bf kHz -==,满足:23.325s kHz f kHz ≤≤ 故最佳采样频率为25kHz,采样信号的频谱图如下图所示 :1.5 判断下面的序列是否是周期的,若是周期的,确定其周期,并绘制一个周期的序列图(1)16()cos()58x n A n ππ=-,A 是常数 解:2251685N wπππ===,所以x(n)是周期的,且最小正周期为5 1285()cos()40n x n A π-= 绘图:方法一:计算法 当n=0时,1()cos()8x n A π-==0.99A 当n=1时,123()cos()40x n A π==3cos()40A ππ+=0.998A 当n=2时,251()cos()40x n A π==11cos()40A π=0.863A 当n=3时,379()cos()40x n A π==19cos()40A ππ+=0.996A 当n=4时,507()cos()40x n A π==27cos()40A π=0.999A 当n=5时,635()cos()40x n A π==6405cos()40A π-= 方法二:Matlab 法> xn=cos(16*n*pi/5-pi/8); > plot(n,xn) > n=0:4;> xn=cos(16*n*pi/5-pi/8); > plot(n,xn)> plot(n,xn,'.');grid;1.6对如下差分方程所述系统,试分析其线性特性与时变特性(1)()()2(1)y n x n x n =-- (3)2()4()y n x n =解:(1)121212[()()]()()2{(1)(1)}T ax n bx n ax n bx n ax n bx n +=+--+- 121122[()][()]()2(1)()2(1)T ax n T bx n ax n ax n bx n bx n +=--+--故 1212[()()][()][()]T ax n bx n T ax n T bx n +=+所以y(n)为线性又0000[()]()2(1)()T x n n x n n x n n y n n -=----=- 所以y(n)为时不变(3)2222121122[()()]4()8()()4()T ax n bx n a x n abx n x n b x n +=++ 221212[()][()]4()4()T ax n T bx n ax n bx n +=+ 故1212[()()][()][()]T ax n bx n T ax n T bx n +≠+ y(n)为非线性又2000[()]4()()T x n n x n n y n n -=-=- 故y(n)为时不变1.7试判断如下算法是否是因果的?是否是稳定的?并说明理由。

东南大学仪科数字信号处理作业

东南大学仪科数字信号处理作业
bn= f1 (t)= 2)因为 f2(t)为奇函数,bn=0 a0=
an= f(t)= +
5. 将下列信号在区间(0,1)中展开为指数形式的傅里叶级数: 1)f1(t)=t4 2)f2(t)=e2t 解:
T=1,w1
1)F(n)=
f(t)=
精品
. 2)F(n)=
精品
. f(t)= 6. 已知如图 2-17 所示的信号 f(t),利用微分性质求该信号的傅立叶变换 F(w)。
[ 精品
.
如有侵权请联系告知删除,感谢你们的配合!
精品
1)
ห้องสมุดไป่ตู้
2) (2+2t)f(t-1) 3)
4) f(t)* 5) 解:
1)
6) f(t)sin[w0(t+a)]
2)
3) 4) 5) 6)f(t)sin[w0(t+a)]=f(t)sin(w0t)cos(w0a)+f(t)cos(w0t)sin(w0a)
2. 已知如图 2-15 所示的信号 f(t),求:1)指数形式与三角形式的傅里叶变换级数;2) 精品
0
答案: 由图可知,
T/2
T
0
图 2-18
T/2
T
10. 求下列频谱函数对应的时间函数:
1)
2)sin(w/2) 3)
4)sin(2w)cos(w) 5) 答案: 1)
2)
3) 精品
. 精品
. 4) 5)
11. 已知
如图 2-19 所示,求其傅立叶反变换 f1(t)。
答案:
1
0
1
2
图 2-19
012 -2
1f1t2t2f2t05t1因为f1t为奇函数a00an0bncos?f1t2因为f2t为奇函数bn0a0w1ntdtw1ntdtw1ntdt已知如图217所示的信号ft利用微分性质求该信号的傅立叶变换fw

数字信号处理系统作业

数字信号处理系统作业

数字信号处理系统作业《数字信号处理系统》实践任务报告学号:3110411072姓名:王伟东班级:11级信计2班⼀、问题提出:1.1实验背景:滤波器是⼀种⽤来消除⼲扰杂讯的器件,将输⼊或输出经过过滤⽽得到纯净的直流电。

对特定频率的频点或该频点以外的频率进⾏有效滤除的电路,就是滤波器,其功能就是得到⼀个特定频率或消除⼀个特定频率。

经典滤波器按通频带分类可以分为低通(LP)、⾼通(HP)、带通(BP)、带阻(BS),按处理信号类型可以分为模拟滤波器和数字滤波器。

通过Matlab语⾔可以快捷的设计出有软件组成的数字滤波器,其中FIR滤波器(即有限冲击响应滤波器)最⼤的优点就是在满⾜幅频特性的同时,还可以获得严格的线性相位特性。

1.2实验要求:某⼀数字信号由500Hz、1000Hz和1500Hz三个频率组成,采样频率为4000Hz,请1.设计⼀个低通滤波器,将上述数字信号送⼊滤波器后只留下500Hz频率分量;2.设计⼀个带通滤波器,将上述数字信号送⼊滤波器后只留下1000Hz频率分量;3.设计⼀个⾼通滤波器,将上述数字信号送⼊滤波器后只留下1500Hz频率分量;4.设计⼀个带阻滤波器,将上述数字信号送⼊滤波器后滤除1000Hz频率分量;要求从时域和频域两个⾓度说明滤波前后信号的变化。

1.3实验⽬的:1. 熟悉FIR 滤波器的滤波原理2. 熟悉FIR 滤波器的汇编实现3.熟悉理想采样的性质,了解信号采样前后的频谱变化,加深对采样定理的理解4. 学习⽤MA TLAB 设计滤波器⼆、试验⽅法和过程2.1实验原理:数字滤波器的设计,是对提出的设计要求给出响应的性能指标,再通过计算,使物理可实现的实际滤波器频率响应特性,逼近给出的频率响应特性。

设计完成后,可根据计算结果在MATLAB或DSP上实现。

应⽤Matlab 设计FIR 滤波器,即根据给定的性能指标,设计⼀个H (z 、),使其逼近这⼀指标,进⽽计算并确定滤波器的系数b (n ),再将所设计滤波器的幅频响应、相频响应曲线作为输出,与设计要求⽐较,对设计的滤波器进⾏优化,从⽽达到滤波的⽬的。

数字信号处理大作业

数字信号处理大作业

数字信号处理大作业班级学号姓名陈志豪一.要求本次作业要求对一段音乐进行处理,该音乐包含了蜂鸣噪声,根据该段音乐,我们需处理以下问题:1. 利用matlab软件对audio1211.wav音频信号进行数字信号采样,分别对采样后的信号进行时/频域分析,并提供仿真图和分析说明;2. 设计合理的数字滤波器,滤去音频信号中的蜂鸣音,给出详细设计流程,并提供频域仿真图和分析说明;3. 将数字滤波后的数字信号转换成wav格式音频文件二.分析(1)通过播放所给音乐文件,很明显能听出wav文件中包含蜂鸣噪音,所以我们应该先分析频谱。

在matlab下可以用函数wavread/audioread读入语音信号进行采样,我们可以通过wavread得到声音数据变量x和采样频率fs、采样精度nbits,在读取声音信号之后,利用读出的采样频率作为参数,这段音频读出的采样精度为16,fs为44100hz,所以我们将此后采集时间、fft的参数设置为fs,也就是44100hz。

最后我们通过plot函数绘制出了音频信号与时间的关系图pic1,使用fft函数进行fft处理。

处理后的信号频谱pic2,如下所示图1.音频信号与时间的关系图从图1横坐标我们看到t在9-10s之间截止,与我们在音乐播放器中显示的时间一致。

图2.fft之后得到的频域分析结果图3.噪声读取图2为运用fft后得到的处理结果,可以从中读取到,在293.7hz、4671hz附近幅值突然增大,可以确定为噪声干扰。

所以我们应该针对频率附近进行滤波。

如果针对性进行滤波处理,应该使用低通滤波器进行处理,去除这部分的噪音。

之后需要选定滤波器并进行程序设计,在4671hz附近进行滤波,去除蜂鸣杂音。

(2)我们需要对蜂鸣音进行除去,自然需要用到滤波器。

所以第二步我们需要设计滤波器并给出详细流程。

在第一问的频谱分析中,通过FFT我们已经知道噪音所在,所以我们需要针对这个问题设计参数。

在这里我们选用巴特沃斯低通滤波器进行处理,我们需要设定好的参数有通带边界频率、阻带边界频率、通带最大衰减和通过阻带的最小分贝数(由buttord 在matlab定义得)。

23春-数字信号处理-在线作业

23春-数字信号处理-在线作业

23春-数字信号处理-在线作业交卷时间2023-05-14 12:28:41一、单选题(每题3分,共20道小题,总分值60分)1.由模拟信号进行采样得到时域离散信号时,同样要满足定理,(3分)采样位移反折对称正确答案A您的答案是A回答正确展开2.序列x(n)的部分x e(n)对应着X(e jω)的实部X R(e jω)。

(3分)对称共轭对称反对称共轭反对称正确答案B您的答案是D回答错误展开3.栅栏效应的存在,有可能漏掉的频谱分量。

(3分)大小高低正确答案A您的答案是A回答正确展开4.如果信号的自变量和函数值均取离散值,则称为。

(3分)模拟信号数字信号离散信号抽样信号正确答案B您的答案是C回答错误展开5.采用按时间抽取的基-2 FFT算法计算N=8点DFT,需要计算______次复数乘法(3分)8165664正确答案D您的答案是D回答正确展开6.因果(可实现)系统其系统函数H(z)的收敛域一定包含点。

(3分)12∞正确答案D您的答案是D回答正确展开7.如果信号的自变量和函数值都取连续值,则称这种信号为或者称为时域连续信号,例如语言信号、温度信号等(3分)模拟信号数字信号离散信号抽样信号正确答案A您的答案是A回答正确展开8.由傅里叶变换理论知道,若信号的频谱有限宽,则其持续时间必然为。

(3分)有限长无限长不确定正确答案B您的答案是B回答正确展开9.序列x(n)的部分x o(n)对应着X(e jω)的虚部(包括j)。

(3分)对称共轭对称反对称共轭反对称正确答案D您的答案是D回答正确展开10.对连续信号进行谱分析时,首先要对其采样,变成时域后才能用DFT(FFT)进行谱分析。

(3分)模拟信号数字信号离散信号抽样信号正确答案C您的答案是A回答错误展开11.维持Fs不变,为提高可以增加采样点数N。

(3分)频率周期频率分辨率数字分辨率正确答案C您的答案是C回答正确展开12.离散序列x(n)只在n为时有意义。

(3分)自然数整数实数复数正确答案B您的答案是B回答正确展开采用按时间抽取的基-2 FFT算法计算N=8点DFT,需要计算______次复数加法(3分)8165664正确答案C您的答案是D回答错误展开14.所谓信号的谱分析,就是计算信号的。

数字信号处理的实际应用

数字信号处理的实际应用

数字信号处理的实际应用数字信号处理(DSP)是通过对数字信号进行采样、量化和编码等处理,基于数字计算技术对信号进行处理和分析的过程。

随着计算机技术的快速发展,数字信号处理在多个领域得到了广泛的应用。

本文将介绍数字信号处理在音频、图像、通信和生物医学领域的具体应用。

一、音频领域1. 数字音频处理数字信号处理在音频领域有着广泛的应用,其中之一就是数字音频处理。

数字音频处理使用数字信号处理技术对音频信号进行去噪、降噪、音频合成和音频压缩等处理。

通过数字滤波器和频谱分析等算法,可以实现音频信号的高质量处理和改变。

2. 数字音频合成数字音频合成是利用数字信号处理技术生成具有各种音色和音效的音频信号。

通过对已有的音频资料进行特征提取和分析,再利用信号合成算法生成新的音频信号。

数字音频合成在音乐创作、语音合成和声音效果设计等方面有广泛的应用。

二、图像领域1. 数字图像处理数字信号处理在图像领域的应用主要体现在数字图像处理。

通过对图像信号进行数字滤波、图像增强、图像压缩和图像识别等处理,可以实现图像的降噪、清晰化、特征提取和图像数据的压缩等功能。

数字图像处理在医学影像、人脸识别和图像检索等领域有重要的应用。

2. 视频编解码数字信号处理技术在视频编解码中发挥了重要的作用。

通过对视频信号的采样和压缩,可以有效地减小视频文件的大小并实现高质量的视频传输。

常见的视频编解码算法有H.264和HEVC等,它们广泛用于数字电视、视频会议和视频监控等领域。

三、通信领域1. 数字信号调制与解调数字信号处理在通信领域的应用主要体现在数字信号的调制和解调。

通过对数字信号进行调制,可以将其转换为适合在传输介质上发送的模拟信号;而解调是将接收到的模拟信号转换回数字信号进行处理和分析。

数字信号调制与解调在移动通信、无线电广播和卫星通信等领域有着广泛应用。

2. 信道编码与解码数字信号处理还在通信领域中应用于信道编码和解码。

通过对数字信号进行编码,可以增加冗余信息以提高信道传输的可靠性和容错性。

数字信号处理作业 (3)

数字信号处理作业 (3)

数字信号处理作业(一) 通信11-2班 陈亚环 2011102040172.4、图p2.1画出了用相加、相乘和延迟这三个基本运算生成的4个运算的示意图。

对每个运算将y[n]表示成x[n]的函数。

(a)、]2[]1[]1[]2[][][21120-----+-+=n y d n y d n x p n x p n x p n y (b )、]2[]1[]1[]2[][][21120-----+-+=n y d n y d n x p n x p n x p n y(c )、]6[]0[]5[)](0[]4[)](0[]3[)](0[]2[)](0[]1[)](0[][]0[][232122132122232112232211231211221113211312232223212122211313121122132312211223112211131211131211232221131211-+-+++-++++++-+++++++-++++++-+++=n x h n x h n x h n x h n x h n x h n x h n y ββββββββββββββββββββββββββββββββββββββββββββββββββββββ(d )、]4[]2[]3[])3[]2[(])5[]1[](1[])6[][](0[][-+-++-+-+-+=n x h n x h h n x n x h n x n x h n y2.39、求下列周期序列的基本周期: (a)、n j a e n x π25.0~][=(b)、)3.06.0cos(][~ππ+=n n x b (c)、)()Re(][5/8/~ππjn n j c e In e n x +=解:(a)、由r N πω20=,当0ω=0.25π时,N=8,r=1 (b )、由r N πω20=,当0ω=0.6π时,N=10,r=3(c )、)125.0c o s ()Re (8/n e n j ππ=,由r N πω20=,当0ω=0.125π时,1N =161r =1)2.0sin()(5/n j e In jn ππ=,由r N πω20=,当0ω=0.2π时,2N =10,2r =1 LCM(1N ,2N )=80Mtlab2.6t=0:0.001:0.85; g1=cos(6*pi*t); g2=cos(14*pi*t); g3=cos(26*pi*t);plot(t/0.85,g1,'-',t/0.85,g2,'--',t/0.85,g3,':'); xlabel('时间');ylabel('幅值'); hold on ; n=0:1:8;gs=cos(0.6*pi*n); plot(n/8.5,gs,'o'); hold off结果如图所示:3.16、求每个序列的DTFT. (a )、1]-u[n a ][ n 1=n x(b )、u [n]a ][ n 2n n x =(c )、1]u[n a ][ n 3+=n x解:已知u[n]a ][ n =n x 的FT 为ωωj -ae -11][ =j e X1]-u[n a a 1]-u[n a ][ 1-n n 1⨯==n x1]u[n a11]u[n a ][ 1n 3+⨯=+=+n a n x )ae -a(1e ][ae -1ae ][FT ][x j -jw3j --jw11ωωωωae X e X n j j -==∴同理为的可以得到由时移定理和线性定理为的到由频域微分定理可以得FT ][ 2n x 2j --jwj -2)ae -1ae )ae -11(][ωωωω(==d d j e X j 3.18、求如下有限长序列的DTFT.(a )、⎩⎨⎧≤≤-=其他,,01][ 1N n N n y(b )、⎩⎨⎧≤≤=其他,,001][y 2N n n解:()()[]jw jw N jw jwNN jw jwNjwnNN j e e e eeeen e -+--------=+++==∑11][y ][Y 12111 ω()[]jwjw N jw jwNjwjjwnNNj ee e eeeen e -+---------=+++==∑111][y ][Y 1222 ωω3.61、连续时间信号)(t x a 由频率为300Hz 、500Hz 、1.2kHz 、2.15kHz 和3.5kHz 的正弦信号的线性组合组成。

数字信号处理课后作业

数字信号处理课后作业

数字信号处理课后作业1、作业分两次完成,第一次就是第三周结束,第5周星期五之前交电子工程学院604房间。

内容就是第一章到第四章的作业2、第二次第8周星期五之前叫过来,内容就是第六章到第七章的作业(必须抄题)3、最后一次作业必须附上50个数字信号处理相关的英文缩写,英文全称和中文意思。

4、附上十款以上去年到今年,大的dsp 厂家生产的dsp 型号及对应的网站地址。

(做了的加平时成绩1分)数字信号处理作业第一章3. 判断下面的序列是否是周期的,若是周期的,确定其周期。

(2))81()(π-=n j en x5. 设系统分别由下面的差分方程描述,x(n)与y(n)分别表示系统输入和输出,判断系统是否是线性非时变的。

)、(为整常数300) ,n x(n-n y(n)=)s i n ()()(8n n x n y ω=)、(6. 给定下述系统的差分方程,判断系统是否是因果稳定系统,并说明理由。

∑-=-=1)(1)(1N k k n x Nn y )、(∑+-==0)()()3(n n n n k k x n y 、)()(5n x en y =)、(7. 设线性时不变系统的单位脉冲响应h(n)和输入序列x(n)如题7图所示,要求画出输出y(n)到波形。

11. 设系统由下面差分方程描述:)1(21)()1(21)(-++-=n x n x n y n y设系统是因果的,利用地推法求系统的单位取样响应。

12. 有一连续信号2/,20),2cos()(π??π==+=Hz f ft t x a 式中,(1)求出)(t x a 的周期;(2)用采样间隔T=0.02s 对)(t x a 进行采样,试写出采样信号)(^t x a 的表达式;(3)画出对应)(^t x a 的时域离散信号(序列)x(n)的波形,并求出x(n)的周期。

第二章 4. 设==其它,01,0,1)(n n x将x(n)以4为周期进行周期延拓,形成周期序列)(~n x ,画出x(n)和)(~n x 的波形,求出)(~n x 的离散傅立叶级数)(~k X 和傅立叶变换。

数字信号处理作业_答案

数字信号处理作业_答案

数字信号处理作业DFT 习题1. 如果)(~n x 是一个周期为N 的周期序列,那么它也是周期为N 2的周期序列。

把)(~n x 看作周期为N 的周期序列,令)(~1k X 表示)(~n x 的离散傅里叶级数之系数,再把)(~n x 看作周期为N 2的周期序列,再令)(~2k X 表示)(~n x 的离散傅里叶级数之系数。

当然,)(~1k X 是周期性的,周期为N ,而)(~2k X 也是周期性的,周期为N 2。

试利用)(~1k X 确定)(~2k X 。

(76-4)2. 研究两个周期序列)(~n x 和)(~n y 。

)(~n x 具有周期N ,而)(~n y 具有周期M 。

序列)(~n w 定义为)()()(~~~n y n x n w +=。

a. 证明)(~n w 是周期性的,周期为MN 。

b. 由于)(~n x 的周期为N ,其离散傅里叶级数之系数)(~k X 的周期也是N 。

类似地,由于)(~n y 的周期为M ,其离散傅里叶级数之系数)(~k Y 的周期也是M 。

)(~n w 的离散傅里叶级数之系数)(~k W 的周期为MN 。

试利用)(~k X 和)(~k Y 求)(~k W 。

(76-5)3. 计算下列各有限长度序列DFT (假设长度为N ):a. )()(n n x δ= b .N n n n n x <<-=000)()(δc .10)(-≤≤=N n an x n(78-7)4. 欲作频谱分析的模拟数据以10千赫速率被取样,且计算了1024个取样的离散傅里叶变换。

试求频谱取样之间的频率间隔,并证明你的回答。

(79 -10)5. 令)(k X 表示N 点序列)(n x 的N 点离散傅里叶变换(a ) 证明如果)(n x 满足关系式:)1()(n N x n x ---=,则0)0(=X 。

(b ) 证明当N 为偶数时,如果)1()(n N x n x --=,则0)2/(=N X 。

数字信号处理大作业

数字信号处理大作业

数字信号处理上机实验学院:电子工程学院班级:021061学号: 02106013姓名:岳震震实验一:信号、系统及系统响应02106013 岳震震一,实验目的(1) 熟悉连续信号经理想采样前后的频谱变化关系,加深对时域采样定理的理解。

(2)熟悉时域离散系统的时域特性。

(3)利用卷积方法观察分析系统的时域特性。

(4)掌握序列傅里叶变换的计算机实现方法,利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。

二,实验原理与方法(1) 时域采样。

(2)LTI系统的输入输出关系。

三,实验内容及步骤(1)认真复习采样理论、离散信号与系统、线性卷积、序列的傅里叶变换及性质等有关内容,阅读本实验原理与方法。

(2)编制实验用主程序及相应子程序。

①信号产生子程序,用于产生实验中要用到的下列信号序列:a .Xa(t)=Ae-at sin(Ω0t)U(t)b.单位脉冲序列:xb(n)=δ(n)c.矩形序列:xc(n)=RN(n),N=10②系统单位脉冲响应序列产生子程序。

本实验要用到两种FIR系统。

a .ha(n)=R10(n);b. hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3)③有限长序列线性卷积子程序用于完成两个给定长度的序列的卷积。

可以直接调用MATLAB语言中的卷积函数conv。

conv用于两个有限长度序列的卷积,它假定两个序列都从n=0开始。

调用格式如下:y=conv(x,h)调通并运行实验程序,完成下述实验内容:①分析采样序列的特性。

a. 取采样频率fs=1 kHz, 即T=1 ms。

b.改变采样频率,fs=300Hz,观察|X(ejω)|的变化,并做记录(打印曲线);进一步降低采样频率,fs=200Hz,观察频谱混叠是否明显存在,说明原因,并记录(打印)这时的|X(ejω)|曲线。

②时域离散信号、系统和系统响应分析。

a.观察信号xb(n)和系统hb(n)的时域和频域特性;利用线性卷积求信号xb(n)通过系统hb(n)的响应y(n),比较所求响应y(n)和hb(n)的时域及频域特性,注意它们之间有无差别,绘图说明,并用所学理论解释所得结果。

数字信号处理实验作业

数字信号处理实验作业

实验5 抽样定理一、实验目的:1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。

2、进一步加深对时域、频域抽样定理的基本原理的理解。

3、观察信号抽样与恢复的图形,掌握采样频率的确定方法和内插公式的编程方法。

二、实验原理:1、时域抽样与信号的重建 (1)对连续信号进行采样例5-1 已知一个连续时间信号sin sin(),1Hz 3ππ=0001f(t)=(2f t)+6f t f ,取最高有限带宽频率f m =5f 0,分别显示原连续时间信号波形和F s >2f m 、F s =2f m 、F s <2f m 三情况下抽样信号的波形。

程序清单如下:%分别取Fs=fm ,Fs=2fm ,Fs=3fm 来研究问题 dt=0.1; f0=1; T0=1/f0; m=5*f0; Tm=1/fm; t=-2:dt:2;f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); subplot(4,1,1); plot(t,f);axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3;fs=i*fm;Ts=1/fs; n=-2:Ts:2;f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); subplot(4,1,i+1);stem(n,f,'filled');axis([min(n),max(n),1.1*min(f),1.1*max(f)]); end程序运行结果如图5-1所示:原连续信号和抽样信号图5-1(2)连续信号和抽样信号的频谱由理论分析可知,信号的频谱图可以很直观地反映出抽样信号能否恢复原模拟信号。

因此,我们对上述三种情况下的时域信号求幅度谱,来进一步分析和验证时域抽样定理。

例5-2编程求解例5-1中连续信号及其三种抽样频率(F s>2f m、F s=2f m、F s<2f m)下的抽样信号的幅度谱。

数字信号处理作业及答案(全)

数字信号处理作业及答案(全)

数字信号处理作业(1)1、画出离散信号的波形 (1))2(3)3(2)(1++-=n n n x δδ (2))2()(2+-=n u n x (3))5()()(3--=n u n u n x(4))()()(214n u n x n ⋅= (5))()25.0sin(3)(5n u n n x ⋅⋅=π2、设x (n )、y (n )分别为系统的输入、输出变量,根据定义确定系统是否为:(1)线性,(2)稳定,(2)因果 ① )()]([ )(2n ax n x T n y == ② b n x n x T n y +==)()]([ )(③ )0()()]([ )(00>-==n n n x n x T n y ④ ∑+-=>=)0()( )(0n n n n m n m x n y3、已知:描述系统的差分方程为 )()1(5- )(n x n y n y =- 且初始条件为: 0)1(=-y 求:系统的单位冲激响应h (n )4、已知:线性时不变系统的单位脉冲响应为 10 , )( )(<<⋅=a n u a n h n 求:该系统的单位阶跃响应。

数字信号处理作业(1)解答1、画出离散信号的波形 (1))2(3)3(2)(1++-=n n n x δδ (2))2()(2+-=n u n x (3))5()()(3--=n u n u n x(4))()()(214n u n x n ⋅= (5))()25.0sin(3)(5n u n n x ⋅⋅=π2、设x (n )、y (n )分别为系统的输入、输出变量,根据定义确定系统是否为:(1)线性,(2)稳定,(3)因果因果:输出只取决于当前和之前的输入。

线性移不变系统的因果的充要条件:h (n )=0 , n < 0稳定系统:有界输入产生有界输出。

线性移不变系统稳定的充要条件:∞<=∑∞-∞=P n h m )(① )()]([ )(2n ax n x T n y ==(非线性,稳定,因果) ② b n x n x T n y +==)()]([ )((非线性,稳定,因果) ③ )0( )()]([ )(00>-==n n n x n x T n y (线性,稳定,因果) ④ )0( )( )(0>=∑+-=nm x n y n n n n m (线性,稳定,非因果)注意:非线性系统的稳定、因果只能按定义判断,不能按线性、移不变系统的h (n )特点判断。

数字信号处理作业3(附答案)

数字信号处理作业3(附答案)

6
7
8
a. b. c. 1 2 3 1 1 3 3 3 | | | 2 3 2 3 2 1 1 1 | 9 1 1 3 1 3 | | 26/25 7 18 7 18 6 24 6 24 ROC: -3<Re(s)<1
With the ROC in b., we have 1 2 3
-1
-0.5
0 t
0.5
1
1.5
2
5
1 1 2 2 5 6 1 1 2 2 3
For each non-empty ROC, X(s) corresponds to a x(t). X(s) has 3 poles: -1, -2 and -3. There are 8 possible ROCs (shown in the following table). And 4 of them are not empty, so there are 4 distinct signals with the expression X(s). -3 C C C C AC AC AC AC -2 C C AC AC C C AC AC -1 ROC C Re(s)>-1 AC -2<Re(s)<-1 C empty AC -3<Re(s)<-2 C empty AC empty C empty AC Re(s)<-3
ECE-S 631 Homework 3 solution 1 The following solutions are based on the assumption that the signals are causal.
பைடு நூலகம்
2 The following solutions are given based on the assumption that the signals are causal.

数字信号处理陈后金版作业

数字信号处理陈后金版作业

数字信号处理陈后金版作业数字信号处理作业:基于MATLAB的有噪声的语音信号分析与处理设计一、题目要求:1) 选择一个语音信号作为分析对象,或录制一段语音信号;2) 对语音信号进行采样,画出采样后语音信号的时域波形和频谱图;3) 利用MATLAB中的正弦函数产生噪声加入到语音信号中,使语音信号被污染,然后进行频谱分析;4) 设计FIR和IIR数字滤波器,并对被噪声污染的语音信号进行滤波,画出滤波前后信号的时域波形和频谱,并对滤波前后的信号进行比较,分析信号的变化;5) 进行信号频谱的搬移,进行变声实验6) 回放语音信号二、题目的实现方法1) 语音信号的采集利用PC 机上的声卡和WINDOWS 操作系统可以进行数字信号的采集。

将话筒输入计算机的语音输入插口上,启动录音机。

按下录音按钮,接着对话筒说话“语音信号处理”,说完后停止录音,屏幕左侧将显示所录声音的长度。

点击放音按钮,可以实现所录音的重现。

以文件名“sound”保存入g :\the sound文件夹中。

可以看到,文件存储器的后缀默认为. wav ,这是WINDOWS 操作系统规定的声音文件存的标准。

2) 语音信号的时频分析Matlab软件平台下,利用wavread函数对语音信号进行采样,记住采样频率和采样点数Wavread 函数调用格式y=wavread(file)%读取file所规定的wav文件,返回采样值放在向量y中。

[y,Ft,nbits]=wavread(file) %采样值放在向量y中,fs表示采样频率(Hz),nbits表示采样位数。

y=wavread(file,N)%读取钱N点的采样值放在向量y中。

y=wavread(file,[N1,N2])%读取从N1到N2点的采样值放在向量y中。

对语音信号shengyin.wav进行采样其程序如下:[y,Ft,nbits]=wavered (shengyin); %把语音信号进行加载入Matlab 仿真软件平台中Ft = 22050nbits = 16首先画出语音信号的时域波形,然后对语音信号进行频谱分析。

西电数字信号处理大作业

西电数字信号处理大作业

实验一、信号的采样clc,clear;dt=0.001;tf=6;t=0:dt:tf;xa=sqrt(t)+cos(t);T=0.5;n=0:tf/T;x=sqrt(n*T)+cos(n*T);figure(1)subplot(2,1,1)plot(t,xa),grid on ;title('original image')subplot(2,1,2)stem(n*T,x),grid on ,title('digital image')实验二、信号与系统的时域分析差分方程为)()2()1()(21n bx n y a n y a n y +----=,其中8.01-=a ,64.02=a ,866.0=b 。

系统单位脉冲响应)(n ha1=-0.8;a2=0.64;b=0.866;ys=0;xn=[1,zeros(1,49)];B=1;A=[1,a1,a2];xi=filtic(B,A,ys);yn=filter(B,A,xn,xi);n=0:length(yn)-1;subplot(1,1,1);stem(n,yn,'.')title('(a)');xlabel('n');ylabel('y(n)')输入x(n)=cos(n)T=0.1;z=cos(n*T);zn=conv(yn,z); figure(2);n1=1:99;stem(n1,zn,'.')实验三、系统的频域和Z域分析程序代码(画出dtft的幅度和频率谱)clc,clear;n=0:1:7;x=(0.9*exp(j*pi/3)).^n;w=0:pi/200:pi;X=x*exp(-j).^(n'*w);realX=real(X);imagX=imag(X);angX=angle(X);magX=abs(X);subplot(2,2,1);plot(w/pi,magX);grid xlabel('frequency in pi unit');title('magnitude part');subplot(2,2,2);plot(w/pi,realX);grid xlabel('frequency in pi unit');title('real part');subplot(2,2,3);plot(w/pi,imagX);grid xlabel('frequency in pi unit');title('imaginary part');subplot(2,2,4);plot(w/pi,angX);grid xlabel('frequency in pi unit');title('angel part');clc,clear;a=[1,-0.5,0.06];b=[1,1,0];m=0:length(b)-1;l=0:length(a)-1;w=0:pi/500:pi;num=b*exp(-j*m'*w);den=a*exp(-j*l'*w);H=num./den;magH=abs(H);angH=angle(H);H1=freqz(b,a,w);magH1=abs(H1);angH1=angle(H1);subplot(2,2,2);plot(w/pi,angH);grid;xlabel('w(frequency in pi units)');ylabel('Ïàλrad/w');subplot(2,2,1);plot(w/pi,magH);grid;xlabel('w(frequency in pi units)');ylabel('·ù¶È|H|');subplot(2,2,3);plot(w/pi,magH1);grid;xlabel('w(frequency in pi units)');ylabel('·ù¶È|H1|');subplot(2,2,4);plot(w/pi,angH);grid;xlabel('w(frequency in pi units)');ylabel('Ïàλrad/w');axis([0,1,-0.8,0]); figure(2);zplane(b,a);实验四、信号的频谱分析程序代码clc,clear;n=0:7;k=0:7;N=8;w=n*(2*pi)/8;x=(0.9*exp(j*pi/3)).^n;X1=[x zeros(1,8)];X2=[X1 zeros(1,16)];XK=x*exp(-j*k'*w);k1=0:15;n1=0:15;w1=n1*(2*pi)/16;XK1=X1*exp(-j*k1'*w1);k2=0:31;n2=0:31;w2=n2*(2*pi)/16;XK2=X2*exp(-j*k2'*w2);w3=0:pi/200:2*pi;X=x*exp(-j*n'*w3);magX=abs(X);angX=angle(X);magXK=abs(XK);angXK=angle(XK);magXK1=abs(XK1);angXK1=angle(XK1);magXK2=abs(XK2);angXK2=angle(XK2);subplot(4,2,1);plot(w3/pi,magX);xlabel('w/pi');ylabel('·ù¶È|X|');grid on;subplot(4,2,2);plot(w3/pi,angX);xlabel('w/pi');ylabel('Ïàλrad/pi'); subplot(4,2,3);stem(n,magXK);xlabel('K');ylabel('·ù¶È|XK|');subplot(4,2,4);stem(n,magXK);xlabel('K');ylabel('Ïàλrad/pi'); subplot(4,2,5);stem(n1,magXK1);xlabel('K1');ylabel('·ù¶È|XK1|'); subplot(4,2,6);stem(n1,magXK1);xlabel('K1');ylabel('Ïàλrad/pi'); subplot(4,2,7);stem(n2,magXK2);xlabel('K2');ylabel('·ù¶È|XK2|'); subplot(4,2,8);stem(n2,magXK2);xlabel('K2');ylabel('Ïàλrad/pi');实验五、IIR数字滤波器设计IIR汉宁窗低通高通低通巴特沃斯通带截止频率wp=0.2pi 通带最大衰减R=1dB阻带截止频率wp=0.35pi 阻带最大衰减R=10dBclc,clear;Wp=0.2;Ws=0.35;Rp=1;Rs=100;[N,Wc]=buttord(Wp,Ws,Rp,Rs);[Bz,Az]=butter(N,Wc)w=0:0.1:pi;[H,w1]=freqz(Bz,Az,w);;ang=angle(H);H=20*log10(abs(H))subplot(4,2,1); plot(w/pi,H) ;gridon ;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|/dB');title('µÍͨÂ˲¨Æ÷')subplot(4,2,2);plot(w/pi,ang);gridon ;xlabel('\omega/\pi');ylabel('Phase/dB')[Bz1,Az1]=butter(N,Wc,'high')w=0:0.1:pi;[H1,w2]=freqz(Bz1,Az1,w);ang1=angle(H1);H1=20*log10(abs(H1))subplot(4,2,3); plot(w/pi,H1) ;gridon ;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|/dB');title('¸ßͨÂ˲¨Æ÷')subplot(4,2,4);plot(w/pi,ang1);gridon ;xlabel('\omega/\pi');ylabel('Phase/dB')Wp1=[0.2 0.8];Ws1=[0.35 0.65];[N2,Wc1]=buttord(Wp1,Ws1,Rp,Rs);[Bz2,Az2]=butter(N2,Wc1,'stop')w=0:0.1:pi;[H2,w3]=freqz(Bz2,Az2,w);ang2=angle(H2);H2=20*log10(abs(H2))subplot(4,2,5); plot(w/pi,H2) ;gridon ;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|/dB');title('´ø×èÂ˲¨Æ÷')subplot(4,2,6);plot(w/pi,ang2);gridon ;xlabel('\omega/\pi');ylabel('Phase/dB')Wp1=[0.2 0.8];Ws1=[0.35 0.65];[N2,Wc1]=buttord(Wp1,Ws1,Rp,Rs);[Bz3,Az3]=butter(N2,Wc1)w=0:0.1:pi;[H3,w4]=freqz(Bz3,Az3,w);ang3=angle(H3);H3=20*log10(abs(H3))subplot(4,2,7); plot(w/pi,H3) ;gridon ;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|/dB');title('´øͨÂ˲¨Æ÷')subplot(4,2,8);plot(w/pi,ang3);gridon ;xlabel('\omega/\pi');ylabel('Phase/dB')切比雪夫1型通带截止频率wp=0.7pi 通带最大衰减R=1dB阻带截止频率wp=0.5pi 阻带最大衰减R=40dBclc,clear;Wp=0.7;Ws=0.5;Rp=1;Rs=40;[N,Wpo]=cheb1ord(Wp,Ws,Rp,Rs);[Bz,Az]=cheby1(N,Rp,Wpo)w=0:0.1:pi;[H,w1]=freqz(Bz,Az,w);ang=angle(H);H=20*log10(abs(H))subplot(4,2,1); plot(w/pi,H) ;gridon ;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|/dB');title('µÍͨÂ˲¨Æ÷')subplot(4,2,2);plot(w/pi,ang);gridon ;xlabel('\omega/\pi');ylabel('Phase/dB')[Bz1,Az1]=cheby1(N,Rp,Wpo,'high');w=0:0.1:pi;[H1,w2]=freqz(Bz1,Az1,w);ang1=angle(H1);H1=20*log10(abs(H1))subplot(4,2,3); plot(w/pi,H1) ;gridon ;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|/dB');title('¸ßͨÂ˲¨Æ÷')subplot(4,2,4);plot(w/pi,ang1);gridon ;xlabel('\omega/\pi');ylabel('Phase/dB')Wp1=[0.2 0.8];Ws1=[0.35 0.65];[N2,Wpo1]=cheb1ord(Wp1,Ws1,Rp,Rs);[Bz2,Az2]=cheby1(N2,Rp,Wpo1,'stop')w=0:0.1:pi;[H2,w3]=freqz(Bz2,Az2,w);ang2=angle(H2);H2=20*log10(abs(H2))subplot(4,2,5); plot(w/pi,H2) ;gridon ;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|/dB');title('´ø×èÂ˲¨Æ÷')subplot(4,2,6);plot(w/pi,ang2);gridon ;xlabel('\omega/\pi');ylabel('Phase/dB')Wp1=[0.2 0.8];Ws1=[0.35 0.65];[N2,Wpo1]=cheb1ord(Wp1,Ws1,Rp,Rs);[Bz3,Az3]=cheby1(N2,Rp,Wpo1)w=0:0.1:pi;[H3,w4]=freqz(Bz3,Az3,w);ang3=angle(H3);H3=20*log10(abs(H3))subplot(4,2,7); plot(w/pi,H3) ;gridon ;xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|/dB');title('´øͨÂ˲¨Æ÷')subplot(4,2,8);plot(w/pi,ang3);gridon ;xlabel('\omega/\pi');ylabel('Phase/dB')实验六、FIR数字滤波器设计FIR汉宁窗低通高通低通% 采用Hamming窗设计一个带阻FIR滤波器阻带:0~0.5pi,阻带最小衰减Rs=40dB;通带:0.5~pi,通带最大衰减:Rp=1dB。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

数字信号处理一.数字信号处理的综述数字信号处理的目的是对真实世界的连续模拟信号进行测量或滤波。

因此在进行数字信号处理之前需要将信号从模拟域转换到数字域,这通常通过模数转换器实现。

而数字信号处理的输出经常也要变换到模拟域,这是通过数模转换器实现的。

数字信号处理的算法需要利用计算机或专用处理设备如数字信号处理器(DSP)和专用集成电路(ASIC)等。

数字信号处理技术及设备具有灵活、精确、抗干扰强、设备尺寸小、造价低、速度快等突出优点,这些都是模拟信号处理技术与设备所无法比拟的。

数字信号处理的核心算法是离散傅里叶变换(DFT),是DFT使信号在数字域和频域都实现了离散化,从而可以用通用计算机处理离散信号。

而使数字信号处理从理论走向实用的是快速傅里叶变换(FFT),FFT的出现大大减少了DFT的运算量,使实时的数字信号处理成为可能、极大促进了该学科的发展。

二.数字信号处理的应用广义来说,数字信号处理是研究用数字方法对信号进行分析、变换、滤波、检测、调制、解调以及快速算法的一门技术学科。

但很多人认为:数字信号处理主要是研究有关数字滤波技术、离散变换快速算法和谱分析方法。

随着数字电路与系统技术以及计算机技术的发展,数字信号处理技术也相应地得到发展,其应用领域十分广泛。

⑴数字滤波器的实用型式很多,大略可分为有限冲激响应型和无限冲激响应型两类,可用硬件和软件两种方式实现。

在硬件实现方式中,它由加法器、乘法器等单元所组成,这与电阻器、电感器和电容器所构成的模拟滤波器完全不同。

数字信号处理系统很容易用数字集成电路制成,显示出体积小、稳定性高、可程控等优点。

数字滤波器也可以用软件实现。

软件实现方法是借助于通用数字计算机按滤波器的设计算法编出程序进行数字滤波计算。

⑵在频域中描述信号特性的一种分析方法,不仅可用于确定性信号,也可用于随机性信号。

所谓确定性信号可用既定的时间函数来表示,它在任何时刻的值是确定的;随机信号则不具有这样的特性,它在某一时刻的值是随机的。

因此,随机信号处理只能根据随机过程理论,利用统计方法来进行分析和处理,如经常利用均值、均方值、方差、相关函数、功率谱密度函数等统计量来描述随机过程的特征或随机信号的特性。

实际上,经常遇到的随机过程多是平稳随机过程而且是各态历经的,因而它的样本函数集平均可以根据某一个样本函数的时间平均来确定。

平稳随机信号本身虽仍是不确定的,但它的相关函数却是确定的。

在均值为零时,它的相关函数的傅里叶变换或Z变换恰恰可以表示为随机信号的功率谱密度函数,一般简称为功率谱。

这一特性十分重要,这样就可以利用快速变换算法进行计算和处理。

在实际中观测到的数据是有限的。

这就需要利用一些估计的方法,根据有限的实测数据估计出整个信号的功率谱。

针对不同的要求,如减小谱分析的偏差,减小对噪声的灵敏程度,提高谱分辨率等。

已提出许多不同的谱估计方法。

在线性估计方法中,有周期图法,相关法和协方差法;在非线性估计方法中,有最大似然法,最大熵法,自回归滑动平均信号模型法等。

谱分析和谱估计仍在研究和发展中。

数字信号处理的应用领域十分广泛。

就所获取信号的来源而言,有通信信号的处理,雷达信号的处理,遥感信号的处理,控制信号的处理,生物医学信号的处理,地球物理信号的处理,振动信号的处理等。

若以所处理信号的特点来讲,又可分为语音信号处理,图像信号处理,一维信号处理和多维信号处理等。

⑶数字信号处理系统,无论哪方面的应用,首先须经过信息的获取或数据的采集过程得到所需的原始信号,如果原始信号是连续信号,还须经过抽样过程使之成为离散信号,再经过模数转换得到能为数字计算机或处理器所接受的二进制数字信号。

如果所收集到的数据已是离散数据,则只须经过模数转换即可得到二进制数码。

数字信号处理器的功能是将从原始信号抽样转换得来的数字信号按照一定的要求,例如滤波的要求,加以适当的处理,即得到所需的数字输出信号。

经过数模转换先将数字输出信号转换为离散信号,再经过保持电路将离散信号连接起来成为模拟输出信号,这样的处理系统适用于各种数字信号处理的应用,只不过专用处理器或所用软件有所不同而已。

⑷语音信号处理是信号处理中的重要分支之一。

它包括的主要方面有:语音的识别,语言的理解,语音的合成,语音的增强,语音的数据压缩等。

各种应用均有其特殊问题。

语音识别是将待识别的语音信号的特征参数即时地提取出来,与已知的语音样本进行匹配,从而判定出待识别语音信号的音素属性。

关于语音识别方法,有统计模式语音识别,结构和语句模式语音识别,利用这些方法可以得到共振峰频率、音调、嗓音、噪声等重要参数,语音理解是人和计算机用自然语言对话的理论和技术基础。

语音合成的主要目的是使计算机能够讲话。

为此,首先需要研究清楚在发音时语音特征参数随时间的变化规律,然后利用适当的方法模拟发音的过程,合成为语言。

其他有关语言处理问题也各有其特点。

语音信号处理是发展智能计算机和智能机器人的基础,是制造声码器的依据。

语音信号处理是迅速发展中的一项信号处理技术。

⑸图像信号处理的应用已渗透到各个科学技术领域。

譬如,图像处理技术可用于研究粒子的运动轨迹、生物细胞的结构、地貌的状态、气象云图的分析、宇宙星体的构成等。

在图像处理的实际应用中,获得较大成果的有遥感图像处理技术、断层成像技术、计算机视觉技术和景物分析技术等。

根据图像信号处理的应用特点,处理技术大体可分为图像增强、恢复、分割、识别、编码和重建等几个方面。

这些处理技术各具特点,且正在迅速发展中。

⑹机械振动信号的分析与处理技术已应用于汽车、飞机、船只、机械设备、房屋建筑、水坝设计等方面的研究和生产中。

振动信号处理的基本原理是在测试体上加一激振力,做为输入信号。

在测量点上监测输出信号。

输出信号与输入信号之比称为由测试体所构成的系统的传递函数(或称转移函数)。

根据得到的传递函数进行所谓模态参数识别,从而计算出系统的模态刚度、模态阻尼等主要参数。

这样就建立起系统的数学模型。

进而可以做出结构的动态优化设计。

这些工作均可利用数字处理器来进行。

这种分析和处理方法一般称为模态分析。

实质上,它就是信号处理在振动工程中所采用的一种特殊方法。

⑺地球物理信号处理,为了勘探地下深处所储藏的石油和天然气以及其他矿藏,通常采用地震勘探方法来探测地层结构和岩性。

这种方法的基本原理是在一选定的地点施加人为的激震,如用爆炸方法产生一振动波向地下传播,遇到地层分界面即产生反射波,在距离振源一定远的地方放置一列感受器,接收到达地面的反射波。

从反射波的延迟时间和强度来判断地层的深度和结构。

感受器所接收到的地震记录是比较复杂的,需要处理才能进行地质解释。

处理的方法很多,有反褶积法,同态滤波法等,这是一个尚在努力研究的问题。

数字信号处理在其他方面还有多种用途,如雷达信号处理、地学信号处理等,它们虽各有其特殊要求,但所利用的基本技术大致相同。

在这些方面,数字信号处理技术起着主要的作用。

三.数字信号处理的发展趋势数字信号处理(DigitalSignalProcessing,简称DSP)是一门涉及许多学科而又广泛应用于许多领域的新兴学科,他的发展趋势如下:⑴技术发展趋势:①数字信号处理器的内核结构进一步改善,多通道结构和单指令多重数据(SIMD)、特大指令字组(VLIM)将在新的高性能处理器中将占主导地位,如Analog Devices的ADSP-2116x。

②DSP 和微处理器的融合,微处理器是低成本的,主要执行智能定向控制任务的通用处理器能很好执行智能控制任务,但是数字信号处理功能很差。

而DSP的功能正好与之相反。

在许多应用中均需要同时具有智能控制和数字信号处理两种功能,如数字蜂窝电话就需要监测和声音处理功能。

因此,把DSP和微处理器结合起来,用单一芯片的处理器实现这两种功能,将加速个人通信机、智能电话、无线网络产品的开发,同时简化设计,减小PCB 体积,降低功耗和整个系统的成本。

例如,有多个处理器的Motorola公司的DSP5665x,有协处理器功能的Massan公司FILU-200,把MCU功能扩展成DSP和MCU功能的TI公司的TMS320C27xx以及Hitachi公司的SH-DSP,都是DSP和MCU融合在一起的产品。

互联网和多媒体的应用需要将进一步加速这一融合过程③DSP 和高档CPU的融合,大多数高档GPP如Pentium 和PowerPC都是SIMD指令组的超标量结构,速度很快。

LSI Logic 公司的LSI401Z采用高档CPU的分支预示和动态缓冲技术,结构规范,利于编程,不用担心指令排队,使得性能大幅度提高。

Intel公司涉足数字信号处理器领域将会加速这种融合。

④DSP 和SOC的融合:SOC(System-On-Chip)是指把一个系统集成在一块芯片上。

这个系统包括DSP 和系统接口软件等。

比如Virata公司购买了LSI Logic公司的ZSP400处理器内核使用许可证,将其与系统软件如USB、10BASET、以太网、UART、GPIO、HDLC 等一起集成在芯片上,应用在xDSL上,得到了很好的经济效益。

因此,SOC芯片近几年销售很好,由1998年的1.6亿片猛增至1999年的3.45亿片。

1999年,约39%的SOC产品应用于通讯系统。

今后几年,SOC将以每年31%的平均速度增长,到2004年将达到13亿片。

毋庸置疑,SOC将成为市场中越来越耀眼的明星。

⑤DSP 和FPGA的融合,FPGA是现场编程门阵列器件。

它和DSP集成在一块芯片上,可实现宽带信号处理,大大提高信号处理速度。

据报道,Xilinx 公司的Virtex-II FPGA对快速傅立叶变换(FFT)的处理可提高30倍以上。

它的芯片中有自由的FPGA可供编程。

Xilinx 公司开发出一种称作Turbo卷积编译码器的高性能内核。

设计者可以在FPGA中集成一个或多个Turbo内核,它支持多路大数据流,以满足第三代(3G)WCDMA无线基站和手机的需要,同时大大节省开发时间,使功能的增加或性能的改善非常容易。

因此在无线通信、多媒体等领域将有广泛应用。

⑵市场发展趋势:很明然,在可预见的一段时间内,无线应用仍将是可编程DSP市场的驱动引擎。

不过,嵌入式DSP市场是一个更大的市场,我们将在后续的市场报告中予以讨论。

DSP技术(无论何种形式)是你能够访问窄带、宽带或是无线互联网的唯一手段,它还是新兴的分组(IP)电话市场的关键。

没有DSP就没有对互联网的访问,没有多媒体,也没有无线通信。

因此,尽管遭遇了短期的市场挫折,DSP仍将是整个半导体工业的技术驱动力。

而国内发展DSP的厂商并不多,而主要的应用产品是DVD与无线电话等,因此国内DSP的产值并不高;而在产品应用上,目前重要的DSP应用产品,如行动电话、调制解调器、HDD等个人计算机与通讯领域应用产品,都是采用国际大厂的DSP solution,因此国内厂商尚无插足的余地。

相关文档
最新文档