数字信号处理 课程简介

数字信号处理（Digital Signal Processing）

课程编号：01115161

课程性质：专业方向任选课

开设学期及学时分配：第六期；每学期64学时

适用专业及层次：电子信息工程本科

先行课程：《数字电子技术》、《信号与系统》

后继课程：无

课程目的、内容与要求：

本课程是信息工程本科专业必修课，它是在学生学完了高等数学、概率论、线性代数、复变函数、信号与系统等课程后，进一步为学习专业知识打基础的课程。本课程将通过讲课、练习使学生建立“数字信号处理”的基本概念，掌握数字信号处理基本分析方法和分析工具，为从事通信、信息或信号处理等方面的研究工作打下基础。

教材：吴镇扬编，《数字信号处理》，高等教育出版社，2004年9月第一版。

推荐参考书：

1. 姚天任,江太辉编，《数字信号处理》（第二版），华中科技大学出版社，2000年版。

2. 程佩青著，《数字信号处理教程》（第二版），清华大学出版社出版，2001年版。

3. 丁玉美,高西全编著，《数字信号处理》，西安电子科技大学出版社，2001年版。

4. 胡广书编，《数字信号处理——理论、算法与实现》，清华大学出版社，2004年版。

授课教师：

1.主讲教师要具有中级及以上专业技术职称和硕士研究生及以上学历。

2.能履行教师职责，爱岗敬业，为人师表，具有良好的师德教风和较高的业务水平。

3.本课程要求教师应具备数据结构与算法设计、编译原理和软件工程等方面的理论基础，熟悉傅立叶变换（FFT）与数字滤波器并具备分析和运用MATLAB编程能力。

教学与实验设施：

本课程在多媒体教室开展，多媒体要满足课程教学需要，能同时运行office的课件和MATLAB软件。

教学方法与手段：

本课程教学方法要灵活，可用多媒体课件与板书相结合的教学形式，有些内容可以通过实物或图片演示。教学要充分发挥学生主体性，与学生建立起平等、民主和对话的师生关系，培养学生发现问题、分析问题、解决问题的能力和探究意识，使学生掌握本课程的核心内容外，指导学生对相关外延知识的获取能力。

课程考核与评价：

本课程的考核应该包括平时成绩、期末考试和期中成绩3个部分，实行百分制。其中平时成绩可以通过个人作业、学习态度、到课率及小组讨论等方式进行评定，期末考试可以采用开卷或闭卷形式，重点考查学生对数字信号处理的基本概念、基本理论、基本方法的理解和掌握程度。

自控原理选择题

1.反馈控制系统又称为（B ） A.开环控制系统 B ．闭环控制系统 B.扰动顺馈补尝 D ．输入顺馈补偿 2.位置随动系统的主反馈环节是（A ） A ．电压负反馈 B ．电流负反馈 C ．转速负反馈 D ．位置负反馈 3.如果典型二阶系统的单位阶跃响应为减幅振荡(又称阻尼振荡)，则其阻尼比（C ） Aξ<0 Bξ=0 C0<ξ<1 D ξ≥1 4.G(s)= 1/[(S+1)(S+2)(S+3)(S+4)]环 节的对数相频特性的高频渐近线斜率为（D ） A -20d B B-40dB C-60dB D-80dB 5.某自控系统的开环传递函数G(s)= 1/[(S+1)(S+2)] ，则此系统为（ A ） A ．稳定系统 B ．不稳定系统 C ．稳 定边界系统 D ．条件稳定系统 6．若一系统的特征方程式为(s+1)2(s －2)2+3＝0，则此系统是（C ） A 稳定的 B 临界稳定的 C 不稳定 D ．条件稳定的 7.下列性能指标中的(D )为系统的稳态指标。 A.σP B.t s C.N D.e ss 8.下列系统中属于开环控制的为( C) A.自动跟踪雷达 B.数控加工中心 C.普通车床 D.家用空调器 9.RLC 串联电路的系统应为(D)环节。 A 比例 B.惯性 C.积分 D.振荡 10.输出信号与输入信号的相位差随频率变化的关系是(B )。 A.幅频特性 B.相频特性 C.传递函数 D.频率响应函数 1.奈奎斯特图分析闭环控制系统的 （A ）A.稳态性能 B.动态性能 C.稳态和动态性能 D.抗扰性能 2.有一线性系统，其输入分别为u 1(t)和u 2(t)时，输出分别为y 1(t)和y 2(t)。 当输入为a 1u 1(t)+a 2u 2(t)时(a 1,a 2为常 数)，输出应为（B ） A.a 1y 1(t)+y 2(t) B.a 1y 1(t)+a 2y 2(t) C.a 1y 1(t)-a 2y 2(t) D.y 1(t)+a 2y 2(t) 3.某串联校正装置的传递函数为 G c (S)=K S S T 1T 1+β+(0<β<1)，则该装置 是（ A ） A.超前校正装置 B.滞后校正装置 C.滞后超前 D.超前滞后校正装置 4.1型系统开环对数幅频渐近特性的 低频段斜率为（ B ） A.-40(dB/dec) B.-20(dB/dec) C.0(dB/dec) D.+20(dB/dec) 5.开环传递函数G(s)H(s)=) p s )(p s () z s (K 211+++，其中p 2>z 1>p 1>0，则实轴上的根轨迹（A ） A.（-∞，-p 2],[-z 1,-p 1] B.(- ∞,-p 2] C.[-p 1,+ ∞] D.[-z 1,-p 1] 6.设系统的传递函数为 G(s)=1 s 5s 2512 ++，则系统的阻尼比为 C. 21 7.设单位负反馈控制系统的开环传递函数G o (s)=) a s (s K +，其中K>0,a>0，则闭环控制系统的稳定性与（ D ） A.K 值的大小有关 B .a 值的大小有关 C.a 和K 值的有关 D.a 和K 值的无关 8. 在伯德图中反映系统动态特性的是 B. 中频段 9. 设开环系统的频率特性G(j ω)=2 ) j 1(1 ω+，当ω=1rad/s 时，其频率特性幅值G(1)=(D ) A. 1 B. 2 C.21 D. 41 10. 开环传递函数为G(s)H(s)=) 3s (s K 3+,则实轴上的根轨迹为( D )。 A.［-3,∞］ B. ［0,∞］ C. (-∞,-3) D. ［-3,0］ 1.实验中可以从( D )获取频率特性。A.稳定的线性和非线性系统 B. 不稳定的线性和非线性系统 C.不 稳定的线性系统 D. 稳定的线性系统2.传递函数的概念适用于（D ）系统。 A ．线性、非线性 B. 线性 非时变 C ．非线性定常 D. 线性定常 3.系统的动态性能包括（ B ）。 A 稳定性平稳性 B.平稳性快速性 C 快速性稳定性 D.稳定性准确性 4. 确定系统根轨迹的充要条件是C A 根轨迹的模方程B.根轨迹的相方程 C 根轨迹增益 D 根轨迹方程的阶次 5 ．正弦信号作用于线性系统所产生的频率响应是（ A ） A ．输出响应的稳态分量 B. 输出响应的暂态分 量 C ．输出响应的零输入分量 D. 输出响应的零状态分量 6．系统的传递函数完全决定于系统 的 ( C )。A ．输入信号 B.输出信号C.结构和参数D.扰动信号7．控制系统的相位稳定裕量反咉 了系统的 ( B )。A ．稳定性 B.稳态性能C.快速性 D.动态性能8．一般来说，系统增加积分环节，系统的稳定性将（ B ）。 A ．变好 B.变坏 C.不变 D.可能9．系统开环对数幅频特性L(ω)中频段主要参数的大小对系统的 （ D ）性能无影响。A.动态 B. 稳态 C. 相对稳定性 D. 响应的快速性10．反馈控制系统又称为（ B ） A ．开环控制系统 B ．闭环控制系统 C 扰动顺馈补偿 D 输入顺馈补偿 1.单位斜坡函数f(t)=t 的拉氏变换式F(s)=（ D ） A.s B.1 C ．S 2 D ． 1/S 2 2．单位抛物线输入函数r(t)的数学表达式是r(t)＝（D ） A ．at 2 B ．1/2 Rt 2 C ．t 2 3．当二阶系统特征方程的根为具有负实部的复数根时，系统的阻尼比C A ζ<0 B ζ＝0 C0<ζ<1 D ζ≥1 4．已知单位反馈控制系统在阶跃函数作用下，稳态误差e ss 为常数，则此系统为B ．I 型系统 5．设某环节的传递函数为G(s)＝1 21 +s ，当ω＝0.5rad ／s 时，其频率特性相位移θ(0.5)=(A) A ．－ 4π B ．－6π C ．6π D 4 π 6超前校正装置的最大超前相角趋D ．90°

测控技术与仪器专业简介

测控技术与仪器专业 业务培养目标： 本专业培养具备精密仪器设计制造以及测量与控制方面基础知识与应用能力，能在国民经济各部门从事测量与控制领域内有关技术、仪器与系统的设计制造、科技开发、应用研究、运行管理等方面的高级工程技术人才。 修业年限：四年 授予学位：工学学士 业务培养要求： 本专业学生主要学习精密仪器的光学、机械与电子学基础理论，测量与控制理论和有关测控仪器的设计方法，受到现代测控技术和仪器应用的训练，具有本专业测控技术及仪器系统的应用及设计开发能力。 专业方向介绍 测控技术及仪器专业是仪器科学与技术和控制科学与技术交叉融合而形成的综合性学科。 设2个专业方向。 方向一：检测技术与自动化装置方向； 方向二：测试计量技术及仪器方向。 方向一以集电子技术、先进控制理论、计算机控制技术、自动检测技术、光电技术以及网络技术于一体为特色，以生产过程的机电装备运行状态及其信息为研究对象。本方向旨在培养基础理论扎实、实践能力强、知识面广，外语综合能力和计算机应用能力较强，人文社会科学综合素质较高，具有开拓创性意识，能够从事工业过程控制理论与装备、计算机辅助测试系统、信息处理与状态识别等领域的研究开发、设计制造和运行管理的复合型高级工程技术人才。 方向二以光—机—电—仪器—计算机技术一体化为特色，以传感器技术、信息获取与处理技术、自动化精密机械以及智能仪器仪表为主要研究对象。本方向旨在培养基础理论扎实、实践能力强、知识面广，外语综合能力和计算机应用能力较强，人文社会科学综合素质较高，具有开拓创性意识，能够从事测控仪器、信息技术以及测试计量技术等方面的研究开发、设计制造和运行管理方面的复合型高级工程技术人才。 业务能力 方向一的毕业生应具有较扎实的自然科学基础，较好的人文和社会科学基础及较强的英语与计算机应用能力以及较强的创新意识；系统地掌握检测技术与自动化装置专业方向的基本理论与技术，主要包括电工电子技术、自动检测技术、工程光学、测控仪器电路、工业过程控制、微机控制技术等基本理论基础；掌握光、机、电、计算机控制相结合的现代测控技术和实验研究技能；具备综合运用专业知识解决生产实际问题的初步能力。 方向二的毕业生应具有较扎实的自然科学基础，较好的人文和社会科学基础以及较强的英语和计算机应用能力、较强的创新意识；系统地掌握本专业所需的基本理论和基础知识，主要包括电子技术、工程光学、精密机械学、传感器技术、控制工程等基础知识；掌握光、机、电、计算机相结合的现代测控技术和实验技能，综合运用专业知识解决生产实际问题的初步能力。

数字信号处理习题集(附答案)

第一章数字信号处理概述 简答题： 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==Ω= π

《数字信号处理》课程研究性学习报告解读

《数字信号处理》课程研究性学习报告 指导教师薛健 时间2014.6

【目的】 (1) 掌握IIR 和FIR 数字滤波器的设计和应用； (2) 掌握多速率信号处理中的基本概念和方法 ； (3) 学会用Matlab 计算小波分解和重建。 (4)了解小波压缩和去噪的基本原理和方法。 【研讨题目】 一、 (1)播放音频信号 yourn.wav ，确定信号的抽样频率，计算信号的频谱，确定噪声信号的频率范围； (2)设计IIR 数字滤波器，滤除音频信号中的噪声。通过实验研究s P ,ΩΩ，s P ,A A 的选择对滤波效果及滤波器阶数的影响，给出滤波器指标选择的基本原则，确定你认为最合适的滤波器指标。 （3）设计FIR 数字滤波器，滤除音频信号中的噪声。与（2）中的IIR 数字滤波器，从滤波效果、幅度响应、相位响应、滤波器阶数等方面进行比较。 【设计步骤】 【仿真结果】

【结果分析】 由频谱知噪声频率大于3800Hz。FIR和IIR都可以实现滤波，但从听觉上讲，人对于听觉不如对图像（视觉）明感，没必要要求线性相位，因此，综合来看选IIR滤波器好一点，因为在同等要求下，IIR滤波器阶数可以做的很低而FIR滤波器阶数太高，自身线性相位的良好特性在此处用处不大。【自主学习内容】 MATLAB滤波器设计 【阅读文献】 老师课件，教材 【发现问题】(专题研讨或相关知识点学习中发现的问题)： 过渡带的宽度会影响滤波器阶数N 【问题探究】 通过实验，但过渡带越宽时，N越小，滤波器阶数越低，过渡带越窄反之。这与理论相符合。 【仿真程序】 信号初步处理部分： [x1,Fs,bits] = wavread('yourn.wav'); sound(x1,Fs); y1=fft(x1,1024); f=Fs*(0:511)/1024; figure(1) plot(x1) title('原始语音信号时域图谱'); xlabel('time n'); ylabel('magnitude n'); figure(2) freqz(x1) title('频率响应图') figure(3) subplot(2,1,1); plot(abs(y1(1:512))) title('原始语音信号FFT频谱') subplot(2,1,2); plot(f,abs(y1(1:512))); title(‘原始语音信号频谱') xlabel('Hz'); ylabel('magnitude'); IIR： fp=2500;fs=3500; wp = 2*pi*fp/FS; ws = 2*pi*fs/FS; Rp=1; Rs=15;

数字信号处理实验一

一、实验目的 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;

浙江大学845自动控制原理考研真题试卷

紧急通知 本资料由浙江大学控制科学与工程学院16届专业课129分学长，也就是我本人亲自整理编排而成。大家可以叫我学长，年龄比我大的辞职考的可以叫我小弟。资料不同于市面上那些看起来非常诱人实则是粗制烂造的资料，而是以一个考过845自控的过来人的经验，完全从学生的体验出发，做到资料最全，资料最好，资料最精致。全套资料包括葵花宝典一到葵花宝典九共九本资料，每本资料都是我精心编辑整理的，并做了精美的封面，一共650页完美打印发给大家，大家把这650从头到尾肯透了，再做下我推荐的几本资料书（16年有一道15分的大题就是上面的类似题，第三问很多高手都没做出来，注意不是周春晖那本哈），可以说完全没问题了。这是其它卖家不可能做到的。同时赠送845自控全套电子资料。葵花宝典一完全由我本人原创，里面包含了考浙大845自动控制原理的全部问题，比如考多少分比较保险，怎么复习，有哪些好的资料书，最近几年考题变化及应对策略，浙大常考题型，招生名额，复试资料，导师联系，公共课复习用书及方法以及845近年命题风格分析等一系列问题，全是我的心得和经验，方法，技巧等，说句心里话，我自己都觉得这些资料非常宝贵，能帮助学弟学妹们少走很多弯路。 注意：前面是一些关于我的故事，有些地方可能对你有用，如果不感兴趣，可以直接拉到后面去看，资料清单和图片都在后面。 学长自我介绍 学长姓邓，名某某，男，本科于14年毕业于四川大学电气信息学院自动化专业，考浙大控制考了3次，14年大三时第一次考浙大控制总分没过线。当时我们学校有三个同学征战浙大控制科学与工程，结果全军覆没，只有我一人过了300分，由此可见考浙大控制还是很有难度的，其中一个难点就是专业课的信息和专业课的命题走向的获取，当时我们都不是很清楚，蒙着头自己学，去图书馆借了很多自动控制原理的资料书来看，我自我感觉学得还不错，当时我一个同学考电子科大的自动化，经常跑来问我自控的问题，我基本都能给他解答出来，他说我好牛逼，觉对没有问题，然而最后的结果是他考电子科大自动控制原理137，而我只考了96分。后面我分析了一下，为什么会出现这样的情况？最重要的就是我们对浙大的出题风格不是很了解，不知道它的命题方向和爱考的地方，方向都错了，怎么可能得高分？虽然我把11年以前的真题都做了，但是浙大12年以后的命题风格和以前有所不同，所以还是无济于事。因此即使你的自控基础知识扎实，也未必能够得到高分，这里面有很多方法和技巧，都是我从后面的考试中慢慢总结出来的。 由于不甘心就这么与浙大失之交臂，所以决定二战，但是又不想向家里要钱了，因为学长家在贵州农村，经济条件不是很好。于是我选择平时晚上去给别人做家教，周末去给培训机构上课。这样的好处是我有大把的白天用来复习，只是晚上出去干干活。这个事就说到这里，不是主题。15年专业课考了113，一个中等的分数，本来可以考130，但是为什么没有考到，这些原因我都在葵花宝典一中给大家分析了，希望大家能我的身上汲取经验，别步我的后尘。但是15年死在英语不过线上，差3分，这是我怎么也没有想到的，学长英语虽然不能说特别好，但是最起码四六级大一就过了，高考英语还是我们小县城的单科第一名（山中无老虎），第一年也考了65分。这是我怎么也没有想到的，所以有的时候感觉

信号处理-习题(答案)

数字信号处理习题解答 第二章 数据采集技术基础 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 ======，， 若由理想内插函数将此采样信号恢复成模拟信号，则恢复后的模拟信号

数字信号处理课程实验报告4

数字信号处理课程实验报告 实验名称FIR数字滤 班级姓名 波器设计 教师姓名实验地点实验日期 一、实验内容 1、设计一个最小阶次的低通FIR数字滤波器，性能指标为：通带0Hz~1500Hz，阻带截 止频率2000Hz，通带波动不大于1%，阻带波动不大于1%，采样频率为8000Hz； 2、用一个仿真信号来验证滤波器的正确性（注意：要满足幅度要求和线性相位特性）。 二、实验目的 1、利用学习到的数字信号处理知识解决实际问题； 2、了解线性相位滤波器的特殊结构； 3、熟悉FIR数字滤波器的设计方法。 三、涉及实验的相关情况介绍（包含使用软件或实验设备等情况） 计算机一台(安装MATLAB6.5版本或以上版本) 四、实验记录（以下1~5项必须完成，第6项为选择性试做） 1.原理基础 令希望设计的滤波器的传输函数是H（ejw，hd（n）是与其对应的单位脉冲响应。一般情况下，由Hd（ejw）求出hd（n），然后由Z变换求出滤波器的系统函数。但是通常Hd（ejw）在边界频率处有不连续点，这使得hd（n）是无限长的非因果序列，所以实际是不能实现的。为了构造一个长度为N的线性相位滤波器，可以将hd（n）截取一段来近似，并且根据线性相位的特点，需要保证截取后的序列关于（N-1）/2对称。设截取的一段为h(n)，则 Wr（n）称为矩形窗函数。 当hd（n的对称中心点取值为（N-1)/2时，就可以保证所设计的滤波器具有线性相位。 2 实验流程

1.信号的谱分析 2.信号的采样 3.信号的恢复 3源程序代码 clc; clear all; close all; fs=700;%采样频率 f=[30 40];%截止频率 a=[1 0]; dev=[0.01 0.1]; % dev纹波 [n,fo,ao,w]=remezord(f,a,dev,fs);%n滤波器阶数fo过渡带起止频率ao频带内幅度————firpmord b=remez(n,fo,ao,w);%firpm b=b.*blackman(length(b))'; b=b; a=1; figure(1) % [H,W]=freqz(b,1,1024,Fs); % plot(W,20*log10(abs(H))); freqz(b,1,1024,fs);grid title('滤波器') grid %%%%%%%%%%%%%%%% fc=28; fcl1=50; fcl2=100; fcl3=150; N=1024; n=1:N; % x=2*cos(2*pi*fc/fs*n)+j*2*sin(2*pi*fc/fs*n)+cos(2*pi*fcl/fs*n)+j*sin(2*pi*fcl/fs*n)+1*r and(1,N); xc=2*cos(2*pi*fc/fs*n); x=2*cos(2*pi*fc/fs*n)+2*cos(2*pi*fcl1/fs*n)+2*cos(2*pi*fcl2/fs*n)+0.1*rand(1,N); % x=2*cos(2*pi*fc/fs*n); xfft=abs(fft(x,N));

数字信号处理实验一

实验一 离散时间信号分析 班级 信息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）完成实验报告。 五、实验结果 六、实验总结

数字信号处理基础书后题答案中文版

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

数字信号处理实验报告(实验1_4)

实验一 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;

基于SPI Flash的Sharc系列DSP程序加载

技术创新 DSP 开发与应用 您的论文得到两院院士关注 基于SPI Flash 的Sharc 系列DSP 程序加载 The Procedures-load of Sharc Series of DSP Based on SPI Flash (四川大学) 雷春梅周群刘谋云应俊 LEI Chun-mei ZHOU Qun LIU Mou-yun YING Jun 摘要:DSP 是针对数字信号处理需要而设计的一种具有特殊结构的微处理器。它需要通过程序加载来进行信号处理,本文通过具体实例介绍了sharc 系列的DSP 通过其SPI 口的程序加载,同时也介绍了相关系统的硬件设计和软件的调试。该技术可靠性高,使用灵活方便,具有很强的实用性。关键词:程序加载;SPI 接口;数字信号处理;调试中图分类号:TP319文献标识码:B Abstract:DSP is designed to meet the need of Digital Signal Processing as a specific structuralized microprocessor,in which the signal processes through the procedures-load.In this paper,we introduce the course of loading the procedures via SPI port of DSP,which belongs to the sharc line,and also the relevant system ’s design of hardware and the debugging of software is introduced.With flexible and convenient,this method is of high reliability and highly practical.Key words:Procedures-load;SPI port;Digital Signal Processing;degug 文章编号:1008-0570(2010)04-2-0115-02 1引言 DSP 芯片主要用于数字信号处理。当前占据市场大部分份额的DSP 生产厂商是美国的德州仪器与模拟器件两家公司。在国内,尤以德州仪器生产的芯片使用更加广泛。 使用DSP 进行信号处理,就涉及到程序加载问题,不同厂家 生产的不同系列DSP 的程序加载方式也大同小异。 现在应用得比较广泛的是从flash 中加载程序。因为flash 是一种可在线编程的存储器,而且具有掉电后数据不丢失、功耗低、容量大、读取快、价格低、操作简单等优点。Flash 存储器芯片又分为并口 flash 和串口flash 。 对于并口flash,操作更为简单,数据并行传输,但占用较多的外围接口资源,因此这种程序加载方式比较适合芯片接口资源比较丰富的系统设计;而串口flash,虽然对其进行编程等操作相对复杂一些,但在接口资源比较紧张时却具有其他编程方式无可比拟的优势。 对flash 芯片的烧写,一般有两种方法:编程器烧写和在线编程。编程器烧写,即是通过硬件编程器,将工程文件经过编译链接后得到的输出文件写入flash 芯片中,这种程序烧写方式实现起来比较简单,但需要额外的硬件编程器,因而一定程度上提高了开发成本。另一方面,表贴式封装的flash 芯片不可能使用编程器烧写,只能采用在线编程。本文将通过一个AD 公司生产的sharc 系列DSP 通过其SPI 口的程序加载的具体实例来详细介绍。后面将通过硬件设计和软件架构来介绍其硬件和软件结构并总结了一些调试经验。 2硬件设计 2.1flash 芯片M25P20-VMN6简介 m25p20是美国ST 公司生产的一款flash 存储器芯片。该芯片具有与工业标准SPI 接口兼容的外部引脚。该芯片具有 2Mbit 的存储空间,即256Kbytes 。存储空间被分成四段,每段56Kbytes,并且每段存储空间由256页组成,每页大小为 256Bytes 。 对该芯片进行编程,一次性连续写入最多256Bytes 。该芯片具有一个状态寄存器,通过对其编程,可以实现芯片存储空间的保护,以及芯片写使能等操作。另外,状态寄存器中也具有 BUSY 标志位,通过该标志位,可以查询芯片当前状态。 除了读写状态寄存器以外,该存储器芯片还具有读、写、整片擦除、段擦除等命令。其器件编程的指令时序见参考文献[2]。 2.2flash 芯片与DSP 芯片的硬件连接 使用的flash M25P20芯片具有与工业标准SPI 接口兼容的外部引脚,我们使用的Sharc 系列DSP 芯片也具有同样的SPI 接口,因此,只需要将两个芯片SPI 接口的引脚对应连接即可。这就大大的节省了DSP 芯片的外部接口资源。由于Sharc 系列DSP 芯片的外设接口资源相当的丰富,而外部引脚数却极其的有限,因此,该系列DSP 芯片的引脚复用情况相当严重。 图1 Sharc 系列DSP 的外设接口分为两个部分:DAI 与DPI 。每个部分都含有各种类型的外设接口,比如,DAI 中有PCG(精确时钟产生器),DPI 中包括Timer(定时器)、UART(串口)等。 我们使用的SPI 接口就包含在DPI 中。DPI,即数字外设接口,其占有的14个外部引脚可以使用SRU 宏分配给它包括的任意一个外设接口。Sharc 系列DSP 包含两个SPI 接口,分别是DPI 中的SPI 和DAI 中的SPI B 。在我们的硬件设计中,使用DPI 中的SPI 。SPI 通信具有主从两种工作模式,在实例中,DSP 的SPI 接口工作于主模式,flash 芯片工作于从模式。对于sharc 系列DSP,根据 雷春梅:硕士研究生

数字信号处理基础书后题答案中文版

数字信号处理基础书后题答案中文版

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 幅度 频

数字信号处理实验答案完整版

数字信号处理实验答案 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),绘出四个周期。

电气工程及其自动化考研总况

电气工程及其自动化考 研总况 IMB standardization office【IMB 5AB- IMBK 08- IMB 2C】

电气工程及其自动化考研总况 一、全国电气工程及其自动化专业学校排名 1.清华大学 2.西安交通大学 3.华中科技大学 4.浙江大学 5.重庆大学 6.天津大学 7.哈尔滨工业大学 8.上海交通大学 9.华北电力大学10.东南大学11.西南交通大学12.沈阳工业大学13.中国矿业大学14.华南理工大学15.南京航空航天大学16.北京交通大学17.武汉大学18.哈尔滨理工大学19.四川大学20.河海大学21.哈尔滨工程大学22.郑州大学23.广西大学24.陕西科技大学 二，电气工程与自动化专业 （1）业务培养目标： 业务培养目标：本专业培养在工业与电气工程有关的运动控制、工业过程控制、电气工程、电力电子技术、检测与自动化仪表、电子与计算机技术等领域从事工程设计、系统分析、系统运行、研制开发、经济管理等方面的高级工程技术人才。 业务培养要求：本专业学生主要学习电工技术、电子技术、自动控制理论、信息处理、计算机技术与应用等较宽广领域的工程技术基础和一定的专业知识。学生受到电工电子、信息控制及计算机技术方面的基本训练，具有工业过程控制与分析，解决强弱电并举的宽口径专业的技术问题的能力。

（2）主干课程： 主干学科：电气工程、控制科学与工程、计算机科学与技术 主要课程：电路原理、电子技术基础、计算机原理及应用、计算机软件基础、控制理论、电机与拖动、电力电子技术、信号分析与处理、电力拖动控制系统、工业过程控制与自动化仪表等。高年级可根据社会需要设置柔性的专业方向模块课及选修课。 主要实践性教学环节：包括电路与电子基础实验、电子工艺实习、金工实习、专业综合实验、计算机上机实践、课程设计、生产实习、毕业设计。 主要实验：运动控制实验、自动控制实验、计算机控制实验、检测仪表实验、电力电子实验等 （3）修业年限： 四年 （4）授予学位： 工学学士 （5）相近专业： 微电子学自动化电子信息工程通信工程计算机科学与技术电子科学与技术生物医学工程电气工程与自动化信息工程信息科学技术软件工程