DSP实习实验报告

合集下载

dsp实验报告总结doc

dsp实验报告总结doc
3.1 硬件总体结构
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原理与应用实验报告总结DSP(Digital Signal Processing)数字信号处理是利用数字技术对信号进行处理和分析的一种方法。

在本次实验中,我们探索了DSP的原理和应用,并进行了一系列实验以验证其在实际应用中的效果。

以下是对实验结果的总结与分析。

实验一:数字滤波器设计与性能测试在本实验中,我们设计了数字滤波器,并通过性能测试来评估其滤波效果。

通过对不同类型的滤波器进行设计和实现,我们了解到数字滤波器在信号处理中的重要性和应用。

实验二:数字信号调制与解调本实验旨在通过数字信号调制与解调的过程,了解数字信号的传输原理与方法。

通过模拟调制与解调过程,我们成功实现了数字信号的传输与还原,验证了调制与解调的可行性。

实验三:数字信号的傅里叶变换与频谱分析傅里叶变换是一种重要的信号分析方法,可以将信号从时域转换到频域,揭示信号的频谱特性。

本实验中,我们学习了傅里叶变换的原理,并通过实验掌握了频谱分析的方法与技巧。

实验四:数字信号的陷波滤波与去噪处理陷波滤波是一种常用的去除特定频率噪声的方法,本实验中我们学习了数字信号的陷波滤波原理,并通过实验验证了其在去噪处理中的有效性。

实验五:DSP在音频处理中的应用音频处理是DSP的一个重要应用领域,本实验中我们探索了DSP在音频处理中的应用。

通过实验,我们成功实现了音频信号的降噪、均衡和混响处理,并对其效果进行了评估。

实验六:DSP在图像处理中的应用图像处理是另一个重要的DSP应用领域,本实验中我们了解了DSP在图像处理中的一些基本原理和方法。

通过实验,我们实现了图像的滤波、边缘检测和图像增强等处理,并观察到了不同算法对图像质量的影响。

通过以上一系列实验,我们深入了解了DSP的原理与应用,并对不同领域下的信号处理方法有了更深刻的认识。

本次实验不仅加深了我们对数字信号处理的理解,也为日后在相关领域的研究与实践提供了基础。

通过实验的结果和总结,我们可以得出结论:DSP作为一种数字信号处理的方法,具有广泛的应用前景和重要的实际意义。

DSP实验报告6

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实验报告+心得体会

DSP实验报告+心得体会
// For this example use the following configuration:
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实验报告――离散信号的产生及运算(精)

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实验报告_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实验报告(综合)

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实验报告(一)

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倍,则采样以后不会发生频谱混淆现象。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 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分贝,故可以选择汉宁窗、海宁窗、布拉克曼窗,我们选择的是汉宁窗。

相关文档
最新文档