北京理工大学DSP实验一
DSP实验报告
实验一 采样、系统性质及滤波一、实验目的和任务（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 两信号波形不同。
dsp实验注释
dsp实验注释实验一：简单指令程序运行实验代码（含注释）：.mmregs.global _main_main:stm #3000h,sp ;sp为堆栈指针寄存器，stm为存储器映射寄存器寻址ssbx xf ;xf赋值为1，灯亮call delay ;调用delay函数，延迟0.5秒rsbx xf ;xf赋值为0，灯灭call delay ;调用delay函数b _main ;无条件调用_main函数nopnop;delay .5 second ;延时5秒delay:stm 270fh,ar3 ;ar3赋值207fh，十进制为9999dloop1:stm 0f9h,ar4 ;ar4赋值0f9h，十进制为249dloop2:banz loop2,*ar4- ;若不为0，ar4减1banz loop1,*ar3- ;若不为0，ar3减1，共进行10000*250次跳转retnopnop.end实验操作：可见XF灯以一定频率闪烁；单击“Halt”暂停程序运行，则XF 灯停止闪烁，如再单击“Run”，则“XF”灯又开始闪烁。
实验二：资料存储实验本实验程序将对0x1000开始的8个地址空间，填写入0xAAAA 的数值，然后读出，并存储到0X1008开始的8个地址空间。
在CCS 中可以观察DATA内存空间地址0X1000~0X100F值的变化。
代码（含注释）：.mmregs.global _main_main:;store data ;存储数据stm 1000h,ar1 ;将立即数1000h送入辅助寄存器ar1 (内存地址) rpt #07h ;循环执行下一条指令，8次st 0aaaah,*ar1+ ;将立即数0xaaaah赋给辅助寄存器ar1的1000h地址内，;之后ar1的地址加1，8次；ar1的地址变1008，内存0x1000-- ox1007中的数据均为0xaaaah;read data then re-store ;读入数据重新存储stm 7h,ar3 ;令辅助寄存器ar3的初值为07hstm 1000h,ar1 ;重新将立即数1000h送入辅助寄存器ar1stm 1008h,ar2 ;将立即数1008h送入辅助寄存器ar2loop: ;进入循环ld *ar1+,t ;将辅助寄存器ar1的值0xaaaah存入T寄存器中，且ar1的地址加1 st t,*ar2+ ;将T寄存器内容0xaaaah存入辅助寄存器ar2 ，并且ar2的地址加1 banz l oop,*ar3- ;寄存器ar3的值不为0时，执行循环loop，ar3的地址值减1直至为0时退出循环here:b here.end ;结束（地址0x1000-- 0x100F，程序完成对16个内存单元赋值）实验操作：1.用“View”下拉菜单中的“Memory”查看内存单元；2.输入要查看的内存单元地址，本实验要查看0x1000H~0x100FH单元的数值变化，输入地址0x1000H；3.查看0x1000H~0x100FH单元的初始值，单击“Run”运行程序，也可以“单步”运行程序；4.单击“Halt”暂停程序运行；5.查看0x1000H~0x100FH单元内数值的变化；6.关闭各窗口，本实验完毕。
北理工数字信号处理实验报告
本科实验报告实验名称：数字信号处理实验实验1 利用DFT 分析信号频谱一、实验目的1.加深对DFT 原理的理解。
2.应用DFT 分析信号频谱。
3.深刻理解利用DFT 分析信号频谱的原理，分析现实过程现象及解决办法。
二、实验原理1、DFT 和DTFT 的关系有限长序列()x n 的离散时间傅里叶变换()j X e ω在频率区间(02)ωπ≤≤的N个等分点{(0),(1),(),(1)}x x x k x N-……上的N 个取样值可以由下式表示：2120()|()()01(21)N jkn j Nk k X e x n eX k k N πωωπ--====≤≤--∑由上式可知，序列()x n 的N 点DFT ()X k ，实际上就是()x n 序列的DTFT 在N 个等间隔频率点{(0),(1),(),(1)}X X X k X N -……上样本()X k 。
2、利用DFT 求DTFT方法1：由()X k 恢复出()j X eω的方法如图2.1所示：图 2.1.由 N 点DFT 恢复频谱DTFT 的流程由图2.1所示流程图可知：01()()()(22)j j nkn j n N n n k X e x n eX k W e N ωωω∞∞∞---=-∞=-∞=⎡⎤==-⎢⎥⎣⎦∑∑∑ 由式2-2可以得到12()()()(23)Nj k kx e X k N ωπφω==--∑其中()x φ为内插函数12sin()2()(24)sin()2N j N e N ωωφωω--=•-方法2：然而在实际MATLAB 计算中，上诉插值公式不见得是最好的方法。
由于DFT 是DTFT 的取样值，其相邻的两个频率样本点的间距为2Nπ,所以如果我们增加数据的长度N ，使得得到的DFT 谱线就更加精细，其包络就越接近DTFT 的结果，这样可以利用DFT 来近似计算DTFT 。
如果没有更多的数据，可以通过补零来增加数据长度。
DSP智能车制作实验指导书
智能车控制大赛实验指导书北京理工大学.自动化学院2009-5-11制作人：张鑫、黄晓江、韩立博目录前言 ----------------------------------------------------------------------------------------------- 3智能车比赛介绍 -------------------------------------------------------------------------------- 4组队参赛与分工 -------------------------------------------------------------------------------- 5比赛器材与工具使用 -------------------------------------------------------------------------- 6智能车硬件制作 ------------------------------------------------------------------------------ 11 硬件设计中问题解决方案 ------------------------------------------------------------------ 21 DSP系统使用方法 --------------------------------------------------------------------------- 23一、组成和特点 ------------------------------------------------------------------------- 23二、使用注意事项-------------------------------------------------------------------------- 32三、编程指南-------------------------------------------------------------------------------- 33 软件调试方法和常见问题 ------------------------------------------------------------------ 44 智能车实战注意 ------------------------------------------------------------------------------ 54 技术文档制作要求 --------------------------------------------------------------------------- 55 参赛心得 --------------------------------------------------------------------------------------- 56 智能车传承与后记 --------------------------------------------------------------------------- 58前言学习制作和设计智能车是目前自动化类学生入门最简单的小型系统设计案例。
北京理工大学通信电路软件实验
本科实验报告实验名称：通信电路与系统软件实验实验一简单基带传输系统分析一、分析内容构造一个简单示意性基带传输系统。
以双极性PN码发生器模拟一个数据信源，码速率为100bit/s，低通型信道噪声为加性高斯噪声（标准差＝0.3v）。
要求：1．观测接收输入和滤波输出的时域波形；2．观测接收滤波器输出的眼图。
二、分析目的掌握观察系统时域波形，特别是眼图的操作方法。
三、系统组成及原理简单的基带传输系统原理框图如图2-1-1所示，该系统并不是无码间干扰设计的，为使基带信号能量更为集中，形成滤波器采用高斯滤波器。
图2-1-1 简单基带传输系统组成框图四、实验步骤及实验结果第1步：进入SystemView系统视窗，设置“时间窗”参数如下：①运行时间：Start Time: 0秒；Stop Time: 0.5秒；②采样频率：Sample Rate：10000Hz。
第2步：调用图符块创建如图2-1-2所示的仿真分析系统：图2-1-2 创建的简单基带传输仿真系统系统中各图符块的设置如表2-1-1所示：其中，Token1为高斯脉冲形成滤波器；Token3为高斯噪声产生器，设标准偏差Std Deviation=0.3v，均值Mean=0v；Token4为模拟低通滤波器，来自选操作库中的“LinearSys”图符按钮，在设置参数时，将出现一个设置对话框，在“Design”栏中单击Analog…按钮，进一步单击“Filter PassBand”栏中Lowpass 按钮，选择Butterworth型滤波器，设置滤波器极点数目：No.of Poles=5（5阶），设置滤波器截止频率：LoCuttoff=200 Hz。
第3步：单击运行按钮，运算结束后按“分析窗”按钮，进入分析窗后，单击“绘制新图”按钮，则Sink9～Sink12显示活动窗口分别显示出“PN码输出”、“信道输入”、“信道输出”和“判决比较输出”时域波形，如图2-1-2所示：图2-1-2（a）代表信源的PN码输出波形图2-1-2（b）经高斯脉冲形成滤波器后的码序列波形图2-1-2（c）信道输出的接收波形图2-1-2（d）判决比较输出波形第4步：观察信源PN码和波形形成输出的功率谱。
《DSP技术与实验》实验指导书
2011秋季学期“DSP 原理与应用”实验教学资料实验一 CCS 集成环境的熟悉 ---2学时 实验二 算术运算与数据操作的汇编语言程序设计---2学时 实验三 信号发生器设计 ---2学时 实验四 IIR 滤波器的DSP 实现 ---2学时 实验五 FIR 滤波器的DSP 实现 ---2学时 实验六 FFT 的DSP 实现 ---2学时 实验七 综合实验暨实验考查 ---4学时实验报告内容：用每个实验要求的思考题代替实验报告，提交电子版，在实验结束时一起交。
实验一 CCS 集成开发环境的熟悉一、实验内容及步骤通过编写程序计算y=x1+x2+x3+x4+x5练习CCS 2 (…C5000)的使用。
1．在D 盘或E 盘建立一个自己的文件夹, 如e:\05136101,该文件夹用英文或数字命名。
后面所有与实验内容有关的工程与程序都在此文件夹下保存或操作。
2．运行CCS 2(‘C5000)进入CCS 集成开发环境。
以下操作如无特殊声明一般都在此环境下操作。
3．新建工程。
点击主菜单Project-New ，会出现图1-1所示的Project Creation 对话框。
参照图1-1按如下顺序进行设置或操作： Location ：e:\05136101 Project ：example2图1-1点击完成按钮即可。
4．新建并保存主程序文件、链接命令文件和中断向量表文件。
本次实验三个文件的参考名称为：主程序文件example2.asm 、链接命令文件example2.cmd 、中断向量表文件vectors.asm 。
点击菜单File-New-Source File 打开文本编辑窗口Untitled1，如图1-2。
用菜单File-Save 功能将其保存到e:\05136101\example2中，其中文件名称为example2，保存类型为*.asm ，如图1-3。
用同样的方法新建并保存文件example2.cmd 和vectors.asm 。
DSP技术 实验一报告模板
《DSP技术》课程实验报告学生姓名：所在班级：指导教师：记分及评价：一、实验名称：数据存储实验二、实验目的掌握TMS320F2812的内部结构；熟悉CCS软件的使用三、实验内容修改例子程序，实现以下功能：（1）往0x003F9020地址开始的八个存储单元依次写入0，1，2，……，7八个数；读取0x003F9020地址开始的八个存储单元内容依次写入0x003F9028地址开始的八个存储单元内；（2）从0x003F9020地址开始的八个存储单元的内容依次与0x003F9028地址开始的八个存储单元内容相乘，运算结果依次存入0x003F9030地址开始的八个存储单元内；从0x003F9020地址开始的八个存储单元的内容依次与0x003F9028地址开始的八个存储单元内容相加，运算结果依次存入0x003F9038地址开始的八个存储单元内；（3）将0x003F9020地址开始的32个存储单元的数据从小到大排序，排序后的结果依次存入0x003F9040地址开始的32个存储单元内四、实验程序与结果分析程序：void main(void){int i,n;volatile unsigned int *room = (volatile unsigned int *)0x3f9020;volatile unsigned int *room2= (volatile unsigned int *)0x3f9030;volatile unsigned int *room3= (volatile unsigned int *)0x3f9040;volatile unsigned int *room4= (volatile unsigned int *)0x3f9090;volatile unsigned int *room5= (volatile unsigned int *)0x3f90A0;volatile unsigned int *temp= (volatile unsigned int *)0x3f9090;unsigned int max;// Initialize System Control:// PLL, WatchDog, enable Peripheral ClocksInitSysCtrl();// Disable CPU interruptsDINT;// Disable CPU interrupts and clear all CPU interrupt flags: IER = 0x0000;IFR = 0x0000;/*对相应的内存空间赋值*/for(i=0;i<8;i++){* room=i;*room4=i;*(room+8)=i;*(room4+8)=i;*room2=(* room)*(*(room+8));*room5=(* room)*(*(room+8));*(room2+8)=(* room)+(*(room+8));*(room5+8)=(* room)+(*(room+8));room++;room2++;room4++;room5++;}/*从小到大排序并将32个数存在room3中*/for(i=0;i<31;i++){for(n=i+1;n<32;n++){if(temp[i]>=temp[n]){max= temp[i];temp[i]=temp[n];temp[n]=max;}}room3[i]=temp[i];}room3[i]=temp[31];}结果：五、小结通过本次对DSP的初次编程，感觉对C语言要求较高，进一步了解了DSP的内存分布情况，基本掌握了CSS的使用。
DSP实验报告模版
我们做的dsp实验是实验一，实验二，实验三！模板上只有1 3 ！2要自己仿照实验一写！不要弄错了！实验题目在另外一个共享表格里！实验一数据存储实验一、实验目的1. 掌握ccs的使用2. 掌握 tms320c54x 程序空间的分配；3. 掌握 tms320c54x 数据空间的分配；4. 能够熟练运用tms320c54x 数据空间的指令。
二、实验设备计算机，ccs 3.1版软件，dsp仿真器，e300实验箱，dsp-54xp cpu板。
三、实验步骤与内容1. 在进行 dsp实验之前，需先连接好仿真器、实验箱及计算机，连接方法如下所示：2. e300 底板的开关sw4 的第1位置on，其余位置off。
其余开关设置为off。
sw5全部置on；其余开关不做设置要求3. 上电复位在硬件安装完成后，确认安装正确、各实验部件及电源连接无误后，启动计算机，接通仿真器电源，此时，仿真器上的“红色指示灯”应点亮，否则dsp开发系统与计算机连接存在问题。
4. 运行ccs程序1) 待计算机启动成功后，实验箱220v电源置“on”，实验箱上电 2) 启动ccs3.1，进入ccs界面后，点击“debug—connect”3) 此时仿真器上的“绿色指示灯”应点亮，ccs正常启动，表明系统连接正常；否则仿真器的连接、jtag 接口或ccs 相关设置存在问题，这时需掉电检查仿真器的连接、jtag 接口连接是否正确，或检查ccs相关设置是否存在问题。
5. 成功运行ccs 程序后，首先应熟悉ccs的用户界面；6. 学会在ccs环境下创建工程文件、添加程序文件、编写程序、编译、装载、调试，学习如何使用观察窗口等。
7. 用“project\open”打开“c:\ti5000\myprojects\01_mem\ mem.pjt”.编译并装载“\ 01_mem\debug\mem.out”8.用“edit”下拉菜单中的“memory/fill”编辑内存单元，参数设置如下图：单击“ok”此时以0x1000 为起始地址的16个内存单元被修改成：0x00099.用“view”下拉菜单“memory”观察内存单元变化，输入要查看的内存单元地址，本实验要查看0x1000h～0x100fh 单元的数值变化，输入地址0x1000h；单击“ok”如下图所示：10. 点击“debug\go main”进入主程序，在程序中“加软件断点1”和“加软件断点2”处施加软件断点。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
本科实验报告实验名称：实验一利用DFT对信号进行频谱分析实验一 利用DFT 对信号进行频谱分析一、 实验目的利用DFT 对信号进行频谱分析，研究不同数据长度、补零、加窗等对频率分辨率的影响。
二、 实验原理1.连续周期信号相对于离散周期信号，连续非周期信号相对于离散非周期信号，都可以通过时域抽样定理建立相互关系。
因此，在离散信号的DFT 分析方法基础上，增加时域抽样的步骤，就可以实现连续信号的DFT 分析。
2.利用DFT 计算连续周期信号X （t ）的频谱的分析步骤：(1) 确定周期信号的基本周期T0；(2) 计算一个周期内的抽样点数N 。
若周期信号的最高次谐频为p 次谐波pw0 ，则频谱中有2p+1根谱线；若周期信号的频谱无限宽，则认为集中信号90%以上(或根据工程允许而定)能量的前(p+1)次谐波为近似的频谱范围，其余谐波忽略不计。
取N>=2p+1；(3) 对连续周期信号以抽样间隔T= T0 /N 进行抽样，得到x[k] ；(4) 利用FFT 函数对x[k]作N 点FFT 运算，得到X[m]；(5) 最后求得连续周期信号的频谱为X(nw0)=X[m]/N 。
三、 实验内容与分析对123πππ+x(t ）=cos(2ft)+cos(2f t)cos(2f t) ，其中123f 2,f 2.5,f 3,t ms kHz kHz kHz ===的单位为 进行频谱分析：1.由f 3,f 2h s h kHz f =≥ ，选择10kHz 的采样频率对x （t ）进行时域采样。
2．由x(t)可知，min f 0.5k Hz ∆= ，下面取不同数据长度对x （t ）进行频谱分析：（1） 对连续时间信号采样后取L1=10、L2=20两个数据长度；（2）采样点数L1=10时，分别取N1=32、N2=64；分析：L1=10时，min f 10.5s f kHz f kHz L∆==>∆=，因此截断对应的10点样值序列x1(n)的频率分辨率较差。
对此信号进行N 点的DFT 后，频率分辨率分别为 N=32时，min f (Hz)0.3125(kHz)f s bin f N ∆==<∆ N=64时，min f (Hz)0.15625(kHz)f s bin f N∆==<∆ 其数值都小于原模拟信号x （t ）的最小频率间隔，但由于数据宽度的限制DFT 结果并不能区分实际信号的三个频率分量。
（3）采样点数L2=20时，分别取N1=32、N2=64；分析：L2=20时，min f 0.5=0.5s f kHz f kHz L∆==∆=，因此截断对应的20点样值序列x2(n)恰能区分x(t)的三个频率分量。
对此信号进行N 点的DFT 后，频率分辨率分别为 N=32时，min f (Hz)0.3125(kHz)f s bin f N ∆==<∆ N=64时，min f (Hz)0.15625(kHz)f s bin f N∆==<∆ 其数值都小于原模拟信号x （t ）的最小频率间隔，选用大点数的N 的DFT ，频率分辨率将提高，可以减小频谱误差。
（4）在L2=20 N1=32的数据长度下进行补零操作分析：补零可以克服栅栏效应，使N为2的整数倍便于FFT 计算。
（5）对在L=20，N=1024的采样信号加窗函数分析：窗函数加权可抑制由信号截取产生的频谱泄露问题，用不同的窗函数对单位取样响应进行加权可以得到等效滤波器不同的振幅频率特性。
四、实验总结与思考(1)提高频谱分辨率的方法：答：频率分辨率为将信号中两个靠的很近的谱峰区分开的能力，f 1s pF N T == 其中N 为采样点数，f s 为采样频率，P T 为采样间隔。
信号长度越长，频率分辨率越好。
所以提高真实信号的长度有助于提高频率分辨率。
(2)常用的窗函数的优缺点：答：由上图可以看出，矩形窗有最窄的主瓣，但是旁瓣泄露严重。
汉宁窗和汉明窗虽主瓣较宽，但是旁瓣泄露少，是常选用的窗函数。
(3)序列补零对频谱分析结果的影响。
答：在序列后补零直接的影响就是增加了序列的长度。
但是却提高了频谱分析的精度。
因为序列补零后，序列长度增加了，由于抽样频率没有改变，因此频谱图中谱线之间的间隔变小，从而显示出了更多的细节，提高了频谱分析精度。
但需要注意的是，实际信号长度并没有增加，频率分辨率仍保持不变。
五、实验代码clear all;freq1 = 2;freq2=2.5;freq3=3.;fs=10;omga1=2*pi*freq1;omga2=2*pi*freq2;omga3=2*pi*freq3;L1=10;L2=20;N1=32;N2=1024;dt=1./fs;ddt=dt./100;f=[freq1 freq2 freq3];z=[0 0 0];zz=[1 1 1];t1=0:ddt:(L1-1)*dt;tt1=0:ddt:(L2-1)*dt;xt1=cos(omga1*t1)+cos(omga2*t1)+cos(omga3*t1);xt2=cos(omga1*tt1)+cos(omga2*tt1)+cos(omga3*tt1);k1=0:dt:(L1-1)*dt;kk1=0:dt:(L2-1)*dt;x1=cos(omga1*k1)+cos(omga2*k1)+cos(omga3*k1);x2=cos(omga1*kk1)+cos(omga2*kk1)+cos(omga3*kk1);y1=fft(x1,N1);i1=0:N1-1;y12=fft(x1,N2);i2=0:N2-1;y2=fft(x2,N1);y22=fft(x2,N2);%连续时间信号采样后，不同数据长度时域波形subplot(211);plot(10*t1,xt1,'--');hold on;stem(0:L1-1,x1,'b');hold offtitle('抽样信号x1（t）L1=10');xlabel('n');ylabel('x1(n)');axis([0 L2 min(x1) max(x1)]);subplot(212)plot(10*tt1,xt2,'--');hold on;stem(0:L2-1,x2,'b');hold offtitle('抽样信号x2（t）L2=20');xlabel('n');ylabel('x2(n)');axis([0 L2 min(x2) max(x2)]);%L2=10 N1=32 N2=64subplot(211);stem(i1,abs(y1));holdon;stem(f*N1/fs,z,'xr');hold offtitle('x1(n),L1=10,N1=32');xlabel('k');ylabel('|x1(k)| ');axis([0 max(i1) min(abs(y1)) max(abs(y1))]);subplot(212)stem(i2,abs(y12));hold on;stem(f*N2/fs,z,'xr');hold off title('x2(n),L1=10,N2=64');xlabel('k');ylabel('|x2(k)| ');axis([0 max(i2) min(abs(y12)) max(abs(y12))]);%L2=20 N1=32 N2=64subplot(211);stem(i1,abs(y2));holdon;stem(f*N1/fs,z,'xr');hold offtitle('x1(n),L2=20,N1=32');xlabel('k');ylabel('|x1(k)| ');axis([0 max(i1) min(abs(y2)) max(abs(y2))]);subplot(212)stem(i2,abs(y22));hold on;stem(f*N2/fs,z,'xr');hold off title('x2(n),L2=20,N2=64');xlabel('k');ylabel('|x2(k)| ');axis([0 max(i2) min(abs(y22)) max(abs(y22))]);%补零操作y2_zero1=fft(x2,64);y2_zero2=fft(x2,128);y2_zero3=fft( x2,1024);subplot(411);stem(abs(y2));title('原谱线');axis([0 max(i1) min(abs(y2)) max(abs(y2)+2)]);subplot(412);stem(abs(y2_zero1));title('补32个零'); axis([0 max(i1) min(abs(y2)) max(abs(y2)+2)]);subplot(413);stem(abs(y2_zero2));title('补96个零'); axis([0 max(i1) min(abs(y2)) max(abs(y2)+2)]);subplot(414);plot(abs(y2_zero3));title('补1024个零');%加窗函数w_ham=0.54-0.46*cos(2*pi*kk1/(N2+1));x2_ham=x2.*w_ham;y22_ham=fft(x2_ham,N2);w_bar=1.0-2*abs(kk1)/(N2+1);x2_bar=x2.*w_bar;y22_bar=fft(x2_bar,N2);w_han=(cos(kk1*pi/(N2+1))).^2;x2_han=x2.*w_han;y22_han=fft(x2_han,N2);subplot(311);plot(abs(y22));title('原谱线');subplot(312);plot(abs(y22_han));title('加汉宁窗'); subplot(313);plot(abs(y22_ham));title('加海明窗');。