当前位置:文档之家› 时域离散信号的产生与运算

时域离散信号的产生与运算

时域离散信号的产生与运算
时域离散信号的产生与运算

典型时域离散序列的产生与简单运算

1. 单位冲激序列

程序1:

function [x,n]=impseq(n0,n1,n2)

% generates x(n)=delta(n-n0); n1<=n<=n2

n=[n1:n2];

x=[(n-n0)==0];

调用:[x,n]=impseq(0,-3,4);

stem(n,x)

程序2:

n1=-3;n2=4;n0=0;

n=n1:n2;

x=[n==n0];

stem(n,x,'filled');

axis([n1,n2,0,1.1*max(x)]);

xlabel('时间(n)');ylabel('幅度x(n)');

title('单位脉冲序列');

2. 单位阶跃序列

程序:

n1=-3;n2=4;n0=0;

n=n1:n2;

x=[n>=n0];

stem(n,x,'filled');

axis([n1,n2,0,1.1*max(x)]);

xlabel('时间(n)');ylabel('幅度x(n)');

title('单位阶跃序列');

3. 矩形序列

程序:

10()00n n n δ=?=?≠?1≥0()00n u n n ?=?

n=[-10:10];

xn1=[(n-0)>=0];

xn2=[(n-4)>=0]; %定义两个阶跃序列;

xn=xn1-xn2; 两个阶跃序列之差得到矩形序列;

stem(n,xn,'.');

xlabel('时间(n)');ylabel('幅度x(n)');

title(‘矩形序列');

4. 正弦序列

程序:

n=0:20;

xn=sin(pi/4*n);

stem(n,xn,'.');

xlabel('时间(n)');ylabel('幅度x(n)');

title(‘正弦序列');

5. 指数序列

程序:

n=[0:20];

x=(0.9).^n;

stem(n,x);

xlabel('时间(n)');ylabel('幅度x(n)');

title(‘指数序列');

6. 对conv 进行简单的扩展conv_m ,可以完成任意位置序列的卷积.

对于有限长序列x (n ),h (n ),它们分别的区域为[n xb,n xe]和[n hb,n he],则卷积后的区域为

[n xb+n hb,n xe+n he]

程序:

function[y,ny]=conv_m(x,nx,h,nh)

nyb=nx(1)+nh(1);

nye=nx(length(x))+nh(length(h));

ny=[nyb:nye];

y=conv(x,h);

调用:

x=[3,11,7,0,-1,4,2];

h=[2,3,0,-5,2,1];

nx=[-3:3];

nh=[-1:4];

[y,ny]=conv_m(x,nx,h,nh) ()sin()x n A n ωθ=+n

a n x =)(

实验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文档可自由复制编辑

用MATL新编实现常用的离散时间信号及其时域运算

用M A T L新编实现常用的离散时间信号及其 时域运算 公司内部编号:(GOOD-TMMT-MMUT-UUPTY-UUYY-DTTI-

实验四用MATLAB实现常用的离散时间信号及其时域运算 —— 摘要:在MATLAB中,只能用向量来表示离散时间信号。与连续信号不同,离散时 间信号无法用符号运算来表示。用适当的MATLAB语句表示出信号后,就可以利用MATLAB的绘图命令stem来绘出直观的信号波形图,stem是专门用于绘制离散时 间信号的。在MATLAB中离散序列的时域运算和变换不能用符号运算来实现,而必 须用向量表示的方法,即在MATLAB中离散序列的相加、相乘需表示成两个向量的 相加、相乘,因而参加运算的两序列向量必须有相同的维数。 一、实验目的:(1)学习MATLAB语言及其常用指令; (2)学习和掌握用MATLAB语言产生离散时间信号的编程方法; (3)通过编程绘制出离散时间信号的波形,加深理解信号的时域运 算。 二、实验内容:(1)运用MATLAB的绘图指令绘制离散时间信号; (2)用MATLAB语言实现离散时间信号的时域运算。 三、实验原理:(1)单位阶跃序列和单位样值序列。 离散时间信号只在某些离散的瞬时给出信号的值,因此,它是时间上不连续的 序列。单位阶跃序列和单位样值序列在离散时间信号与系统的分析中是两个非 常典型的序列,分别记为u(n)和δ(n)。它们的定义分别如下: 1 n≥0 1 n≥0 u(n)= δ(n)= 0 n<0 0 n≠0

若单位阶跃序列的起始点为n0,单位样值序列出现在n0时刻,则表达式分别为: 1 n≥n0 1 n=n0 u(n-n0)= δ(n-n0)= 0 n

随机信号处理实验报告讲解

随机信号处理实验报告

目录 一、实验要求: (3) 二、实验原理: (3) 2.1 随机信号的分析方法 (3) 2.2 随机过程的频谱 (3) 2.3 随机过程的相关函数和功率谱 (4) (1)随机信号的相关函数: (4) (2)随机信号的功率谱 (4) 三、实验步骤与分析 (5) 3.1实验方案 (5) 3.2实验步骤与分析 (5) 任务一:(s1 变量)求噪声下正弦信号的振幅和频率 (5) 任务二:(s1 变量)求噪声下正弦信号的相位 (8) 任务三:(s1 变量)求信号自相关函数和功率谱 (11) 任务四:(s变量)求噪声下信号的振幅和频率 (14) 任务五:(s变量)求信号的自相关函数和功率谱 (17) 3.3实验结果与误差分析 (19) (1)实验结果 (19) (2)结果验证 (19) (3)误差分析 (21) 四、实验总结和感悟 (22) 1、实验总结 (22) 2、实验感悟 (23) 五、附低通滤波器的Matlab程序 (23)

一、实验要求: (学号末尾3,7)两个数据文件,第一个文件数据中只包含一个正弦波,通过MA TLAB 仿真计算信号频谱和功率谱来估计该信号的幅度,功率,频率和相位?对第二个文件数据估计其中正弦波的幅度,功率和频率?写出报告,包含理论分析,仿真程序及说明,误差精度分析等。第一文件调用格式load FileDat01_1 s1,数据在变量s1中;第二文件调用格式load FileDat01_2 s ,数据在变量s 中。 二、实验原理: 2.1 随机信号的分析方法 在信号与系统中,我们把信号分为确知信号和随机信号。其中随机信号无确定的变化规律,需要用统计特新进行分析。这里我们引入随机过程的概念,所谓随机过程就是随机变量的集合,每个随机变量都是随机过程的一个取样序列。 随机过程的统计特性一般采用随机过程的分布函数和概率密度来描述,他们能够对随机过程作完整的描述。但由于在实践中难以求得,在工程技术中,一般采用描述随机过程的主要平均统计特性的几个函数,包括均值、方差、相关函数、频谱及功率谱密度等来描述它们。 2.2 随机过程的频谱 信号频谱分析是采用傅立叶变换将时域信号x(t)变换为频域信号X(f),从而帮助人们从另一个角度来了解信号的特征。时域信号x(t)的傅氏变换为: ()()2j ft X f x t e dt π+∞ --∞ =? 信号的时域描述只能反映信号的幅值随时间的变化情况,除只有一个频率分量的简谐波外,一般很难明确揭示信号的频率组成和各频率分量的大小。信号的频谱X(f)代表了信号在不同频率分量处信号成分的大小,它能够提供比时域信号波形更直观,丰富的信息。 在实际的控制系统中能够得到的是连续信号x(t)的离散采样值x(nT),因此需要利用离散信号x(nT)来计算信号x(t)的频谱。

离散信号的产生及运算

离散信号的产生及运算

实验一离散信号的产生及运算 一.实验目的: 1.复习和巩固数字信号处理中离散信号的产生和运算 2.学习和掌握用MATLAB 产生离散信号的方法 3.学习和掌握用MATLAB 对离散信号进行运算 二.实验原理 1.用MATLAB 函数产生离散信号 信号是数字信号处理的最基本内容。没有信号,数字信号处理就没了工作对象。MATLAB7.0 内部提供了大量的函数,用来产生常用的信号波形。例如,三角函数(sin,cos), 指数函数(exp),锯齿波函数(sawtooth), 随机数函数(rand)等。 ⑴产生被噪声污染的正弦信号 用随机数函数产生污染的正弦信号。 ⑵产生单位脉冲序列和单位阶跃序列 按定义,单位脉冲序列为 0 0 0 1, ( ) 0,n n n n n n 单位阶跃序列为 。 0 0 0 1, ( ) 0,n n u n n n n ⑶矩形脉冲信号:在MATLAB 中用rectpuls 函数来表示,其调用形式为: y=rectpuls(t,width),用以产生一个幅值为1,宽度为width,相对于t=0 点左右对称的矩形波信号,该函数的横坐标范围由向量t 决定,是以t=0 为中心向左右各展开width/2 的范围,width 的默认值为1。 例:以t=2T(即t-2×T=0)为对称中心的矩形脉冲信号的MATLAB 源程序如下:(取T=1) t=0:0.001:4; T=1; ft=rectpuls(t-2*T,2*T); plot(t,ft);grid on; axis([0 4 –0.5 1.5]); ⑷周期性矩形波(方波)信号在MATLAB 中用square 函数来表示,其调用形式为:y=square(t,DUTY),用以产生一个周期为2π、幅值为±1 的周期性方波信号,其中的DUTY 参数表示占空比,即在信号的一个周期中正值所占的百分比。例如频率为30Hz 的周期性方波信号的MATLAB 参考程序如下: t=-0.0625:0.0001:0.0625; y=square(2*pi*30*t,75); plot(t,y);axis([-0.0625 0.0625 –1.5 1.5]);grid on ; 2.MATLAB 中信号的运算 乘法和加法:离散信号之间的乘法和加法,是指它的同序号的序列值逐项对应相乘和相加。

离散信号与系统时域分析

目录 第1章设计任务及要求 (1) 1.1课程设计内容 (1) 1.2课程设计要求 (1) 第2章设计原理 (2) 2.1离散信号与系统的时域分析设计 (2) 2.1.1描写系统特性的方法介绍 (2) 2.1.2系统的时域特性 (2) 第3章设计实现 (3) 3.1实验内容与方法 (3) 3.1.1实验内容 (3) 第4章设计结果及分析 (3) 4.1程序设计结果及分析 (4) 总结 (7) 参考文献: (7) 附录: (8)

第1章 设计任务及要求 1.1课程设计内容 编制Matlab 程序,完成以下功能,产生系统输入信号;根据系统差分方程求解单位脉冲响应序列;根据输入信号求解输出响应;用实验方法检查系统是否稳定;绘制相关信号的波形。具体要求如下: (1) 给定一个低通滤波器的差分方程为 ()0.05()0.05(1)0.9(1)y n x n x n y n =+-+- 输入信号分别为182()=()()()x n R n x n u n =, ① 分别求出系统响应,并画出其波形。 ② 求出系统的单位脉冲响应,画出其波形。 (2) 给定系统的单位脉冲响应为1102()=()()() 2.5(1) 2.5(2)(3)h n R n h n n n n n δδδδ=+-+-+-,用线性卷积法求18()=()x n R n 分别对系统h1(n)和h2(n)的输出响应,并画出波形。 (3) 给定一谐振器的差分方程为() 1.8237(1)-0.9801(2)()(2)o o y n y n y n b x n b x n =--++-令b0=1/100.49,谐振器的谐振频率为0.4rad 。 1) 用实验方法检查系统是否稳定。输入信号为u(n)时,画出系统输出波形。 2) 给定输入信号为()=sin(0.014)sin(0.4)x n n n +求出系统的输出响应,并画出其波形。 1.2课程设计要求 1. 要求独立完成设计任务。 2. 课程设计说明书封面格式要求见《天津城市建设学院课程设计教学工作规范》附表1 3. 课程设计的说明书要求简洁、通顺,计算正确,图纸表达内容完整、清楚、规范。 4. 简述离散系统时域分析方法和通过实验判断系统稳定性的方法;完成以上设计实验并对结果进行分析和解释;打印程序清单和要求画出的信号波形;写出本次课程设计的收获和体会。 5. 课设说明书要求: 1) 说明题目的设计原理和思路、采用方法及设计流程。 2) 详细介绍运用的理论知识和主要的Matlab 程序。 3) 绘制结果图形并对仿真结果进行详细的分析。

典型连续信号和离散信号时域波形图

一.典型连续信号和离散信号的时域波形。 1.单边指数信号)()(t u Ae t y t α=; 2.单位冲激信号)()(0t t t y +=δ; 3.单位阶跃信号)()(0t t u t y +=; 4.矩形脉冲信号)]()([)(21t t u t t u A t y +-+?=; 5.正弦信号)()sin()(t u t A t y ω?=; 6.单位序列)()(0n n n y +=δ; 7.单位阶跃序列)()(0n n u n y +=; 8.单位矩形序列)()()(21n n u n n u n y +-+=; 9.指数序列)()(n u a A n y n ?=; 10.正弦序列)()sin()(n u n A n y ω?=。

单边指数信号 function zhishu(A,a,t1,t2,dt) t1=0 t2=10 A=1 A=-0.4 dt=0.01 t=t1:dt:t2; y=A*exp(a*t); plot(t,y) axis([t1,t2,0,1.2]) xlabel('t') ylabel('y(t)') title(' 单边指数信号') 单位冲激信号 function chongji(t1,t2,t0) dt=0.01; t1=10; t2=-5; t=t1:dt:t2; n=length(t); x=zeros(1,n); x(1,(-t0-t1)/dt+1)=1/dt; stairs(t,x); axis([t1,t2,0,1.2/dt]) xlabel('t') ylabel('y(t)') title('单位冲激信号')

离散时间信号表与运算

离散时间信号表与运算

————————————————————————————————作者:————————————————————————————————日期:

实验一 离散时间信号的表示与运算 一 实验目的 1、熟悉MATLAB 的绘图函数; 2、掌握单位取样序列、单位阶跃序列、矩形序列和正余弦序列的产生方法; 3、掌握离散时间信号基本运算的MATLAB 实现; 4、掌握离散时间信号线性卷积和运算的MATLAB 实现。 二 实验设备 1、计算机 2、MA TLAB R2007a 仿真软件 三 实验原理 1)序列相加和相乘 设有序列)(1n x 和)(2n x ,它们相加和相乘如下: ) ()()()()()(2121n x n x n x n x n x n x ?=+= 注意,序列相加(相乘)是对应序列值之间的相加(相乘),因此参加运算的两个序列必须具有相同的长度,并且保证位置相对应。如果不相同,在运算前应采用zeros 函数将序列左右补零使其长度相等并且位置相对应。在MATLAB 中,设序列用x1和x2表示,序列相加的语句为:x=x1+x2;然而要注意,序列相乘不能直接用x=x1*x2,该式表示两个矩阵的相乘,而不是对应项的相乘。对应项之间相乘的实现形式是点乘“.*”,实现语句为:x=x1.*x2。 2)序列翻转 设有序列:)()(n x n y -=,在翻转运算中,序列的每个值以n=0为中心进行翻转,需要注意的是翻转过程中序列的样值向量翻转的同时,位置向量翻转并取反。MATLAB 中,翻转运算用fliplr 函数实现。设序列)(n x 用样值向量x 和位置向量nx 表述,翻转后的序列 )(n y 用样值向量y 和位置向量ny 描述。 3)序列的移位 移位序列)(n x 的移位序列可表示为:)()(0n n x n y -=,其中,00>n 时代表序列右移 0n 个单位;00

实验一 时域离散信号与系统变换域分析(2015)资料

实验一 时域离散信号与系统变换域分析 一、实验目的 1.了解时域离散信号的产生及基本运算实现。 2.掌握离散时间傅里叶变换实现及系统分析方法。 3. 熟悉离散时间傅里叶变换性质。 4. 掌握系统Z 域分析方法。 5. 培养学生运用软件分析、处理数字信号的能力。 二、实验设备 1、计算机 2、Matlab7.0以上版本 三、实验内容 1、对于给定的时域离散信号会进行频谱分析,即序列的傅里叶变换及其性质分析。 2、对于离散系统会进行频域分析及Z 域分析。包括频谱特性、零极点画图、稳定性分析。 3、对于差分方程会用程序求解,包括求单位冲击序列响应,零输入响应、零状态响应、全响应,求其系统函数,及其分析。 4、信号时域采样及其频谱分析,序列恢复。 5、扩展部分主要是关于语音信号的读取及其播放。 四、实验原理 1、序列的产生及运算 在Matlab 中自带了cos 、sin 、exp (指数)等函数,利用这些函数可以产生实验所需序列。 序列的运算包括序列的加法、乘法,序列)(n x 的移位)(0n n x -,翻褶)(n x -等。序列的加法或乘法指同序号的序列值逐项对应相加或相乘,但Matlab 中“+”“.*”运算是对序列的值直接进行加或乘,不考虑两序列的序号是否相同,因此编程时考虑其序号的对应。 2、序列的傅里叶变换及其性质 序列的傅里叶变换定义:)(|)(|)()(ω?ωωω j j n n j j e e X e n x e X ==∑∞-∞=-,其幅度特性为|)(|ωj e X , 在Matlab 中采用abs 函数;相位特性为)(ω?,在Matlab 中采用angle 函数。 序列傅里叶变换的性质:

时域离散信号的产生与基本运算

实验一 时域离散信号的产生与基本运算 一、实验目的 1、了解常用的时域离散信号及其特点。 2、掌握MATLAB 产生常用时域离散信号的方法。 3、掌握时域离散信号简单的基本运算方法。 二、实验内容 1、自己设定参数,分别表示并绘制单位抽样序列、单位阶跃序列、正弦序列、 实指数序列、随机序列。 2、自己设定参数,分别表示并绘制信号移位、信号相加、信号相乘、信号翻转、 信号和、信号积、信号能量。 3、已知信号 (1) 描绘)(n x 序列的波形。 (2) 用延迟的单位脉冲序列及其加权和表示)(n x 序列。 (3) 描绘以下序列的波形:)2()(),2(2)(),2(2)(321n x n x n x n x n x n x -=+=-= 三、实现步骤 1、自己设定参数,分别表示并绘制单位抽样序列、单位阶跃序列、正弦序列、 实指数序列、随机序列。 (1)单位抽样序列 程序: x=zeros(1,10);

x(2)=1; stem(x,'filled') axis([0,10,-0.2,1]); title('μ¥??3é?ùDòáD'); -0.20 0.2 0.4 0.6 0.8 图 1 (2)单位阶跃序列 程序: N=10; u=ones(1,N); stem(u,'filled') axis([-10,10,0,1]); title('μ¥???×??DòáD');

00.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 单位阶跃序列 图 2 (3)正弦序列 程序: x=-20:1:20; y=sin(0.2*pi.*x+0.5*pi); stem(x,y,'filled'); axis([-20,20,-2,2]); title('?y?òDòáD');

随机信号的功率谱分析 (DEMO)

信号的功率谱分析 1、功率谱密度函数的定义 对于随机信号)(t x ,由于其任一样本函数都是时间的无限的函数,一般不能满足傅里叶变换的存在条件(即积分?∞ ∞-dt t x )(必须收敛)。如果将样本函数取在一个有限区间]2 ,2[T T -内,如图所示,令在该区间以外的0)(=t x ,则积分?∞ ∞-dt t x )(收敛,满足傅里叶变换条件,变换后用功率谱密度函数表示。 2、功率谱密度函数(又称功率谱)的物理意义 是在频域中对信号能量或功率分布情况的描述。功率谱表示振动能量在频率 域的分解,其应用十分广泛。功率谱的横坐标是频率,纵坐标是实部、虚部的模 的平方。 功率谱密度函数作为随机信号在频域内描述的函数。对于随机信号而言,它 不存在频谱函数,只存在功率谱密度函数(功率大小在频谱中反映为频谱的面 积)。时域中的相关分析为在噪声背景下提取有用信息提供了途径。功率谱分析 则从频域提供相关技术所能提供的信息,它是研究平稳随机过程的重要方法。 3.功率谱密度函数的应用 (1)结构各阶固有频率的测定 工程结构特别是大型结构(如高层楼房、桥梁、 高塔和重要机械设备等)要防止共振引起的破坏,需要测定其固有频率。如果对 结构加以激励(或以大地的脉动信号作为激励信号),即可测定结构的响应(振动信 号),再对响应信号作自功率谱分析,便可由谱图中谱峰确定结构的各阶固有频

率。 (2)利用功率谱的数学特点求取信号传递系统的频率响应函数。 (3)作为工业设备工作状况的分析和故障诊断的依据 根据功率谱图的变化,可以判断机器设备的运转是否正常。同时.还可根据机器设备正常工作和不正常工作时,振动加速度信号的功率谱的差别,查找不正常工作时,功率谱图中额外谱峰产生的原因以及排除故障的方法。 自功率谱密度函数 定义及其物理意义 假如)(t x 是零均值的随机过程,即0=x μ(如果原随机过程是非零均值的,可以进行适当处理使其均值为零)又假设)(t x 中没有周期分量,那么当∞→τ, 0)(→τx R 。这样,自相关函数)(τx R 可满足傅里叶变化的条件∞?∞ ∞- ττd R x )(。 )(τx R 傅里叶变换)(f S x ,ττπτd e R f S j x x 2)()(-∞∞-? =和逆变换 df e f S R j x x πττ2)()(-∞ ∞ -?=定义)(f S x 为)(t x 的自功率谱密度函数,简称自谱或自功率谱。 出此可见, )(f S x 曲线下和频率轴所包围的面积就是信号的平均功率,)(f S x 就是信号的功率密度沿频率轴的分布,故称)(f S x 为自功率谱密度函数。自功率谱密度函数和幅值谱的关系为 2)(1lim )(f X T f S x T x →=利用这一种关系,就可以通过直接对时域信号作傅里叶变换来计算功率谱。自相关分析可以有效地检测出信号中有无周期成分。自功率谱密度也能用来检测信号中的周期成分。周期信号的频谱是脉冲函数,在某特定频率上的能量是无限的。但是在实际处理时,用矩形窗函数对信号进行截断,这相当于在频域用矩形窗函数的频谱sinc 函数和周期信号的频谱——δ函数实行卷积,因此截断后的周期函数的频谱已不再是脉冲函数,原来为无限大的谱线高度变成有限长,谱

实验1 常见离散信号产生和基本运算

实验项目一 常见离散信号产生和基本运算 1.实验目的 (1)掌握MATLAB 最基本的矩阵运算语句; (2)掌握对常用离散信号的理解与运算实现。 2.实验内容 (1)熟悉MATLAB 的使用环境和方法; (2)练习使用基本的向量生成、矩阵运算、绘图等语句; 利用冒号(:)生成向量: X1=[1 2 3 4 5] X2=[1.000 1.500 2.000 2.500] X3=[5 4 3 2 1] 分别生成3*3,3*4的全0矩阵,全1矩阵和随机矩阵; 分别输入矩阵: 123 456789A = 1.0 1.1 1.22.0 2.1 2.23.0 3.1 3.2 B = 分别计算A+B ,A-B ,A+3,A-4,A*B ,A.*B ,C=inv (A ),A/B,A./B ; 分别计算sin(x1),cos(x1),exp(x1),log(x2),sqrt(x2)。 (3)生成以上五种基本离散信号函数; (4)绘出信号zn e n x =)(,当6)12/1(π j z +-=、6)12/1(π j z +=时、 121=z 、62π j z +=、6π j z =时的信号实部和虚部图; (5)绘出信号)1.0*2sin(5.1)(n n x π=的频率是多少?周期是多少?产生一个数字频率为0.9的正弦序列,并显示该信号,说明其周期并绘图。 3.实验技能要求 掌握并能灵活运用MATLAB 语句对离散时间信号进行基本建立和运算。 4.实验操作要求 在实验操作过程中要注意对基本实验仪器的保护。

5.实验场所 魂芯DSP应用实验室 6.实验课后训练 实验课后训练以实验报告为表现形式,在实验报告中要对实验过程中出现的问题进行分析和思考,对所测得的数据进行数据处理,并根据结果进行总结。

实验用MATLAB产生时域离散信号

实验1用M A T L A B产生时域离散信号 一、.实验目的: 1、了解常用时域离散信号及其特点 2、掌握用MATLAB产生时域离散信号的方法 二、实验内容及步骤 1、阅读并上机验证实验原理部分的例题程序,理解每一条语句的含义。 改变例题中的有关参数(如信号的频率、周期、幅度、显示时间的取值范围、采样点数等),观察对信号波形的影响。 2、编写程序,产生以下离散序列: n1=-3;n2=4;n0=0; n=n1:n2; x=[n==n0]; stem(n,x,'filled'); axis([n1,n2,0,*max(x)]); xlabel('时间(n)');ylabel('幅度x(n)'); title('单位脉冲序列'); (2)n1=-5;n2=5;n0=0; n=n1:n2; x=[n>=n0]; stem(n,x,'filled') axis([n1,n2,0,*max(x)]); xlabel('时间(n)');ylabel('幅度x(n)'); title('单位阶跃序列'); n1=20;a=;w=*pi; n=0:n1; x=exp((a+j*w)*n); subplot(2,2,1);plot(n,real(x)); title('复指数信号的实部'); subplot(2,2,3);stem(n,real(x),'filled'); title('复指数序列的实部'); subplot(2,2,2);plot(n,imag(x)); title('复指数信号的虚部'); subplot(2,2,4);stem(n,imag(x),'filled'); title('复指数序列的虚部');

用MATLAB实现常用的离散时间信号及其时域运算

用M A T L A B实现常用的离散时间信号及其时域 运算 WTD standardization office【WTD 5AB- WTDK 08- WTD 2C】

实验四用M A T L A B实现常用的离散时间信号及其时域运算 —— 摘要:在MATLAB中,只能用向量来表示离散时间信号。与连续信号不同,离散时间信号无法用符号运算来表示。用适当的MATLAB语句表示出信号后,就可以利用MATLAB的绘图命令stem来绘出直观的信号波形图,stem是专门用于绘制离散时间信号的。在MATLAB中离散序列的时域运算和变换不能用符号运算来实现,而必须用向量表示的方法,即在MATLAB中离散序列的相加、相乘需表示成两个向量的相加、相乘,因而参加运算的两序列向量必须有相同的维数。 一、实验目的:(1)学习MATLAB语言及其常用指令; (2)学习和掌握用MATLAB语言产生离散时间信号的编程方法; (3)通过编程绘制出离散时间信号的波形,加深理解信号的时域运算。二、实验内容:(1)运用MATLAB的绘图指令绘制离散时间信号; (2)用MATLAB语言实现离散时间信号的时域运算。 三、实验原理:(1)单位阶跃序列和单位样值序列。 离散时间信号只在某些离散的瞬时给出信号的值,因此,它是时间上不连续的序列。单位阶跃序列和单位样值序列在离散时间信号与系统的分析中是两个非常典型的序列,分别记为u(n)和δ(n)。它们的定义分别如下: 1 n≥0 1 n≥0 u(n)= δ(n)= 0 n<0 0 n≠0 若单位阶跃序列的起始点为n0,单位样值序列出现在n0时刻,则表达式分别为:

DSP实验报告――离散信号的产生及运算 精

离散信号的产生及运算报告 一、实验目的: 1、复习和巩固数字信号处理中离散信号的产生和运算 2、学习和掌握用MATLAB产生离散信号的方法 3、学习和掌握用MATLAB对离散信号进行运算 二、实验原理: 1.用MATLAB函数产生离散信号 信号是数字信号处理的最基本内容。没有信号,数字信号处理就没了工作对象。MATLAB7.0内部提供了大量的函数,用来产生常用的信号波形。例如,三角函数(sin,cos),指数函数(exp),锯齿波函数(sawtooth), 随机数函数(rand)等。 1 产生被噪声污染的正弦信号 用随机数函数产生污染的正弦信号。 2 产生单位脉冲序列和单位阶跃序列 按定义,单位脉冲序列为

单位阶跃序列为 。 3 矩形脉冲信号:在MATLAB 中用rectpuls 函数来表示,其调用形式为: y=rectpuls(t,width,用以产生一个幅值为1,宽度为width,相对于t=0 点左右对称的矩形波信号,该函数的横坐标范围由向量t 决定,是以t=0 为中心向左右各展开width/2 的范围,width 的默认值为1。 例:以t=2T(即t-2×T=0为对称中心的矩形脉冲信号的MATLAB 源程序如下:(取T=1) t=0:0.001:4; T=1; ft=rectpuls(t-2*T,2*T; plot(t,ft;grid on; axis([0 4 –0.5 1.5]; 4 周期性矩形波(方波)信号在MATLAB 中用square 函数来表示,其调用形 式为:y=square(t,DUTY,用以产生一个周期为2π、幅值为±1 的周期性方

《离散时间信号的表示及运算》

实验一 离散时间信号的表示及运算 一、实验目的 1.掌握离散时间信号的时域表示; 2.掌握离散时间信号的基本运算; 3.用MA TLAB 表示的常用离散时间信号及其运算; 4.掌握用MA TLAB 描绘二维图形的方法。 二、实验原理 离散时间信号是指在离散时刻才有定义的信号,简称离散信号,或者序列。离散序列通常用)(n x 来表示,自变量必须是整数。离散时间信号的波形绘制在MATLAB 中一般用stem 函数。 对离散时间序列实行基本运算可得到新的序列,这些基本运算主要包括加、减、乘、除、移位、反折等。两个序列的加减乘除是对应离散样点值的加减乘除,因此,可通过MATLAB 的点乘和点除、序列移位和反折来实现。 一些常用序列 1.单位冲激序列(单位抽样))(n δ ?? ?≠==0,00,1)(n n n δ (1) 2.单位阶跃序列)(n u ???=,,01)(n u 00<≥n n (2) 3.矩形序列)(n R N ???=,,01)(n R N 其他10-≤≤N n (3) 4.正弦序列和指数序列 正弦序列 )s i n ()(0?ω+=n A n x (4) 式中:A 为幅度,0ω为数字域的频率,它反映了序列变化的速率,?为起始相位。 实指数序列 )()(n u a n x n = (5)

式中,a 为实数。当1a 时,序列是发散的。a 为负数时,序列是摆动的。 复指数序列 n j e n x )(0)(ωσ+= (6) 它具有实部和虚部,0ω是复正弦的数字域频率。 三、实验内容 1.用Matlab 编制程序分别产生单位抽样序列)(n δ、单位阶跃序列)(n u 、矩形序列)()(5n R n x =、正弦序列)8 sin(2)(n n x π=、复指数序列n j e n x )641()(π+=,并画波形图; 绘制)(n δ波形 绘制n j e n x ][)()2.01.0(π+-=的实部和虚部的波形。

信号、系统及系统响应,离散系统的时域分析实验报告

实验报告 实验二 信号、系统及系统响应,离散系统的时域分析 一、实验目的 (1) 熟悉连续信号经理想采样前后的频谱变换关系,加深对时域采样定理的理 解; (2) 熟悉时域离散系统的时域特性; (3) 利用卷积方法观察分析系统的时域特性; (4) 掌握序列傅里叶变换的计算机实现方法,利用序列的傅里叶变换对连续信 号、离散信号及系统响应进行频域分析。 (5) 熟悉并掌握离散系统的差分方程表示法; (6) 加深对冲激响应和卷积分析方法的理解。 二、实验原理与方法 1、信号、系统及系统响应 采样是连续信号数字处理的第一个关键环节。对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生的变化以及信号信息不丢失的条件,而且可以加深对傅里叶变换、Z 变换和序列傅里叶变换之间关系式的理解。 我们知道,对一个连续信号xa(t)进行理想采样的过程可用(2-1)表示。 ^ ()()() (21) a a x t x t p t =- 其中^ ()a x t 为()a x t 的理想采样,()p t 为周期冲激脉冲,即 ()() (22) n p t t nT δ∞ =-∞= --∑ ^ ()a x t 的傅里叶变换^ ()a X j Ω为 ^ 1()[()] (23) a a s m X j X j m T ∞ =-∞ Ω=Ω-Ω-∑ (2-3)式表明^ ()a X j Ω为()a X j Ω的周期延拓,其延拓周期为采样角频率

(2/)s T πΩ=。其采样前后信号的频谱只有满足采样定理时,才不会发生频率混叠失真。 将(2-2)带入(2-1)式并进行傅里叶变换: ^ ()[()()]j t a a n X j x t t nT e dt δ∞ ∞ -Ω-∞ =-∞ Ω=-∑? [()()]j t a n x t t nT e dt δ∞ ∞ -Ω-∞ =-∞ = -∑? ()(24) j nT a n x nT e ∞ -Ω=-∞ = -∑ 式中()a x nT 就是采样后得到的序列()x n ,即 ()()a x n x nT = ()x n 的傅里叶变换()j X e ω为 ()()(25) j j n n X e x n e ω ω∞ -=-∞ = -∑ 比较(2-5)和(2-4)可知 在数字计算机上观察分析各种序列的频域特性, 通常对X(ej ω)在[0, 2π]上进行M 点采样来观察分析。 对长度为N 的有限长序列x(n), 有 一个时域离散线性非移变系统的输入/输出关系为 上述卷积运算也可以在频域实现 2、离散系统时域分析 ^ ()() (26) j a T X j X e ωω=ΩΩ=-1 ()()(27) 2,0,1,,1k N j n j k n k X e x m e k k M M ωωπ ω--==-= =???-∑()()()()() (28) m y n x n h n x m h n m ∞ =-∞ =*= --∑()()() (29) j j j Y e X e H e ωωω=-式中

实验一 :熟悉MATLAB的随机信号处理工具箱

《随机信号分析与处理》 实验报告 指导教师:廖红华 班级:03104 学号:03104xx 姓名 : 2012-12-30

实验一:熟悉MATLAB的随机信号处理工具箱 一、实验目的 1、熟悉GUI格式的编程及使用。 2、掌握随机信号的简单分析方法 3、熟悉语音信号的播放、波形显示、均值等的分析方法及其编程 二、实验原理及分析 1、语音的录入与打开 在MATLAB中,[y,Fs,bits]=wavread('Blip',[N1 N2]);用于读取语音,采样值放在向量y中,Fs表示采样频率(Hz),bits表示采样位数。[N1 N2]表示读取从N1点到N2点的值。 波形分析: 此图为一段语音信号的501个点的采样图,横轴为点数,纵轴为语音的强度,如果要加强信号,可采用尺寸扩大,此处主要学习了语音信号的调用。 2、信号加噪 信号加噪采用randn()函数产生噪声,然后加载在原始语音信号上, 原始信号+加噪信号图: 原始信号加倍后+噪声信号: 波形分析:

此图红色曲线为原始语音信号,绿色部分为噪声信号,蓝色部分为加噪后的信号,图1由于原始信号太小而与加噪后的信号差异不明显,于是有图2放大后的信号,在这个模块实验过程中两路信号相加时,维数要相同。 3、时域信号的FFT分析 FFT即为快速傅里叶变换,是离散傅里叶变换的快速算法,它是根据离散傅里叶变换的奇、偶、虚、实等特性,对离散傅里叶变换的算法进行改进获得的。在MATLAB的信号处理工具箱中函数FFT的一种调用格式为 其中X是序列,Y是序列的FFT。 波形分析: 大致一看,此图呈现左右对称的,这是因为FFT 是Z 变换和离散序列傅立叶变换上的单位圆上等间隔取点,而傅立叶和Z变换均包含周期为2pi的特性。那么在单位圆上取点,根据三角函数的特性他们相位相差一百八十度只需要在前面加一个负号(sinx)或者直接不用加(cosx),而我们得到的FFT是幅频特性曲线,高低只代表幅度大小,重点来了:我们在单位原上取的点是一个复数(s域或者z域),复数的大小是实部的平方加虚部的平方再开根号,根据刚刚我们推得的三角函数特性,如果相位差180度,也就是一个pi,他们之间的幅度应该是完全一样的! 4、均值 随机变量X 的均值也称为数学期望,它定义为 对于离散型随机变量,假定随机变量X 有N 个可能取值,各个取值的概率为则均值定义为 上式表明,离散型随机变量的均值等于随机变量的取值乘以取值的概率之和,如果取值是等概率的,那么均值就是取值的算术平均值,如果取值不是等概率的,那么均值就是概率加权和,所以,均值也称为统计平均值。

实验一--常见离散信号的MATLAB产生和图形显示

实验一 常见离散信号的MATLAB 产生和图形显示 授课课时:2学时 一、实验目的: (1)熟悉MATLAB 应用环境,常用窗口的功能和使用方法。 (2)掌握MATLAB 在时域内产生常用离散时间信号的方法。 (3)掌握离散信号的基本运算。 (4)掌握简单的绘图命令。 二、实验原理: (一)信号的表示和产生 ① 单位抽样序列 ? ??=01)(n δ 00 ≠=n n 如果)(n δ在时间轴上延迟了k 个单位,得到)(k n -δ即: ???=-0 1)(k n δ ≠=n k n 参考程序: 例1-1:)2010(()(<<-=n n n x )δ clear all n1=-10;n2=20;n0=0;%在起点为n1,终点为n2的范围内,于n0处产生冲激。 n=n1:n2;%生成离散信号的时间序列 x=[n==n0];%生成离散信号x(n) stem(n,x);%绘制脉冲杆图

xlabel(' n');ylabel('x(n)');%横坐标和纵坐标的标注说明。 title('Unit Sample Sequence');%图形上方标注图名 axis([-10 20 0 1.2]);%确定横坐标和纵坐标的取值范围 ② 单位阶跃序列 ???=0 1 )(n u 00<≥n n 例1-2:)202((u )(<<-=n n n x ) clear all n1=-2;n2=20;n0=0; n=n1:n2;%生成离散信号的时间序列 x=[n>=n0];%生成离散信号x(n) stem(n,x,'filled'); xlabel('n');ylabel('x(n)'); title('Unit step Sequence'); axis([-2 20 0 1.2]); ③ 正弦序列 )sin()(?+=wn A n x 例1-3:一正弦信号的频率为1HZ ,振幅值幅度A 为1V ,在窗口显示2个周期的信号波形,并对该信号的一个周期进行32点采样获得离散信号并显示该连续信号和离散信号的波形。 参考程序: clear f=1; A=1;nt=2;

离散时间信号的表示及运算

第2章 离散时间信号的表示及运算 2.1 实验目的 ● 学会运用MATLAB 表示的常用离散时间信号; ● 学会运用MATLAB 实现离散时间信号的基本运算。 2.2 实验原理及实例分析 2.2.1 离散时间信号在MATLAB 中的表示 离散时间信号是指在离散时刻才有定义的信号,简称离散信号,或者序列。离散序列通常用)(n x 来表示,自变量必须是整数。 离散时间信号的波形绘制在MATLAB 中一般用stem 函数。stem 函数的基本用法和plot 函数一样,它绘制的波形图的每个样本点上有一个小圆圈,默认是空心的。如果要实心,需使用参数“fill ”、“filled ”,或者参数“.”。由于MATLAB 中矩阵元素的个数有限,所以MATLAB 只能表示一定时间范围内有限长度的序列;而对于无限序列,也只能在一定时间范围内表示出来。类似于连续时间信号,离散时间信号也有一些典型的离散时间信号。 1. 单位取样序列 单位取样序列)(n δ,也称为单位冲激序列,定义为 )0()0(01)(≠=???=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])

时域离散信号的产生与运算

典型时域离散序列的产生与简单运算 1. 单位冲激序列 程序1: function [x,n]=impseq(n0,n1,n2) % generates x(n)=delta(n-n0); n1<=n<=n2 n=[n1:n2]; x=[(n-n0)==0]; 调用:[x,n]=impseq(0,-3,4); stem(n,x) 程序2: n1=-3;n2=4;n0=0; n=n1:n2; x=[n==n0]; stem(n,x,'filled'); axis([n1,n2,0,1.1*max(x)]); xlabel('时间(n)');ylabel('幅度x(n)'); title('单位脉冲序列'); 2. 单位阶跃序列 程序: n1=-3;n2=4;n0=0; n=n1:n2; x=[n>=n0]; stem(n,x,'filled'); axis([n1,n2,0,1.1*max(x)]); xlabel('时间(n)');ylabel('幅度x(n)'); title('单位阶跃序列'); 3. 矩形序列 程序: 10()00n n n δ=?=?≠?1≥0()00n u n n ?=?

n=[-10:10]; xn1=[(n-0)>=0]; xn2=[(n-4)>=0]; %定义两个阶跃序列; xn=xn1-xn2; 两个阶跃序列之差得到矩形序列; stem(n,xn,'.'); xlabel('时间(n)');ylabel('幅度x(n)'); title(‘矩形序列'); 4. 正弦序列 程序: n=0:20; xn=sin(pi/4*n); stem(n,xn,'.'); xlabel('时间(n)');ylabel('幅度x(n)'); title(‘正弦序列'); 5. 指数序列 程序: n=[0:20]; x=(0.9).^n; stem(n,x); xlabel('时间(n)');ylabel('幅度x(n)'); title(‘指数序列'); 6. 对conv 进行简单的扩展conv_m ,可以完成任意位置序列的卷积. 对于有限长序列x (n ),h (n ),它们分别的区域为[n xb,n xe]和[n hb,n he],则卷积后的区域为 [n xb+n hb,n xe+n he] 程序: function[y,ny]=conv_m(x,nx,h,nh) nyb=nx(1)+nh(1); nye=nx(length(x))+nh(length(h)); ny=[nyb:nye]; y=conv(x,h); 调用: x=[3,11,7,0,-1,4,2]; h=[2,3,0,-5,2,1]; nx=[-3:3]; nh=[-1:4]; [y,ny]=conv_m(x,nx,h,nh) ()sin()x n A n ωθ=+n a n x =)(

相关主题
文本预览
相关文档 最新文档