离散时间信号处理
实验1离散时间信号的产生与运算
数字信号处理 实验报告 班级： 学号： 姓名：word文档可自由复制编辑
实验1离散时间信号的产生与运算 一、实验目的 (1)了解离散时间信号的特点。 (2)掌握在计算机中生成及绘制各种常用离散时间信号序列的方法。 (3)掌握序列的加、减、乘、除和平移、反转、尺度变换等基本运算及计算机的 实现方法。 二、实验原理 信号是随时间变化的物理量，而计算机只能处理离散信号。离散信号是在某些不连续的时间上有信号值，而在其它时间点上没有定义的一类信号。离散信号一般可以由连续信号通过模数转换得到。 常用的离散信号有单位脉冲序列、单位阶跃序列、复指数序列、正弦信号序列、随机序列等。 离散信号的基本运算包括信号的加、减、乘、除。离散信号的时域变换包括信号的平移、反转、尺度变换等。 三、实验内容与方法 1、编写程序，生成如下数字信号：sqrt(2*k)u(k错误！未找到引用源。3)， δ(k+5)。 (1) f(k)=sqrt(2*k)u(k错误！未找到引用源。3) 代码： k=(1:10); n=3; u=[(k-n)>=0]; a=sqrt(2*k); stem(k,a.*u); title('sqrt(2*k)u(k 3)的图像'); xlabel('时间(k)');ylabel('幅值f(k)'); 运行图： word文档可自由复制编辑
(2) f(k)= δ(k+5) 代码： k1=-10;k2=0;k=k1:k2; n=-5; %单位脉冲出现的位置 f=[(k-n)==0]; stem(k,f,'filled');title('δ(k+5)序列的图像') xlabel('时间(k)');ylabel('幅值f(k)'); 运行图： word文档可自由复制编辑
离散时间信号的表示及运算
第2章 离散时间信号的表示及运算 2.1 实验目的 ● 学会运用MATLAB 表示的常用离散时间信号； ● 学会运用MATLAB 实现离散时间信号的基本运算。 2.2 实验原理及实例分析 2.2.1 离散时间信号在MATLAB 中的表示 离散时间信号是指在离散时刻才有定义的信号，简称离散信号，或者序列。离散序列通常用)(n x 来表示，自变量必须是整数。 离散时间信号的波形绘制在MATLAB 中一般用stem 函数。stem 函数的基本用法和plot 函数一样，它绘制的波形图的每个样本点上有一个小圆圈，默认是空心的。如果要实心，需使用参数“fill”、“filled”，或者参数“.”。由于MA TLAB 中矩阵元素的个数有限，所以MA TLAB 只能表示一定时间范围内有限长度的序列；而对于无限序列，也只能在一定时间范围内表示出来。类似于连续时间信号，离散时间信号也有一些典型的离散时间信号。 1. 单位取样序列 单位取样序列)(n δ，也称为单位冲激序列，定义为 ) 0()0(0 1)(≠=?? ?=n n n δ （12-1） 要注意，单位冲激序列不是单位冲激函数的简单离散抽样，它在n =0处是取确定的值1。在MATLAB 中，冲激序列可以通过编写以下的impDT .m 文件来实现，即 function y=impDT(n) y=(n==0); %当参数为0时冲激为1,否则为0 调用该函数时n 必须为整数或整数向量。 【实例2-1】 利用MATLAB 的impDT 函数绘出单位冲激序列的波形图。 解：MATLAB 源程序为 >>n=-3:3; >>x=impDT(n); >>stem(n,x,'fill'),xlabel('n'),grid on >>title('单位冲激序列') >>axis([-3 3 -0.1 1.1])
离散时间信号处理-知识点总结
离散时间信号及系统的DTFT 离散时间信号及系统的z变换 DFT的表达式 连续时间信号机系统的Fourier变换 时域-系统的因果性及稳定性P21、P32、P48 z域-系统的因果性及稳定性P110 抽样时间信号的频域表示P142 抽样离散信号与原连续信号的时域关系P150 连续信号、采样时间信号与离散信号的频谱关系P157 DTFT的对称性质P56 DTFT的理论及性质P59 DTFT变换对P62 DTFT与原连续信号的频谱关系P147 离散Fourier级数DFS性质P550 DFT性质P576 线性循环卷积P576 重叠保留法、相加法P582 窗函数效应P698 时间依赖Fourier变换P714 Decimation in Time P640、P645 Decimation in Frequency P649、P651 z-Transform变换对P104 z-Transform性质P126
LTI的典型单位冲激响应P31 LTI的特征函数及特征根P40、P46 全通系统P274 最小相位系统P280 线性相位系统P291 线性相位系统与最小相位系统的关系P308 FIR滤波器窗函数P469 FIR滤波器最佳逼近P486 降采样频谱P168、P170 升采样频谱P172、P174 随机信号理论Appendix-A 随机信号的自协方差及自相关序列的时域频域性质P65 平稳随机信号的Fourier分析P723 AD噪声分析P193 数字滤波器中的舍入误差噪声P391 有限字长效应P370 系数量化误差P377 FFT有限寄存器长效应P661 极限循环P415
离散时间系统的时域分析
第七章离散时间系统的时域分析 §7-1 概述 一、离散时间信号与离散时间系统 离散时间信号：只在某些离散的时间点上有值的 信号。 离散时间系统：处理离散时间信号的系统。 混合时间系统：既处理离散时间信号，又处理连 续时间信号的系统。 二、连续信号与离散信号 连续信号可以转换成离散信号，从而可以用离散时间系统（或数字信号处理系统）进行处理： 三、离散信号的表示方法：
1、 时间函数：f(k)<——f(kT)，其中k 为序号，相当于时间。 例如：)1.0sin()(k k f = 2、 （有序）数列：将离散信号的数值按顺序排列起来。例如： f(k)={1,0.5,0.25,0.125,……,} 时间函数可以表达任意长（可能是无限长）的离散信号，可以表达单边或双边信号，但是在很多情况下难于得到；数列的方法表示比较简单，直观，但是只能表示有始、有限长度的信号。 四、典型的离散时间信号 1、 单位样值函数：? ??==其它001)(k k δ 下图表示了)(n k ?δ的波形。
这个函数与连续时间信号中的冲激函数 )(t δ相似，也有着与其相似的性质。例如： )()0()()(k f k k f δδ=， )()()()(000k k k f k k k f ?=?δδ。 2、 单位阶跃函数：? ??≥=其它001)(k k ε 这个函数与连续时间信号中的阶跃函数)(t ε相似。用它可以产生（或表示）单边信号（这里称为单边序列）。 3、 单边指数序列：)(k a k ε
比较：单边连续指数信号：)()()(t e t e t a at εε=，其 底一定大于零，不会出现负数。 (a) 0.9a = (d) 0.9a =? (b) 1a = (e) 1a =? (c) 1.1a = (f) 1.1a =?
华中师范大学离散时间信号处理A卷
院（系 ）： 专业： 年级： 学生 姓名： 学号： --- -- -- --- -- -- -- --- -- -- -- --- -- -- -- --- -- -- -- -- -- -- 密 -- -- -- -- --- -- -- -- --- -- -- -- --- -- -- - 封 --- -- -- --- -- -- -- --- -- -- -- -- -- 线 ---- -- -- -- --- -- -- -- --- -- -- -- --- -- -- -- --- -- -- -- --- -- -- -- --
第 1 页共 3 页
------------------------------------------------- 密 ---------------------------------- 封 ----------------------------- 线 ---------------------------------------------------------
第 2 页(共 3 页)
第 3 页(共 3 页) ------------------------------------------------- 密 ---------------------------------- 封 ----------------------------- 线 ---------------------------------------------------------
离散时间信号与离散时间系统
§7-1 概述 一、 离散时间信号与离散时间系统 离散时间信号：只在某些离散的时间点上有值的 信号。 离散时间系统：处理离散时间信号的系统。 混合时间系统：既处理离散时间信号，又处理连 续时间信号的系统。 二、 连续信号与离散信号 连续信号可以转换成离散信号，从而可以用离散时间系统（或数字信号处理系统）进行处理： 三、 离散信号的表示方法： 1、 时间函数：f(k)<——f(kT)，其中k 为序号，相当于时间。 例如：)1.0sin()(k k f = 2、 （有序）数列：将离散信号的数值按顺序排列起来。例如： f(k)={1,0.5,0.25,0.125,……,} 时间函数可以表达任意长（可能是无限长）的离散信号，可以表达单边或双边信号，但是在很多情况下难于得到；数列的方法表示比较简单，直观，但是只能表示有始、有限长度的信号。 四、 典型的离散时间信号 1、 单位样值函数： ?? ?==其它001)(k k δ 下图表示了)(n k -δ的波形。 连续信号 离散信号 数字信号 取样 量化
这个函数与连续时间信号中的冲激函数)(t δ相似，也有着 与其相似的性质。例如： )()0()()(k f k k f δδ=， )()()()(000k k k f k k k f -=-δδ。 2、 单位阶跃函数： ?? ?≥=其它001)(k k ε 这个函数与连续时间信号中的阶跃函数) (t ε相似。用它可以产生（或表示）单边信号（这里称为单边序列）。 3、 单边指数序列： )(k a k ε 比较：单边连续指数信号： )( )()(t e t e t a at εε=，其底一定大于零，不会出现负数。 4、 单边正弦序列：)()cos(0k k A εφω+ (a) 0.9a = (d) 0.9a =- (b) 1a = (e) 1a =- (c) 1.1a = (f) 1.1a =-
离散时间信号的表示及运算
实验一 离散时间信号的表示及运算 一、实验目的 学会运用MATLAB 表示的常用离散时间信号；学会运用MATLAB 实现离散时间信号的基本运算。 二、实验原理 (一) 离散时间信号在MATLAB 中的表示 离散时间信号是指在离散时刻才有定义的信号，简称离散信号，或者序列。离散序列通常用)(n x 来表示，自变量必须是整数。离散时间信号的波形绘制在MATLAB 中一般用stem 函数。stem 函数的基本用法和plot 函数一样，它绘制的波形图的每个样本点上有一个小圆圈，默认是空心的。如果要实心，需使用参数“fill ”、“filled ”，或者参数“.”。由于MATLAB 中矩阵元素的个数有限，所以MATLAB 只能表示一定时间范围内有限长度的序列；而对于无限序列，也只能在一定时间范围内表示出来。类似于连续时间信号，离散时间信号也有一些典型的离散时间信号。(二) 离散时间信号的基本运算 对离散时间序列实行基本运算可得到新的序列，这些基本运算主要包括加、减、乘、除、移位、反折等。两个序列的加减乘除是对应离散样点值的加减乘除，因此，可通过MATLAB 的点乘和点除、序列移位和反折来实现，与连续时间信号处理方法基本一样。三、实验内容（包括代码与产生的图形） 1. 试用MATLAB 命令分别绘出下列各序列的波形图。 （1）()()n u n x n ?? ? ??=21 （2）()()n u n x n 2=
（3）()()n u n x n ?? ? ??-=21 （4）()()()n u n x n 2-= （5）()()121 -=-n u n x n （6）()()n u n x n 1 21-? ? ? ??= (1)、(2) n=-3:8; a=1/2; x=a.^n.*uDT(n); subplot(221); stem(n,x,'fill','r'),xlabel('n'),grid on title('(1)x(n)=(1/2)^{n}*U(n)') axis([-3 8 -0.1 1.1]) n1=-3:8; b=2; x=b.^n1.*uDT(n1); subplot(222); stem(n1,x,'fill','r'),xlabel('n'),grid on title('(2)x(n)=(2)^{n}*U(n)') axis([-3 4.5 -1.5 18]) 分析：(1)该信号为指数衰减序列与阶跃序的乘积，当n<0时,U(n)=0，所以该 信号为零；当n=0时，U(n)=1,n ?? ? ??21=1，该信号为1；当n>0,U(n)=1,该信号呈现
离散时间信号处理-知识点总结
离散时间信号及系统的DTFT X e jw = x [n ]∞ n =?∞e ?jwn x n = 1 2π X (e jw )e jwn dw π ?π 离散时间信号及系统的z 变换 X (z )= x [n ]∞ n =?∞ z ?n DFT 的表达式W N =e ?j(2π/N) X [k ]= x [n ]N?1 n =0W N kn x [n ]= 1 X [k ]N?1 k =0W N ?kn 连续时间信号机系统的Fourier 变换 F w = f (t )e ?jwt dt ∞ ?∞ 时域-系统的因果性及稳定性 P21、P32、P48 z 域-系统的因果性及稳定性 P110 抽样时间信号的频域表示 P142 抽样离散信号与原连续信号的时域关系 P150 连续信号、采样时间信号与离散信号的频谱关系 P157 DTFT 的对称性质 P56 DTFT 的理论及性质 P59 DTFT 变换对 P62 DTFT 与原连续信号的频谱关系 P147 离散Fourier 级数DFS 性质 P550 DFT 性质 P576 线性循环卷积 P576 重叠保留法、相加法 P582 窗函数效应 P698 时间依赖Fourier 变换 P714 Decimation in Time P640、P645 Decimation in Frequency P649、P651 z-Transform 变换对 P104 z-Transform 性质 P126
LTI的典型单位冲激响应P31 LTI的特征函数及特征根P40、P46 全通系统P274 最小相位系统P280 线性相位系统P291 线性相位系统与最小相位系统的关系P308 FIR滤波器窗函数P469 FIR滤波器最佳逼近P486 降采样频谱P168、P170 升采样频谱P172、P174 随机信号理论Appendix-A 随机信号的自协方差及自相关序列的时域频域性质P65 平稳随机信号的Fourier分析P723 AD噪声分析P193 数字滤波器中的舍入误差噪声P391 有限字长效应P370 系数量化误差P377 FFT有限寄存器长效应P661 极限循环P415
离散时间信号与系统
离散时间信号与系统
实验：离散时间信号与系统的时域分析 一、实验目的 1、熟悉和掌握常用的用于信号与系统时域仿真分析的MATLAB函数； 2、掌握离散时间信号的MATLAB产生，掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB编程； 3、牢固掌握系统的单位序列响应的概念，掌握MATLAB描述LTI系统的常用方法及有关函数，并学会利用MATLAB求解LTI系统响应，绘制相应曲线。 基本要求：掌握用MATLAB描述离散时间信号的方法，能够编写MATLAB程序，实现各种信号的时域变换和运算，并且以图形的方式再现各种信号的波形。掌握线性时不变离散系统的时域数学模型用MATLAB描述的方法，掌握线性常系数差分方程的求解编程。 二、实验原理 信号（Signal）一般都是随某一个或某几个独立变量的变化而变化的，例如，温度、压力、
声音，还有股票市场的日收盘指数等，这些信号都是随时间的变化而变化的，还有一些信号，例如在研究地球结构时，地下某处的密度就是随着海拔高度的变化而变化的。一幅图片中的每一个象素点的位置取决于两个坐标轴，即横轴和纵轴，因此，图像信号具有两个或两个以上的独立变量。 在《信号与系统》课程中，我们只关注这种只有一个独立变量（Independent variable）的信号，并且把这个独立变量统称为时间变量（Time variable），不管这个独立变量是否是时间变量。 在自然界中，大多数信号的时间变量都是连续变化的，因此这种信号被称为连续时间信号（Continuous-Time Signals）或模拟信号（Analog Signals），例如前面提到的温度、压力和声音信号就是连续时间信号的例子。但是，还有一些信号的独立时间变量是离散变化的，这种信号称为离散时间信号。前面提到的股票市场的日收盘指数，由于相邻两个交易日的日收盘指数相隔24小时，这意味着日收盘指数的时间变量是不连续的，因此日收盘指数是离散时间信号。 而系统则用于对信号进行运算或处理，或者
数字信号处理课程设计--离散时间信号处理
离散时间信号处理 课程设计
某雷达系统接收机框架如图1所示。接收机输入信号,其中，中心频率=40MHz ，为相位调制函数(调制带宽B=2MHz)，信号谱如图2所示。输出基带信号I[n]、Q[n]数据率(采样率)为2.5MHz 。 图1 接收机原理框架 图2 接收信号参数 1. 请设计 ADC 的采样率fs ，画出x[n]的频谱X(e jw ) 2. 请推导正交解调器的输出S_I[n]和S_Q[n]的表达式，画出频谱示意图；若要将信号谱搬移到零中频，请确定NCO 的频率W 0 3. 请设计一个FIR 线性相位数字LPF 对正交解调的输出进行处理，要求杂散抑制比超过50dBc ，确定滤波器设计指标、给出设计过程和结果。 4. 请确定采样率变换模块的参数(抽取或内插系数)，画出输出信号谱。 f Xc(f)
1.设计ADC的采样率fs，画出x[n]的频谱X(e jw) Adc采样率fs>=2*f0，fs取100MHZ。 φ(t)为相位调制函数(调制带宽B=2MHz)，取φ(t)=cos(2π*1000000*t)作为单音调制信号，则,,,程序如下： clear; clc; clf; n0=20000;%点数 sprate=100000000;%采样频率100M t0=n0/sprate;%总时长 fca=40000000;%载波频率 fsi=1000000;%信号频率 t=linspace(0,t0,n0); Dsb=cos(2*pi*fca*t+cos(fsi*2*pi*t)); subplot(3,1,1); plot(t(1:100),Dsb(1:100)); title('时域'); xlabel('t');ylabel('voltage'); grid on; legend('modulated signal'); subplot(3,1,2); w=linspace(-1,1,n0); plot(w,fftshift(abs(fft(Dsb)))); title('幅频'); xlabel('w');ylabel('|X(ejw)|'); grid on; legend('modulated signal'); subplot(3,1,3); plot(w,fftshift(angle(fft(Dsb)))); title('相频'); xlabel('w');ylabel('∠X(ejw)'); grid on; legend('modulated signal');