连续时间信号的时域分析
一、实验目的
1、 掌握连续时间信号时域运算的基本方法；
2、 掌握相关格式的调用格式及作用；
3、 掌握连续信号的基本运算；
4、 掌握利用计算机进行卷积的运算的原理和方法；
5、 熟悉连续信号卷积运算函数conv 的应用； 二、实验原理
信号的基本运算包括信号的相加(减)和相乘（除。信号的时域变换包括信号的平移、翻转、倒相尺度变换等，由以下公式所描述：
1、 相加（减）：12(t)f (t)f (t)f =±
2、 乘：12f(t)f (t)f (t)=?
3、 延时或平移：0f(t)f(t t )→-，0t 0>时右移，0t 0<时左移
4、 翻转：→f(t)f(-t)
5、 尺度变换：()()f t f at →，1a >时尺度缩小，1a <时尺度放大，0a <时还必
须包含翻转
6、 标量相乘：()()f t af t →
7、 倒相：()()f t f t →-
8、 微分：()
()df t f t dt
→ 9、 积分：()()t
f t f d ττ-∞
→?
10、
卷积：12()()*()f t f t f t =
三、验证性实验
1、 连续信号的相加 >> clear all;
>> t=0:0.0001:3; >> b=3;
>> t0=1;u=stepfun(t,t0); >> n=length(t); >>fori=1:n
u(i)=b*u(i)*(t(i)-t0); end
>> y=sin(2*pi*t); >> f=y+u; >>plot(t,f);
>>xlabel('时间(t)');ylabel('幅值f(t)');title('连续信号的相加');
2、连续信号的相乘
>> clear all;
>>t=0:0.0001:5;
>>b=3;
>>t0=1;u=stepfun(t,t0);
>>n=length(t);
>>for i=1:n
>>u(i)=b*u(i)*(t(i)-t0);
>>end
>>y=sin(2*pi*t);
>> f=y.*u;
>>plot(t,f)
>>xlabel(‘时间(t)’);ylabel(‘幅值f(t)’);title(‘连续信号的相乘’);
3、移位
>> clear all;
>> t=0:0.0001:2;
>> y=sin(2*pi*t);
>> y1=sin(2*pi*(t-0.2));
>>plot(t,y,'-',t,y1,'--')
4、尺度变换
>> clear all;
>>t=0:0.0001:1;
>>a=2;
>>y=sin(2*pi*t);
>>y1=sin(2*a*pi*t);
>>subplot(2,1,1);
>>plot(t,y);
>>ylabel('y(t)');xlabel('t');
>> title('尺度变换');
>>subplot(2,1,2)
>>plot(t,y1);
>>ylabel('y1(t)');xlabel('t');
四、设计性实验
1、 已知信号1f (t)(t 4)[U(t)U(t 4)]=-+--，2(t)sin(2t)f π=，用MATLAB 绘出下列
信号的时域波形。 (a)311(t)f (t)f (t)f =-+ (b)411f (t)[f (t)f (t)]=--+ (c)523f (t)f (t)f (t)=? (d)612f (t)f (t)f (t)=?
程序如下：
>>clear >>syms t
>>f1=sym('(-t+4)*(heaviside(t)-heaviside(t-4))');
>>f2=sin(2*pi*t);
>>f3=f1+subs(f1,t,-t);
>>f4=-f3;
>>f5=f2*f3;
>>f6=f1*f2;
>>ezplot(f1,[-5,5]);title('f1(t)');
>>ezplot(f2,[-5,5]);title('f2(t)');
>>ezplot(f3,[-5,5]);title('f3(t)');
>>ezplot(f4,[-5,5]);title('f4(t)');
>>ezplot(f5,[-5,5]);title('f5(t)');
>>ezplot(f6,[-5,5]);title('f6(t)'); 图像依次如下：
2、 已知信号f(t)的波形如图1-2-2所示。试画出下列各函数对时间t 的波形。
(a) f(-t) (b) f(-t+2) (c) t(-t-2) (d) f(2t) (e) f(t/2) (f) f(t-2) (g) f(-t/2+1)
(h) 1
[(1)]2
d f t dt + (i)
(2)t
f d ττ-∞
-?
程序如下：
>>clear >>syms t
>>f1=sym('(-t/2+1)*(heaviside(t+2)-heaviside(t-2))'); >>f2=subs(f1,t,-t); >>f3=subs(f1,t,-t+2); >>f4=subs(f1,t,-t-2); >>f5=subs(f1,t,2*t); >>f6=subs(f1,t,t/2); >>f7=subs(f1,t,t-2); >>f8=subs(f1,t,-t/2+1);
>>ezplot(f1,[-3,3]);title('f(t)'); >>ezplot(f1,[-3,3]);title('f1(t)'); >>ezplot(f2,[-3,3]);title('f1(t)'); >>ezplot(f3,[-3,7]);title('f2(t)'); >>ezplot(f4,[-6,2]);title('f3(t)'); >>ezplot(f5,[-3,3]);title('f4(t)'); >>ezplot(f6,[-5,5]);title('f5(t)'); >>ezplot(f7,[-2,7]);title('f6(t)'); >>ezplot(f8,[-3,7]);title('f7(t)'); >> f9=diff(subs(f1,t,t/2+1)); >>ezplot(f9,[-5,5]);title('f8(t)'); >> f10=int(f2);
>>ezplot(f10,[-5,5]);title('f9(t)');
图像依次如下：
图1-2-2
测试技术与信号处理实验报告 机械转子底座的振动测量和分析 一、实验目的 1．掌握磁电式速度传感器的工作原理、特点和应用。
2．掌握振动的测量和数据分析。 二、实验内容和要求 先利用光电式转速传感器测量出电机的转速；然后利用磁电式速度传感器测量机械转子底座在该电机转速下的振动速度；对测量出的振动速度信号进行频谱分析；找出振动信号的主频与电机转速之间的关系。 三、实验步骤 1．启动实验程序“机械转子系统的振动测量.exe”; 输入个人信息，也可以启动之后通过单击“修改”按钮修改个人信息。 2．单击“采样设置”按钮，输入采集卡连接磁电速度传感器的采样通道号，批量采样频率（建议设为10KHz）、批量采样点数（建议设为10000）。 3．打开转子电机的电源，单击“单点采样”。 4．旋转调节旋钮改变转子的转速,观察图形区显示的磁电速度传感器采集到的转子底座振动信号；如果振动信号比较小，可适当提高转子的转速。 5．转子转速的测量： (1) 单击“采样设置”按钮，输入采集卡连接光电转速传感器的 采样通道号、批量采样频率（建议值为10KHz）、批量采样点 数（建议值为10000）。 (2) 单击“批量采样”按钮，开始采样；采样完成之后，采集到 的波形信号会显示在图形窗口，系统会自动计算出转子的速度
并显示出来。记录下此时的转子的转速（单位：r/s）。 (3) 再重复步骤(2)测量2次。以三次测量的平均值作为此时转子 的转速。 转速的测量结果 单点采样采集通道6,测量3组数据 6．振动信号的测量和频谱分析： (1) 单击“采样设置”按钮，输入采集卡连接磁电速度传感器的 采样通道号、批量采样频率（建议设为10KHz）、批量采样点 数（建议设为10000）。 (2) 单击“批量采样”按钮，开始采样；采样完成之后，采集到 的波形信号会显示在图形窗口。如果信号不正常，重复点击“批 量采样”按钮 (3) 单击“保存”按钮，将采集到的磁电传感器的信号数据保存 为文本文件。文件必须保存到“C:\ExperiData\”目录下。可单 击“保存设置”更改文件名。 (4) 打开刚保存的文本文件，文件前面几行保存了个人信息、采 样频率、采样通道、保存的数据个数等信息。文件中共有四列 数据，第一列为数据的序号，第二列为磁电传感器检测到的数 据。
北京交通大学信号与系统时域分析
【研讨题目2】 信号与系统时域分析专题研讨 【目的】 1.研究用离散方法近似计算连续信号的卷积积分； 2.通过分析近似计算卷积积分过程中出现的问题，锻炼学生分析问题和解决问题的能力； 【知识点】 信号时域分析，卷积积分，卷积和 【研讨题目】连续信号卷积积分的数值近似计算 两个连续信号的卷积积分定义为 τττd )()()(-= ?∞ ∞ -t h x t y 为了能用数值方法进行计算，需对连续信号进行抽样。记x [k ]=x (k ), h [k ]=h (k ),为 进行数值计算所选定的抽样间隔，可以证明连续信号卷积积分可近似的表示为 (Δ)Δ([][])y k x k h k ≈?* (1) 由式(1)可知，可以利用Matlab 提供的conv 函数近似计算连续信号的卷积积分。 一、(*)理论分析 为了对近似计算的结果进行分析，用解析的方法计算下列卷积积分，推出卷积积分的解析表达式； (1) 时限信号卷积积分 x 1(t )=u (t )-u (t -1)，y 1(t )=x 1(t )*x 1(t ); 卷积结果为：y1(t)= x 1(t )*x 1(t )=r(t)-2*r(t-1)+r(t-2) (2) 分段常数信号卷积积分 x 2(t )= x 1(t )+2 x 1(t -1)+ x 1(t -2)，h 2(t )= x 1(t )- x 1(t -1)， y 2(t )=x 2(t )*h 2(t ); 卷积结果为：y2(t)= x 2(t )*h 2(t ) =y1(t)+y1(t-1)-y1(t-2)-y1(t-3) =r(t)-r(t-1)-2*r(t-2)+2*r(t-3)+r(t-4)-r(t-5) (3) 非时限信号卷积积分 x 3(t )=u (t )，h 3(t )=e -t u (t ), y 3(t )=x 3(t )*h 3(t ) 卷积结果为：y3= x 3(t )*h 3(t ) =[1-exp(-t)]*u(t) 二、(*)时限信号卷积积分的近似计算 取不同的△值，用Matlab 函数conv 近似计算卷积积分y 1(t )并画出其波形，讨论的取值对计算结果的影响。
实验1 常用连续时间信号的实现 一、实验目的 （1）了解连续信号的特点； （2）熟悉MATLAB的基本使用方法； （3）典型信号的MATLAB表示方法； （4）熟悉MATLAB Plot函数等应用。 二、实验原理 信号的定义 信号是随时间变化的物理量。信息的本质是时间的函数。 信号的描述 1、时域法 时域法是将信号表示成时间的函数f(t)对信号进行描述的方 法。 2、频域法（变换域） 频域法是通过正交变换，将信号表示成其它变量的函数来对信 号进行描述的方法。 信号的分类 确定性信号、连续信号、周期信号、能量信号、奇、偶信号。
三、MATLAB基本使用方法 1、MATLAB运行平台如图所示：在“>>”提示符后输入命令 并回车，执行结果会显示在屏幕上。 例如：输入y=3+5回车，显示结果如图所示。 另外在输入命令时，有时我们希望有些中间过程的结
果不显示在屏幕上，而只显示最后的结果，这时我们需要在不显示执行结果的命令后加上“；”，该命令执行的结果不回显在屏幕上。例如：计算z=x+y，其中x=2+1，y=3+5，我只想看到z的值，输入命令格式如图所示，我们可以看到x，y的结果没有在屏幕上显示，只显示z的值。 2、MATLAB中矩阵的输入方法： 矩阵的输入方法有两种，第一种方式如图所示； 第二种方式如图所示
3、M文件的使用： 在处理一些包含多条命令的问题时，如果在MATLAB的命令窗口中进行处理，当出现错误时不好修改，这时我们需要借助MATLAB提供的M文件方式来处理。M文件类似于批处理件，单击MATLAB菜单中的“file”选项，从下拉列表中选取“new” 选项，从其下拉列表中选取“m_file”，即可打开M文件的编辑窗口如图所示： 另外，M文件还可以编写函数。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 连续时间信号与系统的基本表示与 分析 实验目的： 1. 了解并掌握用Matlab 分析信号与系统时用到的主要基本知识； 2. 掌握用matlab 表示常见的连续信号，并进行信号的基本运算 3. 熟悉使用Matlab 分析系统的方法。 实验内容： 1. 画出下列连续信号的波形 （1）()() 2()t x t e u t -=- （2）[](1) ()(1)(2)t x t e u t u t --=--- （3）()cos(10)()t x t e t u t π-= 2. 判断系统是否为线性系统，是否具有时不变特性。 （1）()(2)y t x t = （2）()2 ()y t x t = 指导资料： 1.1MATLAB 的基础知识 1.1.1数据的表示和运算 MATLAB 中，变量可以通过变大时直接赋值，例如直接输入 >>a=2+2 得到的结果为 a=4 如果输入的表达式后面加上分号“;”，那么结果就不会显示出来。由于MA TLAB 的变量名对字母大小写敏感，因此“a ”和“A ”是两个不同的变量名。 MATLAB 主要用到以下数值运算符： + 加 - 减 * 乘 / 除 ^ 乘方(幂)
' （矩阵）转置 这些符号可以对数值或已经定义过的变量进行运算，并给变量直接赋值。例如，假设变量“a”在上面已经定义过，则 >>b=2*a 得到的结果为 b=8 MATLAB中有一些预定义的变量可以直接使用。信号与系统中常用的变量有： i或j 1 pi π(3.1415926……) 在信号与系统中，常用以下函数进行计算和对变量的赋值： abs 数值的大小（实数的绝对值） angel 复数的角度，以弧度表示 real 求复数的实部 imag 求复数的虚部 cos 余弦函数，假设角度是弧度值 sin 正弦函数，假设角度是弧度制 exp 指数海曙 sqrt 求平方根 例如： >>y=2*(1+4*i) y=2.000+8.000i >>c=abs(y) c=8.2462 >>d=angle(y) d=1.3258 1.1.2矩阵的表示和运算 MATLAB是基于矩阵和向量的代数运算，甚至标量也可以看做是1×1 的矩阵，因此，MATLAB中对矩阵和向量的操作比较简单。 向量可以用两种方法定义。第一种是指定元素建立向量： v = [1 3 5 7]； 这个命令创建了一个1×4 的行向量，元素为1,3,5和7.可以用逗号代替空格来分隔元素： v = [1,3,5,7]； 如果要增加向量的元素，可以表示为 v(5) = 8 得到的向量为v = [1 3 5 7 8]。前面定义过的向量还可以用来定义新的向量，例如，前面已经定义过的向量v ，再定义向量a和b： a = [9 10]; b = [ v a]; 得到向量b为
测控1005班齐伟0121004931725 （18号）实验一差分方程、卷积、z变换 一、实验目的 通过该实验熟悉 matlab软件的基本操作指令，掌握matlab软件的使用方法，掌握数字信号处理中的基本原理、方法以及matlab函数的调用。 二、实验设备 1、微型计算机1台； 2、matlab软件1套 三、实验原理 Matlab 软件是由mathworks公司于1984年推出的一套科学计算软件，分为总包和若干个工具箱，其中包含用于信号分析与处理的sptool工具箱和用于滤波器设计的fdatool工具箱。它具有强大的矩阵计算和数据可视化能力，是广泛应用于信号分析与处理中的功能强大且使用简单方便的成熟软件。Matlab软件中已有大量的关于数字信号处理的运算函数可供调用，本实验主要是针对数字信号处理中的差分方程、卷积、z变换等基本运算的matlab函数的熟悉和应用。 差分方程（difference equation)可用来描述线性时不变、因果数字滤波器。用x表示滤波器的输入，用y表示滤波器的输出。 a0y[n]+a1y[n-1]+…+a N y[n-N]=b0x[n]+b1x[n-1]+…+b M x[n-M] (1) ak,bk 为权系数，称为滤波器系数。 N为所需过去输出的个数，M 为所需输入的个数卷积是滤波器另一种实现方法。 y[n]= ∑x[k] h[n-k] = x[n]*h[n] (2) 等式定义了数字卷积，*是卷积运算符。输出y[n] 取决于输入x[n] 和系统的脉冲响应h[n]。 传输函数H(z)是滤波器的第三种实现方法。 H(z)=输出/输入= Y(z)/X(z) （3）即分别对滤波器的输入和输出信号求z变换，二者的比值就是数字滤波器的传输函数。 序列x[n]的z变换定义为 X (z)=∑x[n]z-n (4) 把序列x[n] 的z 变换记为Z{x[n]} = X(z)。
实验一信号与系统的时域分析 一、实验目的 1、熟悉与掌握常用的用于信号与系统时域仿真分析的MA TLAB函数; 2、掌握连续时间与离散时间信号的MA TLAB产生,掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MA TLAB编程; 3、牢固掌握系统的单位冲激响应的概念,掌握LTI系统的卷积表达式及其物理意义,掌握卷积的计算方法、卷积的基本性质; 4、掌握利用MA TLAB计算卷积的编程方法,并利用所编写的MA TLAB程序验证卷积的常用基本性质; 掌握MA TLAB描述LTI系统的常用方法及有关函数,并学会利用MATLAB求解LTI系统响应,绘制相应曲线。 基本要求:掌握用MA TLAB描述连续时间信号与离散时间信号的方法,能够编写MATLAB程序,实现各种信号的时域变换与运算,并且以图形的方式再现各种信号的波形。掌握线性时不变连续系统的时域数学模型用MA TLAB描述的方法,掌握卷积运算、线性常系数微分方程的求解编程。 二、实验原理 信号(Signal)一般都就是随某一个或某几个独立变量的变化而变化的,例如,温度、压力、 声音,还有股票市场的日收盘指数等,这些信号都就是随时间的变化而变化的,还有一些信号,例如在研究地球结构时,地下某处的密度就就是随着海拔高度的变化而变化的。一幅图片中的每一个象素点的位置取决于两个坐标轴,即横轴与纵轴,因此,图像信号具有两个或两个以上的独立变量。 在《信号与系统》课程中,我们只关注这种只有一个独立变量(Independent variable)的信号,并且把这个独立变量统称为时间变量(Time variable),不管这个独立变量就是否就是时间变量。 在自然界中,大多数信号的时间变量都就是连续变化的,因此这种信号被称为连续时间信号(Continuous-Time Signals)或模拟信号(Analog Signals),例如前面提到的温度、压力与声音 信号就就是连续时间信号的例子。但就是,还有一些信号的独立时间变量就是离散变化的,这种信号称为离散时间信号。前面提到的股票市场的日收盘指数,由于相邻两个交易日的日收盘指数相隔24小时,这意味着日收盘指数的时间变量就是不连续的,因此日收盘指数就是离散时间信号。 而系统则用于对信号进行运算或处理,或者从信号中提取有用的信息,或者滤出信号中某些无用的成分,如滤波,从而产生人们所希望的新的信号。系统通常就是由若干部件或单元组成的一个整体(Entity)。系统可分为很多不同的类型,例如,根据系统所处理的信号的不同,系统可分为连续时间系统(Continuous-time system)与离散时间系统(Discrete-time system),根
课程设计任务书 学生姓名:专业班级: 指导教师:工作单位: 题目: 连续时间信号傅里叶级数分析及MATLAB实现 初始条件： MATLAB 6.5 要求完成的主要任务： 深入研究连续时间信号傅里叶级数分析的理论知识，利用MA TLAB强大的图形处理功能，符号运算功能以及数值计算功能，实现连续时间周期信号频域分析的仿真波形。 1.用MATLAB实现周期信号的傅里叶级数分解与综合。 2.用MATLAB实现周期信号的单边频谱及双边频谱。 3.用MATLAB实现典型周期信号的频谱。 4.撰写《MATLAB应用实践》课程设计说明书。 时间安排： 学习MATLAB语言的概况第1天 学习MATLAB语言的基本知识第2、3天 学习MATLAB语言的应用环境，调试命令，绘图能力第4、5天 课程设计第6-9天 答辩第10天 指导教师签名：年月日 系主任（或责任教师）签名：年月日
目录 摘要................................................................................................................................................ I Abstract .......................................................................................................................................... II 绪论. (1) 1 MATLAB简介 (2) 1.1 MATLAB语言功能 (2) 1.2 MATLAB语言特点 (2) 2 傅里叶级数基本原理概要 (4) 2.1 周期信号的傅里叶分解 (4) 2.2 三角形式和指数形式傅里叶级数及各系数间的关系 (4) 2.3 周期信号的频谱 (5) 3 用MATLAB实现周期信号的傅立叶级数分解与综合 (6) 3.1 合成波形与原波形之间的关系 (6) 3.2 吉布斯现象 (6) 4 用MATLAB实现周期信号的单边频谱及双边频谱。 (8) 4.1 单边，双边（幅度，相位）频谱及其关系 (8) 4.1.1单边，双边（幅度，相位） (8) 4.1.2 单边，双边频谱关系 (9) 4.2以单边幅度频谱为例，研究脉冲宽度与频谱的关系 (10) 4.3以单边幅度频谱为例，研究脉冲周期与频谱的关系 (11) 5用MATLAB实现典型周期信号的频谱 (13) 5.1 周期方波脉冲频谱的MATLAB实现 (13) 5.2 周期三角波脉冲频谱的MATLAB 实现 (14) 6 小结及心得体会 (17) 参考文献 (18) 附录： (19)
Harbin Institute of Technology 匹配滤波器实验报告 课程名称：信号检测理论 院系：电子与信息工程学院 姓名：高亚豪 学号：14SD05003 授课教师：郑薇 哈尔滨工业大学
1. 实验目的 通过Matlab 编程实现对白噪声条件下的匹配滤波器的仿真，从而加深对匹配滤波器及其实现过程的理解。通过观察输入输出信号波形及频谱图，对匹配处理有一个更加直观的理解，同时验证匹配滤波器具有时间上的适应性。 2. 实验原理 对于一个观测信号()r t ，已知它或是干扰与噪声之和，或是单纯的干扰， 即 0()()()()a u t n t r t n t +?=?? 这里()r t ，()u t ，()n t 都是复包络，其中0a 是信号的复幅度，()u t 是确知的归一化信号的复包络，它们满足如下条件。 2|()|d 1u t t +∞ -∞=? 201||2 a E = 其中E 为信号的能量。()n t 是干扰的均值为0，方差为0N 的白噪声干扰。 使该信号通过一个线性滤波系统，有效地滤除干扰，使输出信号的信噪比在某一时刻0t 达到最大，以便判断信号的有无。该线性系统即为匹配滤波器。 以()h t 代表系统的脉冲响应，则在信号存在的条件下，滤波器的输出为 0000()()()d ()()d ()()d y t r t h a u t h n t h τττττττττ+∞+∞+∞ =-=-+-???
右边的第一项和第二项分别为滤波器输出的信号成分和噪声成分，即 00()()()d x t a u t h τττ+∞ =-? 0 ()()()d t n t h ?τττ+∞ =-? 则输出噪声成分的平均功率（统计平均）为 2 20E[|()|]=E[|()()d |]t n t h ?τττ+∞ -? **00*000200 =E[()(')]()(')d d '=2()(')(')d d ' 2|()|d n t n t h h N h h N h ττττττδττττττττ+∞+∞+∞+∞+∞ ---=?? ?? ? 而信号成分在0t 时刻的峰值功率为 22 20000|()||||()()d |x t a u t h τττ+∞ =-? 输出信号在0t 时刻的总功率为 22000E[|()|]E[|()()|]y t x t t ?=+ 22**0000002200E[|()||()|()()()()] |()|E[|()|] x t t x t t t x t x t t ????=+++=+ 上式中输出噪声成分的期望值为0，即0E[()]0t ?=，因此输出信号的功率 成分中只包含信号功率和噪声功率。 则该滤波器的输出信噪比为 222000022000|||()()d ||()|E[|()|]2|()|d a u t h x t t N h τττρ?ττ+∞ +∞-==?? 根据Schwartz 不等式有
实验一 连续时间信号分析 一、实验目的 （一）掌握使用Matlab 表示连续时间信号 1、学会运用Matlab 表示常用连续时间信号的方法 2、观察并熟悉常用信号的波形和特性 （二）掌握使用Matlab 进行连续时间信号的相关运算 1、学会运用Matlab 进行连续时间信号的时移、反褶和尺度变换 2、学会运用Matlab 进行连续时间信号微分、积分运算 3、学会运用Matlab 进行连续时间信号相加、相乘运算 4、学会运用Matlab 进行连续时间信号卷积运算 二、实验条件 一台电脑、winXP 系统、matlab7.0软件 三、实验内容 1、利用Matlab 命令画出下列连续信号的波形图。 （1）)4/3t (2cos π+ 代码： clear all;close all;clc; K=2;a=3; t=0:0.01:3; ft=K*cos(a*t+pi/4); plot(t,ft),grid on axis([-5,5,-2.2,2.2]) title('2cos(3t+4π)')
-5 -4 -3 -2 -1 1 2 3 4 5 -2-1.5-1-0.500.511.5 22cos(3t+4π) （2） )t (u )e 2(t -- -3 -2-10123 -3 -2 -1 1 2 3 指数信号与阶跃信号的乘积
代码： 函数文件： function f=uCT(t) f=(t>=0); 命令文件： clear all;close all;clc; a=-1; t=-5:0.01:5; ft=(2-exp(a*t)).*uCT(t); %y=2-exp(a*t); %plot(t,y),grid on plot(t,ft),grid on axis([-3,3,-3,3]); title('指数信号与阶跃信号的乘积') （3）)]2()(u )][t (cos 1[--+t u t π
实验一 连续时间信号 §1.1 表示信号的基本MATLAB 函数 目的 学习连续时间信号和离散时间信号在MATLAB 中的表示。 相关知识 1．离散时间信号的表示 通常，信号用一个行向量或一个列向量表示。在MATLAB 中全部向量都从1开始编号，如y(1)是向量y 的第1个元素。如果这些编号与你的应用不能对应，可以创建另外一标号向量与信号编号保持一致。 例如，为了表示离散时间信号?? ?≤≤-=n n n n x 其余 033 2][ 首先利用冒号运算符对][n x 的非零样本定义标号向量，然后再定义向量x ，表示在这些时间编号每一点的信号值 >> n=[-3:3]; >> x=2*n;
如果要在一个更宽的范围内检查信号，就需拓宽n和x。例如如要在5 -n画 ≤ 5≤ 出这个信号，可以拓宽标号向量n，然后将这些附加的元素加到向量x上，如>> n=[-5:5]; >> x=[0 0 x 0 0]; >> stem(n,x);
如果要大大扩展信号的范围，可利用zeros函数。 例如如果想要包括100 ≤ -n，这时可键 5≤ -n的范围，而向量x已扩展到5 ≤ 100≤ 入 >> n=[-100:100]; >> x=[zeros(1,95) x zeros(1,95)];
假设要定义][ ][ 1n n xδ =，]2 [ ] [2+ =n n xδ，可编程如下>> nx1=[0:10]; >> x1=[1 zeros(1,10)]; >> nx2=[-5:5]; >> x2=[zeros(1,3) 1 zeros(1,7)]; >> stem(nx1,x1); >> stem(nx2,x2);
随机信号分析 实验报告 目录 随机信号分析 (1) 实验报告 (1) 理想白噪声和带限白噪声的产生与测试 (2) 一、摘要 (2) 二、实验的背景与目的 (2) 背景： (2) 实验目的: (2) 三、实验原理 (3) 四、实验的设计与结果 (4) 实验设计: (4) 实验结果: (5) 五、实验结论 (12) 六、参考文献 (13) 七、附件 (13) 1
理想白噪声和带限白噪声的产生与测试一、摘要 本文通过利用MATLAB软件仿真来对理想白噪声和带限白噪声进行研究。理想白噪声通过低通滤波器和带通滤波器分别得到低通带限白噪声和帯通带限白噪声。在仿真的过程中我们利用MATLAB工具箱中自带的一些函数来对理想白噪声和带限白噪声的均值、均方值、方差、功率谱密度、自相关函数、频谱以及概率密度进行研究，对对它们进行比较分析并讨论其物理意义。 关键词：理想白噪声带限白噪声均值均方值方差功率谱密度自相关函数、频谱以及概率密度 二、实验的背景与目的 背景： 在词典中噪声有两种定义：定义1：干扰人们休息、学习和工作的声音，引起人的心理和生理变化。定义2：不同频率、不同强度无规则地组合在一起的声音。如电噪声、机械噪声，可引伸为任何不希望有的干扰。第一种定义是人们在日常生活中可以感知的，从感性上很容易理解。而第二种定义则相对抽象一些，大部分应用于机械工程当中。在这一学期的好几门课程中我们都从不同的方面接触到噪声，如何的利用噪声，把噪声的危害减到最小是一个很热门的话题。为了加深对噪声的认识与了解，为后面的学习与工作做准备，我们对噪声进行了一些研究与测试。 实验目的: 了解理想白噪声和带限白噪声的基本概念并能够区分它们，掌握用MATLAB 或c/c++软件仿真和分析理想白噪声和带限白噪声的方法，掌握理想白噪声和带限白噪声的性质。
连续时间信号和系统时域分析及MATLAB实 现[精编版] MATLAB课程设计任务书 姓名：王** 学号：2010******010 题目:
连续时间信号和系统时域分析及MATLAB实现 初始条件： MATLAB 7.5.0 ，Windows XP系统 实验任务： 一、用MATLAB实现常用连续时间信号的时域波形（通过改变参数，分析其时域特性）。 1、单位阶跃信号， 2、单位冲激信号， 3、正弦信号， 4、实指数信号， 5、虚指数信号， 6、复指数信号。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘， 4、微分， 5、积分 三、用MATLAB实现信号的时域变换（参数变化，分析波形变化） 1、反转， 2、使移（超时，延时）， 3、展缩， 4、倒相， 5、综合变化 四、用MATLAB实现信号简单的时域分解 1、信号的交直流分解， 2、信号的奇偶分解 五、用MATLAB实现连续时间系统的卷积积分的仿真波形 给出几个典型例子，对每个例子，要求画出对应波形。 六、用MATLAB实现连续时间系统的冲激响应、阶跃响应的仿真波形。 给出几个典型例子，四种调用格式。 七、利用MATLAB实现连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形。 目录
1.1 MATLAB设计目的 (1) 1.2 MATLAB语言特点 (1) 2常用连续时间信号的时域波形 (1) 2.1单位阶跃信号 (1) 2.2单位冲激信号 (2) 2.3正弦信号 (3) 2.4实指数信号 (4) 2.5虚指数信号 (5) 2.6复指数信号 (6) 3 连续时间信号的时域运算 (7) 3.1相加 (7) 3.2相乘 (8) 3.3数乘 (9) 3.4微分 (10) 3.5积分 (11) 4.1反转 (12) 4.2时移 (13) 4.3展缩 (14) 4.4倒相 (15) 4.5综合变化 (16) 5连续时间信号简单的时域分解 (17) 5.1信号的交直流分解 (17)
《测试信号分析与处理》 实验一差分方程、卷积、z变换 一、实验目的 通过该实验熟悉 matlab软件的基本操作指令，掌握matlab软件的使用方法，掌握数字信号处理中的基本原理、方法以及matlab函数的调用。 二、实验设备 1、微型计算机1台； 2、matlab软件1套 三、实验原理 Matlab 软件是由mathworks公司于1984年推出的一套科学计算软件，分为总包和若干个工具箱，其中包含用于信号分析与处理的sptool工具箱和用于滤波器设计的fdatool工具箱。它具有强大的矩阵计算和数据可视化能力，是广泛应用于信号分析与处理中的功能强大且使用简单方便的成熟软件。Matlab软件中已有大量的关于数字信号处理的运算函数可供调用，本实验主要是针对数字信号处理中的差分方程、卷积、z变换等基本运算的matlab函数的熟悉和应用。 差分方程（difference equation)可用来描述线性时不变、因果数字滤波器。用x表示滤波器的输入，用y表示滤波器的输出。 a0y[n]+a1y[n-1]+…+a N y[n-N]=b0x[n]+b1x[n-1]+…+b M x[n-M] (1) ak,bk 为权系数，称为滤波器系数。
N为所需过去输出的个数，M 为所需输入的个数卷积是滤波器另一种实现方法。 y[n]= ∑x[k] h[n-k] = x[n]*h[n] (2)等式定义了数字卷积，*是卷积运算符。输出y[n] 取决于输入x[n] 和系统的脉冲响应h[n]。 传输函数H(z)是滤波器的第三种实现方法。 H(z)=输出/输入= Y(z)/X(z) （3） 即分别对滤波器的输入和输出信号求z变换，二者的比值就是数字滤波器的传输函数。 序列x[n]的z变换定义为 X (z)=∑x[n]z-n (4)把序列x[n] 的z 变换记为Z{x[n]} = X(z)。 由X(z) 计算x[n] 进行z 的逆变换x[n] = Z-1{X(z)}。 Z 变换是Z-1的幂级数，只有当此级数收敛，Z 变换才有意义，而且同一个Z 变换等式，收敛域不同，可以代表不同序列的Z 变换函数。 这三种数字滤波器的表示方法之间可以进行相互转换。 四、实验步骤 1、熟悉matlab软件基本操作指令。读懂下列matlab程序指令，键入程序并 运行，观察运行结果。 Conv.m% 计算两个序列的线性卷积； %-----------------------------------------------------------------
实验二 连续时间信号的频域分析 一、实验目的 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法； 2、观察截短傅里叶级数而产生的“Gibbs 现象”，了解其特点以及产生的原因； 3、掌握连续时间傅里叶变换的分析方法及其物理意义； 4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质； 5、学习掌握利用Matlab 语言编写计算CTFS 、CTFT 和DTFT 的仿真程序，并能利用这些程序对一些典型信号进行频谱分析，验证CTFT 、DTFT 的若干重要性质。 基本要求：掌握并深刻理傅里叶变换的物理意义，掌握信号的傅里叶变换的计算方法，掌握利用Matlab 编程完成相关的傅里叶变换的计算。 二、原理说明 1、连续时间周期信号的傅里叶级数CTFS 分析 任何一个周期为T 1的正弦周期信号，只要满足狄利克利条件，就可以展开成傅里叶级数。 三角傅里叶级数为： ∑∞ =++=1000)]sin()cos([)(k k k t k b t k a a t x ωω 2.1 或： ∑∞ =++ =1 00)cos()(k k k t k c a t x ?ω 2.2 其中1 02T π ω= ，称为信号的基本频率（Fundamental frequency ），k k b a a ，和，0分别是信号)(t x 的直流分量、余弦分量幅度和正弦分量幅度，k k c ?、为合并同频率项之后各正弦谐波分量的幅度和初相位，它们都是频率0ωk 的函数，绘制出它们与0ωk 之间的图像，称为信号的频谱图（简称“频谱”），k c －0ωk 图像为幅度谱，k ?－0ωk 图像为相位谱。 三角形式傅里叶级数表明，如果一个周期信号x(t)，满足狄里克利条件，就可以被看作是由很多不同频率的互为谐波关系（harmonically related ）的正弦信号所组成，其中每一个不同频率的正弦信号称为正弦谐波分量 (Sinusoid component)，其幅度（amplitude ）为k c 。
河南城建学院实验报告 课程名称：信号处理基础 系：电气与电子工程系 专业：自动化 指导教师：梁成武 姓名：肖邓 学号： 122408155 报告上交时间： 2010年 12月 22日 教师评语： 成绩等级 日期：2010年12月 22日第一次实验内容：MATLAB工具使用和信号的时域分析 一、实验目的 1、掌握MATLAB基本知识和基本应用。 2、掌握用MATLAB进行仿真的简单操作。 3、通过实验体会MATLAB的操作，获取一些初步的经验。 二、实验设备： 装有MATLAB6.5版软件的计算机一台 三、实验原理与步骤 1、打开信号处理实验素材文件夹，仔细阅读里边的第01章节，认识MATLAB桌面的基本结构和学习使用一些基本的命令。 2、通过MATLAB做出指数函数2e t和抽样函数() Sa t的图形并分析。 （1）指数函数y=2e t的MATLAB实验程序：
t=0:0.01:10; y=exp(1/2*t); plot(t,y) （2）抽样函数y=() Sa t的MATLAB实验程序： t=-10:0.01:10; y=sinc(t); plot(t,y) 四、实验结果与数据处理 1、指数函数2e t实验结果如下图（一）、图（二）所示： 图（一）指数函数2e t的MATLAB仿真图 2、抽样函数() Sa t的实验结果： 图（二）抽样函数() Sa t的MATLAB仿真图 五、实验体会与讨论 通过本次试验掌握了MATLAB基本知识和基本应用，掌握了用MATLAB进行仿真的简单操作和一些基本的操作指令。 MATLAB 在处理函数、信号，以及了解他们的特性方面，MATLAB是一款非常方便的软件，当然也可以通过它来做出某些复杂函数或信号的图形。 第二次实验内容：线性系统的时域特性 一、实验目的： 了解线性系统的时域特性，并会用计算机进行简单的分析。 二、实验设备： 装有MATLAB 6.5版软件的计算机一台 三、实验原理与步骤： 连续时间信号在食欲的一些基本运算-——尺度变换、平移、翻转、叠加、相乘、
福建农林大学计算机与信息学院 信息工程类 课程设计报告 课程名称：信号与系统 课程设计题目：连续时间信号的频域分析 姓名： 系：电子信息工程 专业：电子信息工程 年级：2008 学号： 指导教师： 职称： 2011 年 1 月10 日
福建农林大学计算机与信息学院信息工程类 课程设计结果评定
目录 1课程设计的目的 (1) 2课程设计的要求 (1) 3课程设计报告内容.....................................................................1-13 3.1连续信号的设计..................................................................1-11 3.2验证傅里叶变换的调制定理 (11) 3.3周期信号及其频谱 (12) 4总结 (13) 参考文献 (14)
连续时间信号的频域分析 1.课程设计的目的 （1）熟悉MATLAB语言的编程方法及MATLAB指令； （2）掌握连续时间信号的基本概念； （3）掌握门函数、指数信号和抽样信号的表达式和波形； （4）掌握连续时间信号的傅里叶变换及其性质； （5）掌握连续时间信号频谱的概念以及幅度谱、相位谱的表示； （6）掌握利用MATLAB进行信号的傅里叶变换以及时域波形和频谱的表示；（7）通过连续时间信号的频域分析，更深刻地理解了连续时间信号的时域和频域间的关系，加深了对连续时间信号的理解。 2.课程设计的要求 （1）自行设计以下连续信号：门函数、指数信号和抽样信号。要求：（a）画出以上信号的时域波形图； （b）实现以上信号的傅里叶变换，画出以上信号的幅度谱及相位谱，并对相关结果予以理论分析； （c）对其中一个信号进行时移和尺度变换，分别求变换后信号的傅里叶变换，验证傅里叶变换的时移和尺度变换性质。 （2）自行设计信号，验证傅里叶变换的调制定理。 （3）自行设计一个周期信号，绘出该信号的频谱，并观察周期信号频谱的特点。 3.课程设计报告内容 3.1（a）①门函数(矩形脉冲)： MATLAB中矩形脉冲信号用rectpuls函数表示： y=rectpuls (t,width) %width缺省值为1 >> t=-2:0.001:2; T=2; yt=rectpuls (t,T); plot(t,yt); axis([-2,2,0,1.5]); grid on; %显示格线
信号检测法——有无法的实验报告 姓名： 指导老师： 时间：
信号检测法——有无法的实验报告 /// /// /// /// (////////学院长春 130000) 摘要：运用信号检测论的“有无法”检测被试辨别力、判断标准和面对不 同先验概率下的反应偏向，通过做ROC曲线测出被试对信号和噪音的辨别能力。 关键词：信号检测论辨别力判断标准 ROC 曲线 1 引言 信号检测论是现代心理物理学的重要组成部分。它假设人们在对刺激进行感知时，干扰（也即噪音）也总是存在的，即人作为一个接受者对刺激的辨别问题可等效于一个在噪音中检测信号的问题。 在心理学领域，信号检测论所指的信号可以理解为刺激。在信号检测论中，噪音就是对信号检测起干扰作用的所有背景，对信号起着干扰作用的因素都可当做“噪音”。 信号检测论把刺激的判断看成对信号的侦察并作出决策的过程，其中既包括感觉过程也包括决策过程。感觉过程是神经系统对信号或噪音的客观反应，它仅取决于外在的刺激的性质，即信号和噪音之间的客观区别；而决策过程受到主观因素的影响。前者决定了被试的感受性大小，信号检测论多选用辨别力指标d’来作为反映客观感受性的指标；后者则决定被试的决策是偏向于严格还是偏向于宽松，信号检测论用似然比β或报告标准C来对反应倾向进行衡量。并学习绘制ROC曲线。 2 研究方法 2.1 被试 1名在校女大学生， 20岁,身体健康，///////学院11级应用心理学专业。2.2仪器及材料JGW—B心理实验台操作箱，100克，104克，108克和112克的砝码重量各一个。 2.3 方法 1．准备工作： （1）把104克、108克和112克的重量分别和100克的重量比较10次，选出一个在十次比较中7次或8次觉得比100克重的重量作为信号刺激。100克的重量作为噪音。 （2）主试按下三种不同的SN和N出现的概率安排实验顺序。 （2）（3）
信号与系统 课程实验报告 实验四 实验名称 信号的时域分析 及Matlab 实现 系 别 教师姓名 实验地点 5309 实验日期 2011-06-20 学生姓名 学号 一、实验内容 1．预习实验原理。 2．对实验内容编写程序(M 文件)，上机运行。 3．记录并整理实验数据。 二、实验目的 1．掌握用Matlab 分析系统时间响应的方法。 2．掌握用Matlab 分析连续、离散系统的冲激响应的方法。 3．理解系统零、极点分布与系统稳定性关系。 三、涉及实验的相关情况介绍（包含使用软件或实验设备等情况） 计算机一台(安装MATLAB6.5版本或以上版本) 四、实验试做记录（含程序、数据记录及分析） 1．设) 2)(1()(p s p s s s H --= 设①p1=-2，p2=-30； ②p1=-2，p2=3 （1）针对极点参数①②，画出系统零、极点分布图，判断该系统稳定性。 程序： clear num=[1]; den=[1,32,60]; [r,poles,k]=residue(num,den); p=roots(den);
z=roots(num); subplot(2,2,1);plot(real(p),imag(p),'*');hold on; plot(real(z),imag(z),'o');grid on T=0:0.1:10; y1=impulse(num,den,T); subplot(2,2,2);plot(T,y1);grid;title('脉冲响应') 结果： 位于s左半平面，稳定 （2）针对极点参数①②，绘出系统的脉冲响应曲线，并观察t→∞时，脉冲响应变化趋势。 程序： clear num=[1]; den=[1,-1,-6]; [r,poles,k]=residue(num,den); p=roots(den); z=roots(num); subplot(2,2,1);plot(real(p),imag(p),'*');hold on; plot(real(z),imag(z),'o');grid on T=0:0.1:10; y1=impulse(num,den,T); subplot(2,2,2);plot(T,y1);grid;title('脉冲响应') 结果：