成都电子科技大学数字信号处理习题课

合集下载

数字信号处理学习指导与题解 4 khdaw

数字信号处理学习指导与题解 4 khdaw

ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
课后答案网,用心为你服务!
最全最多的课后习题参考答案,尽在课后答案网()!
Khdaw团队一直秉承用心为大家服务的宗旨,以关注学生的学习生活为出发点,
爱校园() 课后答案网() 淘答案()
ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
ww

.k h
后 答 案 网
da w. co m
ww

电子科技大学DSP上机题目

电子科技大学DSP上机题目
1 已知 3 阶椭圆 IIR 数字低通滤波器的性能指标为: 通带截止频率 0.4π , 通带波纹为 0.6dB, 最小阻带衰减为 32dB。设计一个 6 阶全通滤波器对其通带的群延时进行均衡。绘制低通滤 波器和级联滤波器的群延时。 %Progranm 1 % Group-delay equalization of an IIR filter. % [n,d] = ellip(3,0.6,32,0.4); [GdH,w] = grpdelay(n,d,512); plot(w/pi,GdH); grid xlabel('\omega/\pi'); ylabel('Group delay, samples'); title('Original Filter'); F = 0:0.001:0.3; g = grpdelay(n,d,F,2); % Equalize the passband Gd = max(g)-g; % Design the allpass delay equalizer [num,den,tau] = iirgrpdelay(6, F, [0 0.3], Gd); [GdA,w] = grpdelay(num,den,512); figure(2); plot(w/pi,GdH+GdA); grid xlabel('\omega/\pi');ylabel('Group delay, samples'); title('Group Delay Equalized Filter'); 2 设计巴特沃兹模拟低通滤波器,其滤波器的阶数和 3-dB 截止频率由键盘输入,程序能根 据输入的参数,绘制滤波器的增益响应。 % Program 2 % Program to Design Butterworth Lowpass Filter % % Type in the filter order and passband edge frequency N = input('Type in filter order = '); Wn = input('3-dB cutoff angular frequency = '); % Determine the transfer function

课件:数字信号处理,第五章,电子科技大学

课件:数字信号处理,第五章,电子科技大学

Example N-point DFT of x[n] cos nRN [n] 6 (N=12) is wanted.
Solution
X [k ] x[n]W
n 0
N 1
kn N
1 [e n 0 2
11
j
2 n 12
e
j
2 n 12
]e
j
2 kn 12
6, k 1, 11 { 0, orthers

Solution 2 Using (B), we get
1 j 12 n j 12 kn 1 j 12 n j 12 kn X [k ] [ e e e e ] 2 n 0 2
~ 11 2 2 2 2
1 11 j 12 ( k 1) n 1 11 j 12 ( k 11) n e e 2 n 0 2 n 0
• Making use of the identity
( k r ) n W N n 0 N 1

N , for k r N , l an integer 0, otherwise
we get
N / 2, G[k ] N / 2, 0, for k r for k N r otherwise
k
Example x[n] cos n 6 has DFS, the coefficients are wanted. Solution 1
1 x[n] e 2
~ 2 j n 12
~
1 e 2
2 j n 12

1 e 2
j
2 n 12
1 e 2
j
2 (11) n 12

电子科大 数字信号处理实验2_FFT的实现

电子科大 数字信号处理实验2_FFT的实现

电 子 科 技 大 学实 验 报 告学生姓名:Shrimp 学 号: 指导教师:一、实验室名称:数字信号处理实验室 二、实验项目名称:FFT 的实现 三、实验原理:一.FFT 算法思想:1.DFT 的定义:对于有限长离散数字信号{x[n]},0 ≤ n ≤ N-1,其离散谱{x[k]}可以由离散付氏变换(DFT )求得。

DFT 的定义为:21[][]N jnk Nn X k x n eπ--==∑,k=0,1,…N-1通常令2jNN eW π-=,称为旋转因子。

2.直接计算DFT 的问题及FFT 的基本思想:由DFT 的定义可以看出,在x[n]为复数序列的情况下,完全直接运算N 点DFT 需要(N-1)2次复数乘法和N (N-1)次加法。

因此,对于一些相当大的N 值(如1024)来说,直接计算它的DFT 所作的计算量是很大的。

FFT 的基本思想在于,将原有的N 点序列分成两个较短的序列,这些序列的DFT 可以很简单的组合起来得到原序列的DFT 。

例如,若N 为偶数,将原有的N 点序列分成两个(N/2)点序列,那么计算N 点DFT 将只需要约[(N/2)2 ·2]=N 2/2次复数乘法。

即比直接计算少作一半乘法。

因子(N/2)2表示直接计算(N/2)点DFT 所需要的乘法次数,而乘数2代表必须完成两个DFT 。

上述处理方法可以反复使用,即(N/2)点的DFT 计算也可以化成两个(N/4)点的DFT (假定N/2为偶数),从而又少作一半的乘法。

这样一级一级的划分下去一直到最后就划分成两点的FFT 运算的情况。

3.基2按时间抽取(DIT )的FFT 算法思想:设序列长度为2L N =,L 为整数(如果序列长度不满足此条件,通过在后面补零让其满足)。

将长度为2L N =的序列[](0,1,...,1)x n n N =-,先按n 的奇偶分成两组:12[2][][21][]x r x r x r x r =+=,r=0,1,…,N/2-1DFT 化为:1/21/212(21)0/21/21221200/21/211/22/2[]{[]}[][2][21][][][][]N N N nk rk r kNNNn r r N N rk k rk NNN r r N N rk k rk N NN r r X k DFT x n x n Wx r Wx r W x r W Wx r W x r WWx r W ---+===--==--=====++=+=+∑∑∑∑∑∑∑上式中利用了旋转因子的可约性,即:2/2rkrkNN W W =。

《数字信号处理》自测题3答案——电子科技大学

《数字信号处理》自测题3答案——电子科技大学

一.选择与填空:1.若一系统为因果系统,则该系统的单位脉冲响应h(n)应满足0,0)(<≡n n h ,该系统的系统函数应满足 其收敛域为∞≤<-z R x ;若该系统为稳定系统,则该系统的单位脉冲响应h(n)应满足∞<∑∞-∞=n n h )(1=。

2.对下列序列,试求其零极点,Z 变换和收敛域。

a) )()(10n R n x = b) )(5sin )(n u n n x ⎪⎭⎫⎝⎛=π解:a) )1(111)(910110--=--=--z z z z z z X 其收敛域为∞≤<z 0,其零极点分布为: 零点:9个一阶零点9,,1 102 ==k ez k j k π;极点:一个9阶极点;0=zb) ⎪⎪⎭⎫ ⎝⎛-==-∞=--∞=∞=-∑∑∑n n n j n n n j n nz e z e j z nz X 0505021)5sin()(πππ ))(()2.0sin()2.0cos(21)2.0sin(111121552111515πππππππj j j j e z e z z zz z z e z e j ---------=+-=⎪⎪⎪⎭⎫⎝⎛---= 其收敛域为∞≤<z 1,其零极点分布为: 零点:1个一阶零点0=z ;极点:2个一阶极点;2.01πj e z =;2.02πj e z -=3.若对话音信号进行离散时间处理,已知话音信号为实信号,其频率范围为300~3400Hz ,则抽样频率为 >6800Hz 。

4.对连续时间非周期信号的谱分析工具是 傅立叶变换 ,对连续时间周期信号的谱分析工具是 傅立叶级数 ,对周期序列的谱分析工具是 离散傅立叶级数 ,对非周期序列的谱分析工具是 离散傅立叶变换 。

5.若某系统的群延时为常数,则该系统具有的相位特性为 线性相位 。

6.脉冲响应不变法存在频谱混叠的现象,是因为 z s →映射为多对一映射 ;双线性变换法不存在频谱混叠的现象,是因为 z s →映射为一对一映射 ;因此,对于同样指标设计的滤波器,双线性变换法的设计结果 优于 (优于/劣于)脉冲响应不变法的设计结果。

电子科技大学 数字信号处理 第二次编程作业

电子科技大学 数字信号处理  第二次编程作业

第二次编程作业加载信号(EEGdata.txt);分析信号的幅度谱,确定低通FIR数字滤波器的指标;分别利用各种窗函数(Rectangular, Hamming, Kaiser)设计此FIR滤波器;对信号加随机噪声,并用设计的FIR滤波器对含噪声信号进行滤波。

要求:1)给出程序,画出滤波器幅度谱及其增益图(dB);分析设计的滤波器是否达到指标。

2)利用设计的三种滤波器对加载的信号分别进行滤波,对比分析滤波结果。

程序:data = load('EEGdata.txt');n=1:length(data);xn=0.1*rand(length(data),1)-0.05;x=xn+data;[h0,w0]=freqz(data);[h1,w1]=freqz(x);figure(1)plot(w0/pi,abs(h0),w1/pi,abs(h1),'r');legend('Ô-dataÐźÅƵÆ×ͼ','¼ÓÔëÐźÅƵÆ×ͼ');wp = 0.15*pi;ws = 0.2*pi;DB = ws-wp;N = ceil(2*3.11*pi/DB+1);wc = (wp+ws)/2;wc = wc/pi;win = rectwin(N+1);b = fir1(N, wc, 'low', win);figure(2)freqz(b,1);wp = 0.15*pi;ws = 0.2*pi;DB = ws-wp;N = ceil(2*3.11*pi/DB+1);wc = (wp+ws)/2;wc = wc/pi;win = hamming(N+1);b1 = fir1(N, wc, 'low', win);figure(3)freqz(b1,1);fpts = [0.15, 0.2];mag = [1, 0];dev = [0.01, 0.01];[N, Wn, beta, ftype] = kaiserord(fpts, mag, dev)Kwin = kaiser(N+1, beta);b2 = fir1(N, Wn, Kwin);[h,omega] = freqz(b2,1,512);figure(4)freqz(b2,1);y=filter(b, 1, x);figure(5)subplot(2,1,1);plot(n,data,n,x,'g',n,y,'r');title('¾ØÐδ°Â˲¨·ù¶ÈÆ×Ч¹û¶Ô±Èͼ');lege nd('Ô-ÐźÅdata','¼ÓÔëÐźÅx','Â˲¨ºóÐźÅy');[H,W]=freqz(y);subplot(2,1,2);plot(w0/pi,abs(h0),W/pi,abs(H),'r');title('¾ØÐδ°Â˲¨ÆµÆ×Ч¹û¶Ô±Èͼ') ;legend('Ô-ÐźÅdata','Â˲¨ºóÐźÅy');y1=filter(b1, 1, x);figure(6)subplot(2,1,1);plot(n,data,n,x,'g',n,y1,'r');title('Hamming´°Â˲¨·ù¶ÈÆ×Ч¹û¶Ô±Èͼ'); legend('Ô-ÐźÅdata','¼ÓÔëÐźÅx','Â˲¨ºóÐźÅy1');[H,W]=freqz(y1);subplot(2,1,2);plot(w0/pi,abs(h0),W/pi,abs(H),'r');title('Hamming´°Â˲¨ÆµÆ×Ч¹û¶Ô±Èͼ');legend('Ô-ÐźÅdata','Â˲¨ºóÐźÅy1');y2=filtfilt(b2, 1, x); figure(7) subplot(2,1,1);plot(n,data,n,x,'g',n,y2,'r');title('Kaiser´°Â˲¨·ù¶ÈÆ×Ч¹û¶Ô±Èͼ');l egend('Ô-ÐźÅdata','¼ÓÔëÐźÅx','Â˲¨ºóÐźÅy2'); [H,W]=freqz(y2); subplot(2,1,2);plot(w0/pi,abs(h0),W/pi,abs(H),'r');title('Kaiser´°Â˲¨ÆµÆ×Ч¹û¶Ô±Èͼ');legend('Ô-ÐźÅdata','Â˲¨ºóÐźÅy2');figure(8)plot(n,data,n,y,'c',n,y1,'r',n,y2,'g');title('ÈýÖÖ´°º¯ÊýÂ˲¨Ð§¹û¶Ô±Èͼ');legend('Ô-dataÐźÅ','¾ØÐδ°Â˲¨','Hamming´°Â˲¨','Kaiser´°Â˲¨');0.10.20.30.40.50.60.70.80.9102468101214161800.10.20.30.40.50.60.70.80.91-3000-2000-1000Normalized Frequency (⨯π rad/sample)P h a s e (d e g r e e s )00.10.20.30.40.50.60.70.80.91-100-5050Normalized Frequency (⨯π rad/sample)M a g n i t u d e (d B )00.10.20.30.40.50.60.70.80.91-3000-2000-1000Normalized Frequency (⨯π rad/sample)P h a s e (d e g r e e s )0.10.20.30.40.50.60.70.80.91-150-100-50050Normalized Frequency (⨯π rad/sample)M a g n i t u d e (d B )0.10.20.30.40.50.60.70.80.91-2000-1500-1000-500Normalized Frequency (⨯π rad/sample)P h a s e (d e g r e e s )0.10.20.30.40.50.60.70.80.91-150-100-50050Normalized Frequency (⨯π rad/sample)M a g n i t u d e (d B )-0.2-0.100.10.2矩形窗滤波幅度谱效果对比图0.10.20.30.40.50.60.70.80.91051015矩形窗滤波频谱效果对比图-0.2-0.100.10.2Hamming 窗滤波幅度谱效果对比图0.10.20.30.40.50.60.70.80.91051015Hamming 窗滤波频谱效果对比图-0.2-0.100.10.2Kaiser 窗滤波幅度谱效果对比图0.10.20.30.40.50.60.70.80.9105101520Kaiser 窗滤波频谱效果对比图结果分析:滤波效果皆不尽如人意,原因是加载的噪声信号是随机信号,各频率皆有,滤波器通带部分的噪声无法滤掉。

电子科技大学《数字信号处理》自测题

电子科技大学《数字信号处理》自测题

1自测题一及答案一.选择与填空:1.数字信号是______(3)___。

(1)时间上和幅度上均连续的理想信号 (2)时间上离散而幅度上连续的信号(3)时间上离散、幅度上量化的非理想信号 (4)时间上连续而幅度上量化的非理想信号2.对于时间上离散的周期信号,可以___(4)___________。

(1)用傅立叶变换(FT )求其频谱 (2)用离散傅立叶变换(DFT )求其频谱 (3)可展成傅立叶级数(FS ) (4)可展成离散傅立叶级数(DFS )3.已知x(n)为16点序列,其DFT 应为___(5)___,IDFT 应为____(6)___。

(1) ∑==15016)(161)(n nK W n x K X (2) nKn W n x K X 1616)()(∑==(3) ∑==1516)()(n nKWn x K X (4) ∑==1515)()(n nKW n x K X(5) ∑==1516)()(K nK WK X n x (6) ∑=-=15016)(161)(K nKW K X n x (7) ∑==15016)(161)(K nK W K X n x (8) ∑=-=15015)(161)(K nKW K X n x 4.数字信号处理中,时域上表征系统的瞬时输出经常用_____(2)_____。

(1)系统的单位脉冲响应 h(n) (2)离散卷积和(3)系统的传递函数H(Z) (4)系统的频响H(e j ω)5.如果对数字滤波器的相频特性不作要求,仅从滤波效率来看,IIR 数字滤波器的效率比FIR 数字滤波 器的效率要___高___,即在相同的通带、阻带滤波性能指标下,__IIR_____数字滤波器结构要简单些, 其滤波器的阶数N 要__低_____。

6.已知某线性相位FIR 数字滤波器的其中一个零点40πjeZ =,根据零点分布规律,其系统函数:H (Z )=)1)(1(44ππjjeeA ---。

电子科技大学20春《数字信号处理》在线作业1.doc

电子科技大学20春《数字信号处理》在线作业1.doc

1.下面关于IIR滤波器设计说法正确的是()。

A.双线性变换法的优点是数字频率和模拟频率成线性关系B.冲激响应不变法无频率混叠现象C.冲激响应不变法不适合设计高通滤波器D.双线性变换法只适合设计低通、带通滤波器【参考答案】: C2.题目及选项如下图所示A.AB.BC.CD.D【参考答案】: A3.下列关于窗函数设计法的说法中错误的是()。

A.窗函数的长度增加,则主瓣宽度减小,旁瓣宽度减小B.窗函数的旁瓣相对幅度取决于窗函数的形状,与窗函数的长度无关C.为减小旁瓣相对幅度而改变窗函数的形状,通常主瓣的宽度会增加D.对于长度固定的窗,只要选择合适的窗函数就可以使主瓣宽度足够窄,旁瓣幅度足够小【参考答案】: B4.一个线性移不变系统稳定的充分必要条件是其系统函数的收敛域包括()。

A.原点B.单位圆C.虚轴D.实轴【参考答案】: B5.题目及选项如下图所示A.AB.BC.CD.D【参考答案】: A6.已知一连续时间信号的最高截止频率是4kHz,则要从该连续时间信号的均匀离散采样值无失真地恢复原信号,则采样频率为()。

A.2kHzB.4kHzC.6kHzD.8kHz【参考答案】: D7.题目及选项如下图所示A.AB.BC.CD.D【参考答案】: A8.题目及选项如下图所示A.AB.BC.CD.D【参考答案】: A9.对连续时间周期信号的谱分析工具是()。

A.傅里叶变换B.傅里叶级数C.离散傅里叶变换D.离散傅里叶级数【参考答案】: B10.下面说法中正确的是()。

A.连续时间非周期信号的频谱为连续周期函数B.连续时间周期信号的频谱为连续周期函数C.离散时间非周期信号的频谱为离散周期函数D.离散时间周期信号的频谱为离散周期函数【参考答案】: D11.题目及选项如下图所示A.AB.BC.CD.D【参考答案】: A12.要处理一个连续时间信号,对其进行采样的频率为6kHz,要无失真的恢复该连续信号,则该连续信号的最高频率可能是()。

电子科技大学 数字信号处理 上机实验四

电子科技大学 数字信号处理 上机实验四

电子科技大学生命科学与技术学院标准实验报告(实验)课程名称数字信号处理2012-2013-第2学期电子科技大学教务处制表电子科技大学实验报告学生姓名:张华博学号:2011091010004 指导教师:李永杰实验地点:清水河校区实验时间:2013年月日一、实验室名称:科二 504 机房二、实验名称:FIR滤波器的特性及应用三、实验学时:2学时四、实验原理:五、实验目的:(详细填写)1. 熟悉相位信息的特性2. 测试FIR滤波器的滤波特性六、实验内容:(详细填写)1. 相位信息的作用2. 测试线性相位FIR滤波器的滤波特性七、实验器材(设备、元器件):八、实验步骤:九、实验数据及结果分析:(详细填写)(包括程序、图、结果等)1、n=0:99;m=2/100:2/100:2;y=sin(0.1*pi*n)+2*cos(0.2*pi*n);XXX = fft(y);plot(m,abs(XXX));title('信号频谱图');figure(2);plot(m,unwrap(angle(XXX))); title('信号相位谱');x1=cos(0.85*pi*n);x2=0.2*n;Pha1=unwrap(angle(XXX))+x1;Pha2=unwrap(angle(XXX))+x2;figure(3)plot(m,Pha1);title('加随机相位的相位谱');figure(4)plot(m,Pha2);title('加线性相位的相位谱');cc1 = abs(XXX).*exp(i*Pha1);dd1 = ifft(cc1, 'symmetric');cc2 = abs(XXX).*exp(i*Pha2);dd2 = ifft(cc2, 'symmetric');figure(5)plot(n,y,'b',n,dd1,'r',n,dd2,': b');title('fft反变换与原信号对比图'); legend('原信号','加随机相位的信号','加线性相位的信号');0.20.40.60.811.21.41.61.820102030405060708090100信号频谱图00.20.40.60.81 1.2 1.4 1.6 1.82510152025信号相位谱00.20.40.60.81 1.2 1.4 1.6 1.82510152025加随机相位后的相位谱00.20.40.60.81 1.2 1.4 1.6 1.8251015202530354045加线性相位后的相位谱102030405060708090100-3-2-1123fft 反变换与原信号对比图2、h=[-1 2 -3 6 -3 2 -1]; %×Ô¼º¹¹Ôìh(n)ÐòÁÐ a=1; %µ±ÎªFIRÂ˲¨Æ÷ʱ£¬È¡a=1; zplane(h,a); M = 100;[H,w] = freqz(h,a,M); figure(2); plot(w,abs(H)); xlabel('\omega'); ylabel('Amplitude'); figure(3);plot(w,unwrap(angle(H))); xlabel('\omega');ylabel('Phase');h1=[-1 2 -3 -3 2 -1]; %×Ô¼º¹¹Ôìh(n)ÐòÁÐ a1=1; %µ±ÎªFIRÂ˲¨Æ÷ʱ£¬È¡a=1; figure(4); zplane(h1,a1); M = 100;[H1,w1] = freqz(h1,a1,M); figure(5); plot(w1,abs(H1));xlabel('\omega');ylabel('Amplitude');figure(6);plot(w1,unwrap(angle(H1)));xlabel('\omega');ylabel('Phase');h2=[-1 2 -3 0 3 -2 1]; %×Ô¼º¹¹Ôìh(n)ÐòÁÐa2=1; %µ±ÎªFIRÂ˲¨Æ÷ʱ£¬È¡a=1;figure(7);zplane(h2,a2);M = 100;[H2,w2] = freqz(h2,a2,M);figure(8);plot(w2,abs(H2));xlabel('\omega');ylabel('Amplitude');figure(9);plot(w2,unwrap(angle(H2)));xlabel('\omega');ylabel('Phase');h3=[-1 2 -3 3 -2 1]; %×Ô¼º¹¹Ôìh(n)ÐòÁÐa3=1; %µ±ÎªFIRÂ˲¨Æ÷ʱ£¬È¡a=1;figure(10);zplane(h3,a3);M = 100;[H3,w3] = freqz(h3,a3,M);figure(11);plot(w3,abs(H3));xlabel('\omega');ylabel('Amplitude');figure(12);plot(w3,unwrap(angle(H3)));xlabel('\omega');ylabel('Phase');-1.5-1-0.50.511.522.5-1.5-1-0.50.511.5Real PartI m a g i n a r y P a r t00.51 1.52 2.53 3.524681012141618A m p l i t u d e0.511.522.533.5-10-9-8-7-6-5-4-3-2-1P h a s e-2-1.5-1-0.500.511.52-1.5-1-0.500.511.5Real PartI m a g i n a r y P a r t00.51 1.52 2.53 3.5123456789ωA m p l i t u d e00.51 1.52 2.53 3.5-5-4-3-2-101234ωP h a s e-2-1.5-1-0.500.511.52-1.5-1-0.50.511.5Real PartI m a g i n a r y P a r t0.511.522.533.5012345678910A m p l i t u d e00.51 1.52 2.53 3.5-12-10-8-6-4-2P h a s e-1-0.500.51-1-0.8-0.6-0.4-0.200.20.40.60.81Real PartI m a g i n a r y P a r t00.51 1.52 2.53 3.524681012ωA m p l i t u d e00.51 1.52 2.53 3.5-5-4-3-2-11ωP h a s eh=[-1 2 -3 6 -3 2 -1]; a=1;coef = sum(abs(h))/sum(a);h = h/coef; zplane(h,a); M = 50;[H,w] = freqz(h,a,M); figure(2); plot(w/pi,abs(H)); figure(3);plot(w/pi,unwrap(angle(H))); n=0:49;xs=sin(0.95*pi*n); xn=sin(0.5*pi*n); x=xs+xn;y = filter(h, a, x); [h2,w2]=freqz(y) figure(4)plot(n,xs,n,y,'r');legend('Ô-ÐźÅ','Â˲¨ºó'); figure(5)plot(w1/pi,abs(h1),w2/pi,abs(h2),'r');title('ƵÆ×ͼ¶Ô±È');legend('Ô-ÐźÅ','Â˲¨ºó');-1.5-1-0.50.511.522.5-1.5-1-0.50.511.5Real PartI m a g i n a r y P a r t00.10.20.30.40.50.60.70.80.910.10.20.30.40.50.60.70.80.910.10.20.30.40.50.60.70.80.91-10-9-8-7-6-5-4-3-2-15101520253035404550-1.5-1-0.50.511.50.10.20.30.40.50.60.70.80.910510152025频谱图对比十、实验结论:(详细填写)1、给一信号相位谱加一线性相位后,做fft 反变换所得图像为原图平移一定相位;若加一随机相位,做反变换后,不仅相位发生平移,而且幅度也发生变化。

电子科大数字信号处理(dsp)matlab作业

电子科大数字信号处理(dsp)matlab作业

M2.22.30(b)n=[-10:10];>> x=sin(0.8*pi*n+0.8*pi); >> stem(n,x)title('2.30(b)的序列')>> xlabel('n')>> ylabel('x')2.30(c)n=[-10:10];>> x=cos(pi*n/5);>> stem(n,x)>> title('the real part')>> xlabel('n')>> ylabel('re(x[n])')>> title('the real part of 2.30(c)')x=sin(pi*n/10);>> stem(n,x)>> title('the imaginary part of 2.30(c)')>> xlabel('n')>> ylabel('im(x[n])')2.30(d)n=[-20:20];>> x=3*cos(1.3*pi*n)-4*sin(0.5*pi*n+0.5*pi);>> stem(n,x)>> title('2.30(d)的序列')>> xlabel('n')>> ylabel('x[n]')2.30(e)n=[-20:20];>> x=5*cos(1.5*pi*n+0.75*pi)+4*cos(0.6*pi*n)-sin(0.5*pi*n); >> stem(n,x)>> title('2.30(e)的序列')>> xlabel('n')>> ylabel('x[n]')M2.7编写的程序如下% Program M2.7% Computation the discreet time system of exercise 2.40alpha=input('alpha=');y0=1;y1=0.5*(y0+(alpha/y0));while abs(y1-y0)>0.0001y2=0.5*(y1+(alpha/y1));y0=y1;y1=y2;enddisp('y[n]收敛于根号alpha=');y2运行结果如下:alpha=9y[n]收敛于根号alpha=y2 =3.0000alpha=30y[n]收敛于根号alpha=y2 =5.4772经两组数据对比,知已用仿真证明,n趋于无穷时,Y[n]收敛于根号下alpha。

电子科技大学 DSP期末考试试卷

电子科技大学 DSP期末考试试卷

学院 姓名 学号 任课老师 选课号……………密……………封……………线……………以……………内……………答……………题……………无……………效…………电子科技大学二零零七至二零零 八学年第 1 学期期 期末 考试《数字信号处理》 课程考试题 A 卷 ( 120 分钟) 考试形式: 开卷 考试日期 2008年1月 24 日课程成绩构成:平时 20分, 期中 20分, 实验 20 分, 期末 40分 Test of Digital Signal Processing(Give the answer of the following problems in English as possible as you can )Note: Problem 1 to 5 must be answered. You can choose to solute either one or both of problem6 and 7, if both, the higher score will be accounted .1. Which one of the following comments is not correct ? ( )(a) The DFT can be thought of the sampled version of the DTFT . (b) The FFT is the efficient computation algorithm for the DFT .(c) A linear-phase digital filter can be realized by either FIR or IIR systems .(d) If H(z) is the transfer function of an LPF filter , then H(-z) can be thought of an HPF filter in general .2. Two systems have impulse responses as follows :]1[][][1−+=n n n h δδ ,]1[][][2−−=n n n h δδ, respectively. (1)Determine the DTFTs of h1[n] and h2[n].(2)Determine and sketch their magnitude responses . (3)Which type of linear-phase system are they belong to?(4)If input x[n] = 2cos(πn) ,determine the output sequences of these two systems respectively.3. Consider the system transfer function:23()3148z H z z z +=−+(1) Sketch the pole-zero plot of H(z);(2) Discuss the type of inverse Z-transform (left-sided, right-sided, and two-sided sequences) associated with each of all possible ROCs of H(z). (3) Determine the expression of causal impulse response h[n].4. Figure 4(a ) is identical to Figure 4.2 of our textbook, which illustrates the discrete-time processing of a continuous-time signal. The ideal sampler and ideal interpolator work at the same sampling frequency T F . Assume the discrete-time processing block represents an idealbandpass digital filter with frequency response)(ωj e H illustrated in Figure 4(b ).Let )1502cos(5.1)1252sin()1002cos()(t t t t x a ×+×+×=πππ. (1) If Hz F T 100=, determine []x n , []y n and ()a y t . (2) If Hz F T150=, determine []x n , []y n and ()a y t .Figure 4(a))Figure 4(b)ω544−5−5. Let X[k] be an 8-point DFT of a length-8 real sequence x[n]. The 5 samples of X[k] are given by X[0]=0, X[1]=0, X[2]=4, X[3]=0, X[4]=0.(1) Determine and sketch X[k].(2) Determine the value of ∑=702][n n x .(3) If ]2[][8>+<=k X k G , sketch G[k] and determine the expression of g[n]which is the 8-point IDFT of G[k].ing windowed Fourier series method, design a causal bandpass FIR digital filter withfollowing specifications:center frequency 5kHzpass band edges 3.5 kHz and 6.5 kHzstop band edges 1.5 kHz and 8.5 kHzstop band attenuation 50 dB,passband ripple can be ignored.Sampling frequency F T=20kHz.Please determine the expression of the impulse response of the FIR digital filter designed.7. The normalized transfer function of order-2 lowpass Butterworth analog filter is as follows: 121)(2++=s s s H an .By means of bilinear transform, design an order-2 low-pass IIR digital filter with -3dB cutoff frequency at 1000Hz and sampling frequency at 4000Hz. You are expected to determine the expression of transfer function in closed form and sketch an implement block-diagram in any cannonic form.。

电子科技大学数字信号处理复习提纲ppt课件

电子科技大学数字信号处理复习提纲ppt课件

N 1
X (k) DFT[x(n)]N x(n)WNkn n0
k=0, 1, …, N-1
x(n)
IDFT[ X (k)]N
1 N
N 1
X (k)WNkn
k0
2) 隐含周期性
k=0, 1, …, N-1
N 1
N 1
X (k mN ) x(n)WN(kmN)n x(n)WNkn X (k)
y(n) x(m)h(n m) x(n) * h(n) m
• (2) x(n)=x(n)*δ(n) ;x(n-n0)=x(n)*δ(n-n0) • (3)
Xˆ n ( j )
1 T
Xa
k
( j
jks )
2
第2章 时域离散信号和系统的频域分析
• 学习要点 • (1)Z变换的正变换和逆变换定义, 以及收敛域与序列特性之间的
Y (e j ) X (e j )H (e j )
• (5)频域卷积定理
若y(n)=x(n)h(n), 则
Y (e j ) 1 H (e j ) X (e j ) 2π
5
第2章 时域离散信号和系统的频域分析
• (6)共轭对称序列和共轭反对称序列
xe (n)
1 [x(n) 2
x (n)]
xo (n)
3
第2章 时域离散信号和系统的频域分析
• 重要公式
• (1)傅里叶变换的正变换和逆变换的公式
X (e j )
x(n)e jn
n
x(n) 1 π X (e j )e jnd
2 -π
• (2)周期序列的离散傅里叶级数变换对
X~
(k
)
DFS[~x (n)]

电子科大《数字信号处理(DSP)》第1章 数字信号处理概述

电子科大《数字信号处理(DSP)》第1章 数字信号处理概述

数字信号的表达
0.85: 0.81: 0.65: 0.49: 0.39: …
1101 1100 1010 0111 0110 …
1001 1111 0110 1101 0011
数字信号的特点
信号采用抽象数字序列表
0.85: 0.81: 0.65: 0.49: 0.39: … 1101 1100 1010 0111 0110 … 1001 1111 0110 1101 0011
达,与物理量没有直接关系, 与物理量没有直接关系, 在传输、保存和处理过程中, 在传输、保存和处理过程中, 信号精度不受环境因素影响, 信号精度不受环境因素影响, 抗干扰性强。 抗干扰性强。
数字信号的特点
信号采用数字序列表达后, 信号采用数字序列表达后,对模拟信号难以 进行的很多处理能够方便地实现,例如: 进行的很多处理能够方便地实现,例如: 大规模长时间的信号存储、 大规模长时间的信号存储、对信号的乘法调 制和各种编码调制、信号的时间顺序处理、 制和各种编码调制、信号的时间顺序处理、信 号的时间压缩/扩张、复杂标准信号的产生… 号的时间压缩/扩张、复杂标准信号的产生…
模拟信号的精度问题
线路上传输的模拟电压会随距离衰减, 线路上传输的模拟电压会随距离衰减,电 路系统内电容中存储的电压会随时间衰减, 路系统内电容中存储的电压会随时间衰减, 这些信号同样会受环境因素影响。 这些信号同样会受环境因素影响。
数字信号的表达
对模拟电压进行等间隔测量, 对模拟电压进行等间隔测量,将各测量值采 用有限精度的数值表达, 用有限精度的数值表达,体现为顺序排布的数 字序列; 字序列;
模拟信号的精度问题
例:模拟放大器与模拟加法器
电路中的电阻值可能由于工艺因素或环境 温度影响而与设计值产生误差! 温度影响而= − AX R

《数字信号处理》第三版答案(非常详细完整)

《数字信号处理》第三版答案(非常详细完整)

答案很详细,考试前或者平时作业的时候可以好好研究,祝各位考试成功!!电子科技大学微电子与固体电子学陈钢教授著数字信号处理课后答案1.2 教材第一章习题解答1. 用单位脉冲序列()n δ及其加权和表示题1图所示的序列。

解:()(4)2(2)(1)2()(1)2(2)4(3)0.5(4)2(6)x n n n n n n n n n n δδδδδδδδδ=+++-+++-+-+-+-+-2. 给定信号:25,41()6,040,n n x n n +-≤≤-⎧⎪=≤≤⎨⎪⎩其它(1)画出()x n 序列的波形,标上各序列的值;(2)试用延迟单位脉冲序列及其加权和表示()x n 序列; (3)令1()2(2)x n x n =-,试画出1()x n 波形; (4)令2()2(2)x n x n =+,试画出2()x n 波形; (5)令3()2(2)x n x n =-,试画出3()x n 波形。

解:(1)x(n)的波形如题2解图(一)所示。

(2)()3(4)(3)(2)3(1)6()6(1)6(2)6(3)6(4)x n n n n n n n n n n δδδδδδδδδ=-+-+++++++-+-+-+-(3)1()x n 的波形是x(n)的波形右移2位,在乘以2,画出图形如题2解图(二)所示。

(4)2()x n 的波形是x(n)的波形左移2位,在乘以2,画出图形如题2解图(三)所示。

(5)画3()x n 时,先画x(-n)的波形,然后再右移2位,3()x n 波形如 5. 设系统分别用下面的差分方程描述,()x n 与()y n 分别表示系统输入和输出,判断系统是否是线性非时变的。

(1)()()2(1)3(2)y n x n x n x n =+-+-; (3)0()()y n x n n =-,0n 为整常数; (5)2()()y n x n =; (7)0()()nm y n x m ==∑。

电子科技大学《数字信号处理》20春期末考试.doc

电子科技大学《数字信号处理》20春期末考试.doc

1.计算序列x(n)的256点DFT,需要的复数乘法次数为 ( ) 。

A.256B.256×256C.256×255D.128×8【参考答案】: B2.下列关于FIR滤波器的说法中正确的是()A.FIR滤波器容易设计成线性相位特性B.FIR滤波器的脉冲响应长度是无限的C.FIR滤波器的脉冲响应长度是确定的D.对于相同的幅频特性要求,用FIR滤波器实现要比用IIR滤波器实现阶数低【参考答案】: C3.实序列傅里叶变换的实部和虚部分别为( )。

A.偶函数和奇函数B.奇函数和偶函数C.奇函数和奇函数D.偶函数和偶函数【参考答案】: A4.序列x(n)=u(n)的能量为()。

A.1B.9C.11D.∞【参考答案】: D5.要处理一个连续时间信号,对其进行采样的频率为3kHz,要不失真的恢复该连续信号,则该连续信号的最高频率可能是为( )。

A.6kHzB.1、5kHzC.3kHzD.2kHz【参考答案】: B6.以下对有限长单位冲激响应(FIR)滤波器特点的论述中错误的是( )。

A.FIR滤波器容易设计成线性相位特性B.FIR滤波器的单位冲激抽样响应h(n)在有限个n值处不为零 C.系统函数H(z)的极点都在z=0处 D.实现结构只能是非递归结构【参考答案】: D7.下面关于IIR滤波器设计说法正确的是( )。

A.双线性变换法的优点是数字频率和模拟频率成线性关系B.冲激响应不变法无频率混叠现象 C.冲激响应不变法不适合设计高通滤波器 D.双线性变换法只适合设计低通、带通滤波器【参考答案】: C8.下列关于用冲激响应不变法设计IIR滤波器的说法中错误的是( )、A.数字频率与模拟频率之间呈线性关系B.能将稳定的模拟滤波器映射为一个稳定的数字滤波器 C.使用的变换是s平面到z平面的多值映射 D.可以用于设计低通、高通和带阻等各类滤波器【参考答案】: D9.在窗函数法的FIR滤波器设计法中,滤波器过渡带宽度近似等于窗函数幅度谱的主瓣宽度的()倍。

电子科技大学 18秋《数字信号处理》在线作业1满分答案

电子科技大学 18秋《数字信号处理》在线作业1满分答案

18秋《数字信号处理》在线作业1
题目及选项如下图所示
A.A
B.B
C.C
D.D
正确答案:A
下列关于因果稳定系统说法错误的是()。

A.极点可以在单位圆外
B.系统函数的Z变换收敛域包括单位圆
C.因果稳定系统的单位冲激响应是因果序列
D.因果稳定系统的系统函数的Z变换收敛域包括
正确答案:A
题目及选项如下图所示
A.A
B.B
C.C
D.D
正确答案:C
一个线性移不变系统稳定的充分必要条件是其系统函数的收敛域包括()。

A.原点
B.单位圆
C.虚轴
D.实轴
正确答案:B
题目及选项如下图所示
A.A
B.B
C.C
D.D
正确答案:A。

(完整word版)数字信号处理(程佩青)课后习题解答(4)

(完整word版)数字信号处理(程佩青)课后习题解答(4)

第四章 快速傅立叶变换运算需要多少时间。

计算需要多少时间,用,问直拉点的,用它来计算每次复加速度为平均每次复乘需如果一台通用计算机的FFT DFT[x(n)]512s 5 s 50.1μμ 解: 解: ⑴ 直接计算:复乘所需时间: 复加所需时间:⑵用FFT 计算:复乘所需时间:复加所需时间:s N T N 01152.0 512log 105 log 105 2251262261=⨯⨯⨯=⨯⨯=--sT T T sN N T 013824.0 002304.0 512log 512105.0 log 105.0 2126262=+=∴=⨯⨯⨯=⨯⨯⨯=--sT T T sN N T 441536.1 130816.0 )1512(512105.0 )1(105.0 21662=+=∴=-⨯⨯⨯=-⨯⨯⨯=--s N T 31072.1 512105 105 26261=⨯⨯=⨯⨯=--运算一次完成。

点试用一个为了提高运算效率值求今需要从值的点实序列是两个已知IFFT N n y n x k Y k X DFT n y n x N k Y k X ,,)(),()(),(,)(),()(),(.2值的过程。

)(),(完成计算点)可用一次()()(综上所述,构造序列)()()()(可得:)()()(再根据都是实序列,)(),(由原题可知:)()()()(()()(性质:又根据可得序列点作对取序列依据题意解 ]Im[ ]Re[ ][][ ][ ).()( )()()( )()();()( ::n y n x IFFT N k jY k X k Z n z n y n z n x n jy n x n z n y n x n jy n x k Y jIDFT k X IDFT k jY k X IDFT DFT n z IFFT N k Z k jY k X k Z k Y n y k X n x +===+=+=+=++=⇔⇔。

电子科技大学智慧树知到“电子信息工程”《数字信号处理》网课测试题答案4

电子科技大学智慧树知到“电子信息工程”《数字信号处理》网课测试题答案4

电子科技大学智慧树知到“电子信息工程”《数字信号处理》网课测试题答案(图片大小可自由调整)第1卷一.综合考核(共15题)1.下面说法中正确的是()。

A.连续非周期信号的频谱为非周期离散函数B.连续周期信号的频谱为非周期离散函数C.离散非周期信号的频谱为非周期离散函数D.离散周期信号的频谱为非周期离散函数2.对5点有限长序列[1 3 0 5 2]进行向左2点圆周移位后得到序列()。

A.[1 3 0 5 2]B.[5 2 1 3 0]C.[0 5 2 1 3]D.[0 0 1 3 0]3.对于序列的傅立叶变换而言,其信号的特点是()。

A.时域连续非周期,频域连续非周期B.时域离散周期,频域连续非周期C.时域离散非周期,频域连续非周期D.时域离散非周期,频域连续周期4.下面关于IIR滤波器设计说法正确的是()。

A.双线性变换法的优点是数字频率和模拟频率成线性关系B.冲激响应不变法无频率混叠现象C.冲激响应不变法不适合设计高通滤波器D.双线性变换法只适合设计低通、带通滤波器5.数字信号的特征是()。

A.时间连续、幅值量化B.时间离散、幅值量化C.时间离散、幅值连续D.时间连续、幅值连续6.下列系统(其中y(n)为输出序列,x(n)为输入序列)中属于线性系统的是()。

A.y(n)=y(n-1)x(n)B.y(n)=x(n)/x(n+1)C.y(n)=x(n)+1D.y(n)=x(n)-x(n-1)7.在窗函数法的FIR滤波器设计法中,滤波器过渡带宽度近似等于窗函数幅度谱的主瓣宽度的()倍。

A.1/8B.1/4C.1D.48.下列关于窗函数设计法的说法中错误的是()。

A.窗函数的截取长度增加,则主瓣宽度减小,旁瓣宽度减小B.窗函数的旁瓣相对幅度取决于窗函数的形状,与窗函数的截取长度无关C.为减小旁瓣相对幅度而改变窗函数的形状,通常主瓣的宽度会增加D.窗函数法不能用于设计FIR高通滤波器9.下列关于FIR滤波器的说法中正确的是()A.FIR滤波器容易设计成线性相位特性B.FIR滤波器的脉冲响应长度是无限的C.FIR滤波器的脉冲响应长度是确定的D.对于相同的幅频特性要求,用FIR滤波器实现要比用IIR滤波器实现阶数低10.已知某序列Z变换的收敛域为|Z|>3,则该序列为()。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
m =−∞ ∞ ∞
∑∫
m =−∞ ∞ −∞
xa (τ )h(t − τ )δ (τ − mT )dτ
π
sin[ (t − mT )] T 内插函数:h(t − mT ) = π (t − mT ) T
π
信号的抽样值xa (mT )经内插函数得到连续信号xa (t )
3.2 采样的恢复
零阶保持器恢复模拟信号 零阶保持器恢复模拟信号
第一章 信号和信号处理
1.信号 1.信号
信号
信息的载体, 信息的载体,携带信息的物理量 携带信息的物理量/ 的物理量/传递信息的函数 传递信息的函数, 的函数, 变量可以是时间 变量可以是时间、 可以是时间、空间、 空间、温度、 温度、距离、 距离、电压等
信号分类
按信号载体的物理特性: 按信号载体的物理特性:点、光声、 光声、磁、机械、 机械、热等 按信号中自变量的数目: 按信号中自变量的数目:
数字信号处理
模拟信号处理: 模拟信号处理:用电阻、 用电阻、电容、 电容、电感等模拟器件处理 数字信号处理: 数字信号处理:实质是数值运算( 实质是数值运算(加、减、乘、除、延时) 延时)
数字信号处理实现方法
软件、 软件、硬件、 硬件、软硬结合
数字信号处理特点
灵活性好、 灵活性好、高精度、 高精度、高稳定性、 高稳定性、便于开发升级、 便于开发升级、功能强、 功能强、便 于大规模集成
π
解:1)由
f 0 = 50 Hz ,得
采样频率应: f s > 2 f 0 = 100 Hz 采样间隔应为:T < 1 / f s = 0.01s
2)选f s = 200 Hz 则采样间隔为:T = 1/ f s = 0.005s
xa (nT ) = sin(2π f 0 nT + π / 8) = sin(2π f 0 n / f s + π / 8) 50 1 = sin(2π n + π / 8) = sin( π n + π / 8) 200 2 ∞ ˆa (t ) = ∑ xa (nT )δ (t − nT ) ∴x
各种数字信息系统
Enhanced Gaming Webpad Wireless Devices: 802.11, Bluetooth, Others Biometrics Military and Government Cellular, Secure Connectivity Telematics Industry-Specific PDAs Digital Media Processing
2.信号处理 2.信号处理
信号处理 从信号中提取尽可能多的有用信息 增强信号的有用分量 估计信号的特征参数 识别信号的特征 抑制或消除不需要的有害分量 需要对信号: 需要对信号:分析和变换, 分析和变换,扩展、 扩展、压缩、 压缩、滤波参数 估计、 估计、特征识别等加工
2.信号处理 2.信号处理
n =−∞

1 π n = ∑ sin( π n + )δ (t − ) n =−∞ 2 8 200
x(n) = xa (t )
t = nT
1 π = sin( π n + ) 2 8
拓展问题 采样序列为矩形脉冲, 采样序列为矩形脉冲,问采样信号的频谱与原始 模拟信号的频谱关系如何? 模拟信号的频谱关系如何? 能否无失真恢复? 能否无失真恢复? 条件? 条件?
−∞
1 h(t ) = 2π


H ( j Ω ) e j Ωt d Ω
= ∫ [ ∑ xa (τ )δ (τ − mT )]h(t − τ )dτ
−∞ ∞


−∞
=
sin[ (t − mT )] T = ∑ xa (mT )h(t − mT ) = ∑ xa (mT ) π m =−∞ m =−∞ (t − mT ) T
课程说明
作业或实验要求
3 次上机实验, 次上机实验,每次4小时, 小时,大作业 作业: 作业:每章内容结束后的下节课交作业
考核方式
理论考试: 理论考试:闭卷考试 2学时 计分方式: 作业、 实验、 计分方式 :理论考试占70%,作业 、实验 、平时成绩占30%
教材及参考书
教材: 教材: 史林、 史林 、赵树杰 “数字信号处理” 数字信号处理”,科学出版社 参考书: 参考书 : 1.Oppenheim A.V., Schafer R.W. “Digital Signal Processing”, Prentice Hall Inc. 丁玉美“ 数字信号处理” 2.丁玉美 “数字信号处理 ”,西安电子科技大学出版社 陈怀深等译“ 实现” 3. 陈怀深等译 “数字信号处理及其MATLAB实现 ”,电子工业出版社
k =−∞ ∧ ∧ ∞
∑ 2π a δ ( Ω − kΩ )
k s
fs = 1 T
3.1 采样定理
2π Pδ ( jΩ ) = T ∧ 1 X a ( jΩ ) = 2π 1 = 2π
采样信号的频谱 是原信号频谱的 周期延拓
k =−∞
∑ δ (Ω − kΩ )
s ∞

X a ( jΩ ) ∗ Pδ ( jΩ ) 2π ⋅ T
一维信号( 一维信号(如语言信号) 如语言信号) 二维信号( 二维信号(如图像信号) 如图像信号) 多维信号
按信号是否为确定函数: 按信号是否为确定函数:
确定信号, 确定信号,随机信号
1.信号 1.信号
信号分类
按信号中自变量和幅度的取值特点
时间连续、 时间连续、幅度连续 模拟信号 时间离散、 时间离散、幅度连续 离散时间信号 时间连续、 时间连续、幅度离散 时间离散、 时间离散、幅度离散 数字信号
作业
1.2, 1.2,1.4, 1.4,1.5, 1.5,1.6 自己仿真1.9 自己仿真1.9, 1.9,1.10
Medical Devices
3.模拟信号的数字处理方法 3.模拟信号的数字处理方法
时间连续信号 时间采样 幅度量化 数字信号
离散时间信号
问题1 问题1:采样前后信号频谱的变化? 采样前后信号频谱的变化? 问题2 问题2:什么条件下, 什么条件下,可以从采样信号无失真恢复 出原信号? 出原信号 ?


n =−∞
∑ x ( nT ) δ ( t − nT )
a

3.1 采样定理
采样信号的频谱分析 采样信号的频谱分析
X a ( jΩ ) = FT [ xa ( t )] X a ( jΩ ) = FT [ x a ( t )] Pδ ( jΩ ) = FT [ Pδ ( t )] = Ω s = 2π T
X a(jΩ )
(a )

-Ω c
0
Ωc P (jΩ ) δ
采样定理
(b )

-Ω s
0
Ωs
^ X a(jΩ )

(c )
-Ω s Ω s 2 0
Ωs 2 Ω s
采样频率 > 信号最高频谱的 2倍
^ X a(js
Ωs
0
Ω cΩ s
2
3.2 采样的恢复
利用低通滤波器还原满足奈奎斯特采样 利用低通滤波器还原满足奈奎斯特采样定理 奈奎斯特采样定理的 定理的 采样信号。 采样信号。 ˆ ( j Ω) X Ya ( jΩ) a


−∞
X a ( jθ ) ∑ δ ( Ω − k Ω s − θ )dθ
k =−∞
1 ∞ ∞ = ∑ ∫ X a ( jθ ) δ ( Ω − k Ω s − θ )dθ T k =−∞ −∞ 1 ∞ = ∑ X a ( jΩ − jk Ω s ) T k =−∞
3.1 采样定理
采样信号的频谱是模拟信号频谱 采样信号的频谱是模拟信号频谱以采样频率 的频谱是模拟信号频谱以采样频率为 以采样频率为 周期进行周期延拓而成 频谱幅度是原信号频谱幅度的1/T 频谱幅度是原信号频谱幅度的1/T倍 1/T倍 若信号的最高频率 Ωs Ωc > , 2 Ωs 为折叠频率 2 则延拓分量产生频 谱混叠
3.1 采样定理
采样示意
理想采样
当 τ→0
ˆa (t ) = x a (t )δT (t ) x
3.1 采样定理
采样信号表示 采样信号表示
Pδ ( t ) = ∑ δ ( t − nT ) n =−∞ ⇒ ∞ ∧ x a ( t ) = xa ( t ) ⋅ Pδ ( t ) = ∑ xa ( t ) δ ( t − nT ) n =−∞ xa (t ) =
x (n ) 解码 x a(n T) 零阶保持 ′ (t ) xa 平滑滤波 x a(t )
卷积
H ( jΩ ) = ∫ e − jΩt dt
T
sin ( Ω T 2 ) − j ΩT =T e ΩT 2
0
2
例:模拟信号xa (t ) = sin(2π f0t + ),其中f0 = 50 Hz 8 1 )求采样频率应为多少?采样间隔应为多少? 2)若选采样频率f s = 200 Hz,采样间隔为多少? ˆa (t )的表达式; 写出采样信号x ˆa (t )的时域离散信号x(n)的波形. 3)画出对应x
H[jΩ]
理想低通滤波器:
T H ( jΩ) = 0 Ωs Ω< 2 Ωs Ω≥ 2
H(jΩ) T
-Ωs/2
0
Ωs/2

ˆ ( jΩ) ⋅ H ( jΩ) = X ( jΩ) Ya ( jΩ) = X a a
讨论
ˆa (t ) → xa (t ) x
Ωs π sin( t ) sin( t ) T Ω s j Ωt 2 = T = e d Ω = Ωs π 2π ∫−Ωs t t 2 T ∞ ˆa (τ )h(t − τ )dτ 输出:ya (t ) = xa (t ) = ∫ x
相关文档
最新文档