DSP实习实验报告
dsp实验报告总结doc
3.2 DSP总体结构
3.3 电源模块设计
3.4 时钟模块设计
3.5 存储器模块设计
3.6复位模块设计
篇二：DSP实验报告
DSP课程设计 实 验 报 告
语音压缩、存储和回放
学 院：电子信息工程学院电子科学与技术专业 设计人员： 吴莲梅 08214085电子0803班 杨 莹 08214088电子0803班指导老师： 日 期：
（1）A律限制采样值为12比特，A律的压缩可以按照下列公式进行定义：
A|x|11?lnA|x|1
(0?|x|?)?sgn(x)(?|x|?1)F(x)?sgn(x)
1?lnAA1?lnAA
式中，A是压缩参数（在欧洲，A=87.6）x是需要压缩的归一化整数。从线性到A律的压缩转换如下表所示：（压缩后的码字组成：比特0-3表示量化值，比特4-6表示段值，压缩后
一、 设计目的
设计一个功能完备，能够独立运行的精简DSP硬件系统，并设计简单的DSP控制程序。
二、 系统分析
1.1设计要求 硬件要求：
（1）使用TMS320VC5416作为核心芯片。 （2）具有最简单的led控制功能。 （3）具有存放程序的外部Flash芯片。 （4）外部输入+5V电源。 （5）绘制出系统的功能框图。
（6） 仪器仪表--如频谱分析、函数发生、锁相环、地震处理等。 （7） 自动控制--如引擎控制、深空、自动驾驶、机器人控制、磁盘控制。 （8） 医疗--如助听、超声设备、诊断工具、病人监护等。（9） 家用电器--如高保真音响、音乐合成、音调控制、玩具与游戏、数字 电话/电视等 DSP 的发展前景 DSP 的功能越来越强，应用越来越广，达到甚至超过了微控制器的功能，比 微控制器做得更好而且价格更便宜， 许多家电用第二代 DSP 来控制大功率电机就 是一个很好的例子。汽车、个人通信装置、家用电器以及数以百万计的工厂使用 DSP 系统。数码相机、IP 电话和手持电子设备的热销带来了对 DSP 芯片的巨大需 求。而手机、
dsp原理与应用实验报告总结
dsp原理与应用实验报告总结DSP（Digital Signal Processing）数字信号处理是利用数字技术对信号进行处理和分析的一种方法。
在本次实验中，我们探索了DSP的原理和应用，并进行了一系列实验以验证其在实际应用中的效果。
以下是对实验结果的总结与分析。
实验一：数字滤波器设计与性能测试在本实验中，我们设计了数字滤波器，并通过性能测试来评估其滤波效果。
通过对不同类型的滤波器进行设计和实现，我们了解到数字滤波器在信号处理中的重要性和应用。
实验二：数字信号调制与解调本实验旨在通过数字信号调制与解调的过程，了解数字信号的传输原理与方法。
通过模拟调制与解调过程，我们成功实现了数字信号的传输与还原，验证了调制与解调的可行性。
实验三：数字信号的傅里叶变换与频谱分析傅里叶变换是一种重要的信号分析方法，可以将信号从时域转换到频域，揭示信号的频谱特性。
本实验中，我们学习了傅里叶变换的原理，并通过实验掌握了频谱分析的方法与技巧。
实验四：数字信号的陷波滤波与去噪处理陷波滤波是一种常用的去除特定频率噪声的方法，本实验中我们学习了数字信号的陷波滤波原理，并通过实验验证了其在去噪处理中的有效性。
实验五：DSP在音频处理中的应用音频处理是DSP的一个重要应用领域，本实验中我们探索了DSP在音频处理中的应用。
通过实验，我们成功实现了音频信号的降噪、均衡和混响处理，并对其效果进行了评估。
实验六：DSP在图像处理中的应用图像处理是另一个重要的DSP应用领域，本实验中我们了解了DSP在图像处理中的一些基本原理和方法。
通过实验，我们实现了图像的滤波、边缘检测和图像增强等处理，并观察到了不同算法对图像质量的影响。
通过以上一系列实验，我们深入了解了DSP的原理与应用，并对不同领域下的信号处理方法有了更深刻的认识。
本次实验不仅加深了我们对数字信号处理的理解，也为日后在相关领域的研究与实践提供了基础。
通过实验的结果和总结，我们可以得出结论：DSP作为一种数字信号处理的方法，具有广泛的应用前景和重要的实际意义。
DSP实验报告6
DSP第六、七次实验报告1. 实验目的:（1）进一步熟悉Matlab实验环境和语言。
（2）熟悉各种滤波器的结构及Matlab实现语言。
（3）掌握用冲击响应不变法和双线性变换法设计IIR滤波器的方法。
（4）掌握用窗函数法和频率抽样法设计FIR滤波器的方法。
2. 实验内容及总结:1.滤波器结构:（1）IIR滤波器各种结构1、直接型结构例如直接型滤波器系统函数, 则有系数向量a=[1,a1,a2,a3],b=[b0,b1,b2], 利用:Y=filter[b,a,x]求信号x(n)通过此滤波器的输出。
2、由系统函数或差分方程求系统的二阶分式（含一阶分式）的级联结构将例如的系统函数重写为二阶分式节的级联型, 利用:[sos,G]=tf2sos(b,a)3、由二阶分式的级联结构转换成系统函数的直接结构是第二步的逆运算, 调用函数:[b,a] = sos2tf(sos)可以求得系数向量a,b, 从而得到H(z)4、由系统函数求部分分式展开（留数及其极点计算）即求z反变换的部分分式展开法, 利用:[r,p,c]=residuez(b,a)其中极点为p, 留数为r, 直接项系数为c。
5、由r,p,c求系统函数即第4步的逆运算, 利用:[b,a]=residuez(r,p,c)6、由直接型结构转换为并联型结构需开发函数:[C,B,A]=tf2par(b,a)其中, b,a为直接型的系数向量, C,B,A为并联型实系数向量, 基本思想是: 1.反复调用[r,p,c]=residuez(b,a)求出极点及留数；2.利用cplxpair函数把极点、留数对按复共轭极点-留数对, 实极点-留数对的顺序排列；3.开发cplxcomp函数, 保证极点和留数相互对应；4.调用[b,a]=residuez(r,p,c)计算并联二阶节的分子分母。
7、由并联型结构转换成直接型结构开发函数:[b,a]=par2tf(C,B,A)为[C,B,A]=tf2par(b,a)的逆函数。
DSP实验报告+心得体会
Gpio_select();
// Step 3. Clear all interrupts and initialize PIE vector table:
// Disable CPU interrupts
GpioMuxRegs.GPAQUAL.all=var3;// Set GPIO input qualifier values
GpioMuxRegs.GPBQUAL.all=var3;
GpioMuxRegs.GPDQUAL.all=var3;
GpioMuxRegs.GPEQUAL.all=var3;
四、实验步骤
1.把2812模块小板插到大板上；
2.在CCS2000环境中打开本实验工程编译Example_7segled.prj，生成输出文件，通过仿真器把执行代码下载到DSP芯片；
3.运行程序;数码管会显示1～8的数字。
4.参考源代码自行修改程序改变显示样式。
五、实验心得体会
通过本次实验中，基本掌握了2812的指令系统的特点，并能够了解并熟悉74HC573的使用方法，进一步加深了对DSP的认识。同时，通过实验操作DSP的IO操作使用方法，对于DSP的IO操作可以熟悉的运用，学到更多的知识。
void delay_loop(void);
void Gpio_select(void);
// Global variable for this example
short codetab[17]=
{0x4020,0x6cc0,0x5800,0x4840,0x6440,0xC040,0xC000,0x4cc0,
DSP实验报告――离散信号的产生及运算(精)
离散信号的产生(chǎnshēng)及运算报告一、实验(shíyàn)目的：1、复习(fùxí)和巩固数字信号处理中离散信号的产生和运算2、学习和掌握(zhǎngwò)用MATLAB产生离散信号的方法3、学习(xuéxí)和掌握用MATLAB对离散信号进行运算二、实验原理：1．用MATLAB函数产生离散信号信号是数字信号处理的最基本内容。
没有信号，数字信号处理就没了工作对象。
MATLAB7.0内部提供了大量的函数，用来产生常用的信号波形。
例如，三角函数（sin,cos）,指数函数（exp）,锯齿波函数（sawtooth）, 随机数函数（rand）等。
1 产生被噪声污染的正弦信号用随机数函数产生污染的正弦信号。
2 产生单位脉冲序列和单位阶跃序列按定义，单位脉冲序列为单位阶跃序列为。
3 矩形脉冲信号：在MATLAB 中用rectpuls 函数来表示，其调用形式为：y=rectpuls(t,width，用以产生一个幅值为1，宽度为width,相对于t=0 点左右对称的矩形波信号，该函数的横坐标范围(fànwéi)由向量t 决定，是以t=0 为中心向左右各展开width/2 的范围，width 的默认值为1。
例：以t=2T(即t-2×T=0为对称中心的矩形脉冲信号(xìnhào)的MATLAB 源程序如下：（取T=1）t=0:0.001:4;T=1;ft=rectpuls(t-2*T,2*T;plot(t,ft;grid on; axis([0 4 –0.5 1.5];4 周期性矩形波（方波）信号在MATLAB 中用square 函数来表示，其调用形式为：y=square(t,DUTY，用以产生一个周期为2π、幅值为±1 的周期性方波信号，其中的DUTY参数表示占空比，即在信号的一个周期中正值(zhènɡ zhí)所占的百分比。
DSP实验报告_6
实验一: 闪灯实验熟悉DSP 软硬件测试系统实验目的1.了解SHARC 系列高性能数字信号处理器的程序开发过程和编程语言；2.熟悉集成开发工具VisualDSP++, 学会使用VisualDSP++进行SHARC 系列ADSP 的程序开发、编译与调试；3.掌握SHARC 系列ADSP 的程序加载设计和加载过程。
实验内容利用波形产生信号板, 结合FPGA 编程技术和程序编程器, 编写测试ADSP21065L 和FPGA 之间硬件连接的应用程序, 同时完成应用程序的加载和脱机操作, 在信号指示灯“HL2”上产生可调周期的脉冲信号, “点亮”与“熄灭”指示灯HL2。
实验要求通过DSP 编程, 在其FLAG11引脚上模拟如下波形的周期信号:要求:(1) 500H T ms >,500L T ms >. (2) 并用示波器查看波形, 测量信号周期。
实验步骤1. 熟悉电路图, 清楚波形产生电路板ADSP21065L 与可编程FPGA 器件之间的连接关系；2. 编写FPGA 程序。
在FPGA 内部将ADSP21065L 的标志引脚FLAG11（引脚号26）设置为输出, 作为FPGA 的输入信号, 在FPGA 内部编程将该信号直接输出在发FPGA 的37引脚号上, 设置37引脚为输出信号, 驱动板上的HL2 LED 指示灯；3. 启动VisualDsp++4.5,选择project 工程选项菜单, 创建一个名称为Test.dpj 的工程文件, 选择处理器的型号为ADSP-21065L ；4．弹出一个对话框, 选择是否需要加入VDSP kernel ,选择“NO ”；5. 在工程中加入以下参考源文件:\exp1\test(boot)\ boot1.asm 和boot1.ldf 6．编译, 链接调试, 生成可执行文件。
7．运行程序, 可以看到波形发生电路板上的指示灯“HL2”不断闪动。
8. 利用示波器观测系统时钟,并测量产生信号的波形和周期。
DSP实验报告(综合)
实验报告||实验名称 D SP课内系统实验课程名称DSP系统设计||一、实验目的及要求1. 掌握用窗函数法设计FIR数字滤波器的原理和方法。
熟悉线性相位FIR 数字滤波器特性。
了解各种窗函数对滤波器特性的影响。
2. 掌握设计IIR数字滤波器的原理和方法。
熟悉IIR数字滤波器特性。
了解IIR数字滤波器的设计方法。
3.掌握自适应数字滤波器的原理和实现方法。
掌握LMS自适应算法及其实现。
了解自适应数字滤波器的程序设计方法。
4.掌握直方图统计的原理和程序设计；了解各种图像的直方图统计的意义及其在实际中的运用。
5.了解边缘检测的算法和用途，学习利用Sobel算子进行边缘检测的程序设计方法。
6.了解锐化的算法和用途，学习利用拉普拉斯锐化运算的程序设计方法。
7.了解取反的算法和用途，学习设计程序实现图像的取反运算。
8.掌握直方图均衡化增强的原理和程序设计；观察对图像进行直方图均衡化增强的效果。
二、所用仪器、设备计算机，dsp实验系统实验箱，ccs操作环境三、实验原理（简化）FIR：有限冲激响应数字滤波器的基础理论，模拟滤波器原理（巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器、贝塞尔滤波器）。
数字滤波器系数的确定方法。
IIR:无限冲激响应数字滤波器的基础理论。
模拟滤波器原理（巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器、贝塞尔滤波器）。
数字滤波器系数的确定方法。
、自适应滤波：自适应滤波器主要由两部分组成：系数可调的数字滤波器和用来调节或修正滤波器系数的自适应算法。
e(n)=z(n)-y(n)=s(n)+d(n)-y(n)直方图：灰度直方图描述了一幅图像的灰度级内容。
灰度直方图是灰度值的函数，描述的是图像中具有该灰度值的像素的个数，其横坐标表示像素的灰度级别，纵坐标是该灰度出现的频率(像素个数与图像像素总数之比)。
图像边缘化：所谓边缘(或边沿)是指其周围像素灰度有阶跃变化。
经典的边缘提取方法是考察图像的每个像素在某个邻域内灰度的变化，利用边缘临近一阶或二阶方向导数变化规律，用简单的方法检测边缘。
DSP实验报告(一)
实验一 信号系统及系统响应一、实验目的1、 熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对采样定理的理解。
2、 熟悉离散信号和系统的时域特性。
3、 熟悉线性卷积的计算编程方法：利用卷积的方法，观察、分析系统响应的时域特性。
4、 掌握序列傅氏变换的计算机实现方法，利用序列的傅氏变换对离散信号、系统及系统响应进行频域分析。
二、实验原理（一）连续时间信号的采样采样是指按一定的频率从模拟信号抽样获得数字信号。
采样是从连续时间信号到离散时间信号的过渡桥梁。
对一个连续时间信号进行理想采样的过程可以表示为该信号的一个周期冲激脉冲的乘积，即()()()ˆa a x t x t M t =（1）其中连续信号的理想采样，是周期冲激脉冲()()n M t t n T d +=-=-å（2）它也可以用傅立叶级数表示为：1()s jm tn M t eT+W =-=å（3)其中T 为采样周期，Ω是采样角频率。
设是连续时间信号的双边拉氏变换，即有：()()ata a X s x t edt+--=ò（4）此时理想采样信号的拉氏变换为()ˆˆ()()1ˆ()1ˆ()1()s s ataa jm tsta m s jm ta m a s m X s x t e dtxt ee dtTxt e dtT X s jm T+--++W -=--++--W =- -++=--====-W òåòåòåò（5）作为拉氏变换的一种特例，信号理想采样的傅立叶变换1ˆ()[()]aa s m X j X j m T+=-W =W-W å（6）由式(5)和式(6)可知，信号理想采样后的频谱是原信号频谱的周期延拓，其延拓周期等于采样频率。
根据Shannon 取样定理，如果原信号是带限信号，且采样频率高于原信号最高频率分量的2倍，则采样以后不会发生频谱混淆现象。
DSP实验报告5_2
DSP 第五次实验1．实验目的：（1）进一步熟悉matlab 实验环境和语言。
（2）掌握求序列圆周翻褶的MATLAB 方法。
（3）掌握求序列DFT 及IDFT 矩阵的MATLAB 方法。
（4）掌握用MATLAB 求解用圆周卷积计算线性卷积的时域的方法。
（5）掌握用FFT 计算有限长序列的线性卷积和线性相关的方法。
2．实验内容及总结：1.圆周翻褶【例3．27】 已知()[2,3,4,5,6],8X n N ==，求x(n)的8点圆周翻褶序列88(())()x n R n -。
代码：clc;clear allx=[2,3,4,5,6];N=8;x=[x,zeros(1,N-length(x))];nx=0:N-1y=x(mod(-nx,N)+1);subplot(121),stem([0:N-1],x);title('原序列');xlabel('n');ylabel('x(n)');grid;subplot(122),stem([0:N-1],y);title('圆周翻褶序列');xlabel('n');ylabel('x((n))8 R8(n)');grid;结果：总结：对于圆周翻褶(0),0()(())()(),11N Nx ny n x n R nx N n n N==-=⎨-≤≤-MA TLAB可用y=x(mode(-nx,N)+1)求得。
因此，要求X(n)=[2,3,4,5,6]，N=8的8点圆周翻褶序列，要先将x(n)补零到8点长度再求圆周翻褶。
x=[x,zeros(1,N-length(x))];nx=0:N-1 %x补零到8点长y=x(mod(-nx,N)+1); %圆周翻褶从一开始，因此得到8点长%序列，应该再加一2.DFT矩阵，IDFT矩阵【例3.29】已知N=4的DFT矩阵w4，求IDFT矩阵w4I。
DSP实验报告(二)
DSP实验报告(二)实验二应用FFT对信号进行频谱分析一、实验目的1、在理论学习的基础上，通过本次实验，加深对快速傅里叶变换的理解，熟悉FFT算法及其程序的编写。
2、熟悉应用FFT对典型信号进行频谱分析的方法。
3、了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用FFT。
二、实验原理与方法①一个连续信号的频谱可以用它的傅立叶变换表示为+ Xa(jW)=-jWtx(t)edtòa-如果对该信号进行理想采样，可以得到采样序列x(n)=xa(nT)同样可以对该序列进行z变换，其中T为采样周期X(z)=+ x(n)z-n+ -令z为ejw，则序列的傅立叶变换X(ejw)=x(n)ejwn-其中ω为数字频率，它和模拟域频率的关系为w=WT=W/fs式中的是采样频率。
上式说明数字频率是模拟频率对采样率的归一化。
同模拟域的情况相似。
数字频率代表了序列值变化的速率，而序列的傅立叶变换称为序列的频谱。
序列的傅立叶变换和对应的采样信号频谱具有下式的对应关系。
1X(e)=Tjw+ - w-2pXa(j)T即序列的频谱是采样信号频谱的周期延拓。
从式可以看出，只要分析采样序列的谱，就可以得到相应的连续信号的频谱。
注意：这里的信号必须是带限信号，采样也必须满足Nyquist定理。
在各种信号序列中，有限长序列在数字信号处理中占有很重要的地位。
无限长的序列也往往可以用有限长序列来逼近。
有限长的序列可以使用离散傅立叶变换。
当序列的长度是N时，定义离散傅立叶变换为：X(k)=DFT[x(n)]=其中W=e2pj-NN-1n=0WNkn它的反变换定义为：1x(n)=IDFT[X(k)]=N根据式和，则有N-1n=0X(k)WNknX(z)|z=Wnk=NN-1n=0x(n)WNnk=DFT[x(n)]j2pN可以得到X(k)2pk的点，就NN是将单位圆进行N等分以后第k个点。
所以，X(k)是z变换在单位圆上的等距采样，或者说是序列傅立叶变换的等距采样。
dsp实验报告
dsp实验报告实验一:CCS入门实验实验目的:1. 熟悉CCS集成开发环境，掌握工程的生成方法;熟悉SEED-DEC643实验环境; 掌握CCS集成开发环境的调试方法。
2.学习用标准C 语言编写程序;了解TI CCS开发平台下的C 语言程序设计方法和步骤; 熟悉使用软件仿真方式调试程序。
3. 学习用汇编语言编写程序; 了解汇编语言与 C 语言程序的区别和在设置上的不同;了解TMS320C6000 汇编语言程序结果和一些简单的汇编语句用法学习在CCS 环境中调试汇编代码。
4. 在了解纯C 语言程序工程和汇编语言程序工程结构的基础上，学习在C 工程中加入汇编编程的混合编程方法; 了解混合编程的注意事项;理解混合编程的必要性和在什么情况下要采用混合编程5. 熟悉CCS集成开发环境，掌握工程的生成方法; 熟悉SEED-DEC643实验环境;掌握CCS集成开发环境的调试方法。
实验原理:CCS 提供了配置、建立、调试、跟踪和分析程序的工具，它便于实时、嵌入式信号处理程序的编制和测试，它能够加速开发进程，提高工作效率。
CCS 提供了基本的代码生成工具，它们具有一系列的调试、分析能力序。
使用此命令后，要重新装载.out 文件后，再执行程序。
使用 CCS常遇见文件简介1. program.c: C 程序源文件;2. program.asm: 汇编程序源文件;3. filename.h: C 程序的头文件，包含DSP/BIOS API模块的头文件;4. filename.lib: 库文件;5. project.cmd: 连接命令文件;6. program.obj: 由源文件编译或汇编而得的目标文件;7. program.out: 经完整的编译、汇编以及连接后生成可执行文件; 8. program.map: 经完整的编译、汇编以及连接后生成空间分配文件; 9.project.wks: 存储环境设置信息的工作区文件。
P.S(CMD文件中常用的程序段名与含义1. .cinit 存放C程序中的变量初值和常量;2. .const 存放C程序中的字符常量、浮点常量和用const声明的常量;3. .text 存放C程序的代码;4. .bss 为C 程序中的全局和静态变量保留存储空间;5. .far 为C 程序中用far声明的全局和静态变量保留空间;6. .stack 为 C 程序系统堆栈保留存储空间，用于保存返回地址、函数间的参数传递、存储局部变量和保存中间结果;7. .sysmem 用于 C 程序中malloc、calloc 和 realloc 函数动态分配存储空间。
DSP实验报告最终版
小组成员
一、实验说明：
为了保证生产质量、生产效率和生产的安全性，在钢铁企业的定期常规设备检修必不可少，在设备检修的过程中，为了保证检修人员的人身安全和设备的安全，通常需设置检修报警提示，如：检修警报、检修指示灯等。
本次实验设计钢铁企业的滚带传送装置的检修控制系统，需要检修时，首先停止传动皮带的工作，切换到检修警报，检修警报开始工作以保证检修人员和设备的安全。当检修工作完成后，停止警报装置，返回传动电机的正常工作状态。
南-绿
北-红
北-黄
北-绿
两个寄存器的地址均映射到2812DSP的扩展空间，CTRLR，DSP通过对该地址的写操作来修改二个寄存器各位的状态，当寄存器某位取‘1’时，相应提示灯被点亮，取‘0’则熄灭。当写入CTRLR的数据（8位有效值）的高两位为‘00’时，数据的低6位将写入EWR寄存器；高两位为‘01’时，数据的低6位将写入SNR寄存器。
图3 利用开关管对直流电动机进行PWM调速控制的原理图和输入输出电压波形 上图是利用开关管对直流电动机进行PWM调速控制的原理图和输入输出电压波形。图中,当开关管MOSFET的栅极输入高电平时，开关管导通，直流电动机电枢绕组两端有电压Us，t1秒后，栅极输入变为低电平，开关管截止，电动机电枢两端电压为0。T2秒后，栅极输入重新变为高电平，开关管的动作重复前面的过程。这样，对应着输入的高低电平，直流电动机电枢绕组两端的电压波形如图中所示。电动机的电枢绕组两端的电压平均值Uo为：
图中PWM输入对应ICETEK-CTR-A板上P4外扩插座第26引脚的S22信号，DSP将在此引脚上给出PWM信号开控制直流电机转速；图中的DIR输入引脚ICETEK-CTR-A板上P4外扩插座第29引脚的S14信号，DSP将在此引脚上给出高电平或低电平来控制直流电机的方向。从DSP输出的PWM信号和转向信号先经过2个与门和1个非门再与各个开关管的栅极相连。
DSP实验报告一
DSP实验报告一引言本实验旨在通过实际操作，探索数字信号处理（DSP）的基本概念和技术。
DSP是一种通过数字计算来处理连续时间信号的技术，被广泛应用于音频处理、图像处理、通信系统等领域。
本实验将重点介绍数字信号的采样、量化和离散化过程，并通过实际编程实现。
实验过程1. 信号的采样1.1 信号的定义在DSP领域，信号是指随着时间变化的某种物理量，可以是声音、图像等。
我们首先需要定义一个连续的信号，用于采样和处理。
在本次实验中，我们选择了一个简单的正弦信号作为示例：x(t) = A \\sin(2\\pi f t)其中，A表示幅值，f表示频率，t表示时间。
1.2 采样过程为了将连续信号转换为离散信号，我们需要对信号进行采样。
采样是指在一定时间间隔内对连续信号进行测量。
我们可以通过模拟采样器来模拟采样过程。
在本实验中，我们选择了采样频率为100Hz，即每秒采样100次。
使用Python编程实现采样过程：import numpy as np# 信号参数设置A =1f =10# 采样频率设置fs =100# 采样点数设置N =100# 生成时间序列t = np.arange(N) / fs# 生成采样信号x = A * np.sin(2* np.pi * f * t)上述代码中，我们通过调整A和f的值来模拟不同的信号。
生成的信号将存储在x变量中，可以用于后续处理。
2. 信号的量化2.1 量化过程量化是指将连续信号的幅值转换为离散的数值。
在实际应用中，我们通常使用有限位数来表示信号的幅值。
常用的量化方式有线性量化和非线性量化。
在本实验中，我们选择了线性量化方式。
具体的量化过程可以通过下列Python代码实现：import math# 量化位数设置bits =8# 量化步长计算step_size =2* A / (2** bits -1)# 信号的量化x_quantized = np.round(x / step_size) * step_size上述代码中，我们通过调整bits的值来控制量化位数。
dsp实验报告
dsp实验报告DSP实验报告一、引言数字信号处理（Digital Signal Processing，DSP）是一种对数字信号进行处理和分析的技术。
它在许多领域中被广泛应用，如通信、音频处理、图像处理等。
本实验旨在通过实际操作，探索和理解DSP的基本原理和应用。
二、实验目的1. 理解数字信号处理的基本概念和原理；2. 掌握DSP实验平台的使用方法；3. 进行一系列DSP实验，加深对DSP技术的理解。
三、实验器材和软件1. DSP开发板；2. 电脑；3. DSP开发软件。
四、实验内容1. 实验一：信号采集与重构在此实验中，我们将通过DSP开发板采集模拟信号，并将其转换为数字信号进行处理。
首先，我们需要连接信号源和开发板，然后设置采样频率和采样时间。
接下来，我们将对采集到的信号进行重构，还原出原始模拟信号，并进行观察和分析。
2. 实验二：滤波器设计与实现滤波器是DSP中常用的模块，用于去除或增强信号中的特定频率成分。
在此实验中，我们将学习滤波器的设计和实现方法。
首先，我们将选择合适的滤波器类型和参数，然后使用DSP开发软件进行滤波器设计。
最后，我们将将设计好的滤波器加载到DSP开发板上，并进行实时滤波处理。
3. 实验三：频谱分析与频域处理频谱分析是DSP中常用的方法，用于分析信号的频率成分和能量分布。
在此实验中，我们将学习频谱分析的基本原理和方法，并进行实际操作。
我们将采集一个包含多个频率成分的信号，并使用FFT算法进行频谱分析。
然后，我们将对频谱进行处理，如频率选择、频率域滤波等，并观察处理后的效果。
4. 实验四：音频处理与效果实现音频处理是DSP中的重要应用之一。
在此实验中，我们将学习音频信号的处理方法，并实现一些常见的音频效果。
例如，均衡器、混响、合唱等。
我们将使用DSP开发软件进行算法设计，并将设计好的算法加载到DSP开发板上进行实时处理。
五、实验结果与分析通过以上实验，我们成功完成了信号采集与重构、滤波器设计与实现、频谱分析与频域处理以及音频处理与效果实现等一系列实验。
硕士信号处理实验报告(3篇)
第1篇一、实验背景随着信息技术的飞速发展，数字信号处理（DSP）技术已成为通信、图像处理、语音识别等领域的重要工具。
本实验旨在通过一系列实验，加深对数字信号处理基本原理和方法的理解，提高实际应用能力。
二、实验目的1. 理解数字信号处理的基本概念和原理。
2. 掌握常用信号处理算法的MATLAB实现。
3. 培养分析和解决实际问题的能力。
三、实验内容本实验共分为五个部分，具体如下：1. 离散时间信号的基本操作（1）实验目的：熟悉离散时间信号的基本操作，如加法、减法、乘法、除法、延时、翻转等。
（2）实验步骤：- 使用MATLAB生成两个离散时间信号。
- 对信号进行基本操作，如加法、减法、乘法、除法、延时、翻转等。
- 观察并分析操作结果。
2. 离散时间系统的时域分析（1）实验目的：掌握离散时间系统的时域分析方法，如单位脉冲响应、零状态响应、零输入响应等。
（2）实验步骤：- 使用MATLAB设计一个离散时间系统。
- 计算系统的单位脉冲响应、零状态响应和零输入响应。
- 分析系统特性。
（1）实验目的：掌握离散时间信号的频域分析方法，如快速傅里叶变换（FFT）、离散傅里叶变换（DFT）等。
（2）实验步骤：- 使用MATLAB生成一个离散时间信号。
- 对信号进行FFT和DFT变换。
- 分析信号频谱。
4. 数字滤波器的设计与实现（1）实验目的：掌握数字滤波器的设计与实现方法，如巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器等。
（2）实验步骤：- 使用MATLAB设计一个低通滤波器。
- 使用窗函数法实现滤波器。
- 对滤波器进行性能分析。
5. 信号处理在实际应用中的案例分析（1）实验目的：了解信号处理在实际应用中的案例分析，如语音信号处理、图像处理等。
（2）实验步骤：- 选择一个信号处理应用案例。
- 分析案例中使用的信号处理方法。
- 总结案例中的经验和教训。
四、实验结果与分析1. 离散时间信号的基本操作实验结果表明，离散时间信号的基本操作简单易懂，通过MATLAB可以实现各种操作，方便快捷。
DSP实验报告_百度文库(精)
实验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课程的实用性和趣味性。
实验二基本算数运算2.1 实验目的和要求加、减、乘、除是数字信号处理中最基本的算术运算。
DSP实验报告(完美版)
DSP实验报告班级：11050641学号：姓名：指导教师：实验一、二 DSP芯片的开发工具及应用实验1.实验目的（1）熟悉CCS集成开发环境，掌握工程的生成方法；（2）熟悉SEED-DTK DAD实验环境；（3）掌握CCS集成开发环境的调试方法。
2.实验设备DSP实验箱，计算机，CCS软件。
3.实验内容及步骤（1）CCS软件的安装；（2）了解SEED-DTK5416实验环境；（3）打开CCS集成开发环境，进入CCS的操作环境；（4）新建一个工程文件○1在c:\ti\myprojects中建立文件夹volume1(如果CCS安装在其他d:\ti ,则在d:\ti\myprojects中)；○2将c:\ti\tutorial\target\volume1拷贝到c:\ti\myprojects\ volume1；○3从在CCS 中的Project 菜单,选择 New；○4在Project Name域中，键入volume1；○5在Location区域中，浏览步骤1所建立的工作文件夹；○6在Project Type 域中，选择Executable(.out)；○7在Target域中，选择CCS配置的目标，并单击完成。
（5）向工程中添加文件○1从Project/Add Files to Project，选择 volume.c，单击Open（或右击Project View图标，选择Add Files to Project ）；○2选择Project/Add Files to Project，在Files of type对话框中，选择Asm Source Files (*.a*, *.s*)。
选择vectors.asm 和 load.asm, 单击Open；○3选择 Project/Add Files to Project，在Files of type 对话框中选择 Linker Command File (*.cmd)，选择volume.cmd，单击Open。
DSP技术实习报告
DSP技术实习报告相关推荐DSP技术实习报告⼀、实习题⽬1、灰度线性变换2、灰度的对数变换3、锯齿波变换4、图像剪取⼆、实习⽬的：1、巩固和深化数字图像处理技术所涉及的数理基础、基本算法和各种图像处理技术⽅法，学习和掌握图像变换。
2、对学习图像处理的基础知识对其应⽤⼯程实践有⼀定的认识，提⾼学⽣对应⽤软件的使⽤能⼒。
3、通过理论联系实际，综合运⽤所学知识，提⾼学⽣独⽴分析和解决实际问题的能⼒，增强学⽣的⼯程意识，打好专业基础三、实习要求:1、能够根据设计题⽬要求查阅检索有关的⽂献资料，结合题⽬选学有关参考书;2、熟悉计算机图像处理的设计⽅法;3、熟悉图像灰度变换程序的设计⽅法;4、掌握图像灰度变换的仿真⽅法;5、完成图像的灰度变换。
四、系统原理描述:灰度变换将输⼊图像映射为输出图像，输出图像每个像素点的灰度值仅由对应的输⼊像素点的值决定。
它常⽤于改变图像的灰度范围及分布，是图像数字化及图像显⽰的重要⼯具。
灰度变换因其作⽤性质有时也被称为对⽐度增强、对⽐度拉伸或点运算，称之为灰度变换。
灰度变换实际上是灰度到灰度的映射过程。
设输⼊图像为A(x，y)输出图像为B(x，y)，则灰度变换可表⽰为：即灰度变换完全由灰度映射函数决定。
显然灰度变换不会改变图像内像素点之间的空间关系。
1. 灰度线性变换设图像灰度值f(x，y)的可能值域为D，但在⼀定条件下，使得其取值范围缩⼩为且如图 7-11所⽰，这种状态常出现于连续图f(x，y)值的动态范围⼩，或摄影曝光不⾜中，因⽽其对应的直⽅图P (D)如图所⽰，集中于某个较⼩的灰度区间内。
但是通常⼈们希望灰度直⽅图在全部灰度区域内均匀分布，最简单地是把低反差图像进⾏灰度变换得到⾼反差图像，即线性变换，如下式：该式可以使变换后的灰度的取值范围扩⼤到⼀般来说，只有当两个相邻像素的灰度值(亮度值)相差到⼀定程度时，⼈的视觉才能分辨。
若灰度值D仅在较⼩区间内时，则⼈眼可分辨的亮度差的总级数则亦很少，从⽽造成⽬标图像灰度值与背景灰度值相接近，⼈眼⽽⽆法分辨检出。
DSP实验报告(精)
DSP上机实验：DTMF信号的编码一、实验要求设计作业：双音多频（DTMF）信号的编码：把自己的电话号码DTMF编码生成为一个.wav文件。
【wavwrite()】技术指标：根据ITU Q.23建议，DTMF信号的技术指标是：传送/接收率为每秒10个号码，或每个号码100ms。
每个号码传送过程中，信号存在时间至少45ms，且不多于55ms，100ms的其余时间是静音。
在每个频率点上允许有不超过±1.5%的频率误差。
任何超过给定频率±3.5%的信号，均被认为是无效的，拒绝接收。
二、实验思想：1、编码：DTMF拨号键盘由一个4*4行列构成，每列代表一个高频信号，每行代表一个低频信号，每当按下一个键时，产生高、低频率的两个正弦信号，代表一个特定的数字或符号，根据ITU Q.23颁布的国际标准，DTMF传送或接受每个号码的时间为100ms，其中每个号码传送的过程中，信号存在时间至少45ms，其余时间静音。
用一个字符串变量来接受输入的电话号码，并将各个数字和符号的ASCII码用一个4*4矩阵表示，每接收到一个数字就对应两个频率，并产生由两个正弦波叠加的信号，完成DTMF编码，利用matlab提供的fft函数画出其频谱，用sound函数发出声音。
2、解码：采用Goertzel算法来检测DTMF信号，它是用IIR滤波器实现DFT算法的一种特殊方法，在实际DTMF解码中，只需要知道输入信号即DTMF信号的离散傅里叶变换X(k)的幅度信息，忽略相位信息，因为只要能得到8个特定频率点的幅度值，看哪两个频率对应的幅度最大，就能知道对应的是哪个数字，达到解码的目的，可以利用matlab提供的goertzel函数来对信号解码。
三、实验代码：d=input('请键入电话号码： ','s');sum=length(d);total_x=[];sum_x=[];sum_x=[sum_x,zeros(1,800)];for a=1:sumsymbol=abs(d(a));tm=[49,50,51,65;52,53,54,66;55,56,57,67;42,48,35,68]; for p=1:4;for q=1:4;if tm(p,q)==abs(d(a));break,endendif tm(p,q)==abs(d(a));break,endendf1=[697,770,852,941];f2=[1209,1336,1477,1633];n=1:400;x=sin(2*pi*n*f1(p)/8000)+sin(2*pi*n*f2(q)/8000); x=[x,zeros(1,400)]; sum_x=sum_x+x;total_x=[total_x,x];endwavwrite(total_x,'soundwave')sound(total_x);t=(1:8800)/8000;subplot(2,1,1);plot(t,total_x);axis([0,1.2,-2,2]);xlabel('时间/s')title('DTMF信号时域波形')xk=fft(x); mxk=abs(xk);subplot(2,1,2);k=(1:800)*sum*8000/800;plot(k,mxk); xlabel('频率');title('DTMF信号频谱');disp('双频信号已生成并发出')四、实验结果图：请键入电话号码： 130xxxxxxxx。
dsp实验报告 3
实验三PWM实验：一、实验目的了解TMS320F28335 的PWM 模块原理二、实验设备(1)装有Windows 的PC 机一台；(2) XDS510 仿真器一套；(3) YX-F28335 开发板一套；(4) 示波器一台；三、实验步骤(1) 首先按照实验一配置CCS4.1.2 软件并打开；(2) 接着把仿真器的USB 与电脑进行连接，将仿真器的另一端JATG 端插到YX-F28335 开发板的JATG 针处；(3) Target->Launch TI Debug 后，点击Target->Connect Target。
(4) 由于工程已经是一个可烧写的可执行文件，所以直接在CCS 中点击Target->LoadProgram……命令，在文件lab24-PWM 下加载Debug 目录下的.out 可执行文件；(5) 在CCS 菜单栏点击Target->Run，之后用户打开示波器，将示波器的地线接到开发板的地线端，另一端接到YX-F28335 开发板J4的第1 脚。
四、实验原理脉冲宽度调制（PWM）是一种对模拟信号电平进行数字编码的方法，其根据相应载荷的变化来调制晶体管栅极或基极的偏置，来实现开关稳压电源输出晶体管或晶体管导通时间的改变，这种方式能使电源的输出电压在工作条件变化时保持恒定，是利用微处理器的数字输出来对模拟电路进行控制的一种非常有效的技术，广泛应用在从测量、通信到功率控制与变换的许多领域中。
PWM 的控制方法：采样控制理论中有一个重要结论：冲量相等而形状不同的窄脉冲加在具有惯性的环节上时,其效果基本相同。
PWM 控制技术就是以该结论为理论基础,对半导体开关器件的导通和关断进行控制,使输出端得到一系列幅值相等而宽度不相等的脉冲,用这些脉冲来代替正弦波或其他所需要的波形。
按一定的规则对各脉冲的宽度进行调制,既可改变逆变电路输出电压的大小,也可改变输出频率。
五、实验过程1、测量波形（电路连接）并在示波器上显示:的转动情况。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
M A T L A B上机实习报告姓名：班级：指导老师：2009.12.13一、实验课程名称：数字信号处理二、实验目的掌握应用MATLAB 解决数字信号处理问题的方法，即熟悉应用MATLAB 处理数字信号处理问题相关的函数。
三、实验设备安装有MATLAB 软件的计算机四、实验内容编写MATLAB 程序，实现下列题目：实验一：设线性时不变（LTI ）系统的冲激响应为h(n)，输入序列为x(n)1、h(n)=(0.8)n ，0≤n ≤4； x(n)=u(n)-u(n-4)2、h(n)=(0.8)n u(n), x(n)=u(n)-u(n-4)3、h(n)=(0.8)n u(n), x(n)=u(n)求以上三种情况下系统的输出y(n)。
实验二：1、设有离散序列 x(n)=cos(0.48πn)+cos(0.52πn)分析下列三种情况下的幅频特性。
（1） 采集数据长度N=16，分析16点的频谱，并画出幅频特性。
采集数据长度N=16，并补零到64点，分析其频谱，并画出幅频特性。
（2） 采集数据长度N=64，分析64点的频谱，并画出幅频特性。
观察三幅不同的幅频特性图，分析和比较它们的特点及形成原因。
2、实现序列的内插和抽取所对应的离散傅里叶变换。
（选做）()()n R n n n x 128]365.1cos 36[cos ⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛=ππ 求()()n x n x 41=和()⎪⎭⎫ ⎝⎛=42n x n x 对应的128点傅里叶变换，比较三个结果所得到的幅频特性，分析其差别产生的原因。
实验三：1、设计一个Butterworth 数字低通滤波器，设计指标如下：通带截止频率：0.2π，幅度衰减不大于1分贝阻带截止频率：0.3π，幅度衰减大于15分贝2、让不同频率的正弦波通过滤波器，验证滤波器性能。
3、分析不同滤波器的特点和结果。
4、编程设计实现IIR 滤波器。
选取合适窗函数设计一个线性相位FIR低通滤波器，使它满足如下性能指标：通带截止频率：ωp=0.5π,通带截止频率处的衰减不大于3分贝；阻带截止频率：ωs=0.66π，阻带衰减不小于40分贝。
五、实验过程1、熟悉与离散信号频谱分析相关的MATLAB函数（参考老师提供的PPT）；2、通过运行附录2中提供的例题，熟悉用MATLAB解决数字信号处理问题的基本方法；3、根据实验内容中各个题目的要求，编写MATLAB程序代码，调试程序，分析并保存结果。
六、实验中的主要结论、实践技能和心得体会本次DSP实习，即MATLAB软件应用于处理数字信号处理问题。
通过这三次实习我对MATLAB 软件由完全不懂到现在能调用相关函数处理比较简单的问题，同时通过网络资源来帮助自己更好的利用该软件为解决现实中的相关问题服务。
在实习上机的过程中，就是对我们数字信号处理基础知识的巩固，虽然很多比较复杂的问题我们可以交给函数处理，但是对题目的理解及如何调用函数、处理的流程这些能力、知识全是我们课堂所学。
锻炼了我们理论联系实际，用理论武器解决现实问题的能力。
但是，要真正熟练掌握MATLAB软件还需要更多的实践。
七、实验结果对实验练习题编写MATLAB程序并运行，在计算机上输出仿真结果。
附录2 解题程序实验一：1、2小题相似，以下给出第2小题程序，该程序可以用于该类问题的求解，仅需改变相关系数即可。
Matlab程序：x=[ones(1,4)]; %输入为矩形脉冲序列，为一行4列向量N1=length(x);n1=0:N1-1; %序列长度为N1N2=100;n2=0:N2-1;h=0.8.^n2; %冲激响应N=N1+N2-1;n=0:N-1;y=conv(x,h); %调用卷积函数，x、h是参数subplot(311);stem(n1,x);subplot(312);stem(n2,h);subplot(313);stem(n,y);3、第3小题将N1、N2取100即可近似实现。
x=[ones(1,100)]; %输入为矩形脉冲序列，为一行4列向量N1=length(x);n1=0:N1-1; %序列长度为N1N2=100;n2=0:N2-1;h=0.8.^n2; %冲激响应N=N1+N2-1;n=0:N-1;y=conv(x,h); %调用卷积函数，x、h是参数subplot(311);stem(n1,x);subplot(312);stem(n2,h);subplot(313);stem(n,y);运行结果：实验二：1、本题即DFT求解，因无DFT函数，改为调用FFT函数求解。
N1=16;N2=64;n1=0:N1-1;n2=0:N2-1;x=cos(0.48*pi*n)+cos(0.52*pi*n);x1=(1:N1);Y1=fft(x1,N1);magY1=abs(Y1);k1=0:N1-1;w1=2*pi/N1*k1;x2=[x1 zeros(1,N2-N1)];Y2=fft(x2,N2);magY2=abs(Y2);k2=0:N2-1;w2=2*pi*k2/N2;Y3=fft(x,N2);magY3=abs(Y3);k3=0:N2-1;w3=2*pi/N2*k3;subplot(231);stem(n1,x1);title('没有足够点的采样信号');subplot(234);stem(w1/pi,magY1);title('信号的频谱');subplot(232);stem(n2,x2);title('添零信号');subplot(235);stem(w2/pi,magY2);title('高密度频谱');subplot(233);stem(n2,x2);title('有足够采样点的信号');subplot(236);stem(w3/pi,magY3);title('高分辨率频谱');试验结果：实验三：1、本题我参考了网上相似的资料，并根据自己在课堂上所学的Butterworth数字低通滤波器设计方法进行理解，基本掌握该设计方法。
另：该资料中不少函数用的是改进过的函数。
Matlab程序如下：wp=0.2*pi; %通带截止频率ws=0.3*pi; %阻带截止频率Rp=1; % 通带幅度衰减As=15; % 阻带幅度衰减Ripple=10^(-Rp/20);Attn=10^(-As/20);T=1; % 令抽样周期T=1OmegaP=wp/T; % 模拟域通带频率OmegaS=ws/T; % 模拟域阻带频率[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As); %系数按降幂排列[b,a]=impinvar(cs,ds,1/T);[C,B,A]=dir2par(b,a);[dBz,magz,phaz,grd,w]=freqz_m(b,a);subplot(2,2,1);plot(w/pi,magz); %数字滤波器幅度响应title('数字滤波器幅度响应');xlabel('Frequency in pi units');axis([0,1,0,1.2]);set(gca,'XTickMode','manual','XTick',[0,wp/pi,ws/pi,1]);set(gca,'YTickMode','manual','YTick',[Attn,Ripple,1]);grid;subplot(2,2,2);plot(w/pi,phaz/pi); %数字滤波器相位响应title('数字滤波器相位响应');xlabel('Frequency in pi units');ylabel('pi units');axis([0,1,-1.1,1.1]);set(gca,'XTickMode','manual','XTick',[0,wp/pi,ws/pi,1]); set(gca,'YTickMode','manual','YTick',[-1,-0.5,0,0.5,1]); grid;subplot(2,2,3);plot(w/pi,dBz); %数字滤波器幅度响应（dB）title('数字滤波器幅度响应（dB）');xlabel('Frequency in pi units');axis([0,1,-30,5]);set(gca,'XTickMode','manual','XTick',[0,wp/pi,ws/pi,1]); set(gca,'YTickMode','manual','YTick',[-As,-Rp,0]); grid;subplot(2,2,4);plot(w/pi,grd); %数字滤波器群延时title('数字滤波器群延时');xlabel('Frequency in pi units');ylabel('samples');axis([0,1,0,10]);set(gca,'XTickMode','manual','XTick',[0,wp/pi,ws/pi,1]); set(gca,'YTickMode','manual','YTick',[2,4,6,8,10]); grid;运行结果：Butterworth Filter Order= 64、IIR滤波器的Matlab程序如下：Wp=0.2*pi;Wt=0.3*pi;Ap=1;At=15;T=1;fs=1/T;wp=(2/T)*tan(Wp/2); %预畸变wt=(2/T)*tan(Wt/2); %预畸变[n,wc]=buttord(wp,wt,Ap,At,'s');%巴特沃斯滤波器阶数估计[z,p,k]=buttap(n); %巴特沃斯模拟低通滤波器原型[bp,ap]=zp2tf(z,p,k); %零极点增益模型转换为传递函数模型[bs,as]=lp2lp(bp,ap,wp); %改变低通模拟滤波器的截止频率[bz,az]=bilinear(bs,as,fs); %双线性变换[h,w]=freqz(bz,az); %计算数字滤波器的频率响应subplot(211);H=20*log10(abs(h));plot(w,H);title('数字滤波器频率响应');grid;subplot(212);pha=angle(h); %求相角plot(w,pha);title('数字滤波器相位响应');grid;运行结果：实验四：题目中要求阻带衰减不小于40分贝，故可以选择汉宁窗、海宁窗、布拉克曼窗，我们选择的是汉宁窗。