数字信号调制实验报告

数字信号调制实验报告（一）

一、2ASK调制

1、搭建ask调制模型，要求，基带数字信号速率为50bit/s，载波信号频率为300Hz。

（1）画出调制模型，并画出调制后的波形图

（2）画出解调模型。（包络解调和相干解调任选一种）

二、2FSK调制

2、搭建2FSK调制模型，要求：基带数字信号速率为50bit/s，载波信号频率为200Hz，500Hz。

（1）画出调制模型，并画出调制后的波形图

BPSK调制及解调实验报告

实验五BPSK调制及解调实验一、实验目的 1、掌握BPSK调制和解调的基本原理； 2、掌握BPSK数据传输过程，熟悉典型电路； 3、了解数字基带波形时域形成的原理和方法，掌握滚降系数的概念； 4、熟悉BPSK调制载波包络的变化； 5、掌握BPSK载波恢复特点与位定时恢复的基本方法；二、实验器材 1、主控&信号源、9号、13号模块各一块 2、双踪示波器一台 3、连接线若干三、实验原理 1、BPSK调制解调（9号模块）实验原理框 PSK调制及解调实验原理框图 2、BPSK调制解调（9号模块）实验框图说明基带信号的1电平和0电平信号分别与256KHz载波及256KHz反相载波相乘，叠加后得到BPSK调制输出；已调信号送入到13模块载波提取单元得到同步载波；已调信号与相干载波相乘后，经过低通滤波和门限判决后，解调输出原始基带信号。四、实验步骤实验项目一 BPSK调制信号观测（9号模块）概述：BPSK调制实验中，信号是用相位相差180°的载波变换来表征被传递的信息。本项目通过对比观测基带信号波形与调制输出波形来验证BPSK调制原理。 1、关电，按表格所示进行连线。

2、开电，设置主控菜单，选择【主菜单】→【通信原理】→【BPSK/DBPSK数字调制解调】。将9号模块的S1拨为0000，调节信号源模块W3使256 KHz载波信号峰峰值为3V。 3、此时系统初始状态为：PN序列输出频率32KHz。 4、实验操作及波形观测。（1）以9号模块“NRZ-I”为触发，观测“I”；（2）以9号模块“NRZ-Q”为触发，观测“Q”。（3）以9号模块“基带信号”为触发，观测“调制输出”。思考：分析以上观测的波形，分析与ASK有何关系？实验项目二 BPSK解调观测（9号模块）概述：本项目通过对比观测基带信号波形与解调输出波形，观察是否有延时现象，并且验证BPSK解调原理。观测解调中间观测点TP8，深入理解BPSK解调原理。 1、保持实验项目一中的连线。将9号模块的S1拨为“0000”。 2、以9号模块测13号模块的“SIN”，调节13号模块的W1使“SIN”的波形稳定，即恢复出载波。 3、以9号模块的“基带信号”为触发观测“BPSK解调输出”，多次单击13号模块的“复位”按键。观测“BPSK解调输出”的变化。 4、以信号源的CLK为触发，测9号模块LPF-BPSK，观测眼图。思考：“BPSK解调输出”是否存在相位模糊的情况？为什么会有相位模糊的情况？五、实验报告 1、分析实验电路的工作原理，简述其工作过程；输入的基带信号由转换开关转接后分成两路，一路经过差分编码控制256KHz的载频，另一路经倒相去控制256KHz的载频。???解调采用锁相解调，只要在设计锁相环时，使它锁定在FSK的一个载频上此时对应的环路滤波器输出电压为零，而对另一载频失锁，则对应的环路滤波器输出电压不为零，那末在锁相环路滤波器输出端就可以获得原基带信号的信息。? 2、分析BPSK调制解调原理。调制原理是：基带信号先经过差分编码得到相对码，再根据相对码进行绝对调相，即将相对码的1电平和0电平信号分别与256K载波及256K反相载波相乘，叠加后得到DBPSK 调制输出。?

《数字信号处理》实验报告

数字信号处理》实验报告年级：2011 级班级：信通 4 班姓名：朱明贵学号： 111100443 老师：李娟福州大学 2013 年11 月

实验一快速傅里叶变换(FFT)及其应用一、实验目的 1. 在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉MATLAB^的有关函数。 2. 熟悉应用FFT对典型信号进行频谱分析的方法。 3. 了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用FFT。 4. 熟悉应用FFT实现两个序列的线性卷积和相关的方法。二、实验类型演示型三、实验仪器装有MATLA爵言的计算机四、实验原理在各种信号序列中，有限长序列信号处理占有很重要地位，对有限长序列，我们可以使用离散Fouier变换(DFT)。这一变换不但可以很好的反映序列的频谱特性，而且易于用快速算法在计算机上实现，当序列x(n)的长度为N时，它的DFT定义为： JV-1 $生反变换为：如-器冃吋科— 有限长序列的DFT是其Z变换在单位圆上的等距采样，或者说是序列Fourier变换的等距采样，因此可以用于序列的谱分析。 FFT并不是与DFT不同的另一种变换，而是为了减少DFT运算次数的一种快速算法。它是对变换式进行一次次分解，使其成为若干小点数的组合，从而减少运算量。常用的FFT 是以2为基数的，其长度A - o它的效率高，程序简单，使用非常方便，当要变换的序列长度不等于2的整数次方时，为了使用以2为基数的FFT，可以用末位补零的方法，使其长度延长至2的整数次方。 (一)在运用DFT进行频谱分析的过程中可能的产生三种误差 1 .混叠序列的频谱是被采样信号频谱的周期延拓，当采样速率不满足Nyquist定理时，就会发生频谱混叠，使得采样后的信号序列频谱不能真实的反映原信号的频谱。避免混叠现象的唯一方法是保证采样速率足够高，使频谱混叠现象不致出现，即在确定采样频率之前，必须

数字信号处理实验报告

语音信号的数字滤波一、实验目的： 1、掌握使用FFT进行信号谱分析的方法 2、设计数字滤波器对指定的语音信号进行滤波处理二、实验内容设计数字滤波器滤除语音信号中的干扰（4 学时） 1、使用Matlab的fft函数对语音信号进行频谱分析，找出干扰信号的频谱； 2、设计数字滤波器滤除语音信号中的干扰分量，并进行播放对比。三、实验原理通过观察原语音信号的频谱，幅值特别大的地方即为噪声频谱分量，根据对称性，发现有四个频率的正弦波干扰，将它们分别滤掉即可。采用梳状滤波器，经过计算可知，梳状滤波器h[n]={1，A，1}的频响|H(w)|=|A+2cos(w)|，由需要滤掉的频率分量的频响w,即可得到A，进而得到滤波器的系统函数h[n]。而由于是在离散频域内进行滤波，所以令w=(2k*pi/N)即可。对原信号和四次滤波后的信号分别进行FFT变换，可以得到它们的幅度相应。最后，将四次滤波后的声音信号输出。四、matlab代码 clc;clear;close all; [audio_data,fs]=wavread('SunshineSquare.wav'); %读取未处理声音 sound(audio_data,fs); N = length(audio_data); K = 0:2/N:2*(N-1)/N; %K为频率采样点

%sound(audio_data,fs); %进行一次FFT变换 FFT_audio_data=fft(audio_data); mag_FFT_audio_data = abs(FFT_audio_data); %画图 figure(1) %原信号时域 subplot(2,1,1);plot(audio_data);grid; title('未滤波时原信号时域');xlabel('以1/fs为单位的时间');ylabel('采样值'); %FFT幅度相位 subplot(2,1,2);plot(K,mag_FFT_audio_data);grid; title('原信号幅度');xlabel('以pi为单位的频率');ylabel('幅度'); %构造h[n]={1，A，1}的梳状滤波器，计算A=2cosW，妻子W为要滤掉的频率%由原信号频谱可知要分四次滤波，滤掉频响中幅度大的频率分量 %第一次滤波 a = [1,0,0,0];%y[n]的系数 [temp,k]=max(FFT_audio_data); A1=-2*cos(2*pi*k/N); h1=[1,A1,1]; audio_data_h1 = filter(h1,a,audio_data); FFT_audio_data_h1=fft(audio_data_h1);

抽样定理和PCM调制解调实验报告

《通信原理》实验报告实验一：抽样定理和PAM调制解调实验系别：信息科学与工程学院专业班级：通信工程1003班学生姓名：陈威同组学生：杨鑫成绩：指导教师：惠龙飞（实验时间：2012 年 12 月 7 日——2012 年 12 月28日）华中科技大学武昌分校

1、实验目的 1对电路的组成、波形和所测数据的分析，加深理解这种调制方法的优缺点。 2.通过脉冲幅度调制实验，使学生能加深理解脉冲幅度调制的原理。 2、实验器材 1、信号源模块一块 2、①号模块一块 3、60M 双踪示波器一台 4、连接线若干 3、实验原理 3.1基本原理 1、抽样定理图3-1 抽样与恢复 2、脉冲振幅调制（PAM ）所谓脉冲振幅调制，即是脉冲载波的幅度随输入信号变化的一种调制方式。如果脉冲载波是由冲激脉冲组成的，则前面所说的抽样定理，就是脉冲增幅调制的原理。自然抽样平顶抽样 ) (t m ) (t T

图3-3 自然抽样及平顶抽样波形 PAM方式有两种：自然抽样和平顶抽样。自然抽样又称为“曲顶”抽样，(t)的脉冲“顶部”是随m(t)变化的，即在顶部保持了m(t)变已抽样信号m s 化的规律（如图3-3所示）。平顶抽样所得的已抽样信号如图3-3所示，这里每一抽样脉冲的幅度正比于瞬时抽样值，但其形状都相同。在实际中，平顶抽样的PAM信号常常采用保持电路来实现，得到的脉冲为矩形脉冲。四、实验步骤 1、将信号源模块、模块一固定到主机箱上面。双踪示波器，设置CH1通道为同步源。 2、观测PAM自然抽样波形。（1）将信号源上S4设为“1010”，使“CLK1”输出32K时钟。（2）将模块一上K1选到“自然”。（3）关闭电源，连接表3-1 抽样实验接线表 (5)用示波器观测信号源“2K同步正弦波”输出，调节W1改变输出信号幅度，使输出信号峰-峰值在1V左右。在PAMCLK处观察被抽样信号。CH1接PAMCLK（同步源），CH2接“自然抽样输出”（自然抽样PAM信号）。

西电数字信号处理上机实验报告

数字信号处理上机实验报告 14020710021 张吉凯第一次上机实验一：设给定模拟信号()1000t a x t e -=，t 的单位是ms 。 (1) 利用MATLAB 绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分量降低到峰值的3%以下的频谱）。 (2) 用两个不同的采样频率对给定的()a x t 进行采样。 ○1()()15000s a f x t x n =以样本秒采样得到。 ()()11j x n X e ω画出及其频谱。 ○2()()11000s a f x t x n =以样本秒采样得到。 ()() 11j x n X e ω画出及其频谱。比较两种采样率下的信号频谱，并解释。（1）MATLAB 程序： N=10; Fs=5; T s=1/Fs; n=[-N:T s:N]; xn=exp(-abs(n)); w=-4*pi:0.01:4*pi; X=xn*exp(-j*(n'*w)); subplot(211) plot(n,xn); title('x_a(t)时域波形'); xlabel('t/ms');ylabel('x_a(t)'); axis([-10, 10, 0, 1]); subplot(212); plot(w/pi,abs(X)); title('x_a(t)频谱图'); xlabel('\omega/\pi');ylabel('X_a(e^(j\omega))');

ind = find(X >=0.03*max(X))*0.01; eband = (max(ind) -min(ind)); fprintf('等效带宽为%fKHZ\n',eband); 运行结果：等效带宽为12.110000KHZ

数字信号处理实验报告(实验1_4)

实验一 MATLAB 仿真软件的基本操作命令和使用方法实验容 1、帮助命令使用 help 命令，查找 sqrt （开方）函数的使用方法； 2、MATLAB 命令窗口（1）在MATLAB 命令窗口直接输入命令行计算3 1)5.0sin(21+=πy 的值；（2）求多项式 p(x) = x3 + 2x+ 4的根； 3、矩阵运算（1）矩阵的乘法已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B

（2）矩阵的行列式已知A=[1 2 3;4 5 6;7 8 9]，求A （3）矩阵的转置及共轭转置已知A=[1 2 3;4 5 6;7 8 9]，求A' 已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B' （4）特征值、特征向量、特征多项式已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征值、特征向量、特征多项式；

（5）使用冒号选出指定元素已知：A=[1 2 3;4 5 6;7 8 9]；求A 中第3 列前2 个元素；A 中所有列第2，3 行的元素； 4、Matlab 基本编程方法（1）编写命令文件：计算1+2+…+n<2000 时的最大n 值；

（2）编写函数文件：分别用for 和while 循环结构编写程序，求 2 的0 到15 次幂的和。

5、MATLAB基本绘图命令（1）绘制余弦曲线 y=cos(t)，t∈[0，2π]

（2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π] （3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求：（a）线形为点划线、颜色为红色、数据点标记为加号；（b）坐标轴控制：显示围、刻度线、比例、网络线（c）标注控制：坐标轴名称、标题、相应文本； >> clear;

通信原理2DPSK调制与解调实验报告

通信原理课程设计报告

一. 2DPSK基本原理 1.2DPSK信号原理 2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用Φ表示本码元初相与前一码元初相之差，并规定：Φ＝0表示0码，Φ＝π表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图2.1所示。图1.1 2DPSK信号在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。定义?Φ为本码元初相与前一码元初相之差，假设： ?Φ=0→数字信息“0”； ?Φ=π→数字信息“1”。则数字信息序列与2DPSK信号的码元相位关系可举例表示如下：数字信息： 1 0 1 1 0 1 1 1 0 1

DPSK信号相位：0 π π 0 π π 0 π 0 0 π 或：π 0 0 π 0 0 π 0 π π 0 2. 2DPSK信号的调制原理一般来说，2DPSK信号有两种调试方法，即模拟调制法和键控法。2DPSK 信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。图1.2.1 模拟调制法 2DPSK信号的的键控调制法框图如图1.2.2所示，其中码变换的过程为将输入的基带信号差分，即变为它的相对码。选相开关作用为当输入为数字信息“0”时接相位0，当输入数字信息为“1”时接pi。图1.2.2 键控法调制原理图码变换相乘载波 s(t)e o(t)

数字信号处理上机实验代码

文件名：tstem.m（实验一、二需要）程序： f unction tstem(xn,yn) %时域序列绘图函数 %xn:被绘图的信号数据序列，yn:绘图信号的纵坐标名称（字符串）n=0:length(xn)-1; stem(n,xn,'.'); xlabel('n');ylabel('yn'); axis([0,n(end),min(xn),1.2*max(xn)]); 文件名：tplot.m（实验一、四需要）程序： function tplot(xn,T,yn) %时域序列连续曲线绘图函数 %xn:信号数据序列,yn:绘图信号的纵坐标名称（字符串） %T为采样间隔 n=0;length(xn)-1;t=n*T; plot(t,xn); xlabel('t/s');ylabel(yn); axis([0,t(end),min(xn),1.2*max(xn)]); 文件名：myplot.m（实验一、四需要）

%(1)myplot;计算时域离散系统损耗函数并绘制曲线图。function myplot(B,A) %B为系统函数分子多项式系数向量 %A为系统函数分母多项式系数向量 [H,W]=freqz(B,A,1000) m=abs(H); plot(W/pi,20*log10(m/max(m)));grid on; xlabel('\omega/\pi');ylabel('幅度(dB)') axis([0,1,-80,5]);title('损耗函数曲线'); 文件名：mstem.m（实验一、三需要）程序： function mstem(Xk) %mstem(Xk)绘制频域采样序列向量Xk的幅频特性图 M=length(Xk); k=0:M-1;wk=2*k/M;%产生M点DFT对应的采样点频率（关于pi归一化值） stem(wk,abs(Xk),'.');box on;%绘制M点DFT的幅频特性图xlabel('w/\pi');ylabel('幅度'); axis([0,2,0,1.2*max(abs(Xk))]); 文件名：mpplot.m（实验一需要）

数字信号处理实验报告92885

目录实验1 离散时间信号的频域分析-----------------------2 实验2 FFT算法与应用-------------------------------7 实验3 IIR数字滤波器的设计------------------------12 实验4 FIR数字滤波器的设计------------------------17

实验1 离散时间信号的频域分析一．实验目的信号的频域分析是信号处理中一种有效的工具。在离散信号的时域分析中，通常将信号表示成单位采样序列δ（n ）的线性组合，而在频域中，将信号表示成复变量e n j ω-或 e n N j π2-的线性组合。通过这样的表示，可以将时域的离散序列映射到频域以便于进一步的处理。在本实验中，将学习利用MATLAB 计算离散时间信号的DTFT 和DFT,并加深对其相互关系的理解。二、实验原理 (1)DTFT 和DFT 的定义及其相互关系。序列x(n)DTFT 定义为()jw X e = ()n x n e ∞ =∞ ∑ω jn -它是关于自变量ω的复函数，且是以2π为周期的连续函数。 ()jw X e 可以表示为()()()jw jw jw re im X e X e jX e =+,其中，()jw re X e 和()jw im X e 分别是 ()jw X e 实部和虚部；还可以表示为 ()jw X e =()|()|jw j w X e e θ，其中， |()|jw X e 和{} ()arg ()j w X e ωθ=分别是()jw X e 的幅度函数和相位函数；它们都是ω的实函数，也是以2π为周期的周期函数。序列()x n 的N 点DFT 定义为2211 ()()()()N N j k j kn kn N N N N n X k X e x n e x n W π π ---==== ∑∑，()X k 是周期为N 的序列。()j X e ω与()X k 的关系：()X k 是对()j X e ω）在一个周期中的谱的等间隔N 点采样，即 2k |()()|jw w N X k X e π = = ，而()j X e ω 可以通过对()X k 内插获得，即

数字信号处理实验八(上机)报告概论

数字信号处理实验报告实验名称：实验八双线性变换法IIR 数字滤波器设计实验时间： 2014 年 12 月 9 日学号： 201211106134 姓名：孙舸成绩：评语：一、实验目的： 1、掌握用双线性变换法设计低通IIR 数字滤波器的基本原理和算法； 2、掌握用双线性变换法设计高通和带通IIR 数字滤波器的基本原理和算法； 3、进一步了解数字滤波器和模拟滤波器的频率响应特性。二、实验原理与计算方法： 1、双线性变换法设计IIR 低通数字滤波器的基本原理和算法双线性变换法设计数字滤波器，采用了二次映射的方法，就是先将整个s 平面压缩到s 1平面的一个T j T j π π ~-的横形条带范围内，然后再将这个条带映射到z 平面上，就能建立s 平面到z 平面的一一对应关系。对于低通数字滤波器，映射关系为 z z T z z T s ++-= +-=--11211211 (1) 其中T 为抽样周期。用双线性变换法设计低通IIR 数字滤波器的基本步骤，首先根据设计要求确定相应的模拟滤波器的传递函数)(s H a ，再应用(1)式得数字滤波器的传递函数)(z H z z T s a s H z H ++-==112) ()( (2) 通常可以给定的参数为：低通数字滤波器通带边界频率p p f πΩ21=、阻带边界频率 s s f πΩ21=和对应的通带衰减函数p α、阻带衰减函数s α。s 1平面中的模拟角频率1Ω与数字角频率ω的关系为线性关系T 1Ωω=，在计算模拟滤波器的阶数N 、极点s i 和传递函数)(s H a 之前，应作预畸变处理 2 2tan 22tan 2 1T f T T T πΩΩ== (3) 模拟滤波器的阶数N 、极点s i 和传递函数)(s H a 的计算方法与冲激响应不变法相同，

数字信号处理上机实验(第三版)

数字信号处理实验（Matlab）实验一: 系统响应及系统稳定性 %实验1：系统响应及系统稳定性 close all;clear all %======内容1：调用filter解差分方程，由系统对u(n)的响应判断稳定性====== A=[1,-0.9];B=[0.05,0.05]; %系统差分方程系数向量B和A x1n=[1 1 1 1 1 1 1 1 zeros(1,50)]; %产生信号x1(n)=R8(n) x2n=ones(1,128); %产生信号x2(n)=u(n) hn=impz(B,A,58); %求系统单位脉冲响应h(n) subplot(2,2,1);y='h(n)';tstem(hn,y); %调用函数tstem绘图 title('(a)系统单位脉冲响应h(n)');box on y1n=filter(B,A,x1n); %求系统对x1(n)的响应y1(n) subplot(2,2,2);y='y1(n)';tstem(y1n,y); title('(b)系统对R8(n)的响应y1(n)');box on y2n=filter(B,A,x2n); %求系统对x2(n)的响应y2(n) subplot(2,2,4);y='y2(n)';tstem(y2n,y); title('(c)系统对u(n)的响应y2(n)');box on %===内容2：调用conv函数计算卷积============================ x1n=[1 1 1 1 1 1 1 1 ]; %产生信号x1(n)=R8(n) h1n=[ones(1,10) zeros(1,10)]; h2n=[1 2.5 2.5 1 zeros(1,10)];

数字信号处理实验报告 (实验四)

实验四离散时间信号的DTFT 一、实验目的１. 运用MA TLAB 计算离散时间系统的频率响应。２. 运用MA TLAB 验证离散时间傅立叶变换的性质。二、实验原理（一）、计算离散时间系统的DTFT 已知一个离散时间系统∑∑==-= -N k k N k k k n x b k n y a 00)()(，可以用MA TLAB 函数frequz 非常方便地在给定的L 个离散频率点l ωω=处进行计算。由于)(ωj e H 是ω的连续函数，需要尽可能大地选取L 的值（因为严格说，在MA TLAB 中不使用symbolic 工具箱是不能分析模拟信号的，但是当采样时间间隔充分小的时候，可产生平滑的图形），以使得命令plot 产生的图形和真实离散时间傅立叶变换的图形尽可能一致。在MA TLAB 中，freqz 计算出序列{M b b b ,,,10 }和{N a a a ,,,10 }的L 点离散傅立叶变换，然后对其离散傅立叶变换值相除得到L l e H l j ,,2,1),( =ω。为了更加方便快速地运算，应将L 的值选为2的幂，如256或者512。例3.1 运用MA TLAB 画出以下系统的频率响应。 y(n)-0.6y(n-1)=2x(n)+x(n-1) 程序： 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));grid title(‘H(e^{j\omega}的实部’)) xlabel(‘\omega/ \pi ’); ylabel(‘振幅’)； subplot(2,1,1) plot(w/pi,imag(h));grid title(‘H(e^{j\omega}的虚部’)) xlabel(‘\omega/ \pi ’); ylabel(‘振幅’)；（二）、离散时间傅立叶变换DTFT 的性质。 1．时移与频移设 )]([)(n x FT e X j =ω, 那么

数字信号处理上机报告-一

数字信号处理第一次上机实验报告实验一：设给定模拟信号()1000t a x t e -=，的单位是ms 。 (1) 利用MATLAB 绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分量降低到峰值的3%以下的频谱）。 (2) 用两个不同的采样频率对给定的进行采样。 ○1 。 ○2 。比较两种采样率下的信号频谱，并解释。实验一MATLAB 程序：（1） N=10; Fs=5; Ts=1/Fs; n=[-N:Ts:N]; xn=exp(-abs(n)); w=-4*pi:0.01:4*pi; X=xn*exp(-j*(n'*w)); subplot(211) plot(n,xn); title('x_a(t)时域波形'); xlabel('t/ms');ylabel('x_a(t)'); t ()a x t ()()15000s a f x t x n =以样本秒采样得到。()()11j x n X e ω画出及其频谱()()11000s a f x t x n =以样本秒采样得到。()() 11j x n X e ω画出及其频谱

axis([-10, 10, 0, 1]); subplot(212); plot(w/pi,abs(X)); title('x_a(t)频谱图'); xlabel('\omega/\pi');ylabel('X_a(e ^(j\omega))'); ind = find(X >=0.03*max(X))*0.01; eband = (max(ind) -min(ind)); fprintf('等效带宽为 %fKHZ\n',eband); 运行结果：

数字信号处理上机实验答案(全)1

第十章上机实验数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。实验一系统响应及系统稳定性。实验二时域采样与频域采样。实验三用FFT 对信号作频谱分析。实验四 IIR 数字滤波器设计及软件实现。实验五 FIR 数字滤波器设计与软件实现实验六应用实验——数字信号处理在双音多频拨号系统中的应用任课教师根据教学进度，安排学生上机进行实验。建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR 数字滤波器设计及软件实现在。学习完第六章进行；实验五在学习完第七章后进行。实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。 10.1 实验一: 系统响应及系统稳定性 1.实验目的（1）掌握求系统响应的方法。（2）掌握时域离散系统的时域特性。（3）分析、观察及检验系统的稳定性。 2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。在计算机上适合用递推法求差分方程的解，最简单的方法是采用MA TLAB 语言的工具箱函数filter 函数。也可以用MATLAB 语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。或者系统的单位脉冲响应满足绝对可和的条件。系统的稳定性由其差分方程的系数决定。实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否都是有界输出，或者检查系统的单位脉冲响应满足绝对可和的条件。可行的方法是在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳定的[19]。系统的稳态输出是指当∞→n 时，系统的输出。如果系统稳定，信号加入系统后，系统输出的开始一段称为暂态效应，随n 的加大，幅度趋于稳定，达到稳态输出。注意在以下实验中均假设系统的初始状态为零。 3．实验内容及步骤（1）编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用filter 函数或conv 函数求解系统输出响应的主程序。程序中要有绘制信号波形的功能。（2）给定一个低通滤波器的差分方程为

数字信号实验报告 (全)

数字信号处理实验报告实验一：用 FFT 做谱分析一、实验目的 1、进一步加深 DFT 算法原理和基本性质的理解。 2、熟悉 FFT 算法原理和 FFT 子程序的应用。 3、学习用FFT 对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用 FFT 。二、实验原理用FFT 对信号作频谱分析是学习数字信号处理的重要内容。经常需要进行谱分析的信号是模拟信号和时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D 和分析误差。频谱分辨率直接和FFT 的变换区间N 有关，因为FFT 能够实现的频率分辨率是2π/N ≤D 。可以根据此时选择FFT 的变换区间N 。误差主要来自于用FFT 作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N 较大时离散谱的包络才能逼近于连续谱，因此N 要适当选择大一些。周期信号的频谱是离散谱，只有用整数倍周期的长度作FFT ，得到的离散谱才能代表周期信号的频谱。如果不知道信号周期，可以尽量选择信号的观察时间长一些。对模拟信号的频谱时，首先要按照采样定理将其变成时域离散信号。如果是模拟周期信号，也应该选取整数倍周期的长度，经过采样后形成周期序列，按照周期序列的谱分析进行。三、实验内容和步骤对以下典型信号进行谱分析： ?? ? ??≤≤-≤≤-=?? ? ??≤≤-≤≤+==其它n n n n n n x 其它n n n n n n x n R n x ,07 4, 330,4)(, 07 4, 830,1)() ()(3241 4() cos 4 x n n π = 5()cos(/4)cos(/8)x n n n ππ=+ 6() cos8cos16cos20x t t t t πππ=++

数字信号处理实验报告(全)

实验一、离散时间系统及离散卷积 1、单位脉冲响应源程序： function pr1() %定义函数pr1 a=[1,-1,0.9]; %定义差分方程y(n)-y(n-1)+0.9y(n-2)=x(n) b=1; x=impseq(0,-40,140); %调用impseq函数 n=-40:140; %定义n从-40 到140 h=filter(b,a,x); %调用函数给纵座标赋值 figure(1) %绘图figure 1 (冲激响应) stem(n,h); %在图中绘出冲激 title('冲激响应'); %定义标题为:'冲激响应' xlabel('n'); %绘图横座标为n ylabel('h(n)'); %绘图纵座标为h(n) figure(2) %绘图figure 2 [z,p,g]=tf2zp(b,a); %绘出零极点图 zplane(z,p) function [x,n]=impseq(n0,n1,n2)%声明impseq函数 n=[n1:n2]; x=[(n-n0)==0]; 结果： Figure 1: Figure 2:

2、离散系统的幅频、相频的分析源程序： function pr2() b=[0.0181,0.0543,0.0543,0.0181]; a=[1.000,-1.76,1.1829,-0.2781]; m=0:length(b)-1; %m从0 到3 l=0:length(a)-1; %l从0 到3 K=5000; k=1:K; w=pi*k/K; %角频率w H=(b*exp(-j*m'*w))./(a*exp(-j*l'*w));%对系统函数的定义 magH=abs(H); %magH为幅度 angH=angle(H); %angH为相位 figure(1) subplot(2,1,1); %在同一窗口的上半部分绘图 plot(w/pi,magH); %绘制w(pi)-magH的图形 grid; axis([0,1,0,1]); %限制横纵座标从0到1 xlabel('w(pi)'); %x座标为 w(pi) ylabel('|H|'); %y座标为 angle(H) title('幅度，相位响应'); %图的标题为:'幅度，相位响应' subplot(2,1,2); %在同一窗口的下半部分绘图 plot(w/pi,angH); %绘制w(pi)-angH的图形 grid; %为座标添加名称

通信原理2DPSK调制与解调实验报告

通信原理课程设计报告一. 2DPSK基本原理 1.2DPSK信号原理 2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用Φ表示本码元初相与前一码元初相之差，并规定：Φ＝0表示0码，

Φ＝π表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图2.1所示。图1.1 2DPSK信号在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。定义?Φ为本码元初相与前一码元初相之差，假设： ?Φ=0→数字信息“0”； ?Φ=π→数字信息“1”。则数字信息序列与2DPSK信号的码元相位关系可举例表示如下：数字信息： 1 0 1 1 0 1 1 1 0 1 DPSK信号相位：0 π π 0 π π 0 π 0 0 π 或：π 0 0 π 0 0 π 0 π π 0 2. 2DPSK信号的调制原理一般来说，2DPSK信号有两种调试方法，即模拟调制法和键控法。2DPSK 信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。

图1.2.1 模拟调制法 2DPSK信号的的键控调制法框图如图1.2.2所示，其中码变换的过程为将输入的基带信号差分，即变为它的相对码。选相开关作用为当输入为数字信息“0”时接相位0，当输入数字信息为“1”时接pi。图1.2.2 键控法调制原理图 3. 2DPSK信号的解调原理 2DPSK信号最常用的解调方法有两种，一种是极性比较和码变换法，另一种是差分相干解调法。 (1) 2DPSK信号解调的极性比较法它的原理是2DPSK信号先经过带通滤波器，去除调制信号频带以外的在信道中混入的噪声，再与本地载波相乘，去掉调制信号中的载波成分，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判决的到基带信号的差分码，再经过逆差分器，就得到了基带信号。它的原理框图如图1.3.1所示。码变换相乘载波 s(t)e o(t) 相乘器低通滤波器抽样判决器2DPSK 带通滤波器延迟T

数字信号实验报告

北京科技大学《信号系统与信号处理综合实验》实验报告学号：__________ 姓名：_____________________ 专业：____________ 年月日

目录： 1实验一CCS使用实验 2实验二、SEED-DTK6446 Linux开发环境搭建3实验三、Linux平台实验 4二、音频采集回放实验 5三、视频采集回放实验 6OSD图像叠加实验 7图像边缘检测实验

课程实验目的 1．数字信号处理是一门理论与实践并重的课程，在学习理论知识的同时再配合经典DSP实验，可以加深对数字信号处理软、硬件的理解与掌握。 2．接触并了解SEED-DTK6446实验箱，学会通过Linux操作平台，利用SEED-DTK6446实验箱完成一些经典的实验历程，加深对数字信号处理的了解。 3. 学习并掌握SEED-DTK6446 CCS开发环境的搭建，建立好所有编译测试环境，为下面的实验做好准备工作。实验一 CCS使用实验一、实验目的 1.熟悉CCS3.3集成开发环境，掌握工程的生成方法； 2.熟悉SEED-DTK6446实验环境； 3. 学习用标准C 语言编制程序； 4.掌握CCS3.3集成开发环境的调试方法；二、实验内容 1.DSP源文件的建立； 2.DSP程序工程文件的建立； 3. 学习使用CCS3.3集成开发工具的调试工具。三、实验步骤 1.创建源文件：选择File →New →Source File 命令；打开配套光盘\03. Examples of program\01.SEEE-DTK6446 CCS Examples\examples\3.1.1 math。 2.创建工程文件：点击Project-->New,创建新工程；点击Project选择add files to project，添加源程序math.c。 3. 设置编译与连接选项：点击Project选择Build Opitions； 4. 工程编译与调试：点击Project →Build all，对工程进行编译；点击File →load program，在弹出的对话框中载入debug 文件夹下的.out可执行文件；点击debug →Go Main回到C程序的入口；运行程序并观察输出结果。四．实验要求：

哈工大数字信号处理实验报告

实验一：用FFT 作谱分析实验目的： (1) 进一步加深DFT 算法原理和基本性质的理解(因为FFT 只是DFT 的一种快速算法，所以FFT 的运算结果必然满足DFT 的基本性质)。 (2) 熟悉FFT 算法原理和FFT 子程序的应用。 (3) 学习用FFT 对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT 。实验原理： DFT 的运算量：一次完整的DFT 运算总共需要2N 次复数乘法和(1)N N -复数加法运算，因而直接计算DFT 时，乘法次数和加法次数都和2N 成正比，当N 很大时，运算量很客观的。例如，当N=8时，DFT 运算需64位复数乘法，当N=1024时，DFT 运算需1048576次复数乘法。而N 的取值可能会很大，因而寻找运算量的途径是很必要的。 FFT 算法原理：大多数减少离散傅里叶变换运算次数的方法都是基于nk N W 的对称性和周期性。（1）对称性 ()*()k N n kn kn N N N W W W --==

（2）周期性 ()(mod`)()()kn N kn n N k n k N N N N N W W W W ++=== 由此可得 ()()/2 (/2)1 n N k N n k nk N N N N N k N k N N W W W W W W ---+?==?=-??=-? 这样： 1.利用第三个方程的这些特性，DFT 运算中有些项可以合并； 2.利用nk N W 的对称性和周期性，可以将长序列的DFT 分解为短序列的DFT 。前面已经说过，DFT 的运算量是与2N 成正比的，所以N 越小对计算越有利，因而小点数序列的DFT 比大点数序列的DFT 运算量要小。快速傅里叶变换算法正是基于这样的基本思路而发展起来的，她的算法基本上可分成两大类，即按时间抽取法和按频率抽取法。我们最常用的是2M N =的情况，该情况下的变换成为基2快速傅里叶变换。完成一次完整的FFT 计算总共需要 2log 2 N N 次复数乘法运算和2log N N 次复数加法运算。很明显，N 越大，FFT 的优点就越突出。实验步骤 (1) 复习DFT 的定义、性质和用DFT 作谱分析的有关内容。 (2) 复习FFT 算法原理与编程思想，并对照DIT-FFT 运算流图和程序框图，读懂本实验提供的FFT 子程序。 (3) 编制信号产生子程序，产生以下典型信号供谱分析用：

数字信号处理上机实验答案完整版

数字信号处理上机实验答案 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】