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数据存取实验报告
竭诚为您提供优质文档/双击可除dsp数据存取实验报告篇一：Dsp实验一数据存取实验《Dsp技术》课程实验报告学生姓名：所在班级：指导教师：记分及评价：一、实验名称Dsp数据存取实验二、实验目的（1）掌握Tms320F2812程序空间的分配（2）掌握Tms320F2812数据控件的分配三、实验内容（3）往0x003F9020地址开始的八个存储单元依次写入0-8的八个数（4）读取0x003F9020地址开始的八个存储单元内容并写入0x003F9028地址开始的八个存储单元内。
（5）从0x003F9020开始的八个存储单元的内容依次与0x003F9020地址开始的八个存储单元相乘，运算结果存入0x003F9000开始的八个存储单元内。
（6）从0x003F9020开始的八个存储单元的内容依次与0x003F9020地址开始的八个存储内容相加，运算结果依次存入0x003F9038地址开始的八个存储单元内。
四、实验程序与结果分析程序和结果如图1：图1结果voidmain(void){inti;volatileunsignedint*room=(volatileunsignedint*)0x3f 9020;volatileunsignedint*room2=(volatileunsignedint*)0x3 f902F;volatileunsignedint*room3=(volatileunsignedint*)0x3 f9030;volatileunsignedint*room4=(volatileunsignedint*)0x3 f903F;//Initializesystemcontrol://pLL,watchDog,enableperipheralclocksInitsysctrl();//DisablecpuinterruptsDInT;//Disablecpuinterruptsandclearallcpuinterruptflags: IeR=0x0000;IFR=0x0000;/*将0xAAAA写入从数据空间的地址0x3f9020开始的8个单元中*/for(i=0;i {*room=0x0000+i;room++;}/*从0x3f9020开始的8个空间读出数据依次写入从0x3f9028开始的8个单元中*/for(i=0;i {*room2=*(room-1);room--;room2--;}room2++;for(i=0;i {*room3=((*room)*(*room2));//0x003F9028开始的八个存储单元的内容依次与0x003F9030地址开始的八个存储单元相乘room++;room2++;room3++;}room--;room2--;for(i=0;i {*room4=((*room)+(*room2));//0x003F9028开始的八个存储单元的内容依次与0x003F9038地址开始的八个存储内容相加room--;room2--;room4--;}}五、小结通过本次实验，我学会了Tms320F2812的寻址方式，明白了试验箱扩展存储器空间的寻址方法以及ccs修改、填充Dsp内存单元的方法，加深了对于ccs2000软件的应用，为接下来的实验提供良好的帮助。
dsp课程设计实验报告
DSP课程设计实验语音信号的频谱分析:要求首先画出语音信号的时域波形, 然后对语音信号进行频谱分析。
在MATLAB中, 可以利用函数fft对信号进行快速傅立叶变换, 得到信号的频谱特性, 从而加深对频谱特性的理解。
其程序为:>> [y,fs,bits]=wavread('I:\xp.wav',[1024 5120]);>> sound(y,fs,bits);>> Y=fft(y,4096);>> subplot(221);plot(y);title('原始信号波形');>> subplot(212);plot(abs(Y));title('原始信号频谱');程序运行结果为:设计数字滤波器和画出频率响应:根据语音信号的特点给出有关滤波器的性能指标:低通滤波器性能指标, =1000Hz, =1200Hz, =100dB, =1dB；高通滤波器性能指标, =4800Hz, =5000Hz, =100dB, =1dB；带通滤波器性能指标, =1200Hz, =3000Hz, =1000Hz, =3200Hz, =100dB, =1dB；要求学生首先用窗函数法设计上面要求的三种滤波器, 在MATLAB中, 可以利用函数firl 设计FIR滤波器；然后再用双线性变换法设计上面要求的三种滤波器, 在MA TLAB中, 可以利用函数butte、cheby1和ellip设计IIR滤波器；最后, 利用MATLAB中的函数freqz画出各种滤波器的频率响应, 这里以低通滤波器为例来说明设计过程。
低通：用窗函数法设计的低通滤波器的程序如下:>> fp=1000;fc=1200;As=100;Ap=1;fs=22050;>> wc=2*fc/fs;wp=2*fp/fs;>> N=ceil((As-7.95)/(14.36*(wc-wp)/2))+1;>> beta=0.1102*(As-8.7);>> Win=Kaiser(N+1,beta);>>b=firl(N,wc,Win);>>freqz(b,1,512,fs);程序运行结果:这里选用凯泽窗设计, 滤波器的幅度和相位响应满足设计指标, 但滤波器长度（N=708）太长, 实现起来很困难, 主要原因是滤波器指标太苛刻, 因此, 一般不用窗函数法设计这种类型的滤波器。
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实验报告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实验报告总结【精品】
我不应把我的作品全归功于自己的智慧，还应归功于我以外向我提供素材的成千成万的事情和人物！——采于网，整于己，用于民2021年5月12日dsp实验报告总结篇一：dsp课程设计实验报告总结DSP课程设计总结（XX-XX学年第2学期）题目：专业班级：电子1103 学生姓名：万蒙学号：指导教师：设计成绩：XX 年6 月目录一设计目的----------------------------------------------------------------------3 二系统分析----------------------------------------------------------------------3 三硬件设计3.1 硬件总体结构-----------------------------------------------------------3 3.2 DSP模块设计-----------------------------------------------------------4 3.3 电源模块设计----------------------------------------------------------4 3.4 时钟模块设计----------------------------------------------------------5 3.5 存储器模块设计--------------------------------------------------------6 3.6 复位模块设计----------------------------------------------------------6 3.7 JTAG模块设计--------------------------------------------------------7 四软件设计4.1 软件总体流程-----------------------------------------------------74.2 核心模块及实现代码---------------------------------------8五课程设计总结-----------------------------------------------------14一、设计目的设计一个功能完备，能够独立运行的精简DSP硬件系统，并设计简单的DSP控制程序。
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开发板上进行实时处理。
五、实验结果与分析通过以上实验，我们成功完成了信号采集与重构、滤波器设计与实现、频谱分析与频域处理以及音频处理与效果实现等一系列实验。
dsp实验报告
dsp实验报告C54x的浮点数的算术运算（实验⼀）实验⽬的：1）了解TMS320C54x汇编语⾔程序的基本格式，以及汇编、链接的基本过程。
2）初步熟悉软件仿真器Simulator的⽤法。
实验内容：A.基础实验将两个⼩数相乘，分离尾数与指数，进⾏算术运算，最后归⼀化。
乘法运算时遵循指数相加尾数相乘的规则。
1)编写浮点乘法程序ex1.asm，完成X1*X2=0.3×（—0.8）运算。
2) 编写链接命令⽂件ex1.cmd。
B.提⾼实验：⽤汇编语⾔实现：y1=x1*a1-x2*a2y2=12/3+1y3=0.3*(－0.5)+1实验程序框图及清单：此次试验中，在⼀个程序⾥完成基础和提⾼部分，四个计算分别作为四个⼦程序。
具体程序及相关注释如下：.title "lab1.asm".mmregsSTACK .usect "STACK", 10H.bss x1,1 ;定点数操作数1.bss x2,1 ;定点数操作数2.bss e1,1 ;指数1.bss m1,1 ;尾数1.bss e2,1 ;指数2.bss m2,1 ;尾数2.bss ep,1 ;乘积的指数.bss mp,1 ;乘积的尾数.bss product,1 ;定点乘积.bss temp,1 ;暂存单元.bss a1,1.bss a2,1.bss b1,1.bss b2,1 ;提⾼实验y1=a1*b1-a2*b2.bss y1,1.bss c1,1 ;提⾼实验y2=12/3+1.bss c2,1.bss c3,1.bss y2,1.bss d1,1 ;提⾼实验y3=0.3*(－0.5)+1.bss d2,1.bss d3,1.bss y3,1.def start.datatable: .word 3*32768/10.word -8*32768/10.word 4*32768/10 ;.word 2*32768/10 ;提⾼实验y1=a1*b1-a2*b2数据.word 12,3,1 ;提⾼实验y2=12/3+1数据.word 3*32768/10.word -5*32768/10 ;提⾼实验y3=0.3*(－0.5)+1数据.word 1 .textstart: STM #0,SWWSRSTM #STACK+10H,SPMVPD table,@x1MVPD table+1,@x2LD @x1,16,A ;EXP AST T,@e1NORM ASTH A,@m1 ;将操作数x1转换成浮点数存⾄e1和m1中LD @x2,16,A ;EXP AST T,@e2NORM ASTH A,@m2 ;将操作数x2转换成浮点数存⾄e2和m2中CALL MULT ;调⽤浮点乘法⼦程序CALL LAB11 ;y1=x1*a1-x2*a2⼦程序CALL LAB12 ;y2=12/3+1⼦程序CALL LAB13 ;y3=0.3*(－0.5)+1⼦程序done: B doneMULT: SSBX FRCT ;设置⼩数相乘浮点乘法⼦程序SSBX SXM ;设置符号扩展LD @e1,A ;ADD @e2,A ;指数相加STL A,@ep ;LD @m1,T ;MPY @m2,A ;尾数相乘EXP A ;ST T,@temp ;NORM A ;将乘积浮点化STH A,@mp ;浮点乘积尾数存mp单元LD @temp,AADD @ep,ASTL A,@ep ;浮点乘积指数存ep单元NEG ASTL A,@tempLD @temp,TLD @mp,16,ANORM ASTH A,@productRETLAB11: STM #a1,AR1 y1=x1*a1-x2*a2⼦程序RPT #3MVPD table,*AR1+SSBX FRCTSSBX SXMLD @a1,TMPY @b1,ALD @a2,TMAS @b2,ASTH A,@y1RETLAB12: STM #c1,AR1 y2=12/3+1⼦程序RPT #2MVPD table+4,*AR1+RSBX FRCTRSBX SXMLD @c1,BRPT #15SUBC @c2,BADD @c3,BSTL B,@y2RETLAB13: STM #d1,AR1 y3=0.3*(－0.5)+1⼦程序RPT #2MVPD table+7,*AR1+SSBX FRCTSSBX SXMLD @d1,TMPY @d2,ALD A,-1,A ;q.15格式转换成q1.14格式ADD @d3,14,A ;与转换成q.14的1相加STH A,@y3 ;q1.14的结果存⼊y3RET.end实验中遇到的问题及解决⽅法主要是Dos命令的熟悉和语法错误的纠正运⾏结果：A.0x0148是基础部分的计算结果（q.15格式）换算回来是B.0x23d6是y1=x1*a1-x2*a2（0.3*0.4-（-0.8）*0.2）的结果0.28045654296875误差为0.28045654296875-0.28=0.00045654296875C.0x0005是y2=12/3+1的结果D.0x666是y3=0.3*(－0.5)+1的结果-0.150146484375误差为-0.150146484375-（-0.15）=-0.00014648375⼩结通过本次试验熟悉并掌握了部分dos命令以及dsp的⼤部分汇编指令，并能在dos调试环境下调程序，检查变量以验证结果正确与否。
dsp信号处理实验报告
dsp信号处理实验报告DSP信号处理实验报告一、引言数字信号处理（DSP）是一种将连续信号转换为离散信号，并对其进行处理和分析的技术。
在现代通信、音频处理、图像处理等领域中，DSP技术被广泛应用。
本实验旨在通过对DSP信号处理的实践，加深对该技术的理解与应用。
二、实验目的本实验旨在通过对DSP信号处理的实践，掌握以下内容：1. 学习使用DSP芯片进行信号采集和处理；2. 理解离散信号的采样和重构过程；3. 掌握常见的DSP信号处理算法和方法。
三、实验原理1. 信号采集与重构在DSP信号处理中，首先需要对模拟信号进行采样，将连续信号转换为离散信号。
采样过程中需要注意采样频率的选择，以避免混叠现象的发生。
采样完成后，需要对离散信号进行重构，恢复为连续信号。
2. DSP信号处理算法DSP信号处理涉及到多种算法和方法，如滤波、频谱分析、时域分析等。
其中，滤波是一种常见的信号处理方法，可以通过滤波器对信号进行去噪、增强等处理。
频谱分析可以将信号在频域上进行分析，了解信号的频率成分和能量分布。
时域分析则关注信号的时序特征，如幅值、相位等。
四、实验步骤1. 信号采集与重构在实验中，我们使用DSP芯片进行信号采集与重构。
将模拟信号输入DSP芯片的模拟输入端口，通过ADC（模数转换器）将模拟信号转换为数字信号。
然后，通过DAC（数模转换器）将数字信号转换为模拟信号输出。
2. 滤波处理为了演示滤波处理的效果，我们选择了一个含有噪声的信号进行处理。
首先，使用FIR滤波器对信号进行低通滤波，去除高频噪声。
然后，使用IIR滤波器对信号进行高通滤波，增强低频成分。
3. 频谱分析为了对信号的频率成分和能量分布进行分析，我们使用FFT（快速傅里叶变换）算法对信号进行频谱分析。
通过观察频谱图，可以了解信号的频率特性。
4. 时域分析为了对信号的时序特征进行分析，我们使用时域分析方法对信号进行处理。
通过计算信号的均值、方差、峰值等指标，可以了解信号的幅值、相位等特性。
硕士信号处理实验报告(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上机实验：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实验报告完美版
DSP实验报告班级：学号：姓名：指导教师：实验一、二 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。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验一 采样、系统性质及滤波一、实验目的和任务（1）熟悉MA TLAB 的主要操作命令。
（2）学会简单的矩阵输入和数据读写。
（3）掌握简单的绘图命令。
（4）用MATLAB 编程并学会创建函数。
（5）观察采样引起的混叠。
（6）判别离散时间系统的时不变性。
（7）卷积计算二、实验内容A 、观察采样引起的混叠设模拟信号为)3sin()2sin(4)5cos()(t t t t x πππ⋅+=，t 的单位为毫秒(ms)。
1. 设采样频率为3kHz ，确定与)(t x 混叠的采样重建信号)(t x a 。
2. 画出)(t x 和)(t x a 在)(60ms t ≤≤范围内的连续波形。
（因数字计算机无法真正画出连续波形，可用较密的离散点的连线来近似。
）3. 分别用"" 和""⨯在两信号波形上标记出3kHz 采样点。
两信号波形是否相同？采样后的两序列是否相同？实验过程与结果：1）实验程序clear% estimate x(t) and xa(t) with a much higher sampling freq. 'fs1'time_period=6; % unit: msfs1=50; % unit: kHzT1=1/fs1; % unit: msn1=0:fix(time_period/T1);x=cos(5*pi*n1*T1)+4*sin(2*pi*n1*T1).*sin(3*pi*n1*T1);xa=cos(pi*n1*T1);% obtain x(nT) and xa(nT) with given sampling freq. 'fs'fs=3;T=1/fs;n=0:fix(time_period/T);x_sample=cos(5*pi*n*T)+4*sin(2*pi*n*T).*sin(3*pi*n*T);xa_sample=cos(pi*n*T);figure,plot(n1*T1,x,'r',n1*T1,xa,'b',n*T,x_sample,'ro'),hold on, stem(n*T,xa_sample,'b:x')legend('x(t)','xa(t)','x(nT)','xa(nT)'),xlabel('t(ms)')2) 运行结果分析与讨论：)(t x 和)(t x a 两信号波形不同。
采样后的两序列相同。
B 、判别离散时间系统的时不变性。
设输入序列为)(n x ，系统)2()(n x n y =实现对)(n x 的抽取。
1. 设500,...,2,1),1002sin()(==n n n x π。
取延迟量D （例如D ＝30）。
记)()(D n x n x D -=，画出)(n x 、)(n x D 的序列波形。
2. 编程求出系统对)(n x 的响应)(n y 以及对)(n x D 的响应)(n y D3. 画出)(D n y -、)(n y D 的波形。
该系统是否为时不变的？实验过程与结果：1）实验程序clear% plot x(n) and x(n-D)D=30;N=500;n=1:N;x=sin(2*pi/100*n);for n=1:N+D,if (n-D)<=0, xD(n)=0;else xD(n)=x(n-D); endendfigure,subplot(2,1,1),plot(1:N,x,'r:',1:length(xD),xD,'b'),legend('x(n)','xD(n)'),,xlabel('n')% plot y(n) and yD(n)for n=1:fix(N/2)y(n)=x(2*n);endfor n=1:length(y)+D,if (n-D)<=0, y_delay(n)=0;else y_delay(n)=y(n-D); endendfor n=1:fix(length(xD)/2)yD(n)=xD(2*n);endsubplot(2,1,2),plot(1:length(y),y,'r:',1:length(y_delay),y_delay,'r.:',1:length(yD),yD,'b.'), legend('y(n)','y(n-D)','yD(n)'),,xlabel('n')axis([0 530 -1 1])2) 运行结果分析与讨论：该系统为时变系统。
)(D n y -= x(2n-D),)(n y D = x(2n-2D),所以)(D n y -不等于)(n y D 。
C 、利用卷积计算出输入信号通过FIR 滤波器的输出，并观察输出信号的input-on 暂态、input-off 暂态和稳态部分。
（计算卷积可用conv 命令）考虑下面两个滤波器，第一个的单位脉冲响应为⎩⎨⎧≤≤*=其它0140)75.0(25.0n h nn ，另一个的单位脉冲响应为]1,5,10,105,1[51--=，－h ；输入为周期方波，在一个周期内⎩⎨⎧≤≤≤≤=492502401)(x x n x 。
1） 分别画出两个滤波器的输出)(n y )2000(≤≤n 的波形，并与书上p144例4.1.8的两幅图比较是否一致。
2） 计算出图中稳态部分的)(n y 值。
实验过程与结果：1）实验程序clearh1=0.25*0.75.^(0:14);h2=1/5*[1 -5 10 -10 5 -1];N=200;n=0:N-1;x1=[ones(1,25) zeros(1,25)]; % one period of 'x(n)' x=[x1 x1 x1 x1];y1=conv(x,h1);y2=conv(x,h2);figure,subplot(2,1,1),,plot(n,x,'r:',n,y1(1:N),'b'), axis([0 200 -0.5 2.5]),grid on,legend('input','output'),xlabel('n')subplot(2,1,2),plot(n,x,'r:',n,y2(1:N),'b'),axis([0 200 -1.5 2.5]),grid on,legend('input','output'),xlabel('n')2) 运行结果实验二一、实验目的和任务（1）熟悉MA TLAB 的主要操作命令。
（2）学会简单的矩阵输入和数据读写。
（3）掌握简单的绘图命令。
（4）用MATLAB 编程并学会创建函数。
（5）观察序列频谱以及通过系统后波形与频谱的变化（6）根据正准型结构编写样本处理算法二、实验内容A 、观察序列频谱，观察信号通过系统后波形与频谱的变化已知输入信号)(2)(5)(21n x n x n x +=，其中)8cos()(1nn x π=，)4cos()(2nn x π=，1,0-=N n ，N 可取5000点。
1、画出)(),(),(21n x n x n x 的前100点波形2、画出)(n x 的DTFT 频谱)(ωX （πω≤≤0）3、某LTI 系统)8()()(--=n u n u n h ，画出系统的幅度频响|)(|ωH4、求系统)(n h 对)(n x 的响应)(n y (可以自己编程也可利用卷积函数)(⋅conv )。
画出)(n y 的波形,并与))4(arg 4cos(|)4(|2))8(arg 8cos(|)8(|5ππππππH n H H n H +++ 的波形比较（各画100点）；画出)(n y 的幅度谱|)(|ωY ，并与|)()(|ωωH X ⋅比较。
5、问：滤波器)(n h 是什么频响类型的滤波器？你从以上实验中观察到什么？与课本上的什么结论相吻合？实验过程与结果：1）实验程序colordef noneN=5000; % set length of the time seriesx1=cos(pi/8*(0:N-1));x2=cos(pi/4*(0:N-1));x=5*x1+2*x2;m=100; figure,plot(1:m,x1(1:m),1:m,x2(1:m),'g',1:m,x(1:m),'b.-'),legend('x1','x2','x')% --- plot X(w)M=5000; % set samples number in frequency domainw=pi/M*(0:M-1); % discreted frequencyj=sqrt(-1);Xw=zeros(1,M);for k=1:M, Xw(k)=sum(x*(exp(-j*w(k)*(0:N-1)'))); endfigure,plot(w/pi,abs(Xw)),xlabel('\omega(\pi)'),ylabel('|X(\omega)|')% --- plot H(w)h=ones(1,8);for k=1:M, Hw(k)=sum(h*(exp(-j*w(k)*(0:7)'))); endfigure,plot(w/pi,abs(Hw)),xlabel('\omega(\pi)'),ylabel('|H(\omega)|')% --- plot Y(w)y=conv(x,h); % calculate the responseH1=sum(h*(exp(-j*pi/8*(0:7)')));H2=sum(h*(exp(-j*pi/4*(0:7)')));m=100;yy=5*abs(H1)*cos(pi/8*(0:m-1)+angle(H1))+2*abs(H2)*cos(pi/4*(0:m-1)+angle(H2)); figure,plot(y(1:m),'r.-'),hold on,plot(yy,'g'),legend('y','yy')Ny=length(y);for k=1:M, Yw(k)=sum(y*(exp(-j*w(k)*(0:Ny-1)'))); endfigure,plot(w/pi,abs(Yw),'r.-') ,hold on,plot(w/pi,abs(Xw.*Hw)),xlabel('\omega(\pi)'),legend('|Y(\omega)|','|X(\omega)|*|H(\omega)|')2) 运行结果分析与讨论：矩形窗口滤波器，从实验中可以观察到时域卷积与频域相乘的图形基本相同，与课本上结构 基本吻合。