FIR 数字滤波器设计函数 1． fir1 功能：基于窗函数的FIR 数字滤波器设计——标准频率响应。 格式：b=fir1(n,Wn) b=fir1(n,Wn,'ftype') b=fir1(n,Wn,Window) b=fir1(n,Wn,'ftype',Window) 说明：fir1函数以经典方法实现加窗线性相位FIR 滤波器设计，它可设计出标准的低通、带通、高通和带阻滤波器。 b=fir1(n,Wn)可得到n 阶低通FIR 滤波器，滤波器系数包含在b 中，这可表示成： n z n b z b b z b --++???++=)1()2()1()(1 这是一个截止频率为Wn 的Hamming(汉明)加窗线性相位滤波器，0≤Wn ≤1，Wn=1相应于0.5fs 。 当Wn=[W1 W2]时，fir1函数可得到带通滤波器，其通带为W1<ω< W2。 b=fir1(n,Wn,'ftype')可设计高通和带阻滤波器，由ftype 决定： ·当ftype=high 时，设计高通FIR 滤波器； ·当ftype=stop 时，设计带阻FIR 滤波器。 在设计高通和带阻滤波器时，fir1函数总是使用阶为偶数的结构，因此当输入的阶次为奇数时，fir1函数会自动加1。这是因为对奇数阶的滤波器，其在Nyquist 频率处的频率响应为零，因此不适合于构成高通和带阻滤波器。 b=fir1(n,Wn,Window)则利用列矢量Window 中指定的窗函数进行滤波器设计，Window 长度为n+1。如果不指定Window 参数，则fir1函数采用Hamming 窗。 Blackman 布莱克曼窗 Boxcar 矩形窗 Hamming 海明窗 Hann 汉宁窗 Kaiser 凯瑟窗 Triang 三角窗 b=fir1(n,Wn,'ftype',Window)可利用ftype 和Window 参数，设计各种加窗的滤波器。 由fir1函数设计的FIR 滤波器的群延迟为n/2。 例如： n=32;wn=1/4;window=boxcar(n+1) b=fir1(n,wn,window)
一．填空题（本题总分12分，每空1分） 1．累加器A分为三个部分，分别为；；。 1．AG,AH,AL 2．TMS320VC5402型DSP的内部采用条位的多总线结构。 2．8,16 3．TMS320VC5402型DSP采用总线结构对程序存储器和数据存储器进行控制。3．哈佛 4．TMS329VC5402型DSP有个辅助工作寄存器。 4．8个 5．DSP处理器TMS320VC5402中DARAM的容量是字。 5．16K字 6．TI公司的DSP处理器TMS320VC5402PGE100有___________个定时器。 6．2 7．在链接器命令文件中，PAGE １通常指________存储空间。 7．数据 8．C54x的中断系统的中断源分为____ ___中断和____ ____中断。 8．硬件、软件 1．TI公司DSP处理器的软件开发环境是__________________。 1．答：CCS(Code Composer Studio) 2．DSP处理器TMS320VC5402外部有___________根地址线。 2．答：20根 3．直接寻址中从页指针的位置可以偏移寻址个单元。 3．答：128 4．在链接器命令文件中，PAGE 0通常指________存储空间。 4．答：程序 5．C54x系列DSP处理器中，实现时钟频率倍频或分频的部件是_____________。 5．答：锁相环PLL 6．TMS320C54x系列DSP处理器上电复位后，程序从指定存储地址________单元开始工作。6．答：FF80h 7．TMS320C54x系列DSP处理器有_____个通用I/O引脚，分别是_________。 7．答：2个，BIO和XF 8．DSP处理器按数据格式分为两类，分别是_______ __；_____ ___。 8．答：定点DSP和浮点DSP 9．TMS329VC5402型DSP的ST1寄存器中，INTM位的功能是。 9．答：开放/关闭所有可屏蔽中断 10．MS320C54X DSP主机接口HPI是________位并行口。 10．答：8 1.在C54X系列中，按流水线工作方式，分支转移指令的分为哪两种类型：_______；_______。 1.答：无延迟分支转移，延迟分支转移 3.C54x的程序中，“.bss”段主要用于_______________。 3.答：为变量保留存储空间 4.从数据总线的宽度来说，TMS320VC5402PGE100是_______位的DSP处理器。 4.答：16位 7.TMS320VC5402型DSP处理器的内核供电电压________伏。 7.答：1.8v
第一章数字信号处理概述 简答题： 1．在A/D变换之前和D/A变换之后都要让信号通过一个低通滤波器，它们分别起什么作用？ 答：在A/D变化之前为了限制信号的最高频率，使其满足当采样频率一定时，采样频率应大于等于信号最高频率2倍的条件。此滤波器亦称为“抗混叠”滤波器。 在D/A变换之后为了滤除高频延拓谱，以便把抽样保持的阶梯形输出波平滑化，故又称之为“平滑”滤波器。 判断说明题： 2．模拟信号也可以与数字信号一样在计算机上进行数字信号处理，自己要增加一道采样的工序就可以了。 （） 答：错。需要增加采样和量化两道工序。 3．一个模拟信号处理系统总可以转换成功能相同的数字系统，然后基于数字信号处理理论，对信号进行等效的数字处理。（） 答：受采样频率、有限字长效应的约束，与模拟信号处理系统完全等效的数字系统未必一定能找到。因此数字信号处理系统的分析方法是先对抽样信号及系统进行分析，再考虑幅度量化及实现过程中有限字长所造成的影响。故离散时间信号和系统理论是数字信号处
理的理论基础。 第二章 离散时间信号与系统分析基础 一、连续时间信号取样与取样定理 计算题： 1．过滤限带的模拟数据时，常采用数字滤波器，如图所示，图中T 表示采样周期（假设T 足够小，足以防止混叠效应），把从)()(t y t x 到的整个系统等效为一个模拟滤波器。 （a ） 如果kHz T rad n h 101,8)(=π截止于，求整个系统的截止频 率。 （b ） 对于kHz T 201=，重复（a ）的计算。 采样（T） () n h () n x () t x () n y D/A 理想低通T c πω=() t y 解 （a ）因为当0)(8=≥ω πωj e H rad 时，在数 — 模变换中 )(1)(1)(T j X T j X T e Y a a j ωω=Ω= 所以)(n h 得截止频率8πω=c 对应于模拟信号的角频率c Ω为 8 π = ΩT c 因此 Hz T f c c 625161 2==Ω= π
实验报告 课程名称：数字信号处理指导老师：刘英成绩：_________________实验名称： FIR数字滤波器设计与使用同组学生姓名：__________ 一、实验目的和要求 设计和应用FIR低通滤波器。掌握FIR数字滤波器的窗函数设计法，了解设计参数（窗型、窗长）的影响。 二、实验内容和步骤 编写MATLAB程序，完成以下工作。 2-1 设计两个FIR低通滤波器，截止频率 C =0.5。 （1）用矩形窗，窗长N=41。得出第一个滤波器的单位抽样响应序列h 1(n)。记下h 1 (n) 的各个抽样值，显示h 1 (n)的图形（用stem(.)）。求出该滤波器的频率响应（的N 个抽样）H 1(k)，显示|H 1 (k)|的图形（用plot(.)）。 （2）用汉明窗，窗长N=41。得出第二个滤波器的单位抽样响应序列h 2(n)。记下h 2 (n) 的各个抽样值，显示h 2(n)的图形。求出滤波器的频率响应H 2 (k)，显示|H 2 (k)|的 图形。 （3）由图形，比较h 1(n)与h 2 (n)的差异，|H 1 (k)|与|H 2 (k)|的差异。 2-2 产生长度为200点、均值为零的随机信号序列x(n)（用rand(1,200)0.5）。显示x(n)。 求出并显示其幅度谱|X(k)|，观察特征。 2-3 滤波 （1）将x(n)作为输入，经过第一个滤波器后的输出序列记为y 1(n)，其幅度谱记为|Y 1 (k)|。 显示|X(k)|与|Y 1 (k)|，讨论滤波前后信号的频谱特征。 （2）将x(n)作为输入，经过第二个滤波器后的输出序列记为y 2(n)，其幅度谱记为|Y 2 (k)|。 比较|Y 1(k)|与|Y 2 (k)|的图形，讨论不同的窗函数设计出的滤波器的滤波效果。 2-4 设计第三个FIR低通滤波器，截止频率 C =0.5。用矩形窗，窗长N=127。用它对x(n)进行滤波。显示输出信号y
实验1 常见离散信号产生与实现 一、实验目得 1、加深对常用离散信号得理解; 2、熟悉使用MATLAB在时域中产生一些基本得离散时间信号。 二、实验原理 1、单位抽样序列 在MATLAB中可以利用函数实现。 2、单位阶越序列 在MATLAB中可以利用函数实现: 3、正弦序列 在MATLAB中实现过程如下:
4、复指数序列 在MATLAB中实现过程如下: 5、指数序列 在MATLAB中实现过程如下: 三、预习要求 1、预先阅读实验讲义(MATLAB基础介绍); 2、讨论正弦序列、复指数序列得性质。 A.绘出信号,当、时、、 时得信号实部与虚部图;当时呢？此时信号周期为多少？程序如下: titlez1=-1/12+j*pi/6; titlez2=1/12+j*pi/6; z3=1/12; z4=2+j*pi/6; z5=j*pi/6;
n=0:20; x1=exp(titlez1*n); x2=exp(titlez2*n); x3=exp(z3*n); x4=exp(z4*n); x5=exp(z5*n); subplot(5,2,1); stem(n,real(x1)); xlabel('n');ylabel('real(x1)'); title('z1=-1/12+j*pi/6时') subplot(5,2,2); stem(n,imag(x1)); xlabel('n');ylabel('imag(x1)'); title('z1=-1/12+j*pi/6时') subplot(5,2,3); stem(n,real(x2)); xlabel('n');ylabel('real(x2)'); title('z2=1/12+j*pi/6时') subplot(5,2,4); stem(n,imag(x2)); xlabel('n');ylabel('image(x2)'); title('z2=1/12+j*pi/6时') subplot(5,2,5); stem(n,real(x3)); xlabel('n');ylabel('real(x3)'); title('z3=1/12时') subplot(5,2,6); stem(n,imag(x3)); xlabel('n');ylabel('image(x3)'); title('z3=1/12时')
数字信号处理习题解答 第二章 数据采集技术基础 2.1 有一个理想采样系统，其采样角频率Ωs =6π，采样后经理想低通滤波器H a (j Ω)还原，其中 ?? ???≥Ω<Ω=Ωππ 3032 1 )(，，j H a 现有两个输入，x 1(t )=cos2πt ，x 2(t )=cos5πt 。试问输出信号y 1(t )， y 2(t )有无失真？为什么？ 分析：要想时域采样后能不失真地还原出原信号，则采样角频率Ωs 必须大于等于信号谱最高角频率Ωh 的2倍，即满足Ωs ≥2Ωh 。 解：已知采样角频率Ωs =6π，则由香农采样定理，可得 因为x 1(t )=cos2πt ，而频谱中最高角频率ππ π32 621 =< =Ωh ， 所以y 1(t )无失真； 因为x 2(t )=cos5πt ，而频谱中最高角频率ππ π32 652 => =Ωh ， 所以y 2(t )失真。 2.2 设模拟信号x (t )=3cos2000πt +5sin6000πt +10cos12000πt ，求： （1） 该信号的最小采样频率； （2） 若采样频率f s =5000Hz ，其采样后的输出信号； 分析：利用信号的采样定理及采样公式来求解。 ○ 1采样定理 采样后信号不失真的条件为：信号的采样频率f s 不小于其最高频
率f m 的两倍，即 f s ≥2f m ○ 2采样公式 )()()(s nT t nT x t x n x s === 解：（1）在模拟信号中含有的频率成分是 f 1=1000Hz ，f 2=3000Hz ，f 3=6000Hz ∴信号的最高频率f m =6000Hz 由采样定理f s ≥2f m ，得信号的最小采样频率f s =2f m =12kHz （2）由于采样频率f s =5kHz ，则采样后的输出信号 ? ?? ? ????? ??-???? ????? ??=? ??? ????? ??+???? ????? ??-???? ????? ??=? ??? ????? ??++???? ????? ??-+???? ????? ??=? ??? ????? ??+???? ????? ??+???? ????? ??=? ?? ? ??====n n n n n n n n n n n f n x nT x t x n x s s nT t s 522sin 5512cos 13512cos 10522sin 5512cos 35112cos 105212sin 5512cos 3562cos 10532sin 5512cos 3)()()(πππππππππππ 说明：由上式可见，采样后的信号中只出现1kHz 和2kHz 的频率成分， 即 kHz f f f kHz f f f s s 25000200052150001000512211 ======，， 若由理想内插函数将此采样信号恢复成模拟信号，则恢复后的模拟信号
21世纪高等院校电子信息类规划教材 安徽省高等学校“十二五”省级规划教材 数字信号处理与DSP实现技术 课后习题与参考答案 主编：陈帅 副主编：沈晓波
淮南师范学院 2015.11 第1章绪论思考题 1．什么是数字信号？ 2．什么是数字信号处理？ 3．数字信号处理系统的实现方法有哪些？ 4．数字信号处理有哪些应用？ 5．数字信号处理包含哪些内容？ 6．数字信号处理的特点是什么？ 第1章绪论参考答案 1．时间和幅度都离散的信号称为数字信号，即信号的时间取离散的值，幅度也取离散的值。 2．数字信号处理是指在数字领域进行数字信号的加工（变换、运算等），即输入是数字信号，采用数字信号处理方法进行处理，输出仍然是数字信号。 3．数字信号处理系统的实现方法有①通用软件方法实现系统；②专用加速处理机方法；③软硬件结合的嵌入式处理方法；④硬件方法。 4．数字信号处理在通信、计算机网络、雷达、自动控制、地球物理、声学、天文、生物医学、消费电子产品等各个领域均有应用，是信息产业的核心技术之一。比如信源编码、信道编码、多路复用、数据压缩，数字语音、汽车多媒体、MP3/MP4/MP5、数字扫面仪、数字电视机顶盒、医院监视系统、生物指纹系统等。 5．数字信号处理主要包含以下几个方面的内容 ①离散线性时不变系统理论。包括时域、频域、各种变换域。 ②频谱分析。FFT谱分析方法及统计分析方法，也包括有限字长效应谱分析。 ③数字滤波器设计及滤波过程的实现（包括有限字长效应）。 ④时频-信号分析（短时傅氏变换），小波变换，时-频能量分布。 ⑤多维信号处理（压缩与编码及其在多煤体中的应用）。 ⑥非线性信号处理。 ⑦随机信号处理。 ⑧模式识别人工神经网络。 ⑨信号处理单片机(DSP)及各种专用芯片(ASIC)，信号处理系统实现。 6．数字信号处理主要具有4个方面优点：①数字信号精度高；②数字信号处理灵活性强；③数字信号处理可实现模拟信号难以实现的特性；④数字信号处理可以实现多维信号处理。
实验3 FFT算法的应用 一、实验目的 1、加深对离散信号的DFT的理解； 2、在MATLAB中实现FFT算法。 二、实验内容 1、2N点实数序列 N=64。用一个64点的复数FFT程序，一次算出，并绘出的图形。
程序代码如下： N=64; n1=0:2:(2*N-1);n2=1:2:(2*N-1); x1=cos(2*pi*7*n1/N)+1/2*cos(2*pi*9*n1/N); x2=cos(2*pi*7*n2/N)+1/2*cos(2*pi*9*n2/N); X1=fft(x1,N);X2=fft(x2,N); X=X1+X2; Y=abs(X); stem(Y); 10203040506070
2、已知某序列在单位圆上的N=64等分样点的Z 变换为: 。 用N 点IFFT 程序计算出和。 程序代码如下： N=64; k=0:1:N-1; X=1./(1-0.8*exp(-j*2*pi*k/N)); x=ifft(X,N); stem(abs(x)); 0.20.40.60.811.21.4
三、利用至少两种频率估计方法编程求给定信号的频率（要求运算 过程中使用小于等于128点的FFT） 无噪声、20dB、15dB、10dB、5dB、0dB、-5dB 已知信号的采样率Fs= 4000 Hz 方法一：谱峰检测法 程序代码如下： function f=wxk3 load signal; N=128; Fs=4000; f=zeros(1,7); y=abs(fft(x,N)); [m,n]=max(y); f(1)=(n-1)*Fs/N; SNR=[20 15 10 5 0 -5]; for i=1:6 x=awgn(x,SNR(i)); y=abs(fft(x,N)); [m,n]=max(y); f(i+1)=(n-1)*Fs/N; end plot(f);hold;
Chapter 2 Solutions 2.1 最小采样频率为两倍的信号最大频率，即44.1kHz 。 2.2 (a)、由ω = 2πf = 20 rad/sec ，信号的频率为f = 3.18 Hz 。信号的奈奎斯特采样频率为6.37 Hz 。 (b)、3 5000π=ω，所以f = 833.3 Hz ，奈奎斯特采样频率为1666.7 Hz 。 (c)、7 3000π=ω，所以f = 214.3 Hz ，奈奎斯特采样频率为428.6 Hz 。 2.3 (a) 1258000 1f 1T S S ===μs (b)、最大还原频率为采样频率的一半，即4000kHz 。 2.4 ω = 4000 rad/sec ，所以f = 4000/(2π) = 2000/π Hz ，周期T = π/2000 sec 。因此，5个周期为5π/2000 = π/400 sec 。对于这个信号，奈奎斯特采样频率为2(2000/π) = 4000/π Hz 。所以采样频率为f S = 4(4000/π) = 16000/π Hz 。因此5个周期收集的采样点为(16000/π samples/sec )(π/400 sec) = 40。 2.5 ω = 2500π rad/sec ，所以f = 2500π/(2π) = 1250 Hz ，T = 1/1250 sec 。因此，5个周期为5/1250 sec 。对于这个信号，奈奎斯特采样频率为2(1250) = 2500 Hz ，所以采样频率为f S = 7/8(2500) = 2187.5 Hz 。采样点数为(2187.5 点/sec)(5/1250 sec) = 8.75。这意味着在模拟信号的五个周期内只有8个点被采样。事实上，对于这个信号来说，在整数的模拟周期中，是不可能采到整数个点的。 2.6 2.7 信号搬移发生在kf S ± f 处，换句话说，频谱搬移发生在每个采样频率的整数倍 0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 频率/kHz
实验二:FIR 数字滤波器设计与软件实现 一、实验指导 1.实验目的 (1掌握用窗函数法设计 FIR 数字滤波器的原理和方法。 (2掌握用等波纹最佳逼近法设计 FIR 数字滤波器的原理和方法。 (3掌握 FIR 滤波器的快速卷积实现原理。 (4学会调用 MA TLAB 函数设计与实现 FIR 滤波器。 2. 实验内容及步骤 (1认真复习第七章中用窗函数法和等波纹最佳逼近法设计 FIR 数字滤波器的原理; (2调用信号产生函数 xtg 产生具有加性噪声的信号 xt ,并自动显示 xt 及其频谱,如图 1所示;
图 1 具有加性噪声的信号 x(t及其频谱如图 (3请设计低通滤波器,从高频噪声中提取 xt 中的单频调幅信号,要求信号幅频失真小于 0.1dB ,将噪声频谱衰减 60dB 。先观察 xt 的频谱,确定滤波器指标参数。 (4根据滤波器指标选择合适的窗函数,计算窗函数的长度 N ,调用 MATLAB 函数 fir1设计一个 FIR 低通滤波器。并编写程序,调用 MATLAB 快速卷积函数 fftfilt 实现对 xt 的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。 (5 重复 (3 , 滤波器指标不变, 但改用等波纹最佳逼近法, 调用MA TLAB 函数 remezord 和 remez 设计 FIR 数字滤波器。并比较两种设计方法设计的滤波器阶数。 提示:○ 1MA TLAB 函数 fir1的功能及其调用格式请查阅教材; ○ 2采样频率 Fs=1000Hz,采样周期 T=1/Fs;
○ 3根据图 1(b和实验要求,可选择滤波器指标参数:通带截止频率 fp=120Hz,阻带截 至频率 fs=150Hz, 换算成数字频率, 通带截止频率 p 20.24 p f ωπ =T=π, 通带最大衰为 0.1dB , 阻带截至频率 s 20.3 s f ωπ =T=π,阻带最小衰为 60dB 。 3、实验程序框图如图 2所示,供读者参考。 图 2 实验程序框图 4.信号产生函数 xtg 程序清单(见教材 二、滤波器参数及实验程序清单 1、滤波器参数选取 根据实验指导的提示③选择滤波器指标参数: 通带截止频率 fp=120Hz,阻带截至频率 fs=150Hz。代入采样频率 Fs=1000Hz,换算成 数字频率,通带截止频率 p 20.24 p f
数字信号处理是将信号以数字方式表示并处理的理论和技术。数字信号处理与模拟信号处理是信号处理的子集。 简介 简单地说，数字信号处理就是用数值计算的方式对信号进行加工的理论和技术，它的英文原名叫digital signal processing，简称DSP。另外DSP也是digital signal processor的简称，即数字信号处理器，它是集成专用计算机的一种芯片，只有一枚硬币那么大。有时人们也将DSP看作是一门应用技术，称为DSP 技术与应用。 《数字信号处理》这门课介绍的是：将事物的运动变化转变为一串数字，并用计算的方法从中提取有用的信息，以满足我们实际应用的需求。 本定义来自《数字信号处理》杨毅明著，由机械工业出版社2012年发行。 特征和分类 信号(signal)是信息的物理体现形式，或是传递信息的函数，而信息则是信号的具体内容。 模拟信号(analog signal)：指时间连续、幅度连续的信号。 数字信号(digital signal)：时间和幅度上都是离散(量化)的信号。 数字信号可用一序列的数表示，而每个数又可表示为二制码的形式，适合计算机处理。 一维(1-D)信号: 一个自变量的函数。 二维(2-D)信号: 两个自变量的函数。 多维(M-D)信号: 多个自变量的函数。 系统：处理信号的物理设备。或者说，凡是能将信号加以变换以达到人们要求的各种设备。模拟系统与数字系统。 信号处理的内容：滤波、变换、检测、谱分析、估计、压缩、识别等一系列的加工处理。 多数科学和工程中遇到的是模拟信号。以前都是研究模拟信号处理的理论和实现。 模拟信号处理缺点：难以做到高精度，受环境影响较大，可靠性差，且不灵活等。数字系统的优点：体积小、功耗低、精度高、可靠性高、灵活性大、易于大规模集成、可进行二维与多维处理 随着大规模集成电路以及数字计算机的飞速发展，加之从60年代末以来数字信号处理理论和技术的成熟和完善，用数字方法来处理信号，即数字信号处理，已逐渐取代模拟信号处理。 随着信息时代、数字世界的到来，数字信号处理已成为一门极其重要的学科和技术领域。 数字信号处理器 DSP芯片，也称数字信号处理器，是一种特别适合于进行数字信号处理运算的微处理器，其主要应用是实时快速地实现各种数字信号处理算法。根据数字信号处理的要求，DSP芯片一般具有如下主要特点： （1）在一个指令周期内可完成一次乘法和一次加法；
随堂练习 随堂练习提交截止时间：2017-06-25 00:00:00 当前页有10题，你已做10题，已提交10题，其中答对10题。 1.最早的单片机是Intel公司的（），它出现在1976年。 A．8088 B．8048 C．8058 D．8068 答题： A. B. C. D. （已提交） 参考答案：B 问题解析： 2.下面哪种操作系统不属于商用操作系统。（） A. Windows 7 B. Linux C. VxWorks D. WinCE 答题： A. B. C. D. （已提交） 参考答案：B 问题解析： 3.下面哪点不是嵌入式操作系统的特点。（） A. 内核精简 B. 专用性强 C. 功能强大 D. 高实时性 答题： A. B. C. D. （已提交） 参考答案：C 问题解析： 4.迄今为止，（）系列的单片机仍然是最为成功的单片机芯片之一，在各种产品中有着非常广泛的应用。 A．8088 B．8048 C．8058 D．8051 答题： A. B. C. D. （已提交） 参考答案：D 问题解析： 5.嵌入式系统的核心是（）和嵌入式操作系统。 A．嵌入式微处理器 B．存储器 C．单片机 D．内部总线 答题： A. B. C. D. （已提交）
参考答案：A 问题解析： 6.下面哪种嵌入式系统是属于软实时系统。（） A．WinCE B．VxWorks C．Nucleus D．eCOS 答题： A. B. C. D. （已提交） 参考答案：A 问题解析： 7.下面哪种嵌入式系统是属于硬实时系统。（） A．WinCE B．RTLinux C．Windows XP D．eCOS 答题： A. B. C. D. （已提交） 参考答案：D 问题解析： 8. Embedded Visual Studio是（） A．嵌入式硬件开发工具 B．嵌入式软件开发工具 C．嵌入式调试工具 D．嵌入式编译工具 答题： A. B. C. D. （已提交） 参考答案：B 问题解析： 9.下面哪些不属于嵌入式控制系统的典型应用。（） A．智能家居 B. 智能手机 C. 智能控制 D. 智能数据挖掘 答题： A. B. C. D. （已提交） 参考答案：D 问题解析： 10.下面哪个系统属于嵌入式系统（）。 A．“天河一号”计算机系统 B．IBMX200笔记本电脑 C．联想S10上网本 D．Iphone手机 答题： A. B. C. D. （已提交） 参考答案：D 问题解析： 随堂练习提交截止时间：2017-06-25 00:00:00
数字信号处理基础书后题答案中文版
Chapter 2 Solutions 2.1 最小采样频率为两倍的信号最大频率，即44.1kHz 。 2.2 (a)、由ω = 2πf = 20 rad/sec ，信号的频率为f = 3.18 Hz 。信号的奈奎斯特采样频率为6.37 Hz 。 (b)、35000π =ω，所以f = 833.3 Hz ，奈奎斯特采样频率为1666.7 Hz 。 (c)、7 3000π =ω，所以f = 214.3 Hz ，奈奎斯特采样频率为428.6 Hz 。 2.3 (a) 1258000 1f 1T S S === μs (b)、最大还原频率为采样频率的一半，即4000kHz 。 2.4 ω = 4000 rad/sec ，所以f = 4000/(2π) = 2000/π Hz ，周期T = π/2000 sec 。因此，5个周期为5π/2000 = π/400 sec 。对于这个信号，奈奎斯特采样频率为2(2000/π) = 4000/π Hz 。所以采样频率为f S = 4(4000/π) = 16000/π Hz 。因此5个周期收集的采样点为(16000/π samples/sec )(π/400 sec) = 40。 2.5 ω = 2500π rad/sec ，所以f = 2500π/(2π) = 1250 Hz ，T = 1/1250 sec 。因此，5个周期为5/1250 sec 。对于这个信号，奈奎斯特采样频率为2(1250) = 2500 Hz ，所以采样频率为f S = 7/8(2500) = 2187.5 Hz 。采样点数为(2187.5 点/sec)(5/1250 sec) = 8.75。这意味着在模拟信号的五个周期内只有8个点被采样。事实上，对于这个信号来说，在整数的模拟周期中，是不可能采到整数个点的。 2.7 信号搬移发生在kf S ± f 处，换句话说，频谱搬移发生在每个采样频率的整数 倍 -200 200 400 600 800 1000 1200 0.10.20.30.40.50.60.70.80.91 幅度 频
数字信号处理课程设计 《数字信号处理》 课程设计报告 FIR数字滤波器设计及MATLAB实现 专业：通信工程 班级：通信1101班 组次：第9组 姓名及学号： 姓名及学号：
目录 一、设计目的 (3) 二、设计任务 (3) 三、设计原理 (3) 3.1窗函数法 (3) 3.2频率采样法 (4) 3.3最优化设计 (5) 3.3.1等波纹切比雪夫逼近准则 (5) 3.3.2仿真函数 (6) 四、设计过程 (7) 五、收获与体会 (13) 参考文献 (13)
FIR数字滤波器设计及MATLAB实现 一、设计目的 FIR滤波器：有限长单位冲激响应滤波器，是数字信号处理系统中最基 本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性， 同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR 滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。滤波器设 计是根据给定滤波器的频率特性，求得满足该特性的传输函数。 二、设计任务 FIR滤波器设计的任务是选择有限长度的() H e满足一定 h n，使传输函数()jw 的幅度特性和线性相位要求。由于FIR滤波器很容易实现严格的线性相位，所以FIR 数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。 设计过程一般包括以下三个基本问题： （1）根据实际要求确定数字滤波器性能指标； （2）用一个因果稳定的系统函数去逼近这个理想性能指标； （3）用一个有限精度的运算去实现这个传输函数。 三、设计原理 FIR滤波器设计的任务是选择有限长度的() H e满足一定 h n，使传输函数()jw 的幅度特性和线性相位要求。由于FIR滤波器很容易实现严格的线性相位，所以FIR数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。 设计过程一般包括以下三个基本问题： （1）根据实际要求确定数字滤波器性能指标； （2）用一个因果稳定的系统函数去逼近这个理想性能指标； （3）用一个有限精度的运算去实现这个传输函数。 3.1窗函数法 设计FIR数字滤波器的最简单的方法是窗函数法，通常也称之为傅立叶级数法。FIR数字滤波器的设计首先给出要求的理想滤波器的频率响应()jw H e,设计 d
数 字 信 号 处 理 发 展 和 应 用 学院：通信学院 专业：电子信息工程 班级：电信1103 姓名：XXX 学号：XXX
数字信号处理发展和应用 【摘要】数字信号处理（DSP）是广泛应用于许多领域的新兴学科，因其具有可程控、可预见性、精度高、稳定性好、可靠性和可重复性好、易于实现自适应算法、大规模集成等优点，广泛应用于实时信号处理系统中。本文概述了DSP 技术的发展历史，各个领域的应用状况，以及在未来的发展趋势。 【关键词】数字信号处理；数据处理；信息技术；发展趋势 一、数字信号处理（DSP）的发展历史 数字信号处理技术的发展经历了三个阶 段。 70 年代DSP 是基于数字滤波和快速傅立叶变换的经典数字信号处理,其系统由分立的小规模集成电路组成,或在通用计算机上编程来实现DSP 处理功能,当时受到计算机速度和存储量的限制,一般只能脱机处理,主要在医疗电子、生物电子、应用地球物理等低频信号处理方面获得应用。 80 年代DSP 有了快速发展,理论和技术进入到以快速傅立叶变换(FFT) 为主体的现代信号处理阶段,出现了有可编程能力的通用数字信号处理芯片,例如美国德州仪器公司(TI 公司) 的TMS32010 芯片,在全世界推广应用,在雷达、语音通信、地震等领域获得应用,但芯片价格较贵,还不能进入消费领域应用。 90 年代DSP 技术的飞速发展十分惊人,理论和技术发展到以非线性谱估计为代表的更先进的信号处理阶段,能够用高速的DSP 处理技术提取更深层的信息,硬件采用更高速的DSP 芯片,能实时地完成巨大的计算量,以TI 公司推出的TMS320C6X芯片为例,片内有两个高速乘法器、6 个加法器,能以200MHZ频率完成8 段32 位指令操作,每秒可以完成16 亿次操作,并且利用成熟的微电子工艺批量生产,使单个芯片成本得以降低。并推出了C2X、C3X、C5X、C6X 不同应用范围的系列,使新一代的DSP 芯片在移动通信、数字电视和消费电子领域得到广泛应用,数字化的产品性能价格比得到很大提高,占有巨大的市场。 二、数字信号处理（DSP）的主要应用领域 1·DSP在电力系统自动化中日益渗透 1.1数字信号处理（DSP）技术在电力系统模拟量采集和测量中的应用 计算机进入电力系统调度后，引入了EMS／DMS／SCADA的概念，而电力系统数据采集和测量是SCADA的基础部分。传统的模拟量的采集和获得，通过变送器将一次PT和CT的电气量变为直流量，再进行A/D转换送给计算机。应用了交流采样技术以后，经过二次PT、CT的变换后，直接对每周波的多点采样值采用DSP处理算法进行计算，得到电压和电流的有效值和相角，免去了变送器环节。这不仅使得分散布置的分布式RTU很快地发展起来，而且还为变电站自动化提供了功能综合优化的手段。 1.2数字信号处理（DSP）在继电保护中的应用 到目前为止，应用于我国电力系统的微机保护产品采用的CPU大多为单片机，由于受硬件资源及计算功能的限制，其采样能力及采样速度很难令人满意。因此，对非正常运行条件下的系统参数测量，在速度和精度上无法满足要求，一些复杂原理和算法的实现，基于常规CPU的保护产品也都难以胜任。基于DSP 的数据采集和处理系统由于其强大的数学运算能力和特殊设计，都使得它在继
实验1 常见离散信号产生和实现 一、实验目的 1、加深对常用离散信号的理解； 2、熟悉使用MATLAB在时域中产生一些基本的离散时间信号。 二、实验原理 1、单位抽样序列 在MATLAB中可以利用函数实现。 2、单位阶越序列 在MATLAB中可以利用函数实现: 3、正弦序列 在MATLAB中实现过程如下： 4、复指数序列
在MATLAB中实现过程如下： 5、指数序列 在MATLAB中实现过程如下： 三、预习要求 1、预先阅读实验讲义（MATLAB基础介绍）； 2、讨论正弦序列、复指数序列的性质。 A．绘出信号，当、时、、 时的信号实部和虚部图；当时呢？此时信号周期为多少？ B．绘出信号的频率是多少？周期是多少？产生一个数字频率为0.9的正弦序列，并显示该信号，说明其周期。 3、使用帮助功能学习square(方波),sawtooth(锯齿波)和sinc函数，并绘图。 四、实验内容 编制程序产生上述5种信号，长度可输入确定，函数需要的参数可输入确定，并绘出其图形。
实验2 离散系统的时域分析 一、实验目的 1、熟悉并掌握离散系统的差分方程表示法； 2、加深对冲激响应和卷积分析方法的理解。 二、实验原理 在时域中，离散时间系统对输入信号或者延迟信号进行运算处理，生成具有所需特性的输出信号，具体框图如下: 其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号, 记系统单位冲激响应，则系统响应为如下的卷积计算式： 当时，h[n]是有限长度的（），称系统为FIR系统；反之，称系统为IIR系统。 三、预习要求
1、在MATLAB中，熟悉利用函数实现差分方程的仿真； 2、在MATLAB中，熟悉用函数计算卷积，用求系统冲激响应的过程。 四、实验内容 1、以下程序中分别使用conv和filter函数计算h和x的卷积y和y1，运行程序，并分析y和y1是否有差别，为什么要使用x[n]补零后的x1来产生y1；具体分析当h[n]有i个值，x[n]有j个值，使用filter完成卷积功能，需要如何补零？ % Program P2_7 clf; h = [3 2 1 -2 1 0 -4 0 3]; %impulse response x = [1 -2 3 -4 3 2 1]; %input sequence y = conv(h,x); n = 0:14; subplot(2,1,1); stem(n,y); xlabel('Time index n'); ylabel('Amplitude'); title('Output Obtained by Convolution'); grid; x1 = [x zeros(1,8)]; y1 = filter(h,1,x1); subplot(2,1,2); stem(n,y1); xlabel('Time index n'); ylabel('Amplitude'); title('Output Generated by Filtering'); grid; 2、编制程序求解下列两个系统的单位冲激响应和阶跃响应，并绘出其图形。要求分别用 filter、conv、impz三种函数完成。 , ，
==============================绪论============================== 1. A/D 8bit 5V 00000000 0V 00000001 20mV 00000010 40mV 00011101 29mV ==================第一章 时域离散时间信号与系统================== 1. ①写出图示序列的表达式 答：3)1.5δ(n 2)2δ(n 1)δ(n 2δ(n)1)δ(n x(n)-+---+++= ②用δ(n) 表示y (n )={2,7,19,28,29,15} 2. ①求下列周期 ) 5 4sin( )8 sin( )4() 51 cos()3() 54sin()2() 8sin( )1(n n n n n π π π π - ②判断下面的序列是否是周期的; 若是周期的， 确定其周期。 （1）A是常数 8ππn 73Acos x(n)??? ? ??-= （2）)8 1 (j e )(π-=n n x 解： （1） 因为ω= 73π, 所以314 π2=ω, 这是有理数， 因此是周期序列， 周期T =14。 （2） 因为ω= 81, 所以ω π2=16π, 这是无理数， 因此是非周期序列。 ③序列)Acos(nw x(n)0?+=是周期序列的条件是是有理数2π/w 0。
3.加法 乘法 序列{2，3，2，1}与序列{2，3，5，2，1}相加为__{4，6，7，3，1}__，相乘为___{4，9，10，2} 。 移位 翻转：①已知x(n)波形，画出x(-n)的波形图。 ② 尺度变换：已知x(n)波形，画出x(2n)及x(n/2)波形图。 卷积和：①h(n)*求x(n)，其他0 2 n 0n 3，h(n)其他03n 0n/2设x(n) 例、???≤≤-=???≤≤= }2 3 ,4,7,4，23{0,h(n)*答案：x(n)= ②已知x (n )={1,2,4,3},h (n )={2,3,5}, 求y (n )=x (n )*h (n ) x (m )={1,2,4,3},h (m )={2,3,5}，则h (-m )={5,3,2}(Step1:翻转) 解得y (n )={2,7,19,28,29,15} ③(n)x *(n)x 3)，求x(n)u(n u(n)x 2),2δ(n 1)3δ(n δ(n)2、已知x 2121=--=-+-+= }{1,4,6,5,2答案：x(n)= 4. 如果输入信号为 ，求下述系统的输出信号。
FIR 数字滤波器设计与实现 一.摘要：数字滤波器是一种具有频率选择性的离散线性系统，在信号数字处理中有着广泛的应 用。其中FIR 滤波器是一种常用的滤波器，它在保证幅度特性满足技术要求的同时，很容易做到严格的线性相位特性，在语音分析、图像处理、雷达监测等对信号相位要求高的领域有着广泛的应用，能实现IIR 滤波器不能实现的许多功能。 二.关键词：FIR 窗函数系统函数MATLAB 三.内容提要： 数字滤波器的功能就是把输入序列通过一定的运算变换成输出序列，因此数字滤波器的结构系 统中就必须包括一定数量和性能的运算器件和运算单元，而运算器件和运算单元的配置必须由数字滤波器的结构特点和性能特点来决定，因此在进行FIR 数字滤波器的设计之前，有必要介绍和总结FIR 数字滤波器的基本结构和相关特性（包括频响曲线（幅度和相位），单位冲激响应等），在介绍完其基本结构和相关特性后，就进行FIR 数字滤波器的设计和实现。 （一）FIR 滤波器的基本结构 在讨论任何一种滤波器时，都要着重分析其系统函数，FIR 滤波器的系统函数为： n N n z n h z H ∑-==1 0)()(。从该系统函数可看出，FIR 滤波器有以下特点： 1)系统的单位冲激响应h(n)在有限个n 值处不为零； 2)系统函数H(z)在|z|>0处收敛，极点全部在z=0处(稳定系统)； 3)结构上主要是非递归结构，没有输出到输入的反馈，但有些结构中(例如频率抽样结构)也包 含有反馈的递归部分。 1.FIR 滤波器实现的基本结构有： 1） 横截型（卷积型、直接型） a.一般FIR 滤波器的横截型(直接型、卷积型)结构： 若给定差分方程为： 。则可以直接由差分方程得出FIR 滤波器结构如 下图所示： 这就是FIR 滤波器的横截型结构，又称直接型或卷积型结构。 b .线性相位FIR 滤波器的横截型结构