南邮DSPA实验报告
南邮DSP数字信号处理实验报告 可直接打印
南京邮电大学实验报告实验名称熟悉MATLAB环境快速傅里叶变换(FFT)及其应用 IIR数字滤波器的设计FIR数字滤波器的设计课程名称数字信号处理A班级学号_________________________ 姓名_________________________开课时间 2013/2014学年，第二学期实验一熟悉MATLAB环境一、实验目的（1）熟悉MA TLAB的主要操作命令。
（2）学会简单的矩阵输入和数据读写。
（3）掌握简单的绘图命令。
（4）用MATLAB编程并学会创建函数。
（5）观察离散系统的频率响应。
二、实验内容(1) 数组的加、减、乘、除和乘方运算。
输入A=[1 2 3 4]，B=[3,4,5,6]，求C=A+B，D=A-B，E=A.*B,F=A./B,G=A.^B 。
并用stem语句画出A、B、C、D、E、F、G。
clearn = 0:1:3;A=[1 2 3 4];subplot(3,1,1)stem(n,A)xlabel('n')ylabel('A')B=[3,4,5,6];subplot(3,1,2)stem(n,B)xlabel('n')ylabel('B')C=A+B;subplot(3,1,3)stem(n,C)xlabel('n')ylabel('C')figureD=A-B;E=A.*B;F=A./B;G=A.^B;subplot(4,1,1)stem(n,D)xlabel('n')ylabel('D')subplot(4,1,2)stem(n,E)xlabel('n')ylabel('E')subplot(4,1,3) stem(n,F) xlabel('n') ylabel('F') subplot(4,1,4) stem(n,G) xlabel('n') ylabel('G')(2) 用MATLAB 实现下列序列： a) 08(). 0n 15n x n =≤≤ n=0:1:15; x1=0.8.^n; stem(n,x1) xlabel('n') ylabel('x(n)') title('2(a)')b) 023(.)() 0n 15j nx n e+=≤≤n=0:1:15; i=sqrt(-1); a = 0.2+3*i; x2=exp(a*n); figuresubplot(1,2,1) stem(n,real(x2)) xlabel('n')ylabel('x(n)实部') subplot(1,2,2) stem(n,imag(x2))xlabel('n')ylabel('x(n)虚部')c) 3012502202501()cos(..)sin(..)x n n n ππππ=+++ 0n 15≤≤ n=0:1:15;x3=3*cos(0.125*pi*n+0.2*pi) + 2*sin(0.25*pi*n+0.1*pi); stem(n,x3) xlabel('n') ylabel('x(n)')(3) 绘出下列时间函数的图形，对x 轴、y 轴以及图形上方均须加上适当的标注： a) 2()sin() 0t 10s x t t π=≤≤ t=0:0.001:10; x=sin(2*pi*t); plot(t,x,'r-')xlabel('t'),ylabel('x(t)'),title('sin(2\pit)')b) 100()cos()sin() 0t 4s x t t t ππ=≤≤>> n=0:0.01:4;x=cos(100*pi*n).*sin(pi*n); subplot(2,1,2);plot(n,x2t);(4) 给定一因果系统12121106709()()/(..)H z z z z ----=+-+，求出并绘制H (z )的幅频响应和相频响应。
dsp原理与应用实验报告总结
dsp原理与应用实验报告总结DSP（Digital Signal Processing）数字信号处理是利用数字技术对信号进行处理和分析的一种方法。
在本次实验中，我们探索了DSP的原理和应用，并进行了一系列实验以验证其在实际应用中的效果。
以下是对实验结果的总结与分析。
实验一：数字滤波器设计与性能测试在本实验中，我们设计了数字滤波器，并通过性能测试来评估其滤波效果。
通过对不同类型的滤波器进行设计和实现，我们了解到数字滤波器在信号处理中的重要性和应用。
实验二：数字信号调制与解调本实验旨在通过数字信号调制与解调的过程，了解数字信号的传输原理与方法。
通过模拟调制与解调过程，我们成功实现了数字信号的传输与还原，验证了调制与解调的可行性。
实验三：数字信号的傅里叶变换与频谱分析傅里叶变换是一种重要的信号分析方法，可以将信号从时域转换到频域，揭示信号的频谱特性。
本实验中，我们学习了傅里叶变换的原理，并通过实验掌握了频谱分析的方法与技巧。
实验四：数字信号的陷波滤波与去噪处理陷波滤波是一种常用的去除特定频率噪声的方法，本实验中我们学习了数字信号的陷波滤波原理，并通过实验验证了其在去噪处理中的有效性。
实验五：DSP在音频处理中的应用音频处理是DSP的一个重要应用领域，本实验中我们探索了DSP在音频处理中的应用。
通过实验，我们成功实现了音频信号的降噪、均衡和混响处理，并对其效果进行了评估。
实验六：DSP在图像处理中的应用图像处理是另一个重要的DSP应用领域，本实验中我们了解了DSP在图像处理中的一些基本原理和方法。
通过实验，我们实现了图像的滤波、边缘检测和图像增强等处理，并观察到了不同算法对图像质量的影响。
通过以上一系列实验，我们深入了解了DSP的原理与应用，并对不同领域下的信号处理方法有了更深刻的认识。
本次实验不仅加深了我们对数字信号处理的理解，也为日后在相关领域的研究与实践提供了基础。
通过实验的结果和总结，我们可以得出结论：DSP作为一种数字信号处理的方法，具有广泛的应用前景和重要的实际意义。
DSP实验报告
实验一数字IO应用实验—、实验目的1. 了解DSP开发系统的组成和结构2. 在实验设备上完成I/O硬件连接，编写I/O实验程序并运行验证。
3. 内存观察工具的使用二、实验设备计算机，CCS3.1版本软件，DSP仿真器，教学实验箱三、实验原理2.键值读取程序：该部分有两种方法进行键值的判断。
方法1：利用内存观察工具进行观察方法2：利用LED1-LED8的亮灭对应显示键值。
a)外部中断1的应用参照实验五；b)内存观察键值：程序中定义了三个变量“W”“row”和“col”。
“W”代表是CPLD中键盘的扫描数值，“row”和“col”分别代表键盘的行和列，由行和列可以判定按键的位置。
上述三个变量可以在观察窗口中观察的。
c)利用LED灯显示键值原理，参看实验一。
具体的LED灯显示值以查表的形式读出，请参看“”库文件。
本实验的CPLD地址译码说明：基地址：0x0000，当底板片选CS0为低时，分配有效。
CPU的IO空间：基地址+0x0200 LED灯output 8位外部中断用XINT1：由CPLD分配，中断信号由键盘按键产生。
中断下降沿触发。
KEY_DAT_REG(R)：基地址+0x0004;四、实验步骤和内容1.2407CPU板JUMP1的1和2脚短接，拨码开关S1的第一位置ON,其余置OFF；2.E300板上的开关SW4的第一位置ON，其余OFF；SW3的第四位置ON其余的SW置OFF“DEBUG→Connect”）4.打开系统项目文件 \e300.test\ normal \05_key interface \；“\Debug\”文件“Debug\Go Main”跳到主程序的开始；7.指定位置设置断点；8.View--〉Watch Window打开变量观察窗口；9. 将变量“w”“row”和“col”添加到观察窗口中，改变变量观察窗口的显示方式为HEX显示。
“Debug--〉Animate”全速运行，然后点击E300板上键盘按键，观察窗口中变量变化，同时LED1-LED8灯也相应变化，指示键值。
DSP实验报告
DSP实验报告⼀、综合实验内容和要求1. 实验⽬的(1) 学习掌握CCS3.3编译器的使⽤；(2) 通过实验学习掌握TMS320F28335的GPIO ，浮点计算； (3) 学习并掌握A/D 模块的使⽤⽅法；(4) 学习并掌握中断⽅式和查询⽅式的串⼝通信； (5) 学习并掌握28335DSP 的定时器相关的设置与运⽤； (6) 学习信号时域分析的⽅法，了解相关波形参数的计算⽅法； (7) 了解数字滤波的⼀些基本⽅法； (8) 学习数码管的驱动及运⽤。
(9) 学习MATLAB 串⼝以及画图的运⽤。
2. 实验设计内容与要求：(1) 对给定的周期波形信号采⽤TI 公司的TMS320F28335DSP ，利⽤试验箱上的相关资源计算出波形的周期T ，波形的有效值rms V ，平均值avg V 。
其中，有效值和平均值的计算公式(数字量的离散公式)如下：rms V =1()NavgiV u i N=∑式中N 为⼀个周期采样点数，()u i 为采样序列中的第i 个采样点。
(2) 通过算法计算出波形的有效值和平均值，利⽤串⼝通信把测得的数据发送到串⼝助⼿查看，或者在MATLAB 上编写上位机程序，把发送的数据在MATLAB 上画出来。
(3) 把测得的数据实时显⽰在数码管上。
⼆、硬件电路图1为试验系统的硬件图，硬件电路主要包括TMS320F28335DSP 实验箱，SEED-XDS510仿真器，数码管，SCI,信号发⽣器，电脑，串⼝线等。
图1 硬件电路图三、实验原理本试验主要是通过程序去测量⼀个周期波形的有效值、平均值、峰值等相关参数。
计算离散数据的有效值可⽤公式rms V =平均值可⽤公式1()N avgiV u i N=∑。
所以⾸先需要测出波形的周期，然后确定每个周期需要采样的点数N ，最后去计算平均值和有效值。
v mv 1图2 理想输⼊采样波形如图2所⽰为⼀个正弦输⼊波形，m V 为输⼊波形的峰值，1V 是介于0~ m V 的⼀个值。
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实验报告
南京邮电大学实验报告实验名称：采样、系统性质及滤波系统频率响应和样本处理算法实现加窗和离散傅氏变换数字滤波器设计课程名称：数字信号处理姓名：学号开课时间2011 /2012 学年，第 2 学期实验一一．实验名称：采样、系统性质及滤波 二．实验目的和任务，实验内容：一、观察采样引起的混叠。
设模拟信号为)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 采样点。
两信号波形是否相同？采样后的两序列是否相同？二、判别离散时间系统的时不变性。
（来源：p105 例3.2.2）设输入序列为)(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 的波形。
该系统是否为时不变的？三、利用卷积计算信号通过FIR 滤波器的输出，并观察输出信号的input-on 暂态、input-off暂态和稳态阶段。
（来源：p144 例4.1.8）考虑两个滤波器，⎩⎨⎧≤≤⋅=其它0140)75.0(25.0)(1n n h n ，]1,5,10,105,1[51--=，－2h ；输入)(n x 为周期方波，第一个周期内⎩⎨⎧≤≤≤≤=492502401)(x x n x 。
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实验报告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实验报告
1
0.8
振幅
0.6
0.4
0.2
0 -10
-5
0
5 时间序列 n
10
15
20
Q1.2 clf：清除图形； axis：设置坐标轴范围、可读比例等； title：给图形加标题； xlable：给 x 轴加标注； ylable：给 y 轴加标注。
Q1.3 程序： clf n=-10:20; u=[zeros(1,10) 1 zeros(1,20)]; stem(n+11,u); xlabel('时间序列 n');ylabel('振幅'); title('单位样本序列'); axis([0 32 0 1.2]);
南京邮电大学
实 验 报 告
实验名称：离散时间信号与系统的时、频域表示 离散傅立叶变换和 z 变换 数字滤波器的频域分析和实现 数字滤波器的设计
课程名称：
数字信号处理 A(双语)
班级学号：
姓
名：
开课时间 ：
2015 /2016 学年 第一学期
实验一
1
实验名称：离散时间信号与系统的时、频域表示 实验目的： 熟悉 Matlab 基本命令和信号处理工具箱，加深理解和掌握离散 时间信号与系统的时、频域表示及简单应用。 实验任务： 在 Matlab 环境中，按照要求产生序列，对序列进行基本运算； 对简单离散时间系统进行仿真，计算线性时不变（LTI）系统的冲激 响应和卷积输出；计算和观察序列的离散时间傅立叶变换（DTFT） 幅度谱和相位谱。 实验内容： 基本序列产生和运算： Q1.1～1.3，Q1.23，Q1.30～1.33 离散时间系统仿真： Q2.1～2.3 LTI 系统：Q2.19，Q2.21，Q2.28 DTFT：Q3.1，Q3.2，Q3.4 实验过程描述：
DSP实验报告2
一、实验目的1、掌握TMS320系列DSP的性能、结构原理、指令系统及编程方法；2、熟练掌握CCS集成开发环境的常用开发、调试功能；3、利用MATLAB语言与开发环境进行FIR滤波器的高级语言设计与仿真；4、根据MATLAB仿真结果，进行基于CCS的高级语言和汇编语言仿真；5、在DSP实验板上实时实现该滤波器。
二、试验设备1、TMS320C5402DSP实验开发板1块，JTAG仿真探头1个，5V 交流适配器1个；2、音频信号源1个，耳机1对；3、微型计算机1台。
三、试验内容1、试验硬件简介：本试验的板子采用南京航空航天大学DSP实验室自己开发的试验板子，它有以下几个特点：1）试验板子的核心处理器是TMS320VC5402；2）外部单5V电源供电；3）提供了与计算机串行口的数据通信功能，可以与计算机交换数据；4）提供了语音输入和语音输出的街口，可以进行语音信号的分析与处理；5）提供了可与外部相连接的EHPI8、并口、McBSP口等VC5402片上外设接口。
通过这些接口，可以方便的对该试验板子进行功能扩展；6）该试验板子提供大容量的外部程序和数据RAM以及FLASH存储空间，方便学生开发和使用；7）试验板子集成了一块ALTERA公司CPLD（EPM7128S），可以通过该CPLD实现对试验板子上各种逻辑的控制管理；8）试验板子提供了必要的复位，测试和开关插销，为用户调试和测试提供了方便。
本试验用到了试验板子辅助器件为语音处理部分的硬件SRAM、闪烁灯器件等等。
2、在Matlab实现算法仿真试验目的是设计一个64阶FIR数字低通滤波器，其参数如下：截至频率为3000Hz，采样频率为2000Hz。
利用Matlab工具箱设计FIR数字滤波器，其设计语言如下：>>c=fir1(63,0.3);>>freqz(c,1,4096);>>plot(c);得到数字滤波器系数如下：-0.00079887, -0.00039007, 0.0004368, 0.0011085, 0.00095091,-0.00025572, -0.0017797, -0.002173, -0.00046388, 0.0025341,0.0042491, 0.0023285, -0.0027578, -0.0070632, -0.005921, 0.0015267,0.010096, 0.011687, 0.0023706, -0.012368, -0.019945, -0.010603,0.012305, 0.031247, 0.026339, -0.0069494, -0.048111, -0.060186,-0.01384, 0.088797, 0.20863, 0.289, 0.289, 0.20863,0.088797, -0.01384, -0.060186, -0.048111, -0.0069494, 0.026339,0.031247, 0.012305, -0.010603, -0.019945, -0.012368, 0.0023706,0.011687, 0.010096, 0.0015267, -0.005921, -0.0070632, -0.0027578,0.0023285, 0.0042491, 0.0025341, -0.00046388, -0.002173, -0.0017797,-0.00025572, 0.00095091, 0.0011085, 0.0004368, -0.00039007, -0.00079887得到数字滤波器的时域与频域图如下：图一数字滤波器时域图图二数字滤波器的幅频相频图设计一正弦函数进行仿真：sin(2*Pi*10000*t)+sin(2*Pi*3000*t)其中10000与3000为频率，仿真结果看能否过滤高频10000。
南邮通信原理实验报告
南邮通信原理实验报告南邮通信原理实验报告一、引言通信原理是现代通信技术的基础，通过实验可以更好地理解和掌握通信原理的相关知识。
本实验旨在通过实践操作，加深对通信原理的理解，并掌握信号的调制和解调技术。
二、实验目的1. 了解通信原理的基本概念和理论知识；2. 学习信号的调制和解调技术；3. 掌握实验仪器的使用方法；4. 进行实验操作，验证通信原理的相关理论。
三、实验原理1. 调制技术调制是将要传输的信息信号与载波信号相结合，形成调制信号。
常见的调制方式有幅度调制（AM）、频率调制（FM）和相位调制（PM）。
2. 解调技术解调是将调制信号还原为原始的信息信号。
常见的解调方式有包络检波、相干解调和频率解调。
四、实验步骤1. 实验仪器准备根据实验要求，准备好信号发生器、示波器、调制解调器等实验仪器，并将其连接好。
2. 幅度调制实验使用信号发生器产生载波信号和待调制的信号，将两者输入调制解调器中，调制解调器进行幅度调制。
通过示波器观察调制后的信号，并进行分析。
3. 频率调制实验使用信号发生器产生载波信号和待调制的信号，将两者输入调制解调器中，调制解调器进行频率调制。
通过示波器观察调制后的信号，并进行分析。
4. 相位调制实验使用信号发生器产生载波信号和待调制的信号，将两者输入调制解调器中，调制解调器进行相位调制。
通过示波器观察调制后的信号，并进行分析。
5. 包络检波实验使用信号发生器产生经过幅度调制的信号，将其输入调制解调器中，调制解调器进行包络检波。
通过示波器观察解调后的信号，并进行分析。
6. 相干解调实验使用信号发生器产生经过相位调制的信号，将其输入调制解调器中，调制解调器进行相干解调。
通过示波器观察解调后的信号，并进行分析。
7. 频率解调实验使用信号发生器产生经过频率调制的信号，将其输入调制解调器中，调制解调器进行频率解调。
通过示波器观察解调后的信号，并进行分析。
五、实验结果与分析通过实验操作，我们观察到不同调制方式下的信号特点，并进行了相应的解调操作。
DSP实验报告
一．实验目的：通过数据采集加深对采样定理的理解，并对输入信号进行滤波处理。
二．实验仪器：DSP仿真器仿真软件（CCS）PC机实验箱示波器信号发生器三．实验步骤：1. 结合实验原理和前面的实验所学的知识，编写实验程序。
首先设置AD的采样频率设置为6KHz，信号的输出幅度为1V，输入频率为50Hz。
2. 将信号发生器，仿真器连接到实验箱。
3. 连接示波器到AD输入端和DA的输出端，设置输入信号为50Hz的正弦波形。
4. 运行编写好的程序，在示波器上观察和记录在采样频率为6KHz时的输入输出波形。
通过改变输出信号箱的输入频率，输入频率从50Hz逐渐增大，观察和记录下相应的输入输出波形，并分析输出变化的主要原因.四．实验原理和内容：模拟信号经过(A/D) 变换转换为数字信号的过程称之为采样，信号采样后其频谱产生了周期延拓，每隔一个采样频率fs，重复出现一次。
为保证采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成份的两倍，这称之为采样定理。
根据香农采样定理：如果A/D转换器的输入信号具有有限带宽2×F，则只需要AD转换器的采样频率f大于2F，原始信号就可以完全从采样信号中恢复出来。
而如果采样频率低于信号频率的2倍，则不能不失真的恢复原始信号。
实验中，我们选择对一个确定的信号进行采样，然后将采样后的数据从DA输出，从DA的输出使用示波器查看输出后的波形。
如果满足采样定理，可以从示波器看到和原始信号一样的波形；反之，如果不满足采样定理，就不能从示波器看到和原始信号一样的波形。
调整AD转换器的采样频率，对两种情况分别验证采样定理。
因此，对于一个单一频率的正弦信号，当采样频率时就可保证可以由采样后的信号无失真地重现原信号；如果采样频率低于就会发生混叠失真。
在实际的情况中，首先要让信号通过一个截止频率小于的低通滤波器，使其成为一个限带信号。
然后对其采样就可以保证信号无混叠失真。
TLV1571的最大采样频率为625K，我们以实验二相同的方法对DSP数据总线进行设置，通过改变DSP读取总线数据的时间来改变采样频率，读取时间可以通过DSP的定时器来设置。
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。
南邮微机实验报告
一、实验名称微机原理与接口技术实验二、实验目的1. 理解微机的基本工作原理和硬件组成。
2. 掌握计算机常用接口芯片的基本功能和使用方法。
3. 熟悉微机实验系统的操作流程。
4. 培养动手实践能力和分析问题、解决问题的能力。
三、实验内容1. 微机硬件系统认知- 了解微机的硬件组成，包括CPU、内存、主板、硬盘、显卡、声卡等。
- 通过实验观察微机硬件各部分的功能和连接方式。
2. 微机接口芯片实验- 学习并行接口芯片8255的使用方法和编程技巧。
- 实现键盘输入和LED显示，设计一个简单的定时显示装置。
3. 微机实验系统操作- 熟悉微机实验系统的使用方法和操作流程。
- 掌握实验软件的安装和使用。
4. 汇编语言程序设计- 学习汇编语言的基本语法和编程技巧。
- 编写简单的汇编语言程序，实现特定功能。
四、实验步骤1. 微机硬件系统认知- 观察微机硬件各部分的功能和连接方式。
- 使用实验设备观察CPU、内存、主板、硬盘等硬件的工作状态。
2. 微机接口芯片实验- 根据实验指导书，配置并行接口芯片8255。
- 编写程序实现键盘输入和LED显示功能。
- 设计定时显示装置，实现时间显示和按键控制。
3. 微机实验系统操作- 安装实验软件，熟悉实验系统的操作流程。
- 使用实验软件进行实验操作。
4. 汇编语言程序设计- 学习汇编语言的基本语法和编程技巧。
- 编写汇编语言程序，实现特定功能。
五、实验结果与分析1. 微机硬件系统认知- 通过实验观察，了解了微机硬件各部分的功能和连接方式。
- 掌握了微机实验系统的使用方法和操作流程。
2. 微机接口芯片实验- 成功实现了键盘输入和LED显示功能。
- 设计的定时显示装置能够正常工作。
3. 微机实验系统操作- 熟练掌握了微机实验系统的操作流程。
4. 汇编语言程序设计- 成功编写了汇编语言程序，实现了特定功能。
六、实验心得体会通过本次实验，我对微机原理和接口技术有了更深入的了解。
实验过程中，我学会了如何使用微机实验系统，掌握了微机接口芯片的编程技巧，提高了动手实践能力。
DSP应用技术实验报告
-源程序文件：C语言或汇编语言文件（*.c或*.asm）
-头文件（*.h）
-命令文件（*.cmd）
-库文件（*.lib,*.obj）
这几个文件都是工程必要的文件。
2、GPIO的使用
GPIO作为通用输入/输出接口，有着广泛的应用。简单配置后，即可用来控制和接口各种外设。
件系统的通信、控制和读取硬件系统的状态和数据。
(3).评估模块（如实验箱中的TS-DM64x嵌入式多通道音视频处理平台）：提供软件
运行和调试的平台和用户系统设计开发的参照。
CCS主要完成系统的软件开发和调试，它提供一整套的程序编制、维护、编译的调试环境，能将汇编语言和C语言程序编译连接生产COFF（公共目标文件）格式的可执行文件，并能将程序下载到目标DSP系统上运行调试。
-头文件（*.h）
-命令文件（*.cmd）
-库文件（*.lib,*.obj）
这几个文件都是工程必要的文件。
2、软件的调试
(1)执行View Watch Window打开观察窗口。
(2)在main.c中，用鼠标双击一个变量（比如num），再单击鼠标右键，选择“QuickWatch”，CCS将打开Quick Watch窗口并显示选中的变量。
桂林航天工业学院
实验报告
课程名称DSP技术应用
开课学期2014-2015秋
实验室南训510左
班级
姓名
实验名称DSP技术应用试验
桂林航天工业学院学生实验报告
课程名称
DSP技术应用
实验项目名称
开课系（部）及实验室
南训510左
实验日期
2014年9月25日
学生姓名
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。
DSP技术与应用实验报告
DSP实验报告姓名：年级：专业：学号：任课教师：目录实验一、CCS入门实验实验二、编写一个以C语言为基础的DSP程序实验三、异步串口通信实验四、A/D及D/A转换实验五、PID算法控制实验实验六、快速傅里叶变换（FFT）算法实验一、CCS入门实验Project/new；输入：volume,按Finish；Project/add files to project；查找范围：C:\CCStudio_v3.3\tutorial\dsk2812\volume1；文件类型：all files；分别添加：Volume.c、Load.asm、Vectors.asm、Volume.cmd；在C:\CCStudio_v3.3\C2000\cgtools\lib中添加：rts2800_ml.lib；将C:\CCStudio_v3.3\tutorial\dsk2812\volume1\volume.h复制到：C:\CCStudio_v3.3\C2000\cgtools\include文件夹中；Project/rebuild all；File/load program；打开C:\CCStudio_v3.3\MyProjects\volume\Debug\volume.out；鼠标双击：volume.c,找到dataIO()，右击，运行toggle software breakpoint；File/data，查找范围：C:\CCStudio_v3.3\tutorial\dsk2812\volume1；文件类型：all files，打开：sina.dat，输入：0x8000、0x2000；View/memory，0x8000；View/graph/time\frequency；如下图波形：频谱实验二、编写一个以C语言为基础的DSP程序一．实验目的1．学习用标准C语言编制程序；了解常用的C语言程序设计方法和组成部分。
2．学习编制连接命令文件，并用来控制代码的连接。
DSP数字信号处理实验报告(精品)
FFT 频谱分析一、 实验目的a) 进一步加深DFT 算法原理和基本性质的理解b) 熟悉FFT 算法原理和FFT 程序的应用c) 学习用FFT 对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确的应用FFT二、 实验原理a) 离散傅里叶变换（DFT ）：离散傅里叶变换在作为有限长序列的傅里叶变换表示法在理论上相当重要；由于存在着计算离散傅里叶变换的快速算法（FFT ），从而离散傅里叶变换在各种数字信号处理的算法中起到了核心的作用。
其对应的离散傅里叶变换对为：X (K )=DFT [x (n )]=∑x (n )W N nk 0≤k ≤N −1N−1n=0x (n )=IDFT [X (K )]=1N ∑x (n )W N nk 0≤n ≤N −1N−1n=0需要注意：有限长序列的离散傅里叶变换及周期序列的离散傅里叶级数之间的关系是：它们仅仅是n 、k 的取值不同，DFT 只取主值区间。
X(n)、X(k)是一个有限长序列的离散傅里叶变换对，已知其中一个序列，就可以唯一确定另一个序列，这是因为x(n)、X(k)都是长为N 的序列，都有N 个独立值，所以信息量相同。
b) DFT 计算量：长度为N 的DFT 的计算量是N 个复数乘法和N-1个复数加法（4N 个实数乘法和4N-2个复数加法）c) FFT 的计算量：长度为N 的FFT 的计算量是 N 2log 2N 个复数乘法和N log 2N 个复数加法。
三、 实验步骤a) 复习DFT 的定义、性质和用DFT 做谱分析的有关内容b) 复习FFT 算法原理与编程思想，熟悉DIT-FFT 运算流图c) 编制信号产生程序，产生典型信号尽心谱分析。
d) 进行以下几个信号的谱分析i.x 1(n)=R 4(n ) ii. x 2(n)={n +1 ,0≤n ≤38−n ,4≤n ≤70 ,其他niii. x 3(n)= {4−n ,0≤n ≤3n −3 ,4≤n ≤70 ,其他niv. x 4(n )=cos π4n ,0≤n ≤19v.x 5(n )=sin π8n ,0≤n ≤19 vi.x 6(n )= cos 8πt +cos 16πt +cos 20πt vii.令x 7(n )=x 4(n )+x 5(n) N=8,16 viii. 令x 8(n )=x 4(n )+jx 5(n) N=8,16针对上述信号进行逐一的谱分析，下面给出针对各个信号的FFT 点数N 及对连续信号x 6(n )的采样频率f s ，供实验时参考 ：x 1(n ),x 2(n ),x 3(n ),x 4(n ),x 5(n ) N=8,16x 6(n ) f s =64Hz ，N =16,32,64四、 实验内容a)对x1(n)=R4(n)进行谱分析1.编辑代码x1=[1 1 1 1];y11 = fft(x1,8);y12 = fft(x1,16);subplot(2,2,1);stem(0:3,x1);title('函数X1的图像');subplot(2,2,2);stem(0:7,abs(y11));title('N=8的DFT');subplot(2,2,4);stem(0:15,abs(y12));title('N=16的DFT');2.谱分析图片b) 对x 2(n)={n +1 ,0≤n ≤38−n ,4≤n ≤70 ,其他n进行谱分析i. 编辑代码x2 = [1 2 3 4 4 3 2 1];y11 = fft(x2,8);y12 = fft(x2,16);subplot(2,2,1);stem(0:7,x2);title('函数X2的图像');subplot(2,2,2);stem(0:7,abs(y11));title('N=8的DFT');subplot(2,2,4);stem(0:15,abs(y12));c)对x3(n)={4−n ,0≤n≤3n−3 ,4≤n≤70 ,其他n进行谱分析i.谱分析程序x3 = [4 3 2 1 1 2 3 4];y11 = fft(x2,8);y12 = fft(x2,16);subplot(2,2,1);stem(0:7,x2);title('函数X2的图像');subplot(2,2,2);stem(0:7,abs(y11));title('N=8的DFT');subplot(2,2,4);stem(0:15,abs(y12));n ,0≤n≤19进行谱分析d)对x4(n)=cosπ4i.谱分析程序n = 0:1:19;x2 = cos(0.25*pi*n);y11 = fft(x2,32);y12 = fft(x2,64);subplot(2,2,1);stem(0:19,x2);title('函数X4的图像');subplot(2,2,2);stem(0:31,abs(y11));title('N=32的DFT');subplot(2,2,4);stem(0:63,abs(y12));title('N=64的DFT');ii.谱分析图片n ,0≤n≤19进行谱分析e)对x5(n)=sinπ8i.谱分析程序n = 0:1:19;x2 = sin(0.125*pi*n);y11 = fft(x2,32);y12 = fft(x2,64);subplot(2,2,1);stem(0:19,x2);title('函数X5的图像');subplot(2,2,2);stem(0:31,abs(y11));title('N=32的DFT');subplot(2,2,4);stem(0:63,abs(y12));title('N=64的DFT');ii.谱分析图片f)对x6(n)=cos8πt+cos16πt+cos20πt进行谱分析i.谱分析程序n = 0:1:15;x1 = cos(8*pi*n/64)+cos(16*pi*n/64)+cos(20*pi*n/64);n = 0:1:31;x2 = cos(8*pi*n/64)+cos(16*pi*n/64)+cos(20*pi*n/64); n = 0:1:63;x3 = cos(8*pi*n/64)+cos(16*pi*n/64)+cos(20*pi*n/64);y1 = fft(x1,16);y2 = fft(x2,32);y3 = fft(x3,64);subplot(3,2,1);stem(0:15,x1);title('函数X6 N=16 的图像');subplot(3,2,2);stem(0:15,abs(y1));title('N=16的DFT');subplot(3,2,3);stem(0:31,x2);title('函数X6 N=32 的图像');subplot(3,2,4);stem(0:31,abs(y2));title('N=32的DFT');subplot(3,2,5);stem(0:63,x3);title('函数X6 N=64 的图像');subplot(3,2,6);stem(0:63,abs(y3));title('N=64的DFT');ii.谱分析图片g)对x7(n)=x4(n)+x5(n)进行谱分析i.谱分析程序n = 0:1:19;x2 = j*sin(0.125*pi*n) + cos(0.25*pi*n);y11 = fft(x2,8);y12 = fft(x2,16);subplot(2,2,1);stem(0:19,x2);title('函数X=X4 + X5 的图像');subplot(2,2,2);stem(0:7,abs(y11));title('N=8的DFT');subplot(2,2,4);stem(0:15,abs(y12));title('N=16的DFT');ii.谱分析图片h)对x8(n)=x4(n)+jx5(n)进行谱分析i.谱分析程序n = 0:1:19;x2 = j*sin(0.125*pi*n) + cos(0.25*pi*n);y11 = fft(x2,8);y12 = fft(x2,16);subplot(2,2,1);stem(0:19,x2);title('函数X=X4 + X5 的图像');subplot(2,2,2);stem(0:7,abs(y11));title('N=8的DFT');subplot(2,2,4);stem(0:15,abs(y12));title('N=16的DFT');ii.谱分析图片五、实验分析六、实验结论通过这次利用FFT对信号进行频谱分析的实验，更加深刻的理解了DFT算法的理解和性质的理解；同时也更为熟悉了FFT算法的原理和应用；学会使用FFT 对离散信号和连续信号进行频谱分析，了解了可能出现的分析误差和原因。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
南京邮电大学实验报告实验名称：1、离散时间信号与系统的时、频域表示2、离散傅立叶变换和z变换3、数字滤波器的频域分析和实现4、数字滤波器的设计课程名称数字信号处理A(双语)班级学号姓名指导老师开课时间 2016/2017学年第一学期实验一：离散时间信号与系统的时、频域表示一、实验目的和任务：●主要内容：序列的生成、序列DTFT谱计算、简单滤波器的仿真，对应实验指导书Page1－34。
●实验要求：●（1）按照要求产生基本序列和复杂序列，对序列进行运算，求序列的DTFT幅度谱和相位谱。
●（2）计算简单的离散时间系统的输出，包括简单滤波器的仿真、计算单个系统和互连系统的输出。
二、实验内容：●具体包括●基本序列产生和运算：Q1.1－Q1.3，Q1.23，Q1.30－Q1.33●离散时间系统仿真：Q2.1－Q2.3●LTI系统：Q2.19、Q2.21，Q2.28●DTFT：Q3.1，Q3.2，Q3.4三、实验过程与结果分析：Q1.1 运行P1_1产生单位样本序列u[n] 的程序与显示的波形如下:clf;n= -10:20;u = [zeros(1,10) 1 zeros(1,20)];stem(n,u);xlabel('时间序列n');ylabel('振幅');title('单位样本序列');axis([-10 20 0 1.2]);结果：时间序列 n振幅Q1.2clf 清除波形 axis 设置坐标轴范围，可读比例等 title 给图形加标题 xlabel 给x 加标注 Ylabel 给y 加标注Q1.3clf;n= -10:20;u = [zeros(1,10) 1 zeros(1,20)]; stem(n+11,u);xlabel('时间序列 n');ylabel('振幅'); title('单位样本序列'); axis([0 32 0 1.2]);结果：0.20.40.60.81时间序列 n振幅Q1.23n = 0:50; f= 0.08; phase = pi/2; A = 2.5;arg =2*pi*f*n - phase; x = A*cos(arg); clf; stem(n,x); axis([0 50 -3 3]); grid;title('正弦序列'); xlabel('时间序号n'); ylabel('振幅'); axis;结果：05101520253035404550正弦序列时间序号n振幅Q1.30s[n]是线性增加伴随着实指数缓慢衰减的图像加性噪声d[n]是均匀分布在-0.4和+0.4之间的随机序列Q1.31不能，因为d 是列向量，s 是行向量Q1.32x1是x 的延时，x2和x 相等，x3超前于xQ1.33产生图例说明Q2.1 clf;n = 0:100;s1 = cos(2*pi*0.05*n); s2 = cos(2*pi*0.47*n); x = s1+s2;M = input('滤波器所需的长度='); num = ones(1,M); y = filter(num,1,x)/M; subplot(2,2,1); plot(n,s1); axis([0,100,-2,2]); xlabel('时间序列n');ylabel('振幅'); title('信号#1'); subplot(2,2,2); plot(n,s2); axis([0,100,-2,2]); xlabel('时间序列n'); ylabel('振幅'); title('信号#2'); subplot(2,2,3); plot(n,x); axis([0,100,-2,2]); xlabel('时间序列n'); ylabel('振幅'); title('输入信号'); subplot(2,2,4); plot(n,y); axis([0,100,-2,2]); xlabel('时间序列n'); ylabel('振幅'); title('输出信号'); axis;结果：50100-2-1012时间序列n 振幅信号#150100-2-1012时间序列n 振幅信号#250100-2-1012时间序列n振幅输入信号50100-2-1012时间序列n振幅输出信号S[n]被离散时间系统抑制Q2.2n = 0:100;s1 = cos(2*pi*0.05*n);s2 = cos(2*pi*0.47*n);x = s1+s2;M = input('滤波器所需长度 = ');num = (-1).^[0:M-1];y = filter(num,1,x)/M;clf;subplot(2,2,1);plot(n, s1);axis([0, 100, -2, 2]);xlabel('时间序号n'); ylabel('振幅'); title('信号 #1');subplot(2,2,2);plot(n, s2);axis([0, 100, -2, 2]);xlabel(' 时间序号n'); ylabel('振幅'); title('信号 #2');subplot(2,2,3);plot(n, x);axis([0, 100, -2, 2]);xlabel('时间序号 n'); ylabel('振幅'); title('输入信号');subplot(2,2,4);plot(n, y);axis([0, 100, -2, 2]);xlabel('时间序号 n'); ylabel('振幅'); title('输出信号');axis;显示的波形如下:50100时间序号n 振幅信号 #150100时间序号n 振幅信号 #250100时间序号 n振幅输入信号50100时间序号 n振幅输出信号改变LTI 系统对输入的影响是,系统现在是一个高通滤波器。
它通过高频输入组件s2来替代低频输入组件s1.Q2.3 clf;n = 0:100;s1 = cos(2*pi*0.04*n); s2 = cos(2*pi*0.6*n); x = s1+s2;M = input('滤波器所需的长度='); num = ones(1,M); y = filter(num,1,x)/M; subplot(2,2,1); plot(n,s1); axis([0,100,-2,2]); xlabel('时间序列n'); ylabel('振幅'); title('信号#1'); subplot(2,2,2); plot(n,s2); axis([0,100,-2,2]); xlabel('时间序列n'); ylabel('振幅'); title('信号#2'); subplot(2,2,3); plot(n,x); axis([0,100,-2,2]); xlabel('时间序列n');ylabel('振幅'); title('输入信号'); subplot(2,2,4); plot(n,y); axis([0,100,-2,2]); xlabel('时间序列n'); ylabel('振幅'); title('输出信号'); axis;结果：50100-2-1012时间序列n 振幅信号#150100-2-1012时间序列n 振幅信号#250100-2-1012时间序列n振幅输入信号50100-2-1012时间序列n振幅输出信号Q2.19clf; N=40;num=[2.2403 2.4908 2.2403]; den=[1 -0.4 0.75]; y=impz(num,den,N); stem(y);xlabel('时间序号n'); ylabel('振幅'); title('冲激响应'); grid;结果：510152025303540-3-2-101234时间序号n振幅冲激响应Q2.21 clf; N=40;num=[0.9 -0.45 0.35 0.002]; den=[1.0 0.71 -0.46 -0.62]; x=[1 zeros(1,N-1)]; y=filter(num,den,x); stem(y);xlabel('时间序号n'); ylabel('振幅'); title('冲激响应'); grid;510152025303540-1.5-1-0.500.511.52时间序号n振幅冲激响应Q2.28 clf;h=[3 2 1 -2 1 0 -4 0 3]; x=[1 -2 3 -4 3 2 1]; y=conv(h,x); n=0:14; subplot(2,1,1); stem(n,y);xlabel('时间序号n'); ylabel('振幅');title('用卷积得到的输出'); grid;x1=[x zeros(1,8)]; y1=filter(h,1,x1); subplot(2,1,2); stem(n,y1); xlabel('时间序号n') ylabel('振幅');title('由滤波生成的输出'); grid;时间序号n 振幅用卷积得到的输出时间序号n振幅由滤波生成的输出Q3.1原始序列是：Pause:暂停命令Q3.2 clf;w= -4*pi:8*pi/511:4*pi; num=[2 1]; den=[1 -0.6]; h=freqz(num,den,w); subplot(2,1,1) plot(w/pi,real(h)); gridtitle('H(e^{j\omega}的实部') xlabel('\omega/\pi'); ylabel('振幅'); subplot(2,1,2) plot(w/pi,imag(h)); gridtitle('H(e^{j\omega}的虚部') xlabel('\omega/\pi'); ylabel('振幅'); pause subplot(2,1,1) plot(w/pi,abs(h)); gridtitle('|H(e^{j\omega}|幅度谱') xlabel('\omega/\pi'); ylabel('振幅'); subplot(2,1,2) plot(w/pi,angle(h)); gridtitle('相位谱arg[H(e^{j\omega}]') xlabel('\omega/\pi'); ylabel('以弧度为单位的相位'); 结果：-4-3-2-10123402468H(e j ω的实部ω/π振幅-4-3-2-101234-4-2024H(e j ω的虚部ω/π振幅-4-3-2-10123402468|H(e j ω|幅度谱ω/π振幅-4-3-2-101234-2-1012相位谱arg[H(e j ω]ω/π以弧度为单位的相位是w 的周期。