实验6 基于MATLAB 的数字滤波器设计
实验目的：加深对数字滤波器的常用指标和设计过程的理解。
实验原理：低通滤波器的常用指标：
P
P j P for e G ωωδδω
≤+≤≤-,1)(1
π
ωωδω
≤≤≤S S j for e
G ,)(
通带边缘频率：
P
ω，阻带边缘频率：S
ω ，
通带起伏：P δ，通带峰值起伏：
]
)[1(log
2010
dB p p
δα
--=，阻带起伏：s
δ，
最小阻带衰减：
]
)[(log 2010dB s S δα-=。
数字滤波器有IIR 和FIR 两种类型，它们的特点和设计方法不同。 在MATLAB 中，可以用[b ，a]=butter （N,Wn ）等函数辅助设计IIR 数字滤波器, 也可以用b=fir1(N,Wn,’ftype’) 等函数辅助设计FIR 数字滤波器。
实验内容： 利用MATLAB 编程设计一个数字带通滤波器，指标要求如下：
通带边缘频率：10.4P ωπ
=，20.6P ωπ
=，通带峰值起伏：
]
[1dB p
≤α
。 阻带边缘频率：
π
ω3.01=S ，
20.7S ωπ
=，最小阻带衰减：
]
[40dB S ≥α。
分别用IIR 和FIR 两种数字滤波器类型进行设计。
实验要求：给出IIR 数字滤波器参数和FIR 数字滤波器的冲激响应，绘出它们
的幅度和相位频响曲线，讨论它们各自的实现形式和特点。
实验步骤：
①Butterworth 滤波器的设计（IIR ）
>>wp=[0.4*pi,0.6*pi]; >>wr=[0.3*pi,0.7*pi]; >>Ap=1; >>Ar=40;
S
ω-P
ω-P S
ωPassband Stopband
Transition band
Fig 7.1 Typical magnitude
specification for a digital LPF
>>[N,Wn]=buttord(wp/pi,wr/pi,Ap,Ar)
N =
7
Wn =
0.3854 0.6146
>> [b,a]=butter(N,Wn,'bandpass')
b =
Columns 1 through 12
0.0002 0 -0.0014 0 0.0042 0 -0.0071 0 0.0071 0 -0.0042 0
Columns 13 through 15
0.0014 0 -0.0002
a =
Columns 1 through 12
1.0000 0.0000 3.7738 0.0000 6.5614 0.0000 6.6518 0.0000 4.2030 0.0000 1.6437 0.0000
Columns 13 through 15
0.3666 0.0000 0.0359
>> [H,w]=freqz(b,a);
>>mag=abs(H);
>>plot(w/pi,mag);
>>xlabel('角频率(\Omega)');
>>ylabel('幅度|Ha(j\Omega)|');
>> title('数字butterworth带通滤波器幅度响应|Ha(j\Omega)|');
>>phase=angle(H);
>>plot(w/pi,phase);
>>xlabel('角频率(\Omega)');
>>ylabel('相位');
>> title('数字butterworth带通滤波器相位响应曲线');
>>
②FIR滤波器的设计：
>> wp1 = 0.4*pi; wp2 = 0.6*pi;
>> ws1 = 0.3*pi; ws2 = 0.7*pi;
>>tr_width = min((wp1-ws1),(ws2-wp2)) tr_width =
0.3142
>> M = ceil(6.2*pi/tr_width) + 1
M =
63
>> n=[0:1:M-1];
>> wc1 = (ws1+wp1)/2; wc2 = (wp2+ws2)/2; >>wc=[wc1/pi,wc2/pi];
>>window= hanning(M);
>> [h1,w]=freqz(window,1);
>>figure(1);
>>subplot(2,1,1)
>>stem(window);
>>axis([0 60 0 1.2]);
>>grid;
>>xlabel('n');
>> title('Hanning窗函数');
>>subplot(2,1,2)
>>plot(w/pi,20*log(abs(h1)/abs(h1(1)))); >>axis([0 1 -350 0]);
>>grid;
>>xlabel('w/\pi');
>>ylabel('幅度(dB)');
>> title('Hanning窗函数的频谱');
>>hn = fir1(M-1,wc, hanning (M));
>> [h2,w]=freqz(hn,1,512);
>>figure(2);
>>subplot(2,1,1)
>>stem(n,hn);
>>axis([0 60 -0.25 0.25]);
>>grid;
>>xlabel('n');
>>ylabel('h(n)');
>> title('Hanning窗函数的单位脉冲响应'); >>subplot(2,1,2)
>>plot(w/pi,20*log(abs(h2)/abs(h2(1)))); >>grid;
>>xlabel('w/\pi');
>>ylabel('幅度(dB)');
>>figure(3);
>>phase=angle(h1);
>>plot(phase);
>>axis([1 pi -1 0]);
>>xlabel('w/\pi');
>>ylabel('线性相位');
>> title('Hanning窗函数相位特性曲线'); >>
实现形式及特点分析：
1.在本例中，相同的技术指标下，IIR滤波器实现的阶数为N=7，
而FIR滤波器的阶数N=63。因此，相同的技术指标，用IIR滤波器实现的阶数远远小于用FIR滤波器的阶数。这是由于IIR 滤波器存在着输出对输入的反馈。
2.从实验中绘制的相位特性曲线可以看出，FIR滤波器可以得到
严格的线性相位，而IIR滤波器做不到这一点。IIR滤波器的选择性愈好，其相位的非线性愈严重。因而，如果IIR滤波器要得到线性相位，又要满足滤波的技术要求，必须加全通网络进行相位校正，这同样会大大增加滤波器的阶数。因此，从相位特性考虑，FIR滤波器又优于IIR滤波器。
第一章绪论 1.1 本论文的背景 随着信息技术的飞速发展，数字信号处理技术已经发展成为一门关键的技术学科，而DSP芯片的出现则为数字信号处理算法的实现提供了可能，这一方面促进了数字信号处理技术的进一步发展,也使数字信号处理的应用领域得到了极大的拓展。在近20年里，DSP芯片已经在通信和家用电器等领域得到了广泛的应用。 1.1.1 数字信号处理器的发展状况 DSP(Digital Signal Processing)也称数字信号处理器，是一种具有特殊结构的微处理器，是建立在数字信号处理的各种理论和算法基础上，专门完成各种实时数字信息处理的芯片。与单片机相比，DSP有着更适合数字信号处理的优点。芯片部采用程序和数据分开的哈佛结构，具有专门的硬件乘法器，广泛采用流水线操作，具有良好的并行特性，提供特殊的DSP指令，可以快速地实现各种数字信号处理算法[1]。 DSP发展历程大致分为三个阶段：70年代理论先行，80年代产品普及，90年代突飞猛进。在DSP出现之前数字信号处理主要依靠MPU(微处理器)来完成。但MPU 较低的处理速度无法满足高速实时的要求。因此，直到70年代才提出了DSP的理论和算法基础。随着大规模集成电路技术的发展，1982年世界上诞生了首枚通用可编程DSP芯片TI的TMS32010。DSP芯片的问世是个里程碑，它标志着DSP应用系统由大型系统向小型化迈进了一大步。进入80年代后期，随着数字信号处理技术应用围的扩大，要求提高处理速度，到1988年出现了浮点DSP，同时提供了高级语言的编译器，使运算速度进一步提高，其应用围逐步扩大到通信、计算机领域。90年代相继出现了第四代和第五代DSP器件。以DSP作为主要元件，再加上外围设备和特定功能单元综合成的单一芯片，加速了DSP解决方案的发展，同时产品价格降低，运算速度和集成度大幅提高[2]。 进入21世纪，现在DSP向着高速，高系统集成，高性能方向发展。当前的DSP 多数基于RISC(精简指令集计算机)结构，且进入了VLSI(超大规模集成电路)阶段。如TI公司的TMS320C80代表了新一代芯片集成技术，它将4个32位的DSP，1个32位RISC主处理器，1个传输控制器，2个视频控制器和50Kb SRAM集成在一个芯片上。这样的芯片通常称之为MVP(多媒体视频处理器)。它可支持各种图像规格和各种算法，功能相当强。而第六代TMSC6000系列则是目前速度最快，性能最高的DSP芯片，该系列芯片的发展蓝图中有高至5000MIPS，3G FLOPS的处理性能。
长安大学 数字信号处理综合设 计 专业_______电子信息工程_______ 班级__24030602___________ 姓名_______张舒_______ 学号2403060203 指导教师陈玲 日期_______2008-12-27________
一、课程设计目的： 1. 进一步理解数字信号处理的基本概念、基本理论和基本方法； 2.熟悉在Windows环境下语音信号采集的方法； 3.学会用MATLAB软件对信号进行分析和处理； 4.综合运用数字信号处理理论知识，掌握用MATLAB软件设计FIR和IIR数字滤波器的方法； 5. 提高依据所学知识及查阅的课外资料来分析问题解决问题的能力。 二、课程设计内容： 1.语音信号的采集 利用windows下的录音机录制一段自己的话音，时间控制在1秒左右；并对语音信号进行采样，理解采样频率、采样位数等概念。 2.语音信号的频谱分析 利用函数fft对采样后语音信号进行快速傅里叶变换，得到信号的频谱特性。 3.设计数字滤波器 采用窗函数法和双线性变换法设计下列要求的三种滤波器，根据语音信号的特点给出有关滤波器的性能指标： 1）低通滤波器性能指标，fp=1000Hz，fc=1200Hz， As=100dB，Ap=1dB； 2）高通滤波器性能指标，fc=4800Hz，fp=5000Hz ，As=100dB，Ap=1dB； 3）带通滤波器性能指标，fp1=1200Hz，fp2=3000Hz，fc1=1000Hz，fc2=3200Hz，As=100dB，Ap=1dB。 4.对语音信号进行滤波 比较用两种方法设计的各滤波器的性能，然后用性能好的滤波器分别对采集的语音信号进行滤波；并比较滤波前后语音信号的波形及频谱，分析信号的变化。 5.回放语音信号，感觉滤波前后的声音变化。 三、实验原理 （一）基于双线性Z变换法的IIR数字滤波器设计 由于的频率映射关系是根据推导的，所以使jΩ轴每隔2π/Ts便映射到单位圆上一周，利用冲激响应不变法设计数字滤波器时可能会导致上述的频域混叠现象。为了克服这一问题，需要找到由s平面到z平面的另外的映射关系，这种关系应保证： 1) s平面的整个jΩ轴仅映射为z平面单位圆上的一周； 2) 若G(s)是稳定的，由G(s)映射得到的H(z)也应该是稳定的； 3) 这种映射是可逆的，既能由G(s)得到H(z)，也能由H(z)得到G(s)； 4) 如果G(j0)=1，那么。 双线性Z变换满足以上4个条件的映射关系，其变换公式为
一、综合实验内容和目的 1、实验目的 (1) 通过实验学习掌握TMS320F28335的浮点处理； (2) 学习并掌握A/D模块的使用方法； (3) 学习并掌握中断方式和查询方式的相关知识及其相互之间的转换； (4) 学习信号时域分析的方法，了解相关电量参数的计算方法； (5) 了解数字滤波的一些基本方法。 2、实验内容 要求1：对给定的波形信号，采用TMS320F28335的浮点功能计算该信号的以下时域参数：信号的周期T，信号的均方根大小V rms、平均值V avg、峰-峰值V pp。 其中，均方根V rms的计算公式如下： V= rms 式中N为采样点数，()u i为采样序列中的第i个采样点。 要求2：所设计软件需要计算采样的波形周期个数，并控制采样点数大于1个波形周期，且小于3个波形周期大小。 要求3：对采集的数据需要加一定的数字滤波。 二、硬件电路 相关硬件：TMS320F28335DSP实验箱，仿真器。
硬件结构图 三、程序流程图 1、主程序流程图 程序的主流程图2、子程序流程图
参数计算的流程图 四、实验结果和分析 1、实验过程分析 (1) 使用的函数原型声明 对ADC模件相关参数进行定义：ADC时钟预定标，使外设时钟HSPCLK 为25MHz，ADC模块时钟为12.5MHz，采样保持周期为16个ADC时钟。 (2) 定义全局变量 根据程序需要，定义相关变量。主要有：ConversionCount、Voltage[1024]、Voltage1[1024]、Voltage2[1024]、filter_buf[N]、filter_i、Max、Min、T、temp、temp1、temp2、temp3、Num、V、Vav、Vpp、Vrm、fre。这些变量的声明请见报告后所附的源程序。 (3) 编写主函数 完成系统寄存器及GPIO初始化；清除所有中断，初始化PIE向量表，将程
目录 摘 要 (2) Abstract (3) 1 绪论 (4) 1.1概述 (4) 1.2选题的目的、意义 (4) 1.3 选题的背景 (5) 1.4 本文所研究的内容 (6) 2 波形信号发生器的原理及方案选择 (7) 2.1任意波形信号发生器的原理 (7) 2.1.1 直接模拟法 (7) 2.1.2 直接数字法 (7) 2.2 任意波形发生器的设计方案 (9) 2.2.1 查表法 (9) 2.2.2计算法 (9) 2.2.3传统方法 (10) 3 基于DSP 5416的任意波形信号发生器的软件设计 (12) 3.1 TMS320C5416的开发流程 (12) 3.2软件开发环境 (13) 3.3任意波形信号发生器的软件编程 (14) 3.3.1 计算法实现波形输出 (14) 3.3.2 D/A转换 (15) 3.3.3波形控制及软件设计流程图 (16) 3.4参数的设定 (18) 4 基于DSP 5416的任意波形信号发生器的硬件设计 (20) 4.1 TMS320VC5416开发板 (20) 4.2 TMS320VC5416实验箱的连接 (23) 4.3 波形信号发生器的硬件测试过程 (23) 5 任意波形信号发生器展望 (28) 结束语 (29) 致谢 (30) 参考文献 (31)
摘 要 任意波形发生器是信号源的一种，它是具有信号源所具有的特点，更因它高的性能优势而倍受人们青睐。信号源主要给被测电路提供所需要的已知信号（各种波形），然后用其它仪表测量感兴趣的参数。可见信号源在各种实验应用和试验测试处理中，它不是测量仪器，而是根据使用者的要求，作为激励源，仿真各种测试信号，提供给被测电路，以满足测量或各种实际需要。 随着无线电应用领域的扩展，针对广播、电视、雷达、通信的专用信号发生器获得了长足的发展，表现在载波调制方式的多样化，从调幅、调频、调相到脉冲调制。如果采用多台信号发生器获得测量信号显然是很不方便的。因此需要任意波形发生器（Arbitrary Waveform Generator，AWG），使其能够产生任意频率的载频信号和多种载波调制信号。 目前我国已经开始研制任意波形发生器，并取得了可喜的成果。但总的来说，我国任意波形发生器还没有形成真正的产业。并且我国目前在任意波形发生器的种类和性能都与国外同类产品存在较大的差距，因此加紧对这类产品的研制显得迫在眉睫。 本文主要工作分为以下几个方面：首先，介绍研制任意波形信号发生器的目的、意义、背景，以及利用CCS仿真工具用软件实现任意波形信号发生器的的过程 ；之后，对硬件的连接及测试结果作介绍；最后，简要的对任意波形信号发生器的未来作一下展望。 关键词：DSP，任意波形信号发生器,DDS
数字信号处理作业提交日期：2016年7月15日
实验一 维纳滤波器的设计 第一部分 设计一维纳滤波器。 (1)产生三组观测数据，首先根据()(1)()s n as n w n =-+产生信号()s n ，将其加噪(信噪比分别为20,10,6dB dB dB )，得到观测数据123(),(),()x n x n x n 。 (2)估计()i x n ，1,2,3i =的AR 模型参数。假设信号长度为L ，AR 模型阶数为N ，分析实验结果，并讨论改变L ，N 对实验结果的影响。 1 实验原理 滤波技术是信号分析、处理技术的重要分支，无论是信号的获取、传输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传递是至关重要的。信号分析检测与处理的一个十分重要的内容就是从噪声中提取信号，实现这种功能的有效手段之一是设计一种具有最佳线性过滤特性的滤波器，当伴有噪声的信号通过这种滤波器的时候，它可以将信号尽可能精确地重现或对信号做出尽可能精确的估计，而对所伴随噪声进行最大限度地抑制。维纳滤波器就是这种滤波器的典型代表之一。 维纳（Wiener ）是用来解决从噪声中提取信号的一种过滤（或滤波）方法。这种线性滤波问题，可以看做是一种估计问题或一种线性估计问题。 设一线性系统的单位样本响应为()h n ，当输入以随机信号()x n ，且 ()() () x n s n v n =+，其中()s n 表示原始信号，即期望信号。()v n 表示噪声，则输出()y n 为()=()()m y n h m x n m -∑，我们希望信号()x n 经过线性系统()h n 后得到的()y n 尽可能接近 于()s n ，因此称()y n 为估计值，用?()s n 表示。 则维纳滤波器的输入-输出关系可用下面表示。 设误差信号为()e n ,则?()()()e n s n s n =-,显然)(n e 可能是正值，也可能是负值，并且它是一个随机变量。因此，用它的均方误差来表达误差是合理的，所谓均方误差最小即 它的平方的统计期望最小：222?[|()|][|()()|][|()()|]E e n E s n s n E s n y n =-=-=min 。而要使均方误差最小，则需要满足2[|()|]j E e n h ?=0. 进一步导出维纳-霍夫方程为：()()()()*(),0,1,2...xs xx xx i R m h i R m i R m h m m =-==∑ 写成矩阵形式为：xs xx R R h =，可知：1xs xx h R R -=。表明已知期望信号与观测数据的互相关函数以及观测信号的自相关函数时，可以通过矩阵求逆运算，得到维纳滤波器的
摘要 (3) Abstract (4) 第一章绪论 (5) 1.1 选题背景及研究意义 (5) 1.2 国内外发展现状及发展趋势 (6) 1.2.1 电动汽车发展现状及趋势 (6) 1.2 .2 液晶显示技术的发展及其应用前景 (7) 1.2.3数字信号处理器的发展及其应用前景 (8) 1.3本设计研究的主要内容 (9) 第二章系统设计方案 (10) 2.1 DSP软件开发工具简介 (10) 2.1.1 TMS320F2812 (10) 2.2 系统设计概述 (11) 2.3电动汽车几种传感器及其作用 (12) 第三章液晶与液晶显示 (14) 3.1 液晶与液晶显示器件 (14) 3.2 液晶显示器件的基本结构 (15) 3.3 典型的液晶显示器件 (16) 3.3.1 静态驱动段型液晶显示器件 (16) 3.3.2 动态驱动点矩阵型液晶显示器件 (17) 3.4.1 AXG19264 引脚介绍 (19) 3.4.2 图形液晶显示行驱动控制器HD61203U (19) 3.4.3 图形液晶显示列驱动控制器HD61202U (19) 3.4.4 HD61202U 的指令系统 (21) 第四章电动汽车液晶显示系统硬件设计 (24) 4.1 硬件设计分析 (24) 4.1.1 处理器直接访问方式 (24) 4.1.2 处理器间接访问方式 (24) 4.3 硬件设备的调试 (27) 4.4 DSP2812功能模块图 (28) 4.5 DSP与液晶模块硬件接口设计 (29) 4.6 保护电路 (29) 第五章电动汽车液晶显示系统软件设计 (31) 5.1 DSP软件系统开发环境介绍 (31) 2.1.2 CCS的组成 (31) 5.1.2 CCS环境下project的主要成员 (32) 5.2主程序的软件流程图 (34) 5.3软件调试的方法 (34) 第六章总结和展望 (36) 6.1 工作总结 (36) 6.2 后续工作展望 (36) 参考文献 (37) 附录 (38)
淮北煤炭师范学院 2009届学士学位论文 基于MA TLAB的数字滤波器设计 学院、专业物理与电子信息学院 电子信息科学与技术 研究方向基于MATLAB的数字滤波器设计 学生姓名耿博 学号200513432024 指导教师姓名邹锋 指导教师职称讲师 2009 年4 月18
基于MATLAB的数字滤波器设计 耿博 （淮北煤炭师范学院物理与电子信息学院235000） 摘要随着信息时代和数字世界的到来，数字信号处理已成为今一门极其重要的学科和技术领域。数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理应用中，数字滤波器十分重要并已获得广泛应用。 数字滤波是数字信号处理的重要内容，数字滤波器可分为IIR和FIR两大类。对于IIR数字滤波器的设计，需要借助模拟原型滤波器，再将模拟滤波器转化为数字滤波器，文中采用的设计方法是脉冲响应不变法、双向性变换法和完全函数设计法；对于FIR数字滤波器的设计，可以根据所给定的频率特性直接设计，文中采用的设计方法是窗函数法。本文根据IIR滤波器和FIR滤波器的特点，在MATLAB坏境下分别用双线性变换法设计IIR和用窗函数设计FIR数字滤波器，并对采集的语音信号进行分析，最后给出了IIR和FIR对语音滤波的效果。 关键词数字滤波器；IIR ；FIR ；MATLAB
The Design of Digital Filter based on MATLAB Geng Bo School of Physics and Electronics Information, Huaibei Coal Industry Teachers? College, 235000 ABSTRACT Along with the information age and the digital world arrival, the digital signal processing has become a now extremely important discipline and the area of technology.The digital signal processing in the correspondence, the multitudinous domains the pronunciation such as the image, the automatic control, the radar, the military, the aerospace, the medical service and the domestic electric appliances and so on have obtained the widespread application.In the digital signal processing application, the digital filter are extremely important and have obtained the widespread application. The digital filter are the digital signal processing important content, the digital filter may divide into IIR and the FIR two main kinds. As for the IIR digital filter design, we need the help of analog prototype filter, and then transform analog filter into digital filter. In the paper we use the design of the pulse response invariable method, the bilinear method and full function design; as for the FIR filter, we can design it directly based on the giving frequency, in the paper it uses the design of the window function.This article according to the IIR filter and the FIR filter characteristic, uses the bilinearity method of transformation under the MATLAB bad boundary to design IIR and to design the FIR numeral filter separately with the window box number, and carries on the analysis to the gathering pronunciation signal, and finally gives IIR and FIR to the pronunciation filter effect. Keywords Digtial Filter;IIR;FIR;MATLAB
东南大学自动化学院 实验报告 课程名称： DSP技术及课程设计 实验名称：直流无刷电机控制综合实验 院（系）：自动化专业：自动化 姓名：ssb 学号：08011 实验室：304 实验组别： 同组人员：ssb1 ssb2 实验时间：2014年 6 月 5 日评定成绩：审阅教师：
目录 1.实验目的和要求 (3) 1.1 实验目的 (3) 1.2 实验要求 (3) 1.2.1 基本功能 (3) 1.2.2 提高功能 (3) 2.实验设备与器材配置 (3) 3.实验原理 (3) 3.1 直流无刷电动机 (3) 3.2 电机驱动与控制 (5) 3.3 中断模块 (7) 3.3.1 通用定时器介绍及其控制方法 (7) 3.3.2 中断响应过程 (7) 3.4 AD模块 (8) 3.4.1 TMS320F28335A 芯片自带模数转换模块特性 (8) 3.4.2 模数模块介绍 (8) 3.4.3 模数转换的程序控制 (8) 4.实验方案与实验步骤 (8) 4.1 准备实验1：霍尔传感器捕获 (8) 4.1.1 实验目的 (8) 4.1.2 实验内容 (9) 4.1.2.1 准备 (9) 4.1.2.2 霍尔传感器捕获 (9) 4.2 准备实验2：直流无刷电机(BLDC)控制 (10) 4.2.1 程序框架原理 (10) 4.2.1.1 理解程序框架 (10) 4.2.1.2 基于drvlib281x库的PWM波形产生 (11) 4.2.2 根据捕获状态驱动电机运转 (12) 4.2.2.1 目的 (12) 4.2.2.2 分析 (12) 4.3 考核实验：直流无刷电机调速控制系统 (13) 4.3.1 初始化工作 (13) 4.3.2 初始化定时器0.... . (13) 4.3.3初始化IO口 (13) 4.3.4中断模块.... (13) 4.3.5 AD模块 (14) 4.3.6在液晶屏显示 (15) 4.3.7电机控制 (17) 4.3.7.1 控制速度方式选择 (17) 4.3.7.2 控制速度和转向 (18) 4.3.8延时子函数 (19) 4.3.9闭环PID调速 (19)
实验一MATLAB语言的基本使用方法 实验类别：基础性实验 实验目的： （1）了解MATLAB程序设计语言的基本方法，熟悉MATLAB软件运行环境。 （2）掌握创建、保存、打开m文件的方法，掌握设置文件路径的方法。 （3）掌握变量、函数等有关概念，具备初步的将一般数学问题转化为对应计算机模型并进行处理的能力。 （4）掌握二维平面图形的绘制方法，能够使用这些方法进行常用的数据可视化处理。 实验内容和步骤： 1、打开MATLAB，熟悉MATLAB环境。 2、在命令窗口中分别产生3*3全零矩阵，单位矩阵，全1矩阵。 3、学习m文件的建立、保存、打开、运行方法。 4、设有一模拟信号f(t)=1.5sin60πt,取?t=0.001,n=0,1,2,…,N-1进行抽样，得到 序列f(n),编写一个m文件sy1_1.m,分别用stem,plot,subplot等命令绘制32 点序列f(n)(N=32)的图形，给图形加入标注，图注，图例。 5、学习如何利用MATLAB帮助信息。 实验结果及分析： 1）全零矩阵 >> A=zeros(3,3) A = 0 0 0 0 0 0 0 0 0 2）单位矩阵 >> B=eye(3) B = 1 0 0 0 1 0 0 0 1 3）全1矩阵 >> C=ones(3) C = 1 1 1 1 1 1 1 1 1 4)sy1_1.m N=32; n=0:N-1; dt=0.001; t=n*dt; y=1.5*sin(60*pi*t); subplot(2,1,1), plot(t,y); xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('正弦函数'); title('二维图形'); subplot(2,1,2), stem(t,y) xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('序列函数'); title('条状图形'); 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 二维图形 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 条状图形
本科毕业设计论文基于DSP的谱分析仪设计 学生姓名： 班级：电自0913 学号： 指导教师： 所在单位：电气工程学院
答辩日期：2013年6月24日
摘要 随着计算机和微电子技术的飞速发展，基于数字信号处理的频谱分析已经应用到各个领域并发挥着重要作用。频谱分析仪对于信号分析来说是必不可少的，它可以利用频率对信号进行分析。频谱分析仪可应用于诸多领域，如通讯发射机以及干扰信号的测量，频谱的监测，器件的特性分析等，但各行各业对其性能要求也不尽相同。 本课题主要做了以下工作：首先，本文介绍了频谱分析仪的作用、课题背景、现状及发展趋势；然后，设计了以TI公司的定点数字信号处理器（DSP）TMS320VC5402为CPU的开发系统，包括复位电路、时钟电路、存储器扩展、电源模块、AD采样、DA单元、JTAG等的设计；由于CPU采用FFT算法，所以详细介绍了FFT的原理以及其在TMS320VC5402上的实现。最后，简要介绍了用于开发DSP的集成开发环境CCS。 关键词：TMS320VC5402；频谱分析；FFT；功率谱
Abstract With the rapid development of computer and microelectronics technology, spectrum analysis based on digital signal processing (DSP) has been applied to various fields and play an important role. A spectrum analyzer for the signal analysis is indispensable, it can make use of frequency analysis of signals. A spectrum analyzer can be applied to many fields, such as communication transmitter and the interfering signal measurement, spectrum monitoring, device characteristics analysis and so on, but in all walks of life to its performance requirements are also different. This topic mainly done the following work: first of all, this paper introduces the role of a spectrum analyzer, topic background, present situation and development trend; Then, designed by TI company's fixed-point digital signal processor (DSP) TMS320VC5402 as CPU development system, including the reset circuit, clock circuit, memory expansion, a power supply module, AD sampling, DA units, such as JTAG design; Due to the CPU adopts FFT algorithm, so the principle of FFT is introduced and its implementation on TMS320VC5402. Finally, this paper briefly introduces the integrated development environment CCS for the development of DSP. Keywords：TMS320VC5402；Spectrum analyzer；FFT；Power spectrum
数字滤波器的MAT LAB设计与 DSP上的实现 数字滤波器的MATLAB 设计与DSP上的实现 公文易文秘资源网佚名2007-11-15 11:56:42我要投稿添加到百度搜藏 摘要：以窗函数法设计线性相位FIR数字滤波器为例，介绍用MATLAB工具软件设计数字滤波器的方法和在定点DSP上的实现。实现时，先在CCS5000仿真开发，然后将程序加载到TMS320VC5409评估板上实时运行，结果实现了目标要求。文中还讨论了定标、误差、循环寻址等在DSP上实现的关键问题。关键词 摘要：以窗函数法设计线性相位 FIR数字滤波器为例，介绍用 MATLAB工具软件设计数字滤波器的方法和在定点DSP上的实现。实现时，先在 CCS5000仿真开发，然后将程序加载到 TMS320VC5 409评估板上实时运行，结果实现了目标要求。文中还讨论了定标、误差、循环寻址等在DSP上实 现的关键问题。 关键词：数字滤波器MATLAB DSP 引言 随着信息时代和数字世界的到来，数字信号处理已成为今一门极其重要的学科和技术领域。数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应
用。在数字信号处理应用中，数字滤波器十分重要并已获得广泛应用。 1数字滤波器的设计 1.1数字滤波器设计的基本步骤 数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应（IIR ）滤波器和有限长冲激响应（FIR ）滤波器。IIR滤波器的特征是，具有无限持续时间冲激响应。种滤波器一般需要用递归模型来实现，因而有时也称之为递归滤波器。FIR滤波器的冲激响应只能延续一定时间， 在工程实际中可以采用递归的方式实现，也可以采用非递归的方式实现。数字滤波器的设计方法有多种，如双线性变换法、窗函数设计法、插值逼近法和Chebyshev逼近法等等。随着 MATLAB软件尤 其是MATLAB的信号处理工作箱的不断完善，不仅数字滤波器的计算机辅助设计有了可能，而且还可以使设计达到最优化。 数字滤波器设计的基本步骤如下： (1确定指标 在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中，数字滤波器常常被用来实现选频操作。因此，指标的形式一般在频域中给岀幅度和相位响应。幅度指标主要以两种方式给岀。第一种是绝对指标。它提供对幅度响应函数的要求，一般应用于FI R滤波器的设计。第二种指标是相对指标。它以分贝值的形式给岀要求。在工程实际中，这种指标最受欢迎。对于相位响应指标形式，通常希望系统在通频带中人有线性相位。运用线性相位响应指标进行滤波器设计具有如下优点：①只包含实数算法，不涉及复数运算；②不存在延迟失真，只有固定数量的延迟；③长度为N的滤波器(阶数为N-1)，计算量为N/2数量级。因此，本文中滤波器的设计就以线性相位FIR滤波器的设计为例。 (2)逼近
实验0 实验设备安装才CCS调试环境 实验目的： 按照实验讲义操作步骤，打开CCS软件，熟悉软件工作环境，了解整个工作环境内容，有助于提高以后实验的操作性和正确性。 实验步骤： 以演示实验一为例： 1．使用配送的并口电缆线连接好计算机并口与实验箱并口，打开实验箱电源； 2．启动CCS，点击主菜单“Project->Open”在目录“C5000QuickStart\sinewave\”下打开工程文件sinewave.pjt，然后点击主菜单“Project->Build”编译，然后点击主菜单“File->Load Program”装载debug目录下的程序sinewave.out； 3．打开源文件exer3.asm，在注释行“set breakpoint in CCS !!!”语句的NOP处单击右键弹出菜单，选择“Toggle breakpoint”加入红色的断点，如下图所示； 4．点击主菜单“View->Graph->Time/Frequency…”，屏幕会出现图形窗口设置对话框 5．双击Start Address，将其改为y0；双击Acquisition Buffer Size，将其改为1； DSP Data Type设置成16-bit signed integer，如下图所示； 6．点击主菜单“Windows->Tile Horizontally”，排列好窗口，便于观察 7．点击主菜单“Debug->Animate”或按F12键动画运行程序，即可观察到实验结果： 心得体会： 通过对演示实验的练习，让自己更进一步对CCS软件的运行环境、编译过程、装载过程、属性设置、动画演示、实验结果的观察有一个醒目的了解和熟悉的操作方法。熟悉了DSP实验箱基本模块。让我对DSP课程产生了浓厚的学习兴趣，课程学习和实验操作结合为一体的学习体系，使我更好的领悟到DSP课程的实用性和趣味性。
武汉工程大学 数字信号处理实验报告 姓名：周权 学号：1204140228 班级：通信工程02
一、实验设备 计算机，MATLAB语言环境。 二、实验基础理论 1.序列的相关概念 2.常见序列 3.序列的基本运算 4.离散傅里叶变换的相关概念 5.Z变换的相关概念 三、实验内容与步骤 1.离散时间信号（序列）的产生 利用MATLAB语言编程产生和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形表示。 四实验目的 认识常用的各种信号，理解其数字表达式和波形表示，掌握在计算机中生成及绘制数字信号波形的方法，掌握序列的简单运算及计算机实现与作用，理解离散时间傅里叶变换，Z变换及它们的性质和信号的频域分
实验一离散时间信号（序列）的产生 代码一 单位样值 x=2; y=1; stem(x,y); title('单位样值 ') 单位阶跃序列 n0=0; n1=-10; n2=10; n=[n1:n2]; x=[(n-n0)>=0]; stem(n,x); xlabel('n'); ylabel('x{n}'); title('单位阶跃序列');
实指数序列 n=[0:10]; x=(0.5).^n; stem(n,x); xlabel('n'); ylabel('x{n}'); title('实指数序列');
正弦序列 n=[-100:100]; x=2*sin(0.05*pi*n); stem(n,x); xlabel('n'); ylabel('x{n}'); title('正弦序列');
随机序列 n=[1:10]; x=rand(1,10); subplot(221); stem(n,x); xlabel('n'); ylabel('x{n}'); title('随机序列');
深圳大学实验报告课程名称：DSP系统设计 实验项目名称：DSP系统设计实验 学院：机电与控制工程学院 专业：自动化 指导教师：杜建铭 报告人1：. 学号：。班级：3 报告人2：. 学号：。班级：3 报告人3：. 学号：。班级：3 实验时间： 实验报告提交时间： 教务处制
实验一、CCS入门试验 一、实验目的 1. 熟悉CCS集成开发环境，掌握工程的生成方法； 2. 熟悉SEED-DEC2812实验环境； 3. 掌握CCS集成开发环境的调试方法。 二、实验仪器 1．TMS320系列SEED-DTK教学试验箱24套 2. 台式PC机24台 三、实验内容 1．仿真器驱动的安装和配置 2. DSP 源文件的建立； 3. DSP程序工程文件的建立； 4. 学习使用CCS集成开发工具的调试工具。 四、实验准备： 1.将DSP仿真器与计算机连接好； 2.将DSP仿真器的JTAG插头与SEED-DEC2812单元的J1相连接； 3.启动计算机,当计算机启动后，打开SEED-DTK2812的电 源。SEED-DTK_MBoard单元的＋5V，＋3.3V，＋15V，－15V的电源指示灯及SEED-DEC2812的电源指示灯D2是否均亮；若有不亮，请断开电源，检查电源。 五、实验步骤 （一）创建源文件 1.进入CCS环境。
2.打开CCS选择File →New →Source File命令 3.编写源代码并保存 4.保存源程序名为math.c，选择File →Save 5.创建其他源程序（如.cmd）可重复上述步骤。 （二）创建工程文件 1.打开CCS，点击Project-->New,创建一个新工程，其中工程名及路径可任意指定弹 出对话框： 2.在Project中填入工程名，Location中输入工程路径；其余按照默认选项，点击完成 即可完成工程创建； 3.点击Project选择add files to project，添加工程所需文件；
实验一 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;
电子与信息工程学院综合实验课程报告 课题名称基于DSP的温度采集系统
毕业设计（论文）原创性声明和使用授权说明 原创性声明 本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。 作者签名：日期： 指导教师签名：日期： 使用授权说明 本人完全了解大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。 作者签名：日期：
学位论文原创性声明 本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。 作者签名：日期：年月日 学位论文版权使用授权书 本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。 涉密论文按学校规定处理。 作者签名：日期：年月日 导师签名：日期：年月日
%要求设计一butterworth低通数字滤波器，wp=30hz，ws=40hz，rp=0.5，rs=40，fs=100hz。>>wp=30;ws=40;rp=0.5;rs=40;fs=100; >>wp=30*2*pi;ws=40*2*pi; >> [n,wn]=buttord(wp,ws,rp,rs,'s'); >> [z,p,k]=buttap(n); >> [num,den]=zp2tf(z,p,k); >> [num1,den1]=impinvar(num,den); Warning: The output is not correct/robust. Coeffs of B(s)/A(s) are real, but B(z)/A(z) has complex coeffs. Probable cause is rooting of high-order repeated poles in A(s). > In impinvar at 124 >> [num2,den2]=bilinear(num,den,100); >> [h,w]=freqz(num1,den1); >> [h1,w1]=freqz(num2,den2); >>subplot(1,2,1); >>plot(w*fs/(2*pi),abs(h)); >>subplot(1,2,2); >>plot(w1*fs/(2*pi),abs(h1)); >>figure(1); >>subplot(1,2,1); >>zplane(num1,den1); >>subplot(1,2,2); >>zplane(num2,den2);