武汉工程大学信号分析与处理实验文档
信号分析与处理实验报告
实验一1、基本信号的表示及可视化(1)单位冲激信号 (t)程序：t=-1:0.001:1; %定义时间向量for i=1:3; %采用循环语句观察i取不同值时的图形dt=1/(i^4);X=(1/dt)*((t>=(-1/2*dt))-(t>=(1/2*dt)));%计算函数值subplot(1,3,i);%将图像分成三部分同时观察stairs(t,X);title('单位冲激信号δ(t)');end（2）单位阶跃信号程序：t=-0.5:0.001:1;%定义时间变量，间隔为0.001S=stepfun(t,0);%定义单位阶跃信号S1=stepfun(t,0.5);%定义单位阶跃延迟信号figure(1);plot(t,S);axis([-0.5 1 -0.2 1.2]);title('单位阶跃信号')%画出图形figure(2);plot(t,S1);axis([-0.5 1 -0.2 1.2]);title('单位阶跃延迟信号')（3）抽样信号f=sin(t)./t;程序：t=-10:0.6:10; %向量t时间范围t=t1:p:t2,p为时间间隔f=sin(t)./t;plot(t,f,'o'); %显示该信号的时域波形title('f(t)=Sa(t)(时间间隔为0.6s)'); %标题xlabel('t') %横坐标标题axis([-10,10,-0.4,1.1]) %横坐标和纵坐标范围（4）单位样值序列和单位阶跃序列A．单位序列δ(k)B．单位阶跃序列ε(k)程序：n1=-10;n2=10;%输入序列的起始点n=n1:n2;k=length(n);x1=zeros(1,k);x1(1,-n1+1)=1;%产生单位样值序列subplot(2,1,1);%绘图stem(n,x1,'filled');x2=ones(1,k);x2(1,1:-n1)=0;subplot(2,1,2);stem(n,x2,'filled');2、信号的频域分析已知周期方波信号0||2()0||22E t f t T t ττ⎧<⎪⎪=⎨⎪<<⎪⎩，当02T τ=， 04T τ=，08T τ=时，画出其幅度谱和相位谱，观察不同周期下，()f t 的频谱图有何区别。
(完整word版)信号分析与处理实验报告程序+报告内容(非常全)
0,11)N -的输入序列，把它按⎩⎨⎧+=)12()2(r x r x 2N W ，title('傅立叶反变换');3.2.4 实验结果图（workspace里的各数据结果在此不一一列出）由图可知，离散时间信号经过傅里叶变换和傅里叶反变换后，没有发生改变，说明了程序的正确性。
3.2.5 频谱混叠%% 验证主程序clc;clear allTs=0.005; %采样时间间隔N=64; %采样点数T0=N*Ts; %采样时间长度M=N*Ts/T0; %M*T0为采样时间，M=1表示采样为整周期采样,M不等于1会发生频谱泄露M1=N1*Ts1/T1;for n=1:N1xx1(n)=10*sin(2*pi*50*n*Ts1+pi/3); %要进行傅立叶变换的连续时间函数endyy1=fft(xx1,N1); %对时域信号做傅立叶变换deltaf1=1/(T1*M1); %书上（6-78）for n=1:N1/2+1yy2(n)=yy1(n); %书上（6-76）endsubplot(2,1,2)f1=0:deltaf1:N1/2*deltaf1; %stem(f1,abs(yy2))xlabel('f/Hz');ylabel('|X(k)|/(A)');title('傅立叶变换');得到图像：四、实验分析1、比较快速傅里叶变换与离散傅里叶变换计算效率对于N 点DFT ，需进行2N 次复数乘法及(1)N N - 次复数加法，而对于同样点数的FFT ，则只需2log 2NN 次复数乘法和 2Nlog N 次复数加法，大大加快了运算速度。
2、离散傅里叶变换与z 变换的关系：有限长序列()x n 的离散傅里叶变换也可以定义为它的z 变换在z 平面单位圆上N 等分的均匀采样，即()[()]()k Nz W X k DFT x n X z -===同时，z 变换可由离散傅里叶变换通过内插函数构造。
信号分析与处理实验指导书3
二、实验原理 信号的基本运算包括信号的相加和相乘。信号的时域变换包括信号的平移、
反折、倒相以及尺度变换。
三、涉及的 MATLAB 函数 fliplr 函数 功能:实现矩阵行元素的左右翻转。 调用格式: B=fliplr(A):其中 A 指要翻转的矩阵。
四、实验内容与方法
1、常规计算
1)序列的加法
x1=-2:2;
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('连续信号的相加');
if (n1>n2||n0>n2||n0<n1)
error('parameter error');
end;
if (n1<=n2)
for n=1:n2-n1+1
if (n==n0)
x(1,n)=n1-1+n;
x(2,n)=1;
end;
x(1,n)=n1-1+n;
x(2,n)=0;
end;
x(2,n0-n1+1)=1;
4)翻转 信号的翻转就是将信号的波形以纵轴为对称轴翻转 180°，将信号 f(t)中的
自变量 t 替换为-t 即可得到其翻转信号。 clear all t=0:0.02:1;t1=-1:0.02:0; g1=3*t; g2=3*(-t1); grid on; plot(t,g1,'--',t1,g2); xlabel('t');ylabel('g(t)'); title('信号的翻转'); 5)尺度变换 将信号 f(t)中的自变量 t 替换为 at clear all; t=0:0.001:1; a=2; y=sin(2*pi*t); y1=sin(2*a*pi*t); subplot(211) plot(t,y); ylabel('y(t)');xlabel('t'); title('尺度变换');
信号分析与处理实验指导书
实验一 信号频谱的测量一、实验目的1、掌握信号频谱的测量方法，加深对周期信号频谱特点的了解。
2、研究矩形脉冲时域周期和脉宽的变化对频谱结构的影响，了解时域和频域间的关系。
3、学习TH-SG01P 型功率函数信号发生器各旋钮、开关的作用及其使用方法。
4、学习虚拟示波器的使用方法。
二、原理及说明1、周期信号的频谱分为幅度谱、相位谱和功率谱三种，分别是信号各频率分量的振幅，初相和功率按频率由低到高依次排列构成的图形。
通常讲的频谱指幅度谱，它可选频表或波形分析仪逐个频率测试而得，也可用频率谱仪直接显示，现在更多的是应用虚拟示波器的FFT 变换来实现。
2、连续周期信号频谱的特点是离散性、谐波性和幅度总趋势的收敛性，可以通过对正弦波、三角波、方波（或矩形脉冲）频谱的具体测试而得到验证。
（1）、正弦波的频谱特别简单，即本身频率的振幅，如图1-1所示。
图1-1 正弦波及其频谱（2）、宽度为2τ，高度为A 的三角波的频谱，当2T τ=时，2()2k k A A Sa π=，如图1-2所示。
图1-2 三角波及其频谱ω12ω1k ω13ω ω1ωAk A13ω15ωω12ω 24/(5)A π24/(3)A π /2A1k ωω1ω24/A πkA（3）、矩形脉冲的频谱，122k k A A Sa Tωττ⎛⎫=⎪⎝⎭。
当为方波2T τ=时，12k k A A Sa ωτ⎛=⎝图1-3 （4）、周期型矩形脉冲的频谱按122k A Sa Tωττ⎛⎫⎪⎝⎭规律变化，它的第一个零点频率2πτ取决于脉宽τ，谱线的疏密取决于周期T 。
当脉宽τ不变时，在20πτ内谱线会增多而变密；当周期T 不变而脉宽τ减小时，其第一零点频率会增高，从而使20πτ内的谱线增多；谱线高度都会因T 增大或τ减小而降低。
因此，信号的波形和其频谱间是一一对应的，它们不过是对同一信号的两种不同描述方式罢了。
在频域中，常把20πτ的一段频率范围定义为信号的有效频带宽度，对于5T τ≥的矩形脉冲，这种定义就比较精确了。
信号分析与处理实验报告
《信号分析与处理》实验报告华北电力大学前言1．实验总体目标通过实验，巩固掌握课程的讲授内容，使学生对信号分析与线性系统分析的基本理论及分析方法有一个感性认识和更好地理解，使学生在分析问题与解决问题的能力及实践技能方面有所提高。
2．适用专业自动化专业本科生3．先修课程信号分析与处理4．实验课时分配5需要配置微机及MATLAB工具软件。
6．实验总体要求1、掌握信号分解的基本思想及信号在时域、频域和变换域进行分解的基本理论及描述方法，用MATLAB编程语言实现基本信号的表示及可视化，计算和分析信号的频谱；2、掌握在时域、频域和变换域分析LTI系统的方法，及系统在时域、频域和变换域的描述方法，用MATLAB编程语言实现LTI系统的时域分析及频率分析。
3、掌握信号的调制与解调，用MATLAB编程语言仿真分析信号的调制与解调。
⒎ 本实验的重点、难点及教学方法建议实验通过MATLAB编程语言来实现基本信号的表示及可视化，计算分析信号的频谱，实现LTI系统的时域分析及频率分析，并仿真分析信号的调制与解调，使学生对信号分析与线性系统分析的基本理论及分析方法有一个感性认识和更好地理解。
实验的重点及难点是：掌握基本信号的数学表示，信号的频谱特点，计算LTI系统的典型响应，掌握信号的调制与解调。
在这样的理论基础上，学会用MATLAB编程语言来实现对信号与系统响应的可视化及对数字滤波器进行设计。
教学建议：打好理论基础，熟练编程语言。
目录实验一信号的时域与频域分析 3实验二信号的时域与频域处理 4实验三数字滤波器的设计 5实验一一、实验目的1、熟悉MATLAB 平台，高效的数值计算及符号计算功能；2、实现基本信号的表示及可视化计算；3、分析信号的频谱。
二、 实验类型验证型 三、 实验仪器微机，MATLAB 工具软件。
四、 实验原理MATLAB 是功能强大的数学软件，它提供了计算周期连续函数和周期离散序列的频谱的一系列函数。
信号分析与处理实验报告
华北电力大学实验报告||实验名称FFT的软件实现实验（Matlab）IIR数字滤波器的设计课程名称信号分析与处理||专业班级：电气化1308 学生姓名：袁拉麻加学号： 2 成绩：指导教师：杨光实验日期： 2015-12-17快速傅里叶变换实验一、实验目的及要求通过编写程序，深入理解快速傅里叶变换算法（FFT）的含义，完成FFT和IFFT算法的软件实现。
二、实验内容利用时间抽取算法，编写基2点的快速傅立叶变换（FFT）程序；并在FFT程序基础上编写快速傅里叶反变换（IFFT）的程序。
三：实验要求1、FFT和IFFT子程序相对独立、具有一般性，并加详细注释；2、验证例6-4，并能得到正确结果。
3、理解应用离散傅里叶变换（DFT）分析连续时间信号频谱的数学物理基础。
四、实验原理：a.算法原理1、程序输入序列的元素数目必须为2的整数次幂，即N=2M,整个运算需要M 级蝶形运算；2、输入序列应该按二进制的码位倒置排列，输出序列按自然序列排列；3、每个蝶形运算的输出数据军官占用其他输入数据的存储单元，实现“即位运算”；4、每一级包括N/2个基本蝶形运算，共有M*N/2个基本蝶形运算；5、第L级中有N/2L个群，群与群的间隔为2L。
6、处于同一级的各个群的系数W分布相同，第L级的群中有2L-1个系数；7、处于第L级的群的系数是(p=1,2,3，…….,2L-1)而对于第L级的蝶形运算，两个输入数据的间隔为2L-1。
b.码位倒置程序流程图开始检测A序列长度nk=0j=1x1(j)=bitget(k,j);j=j+1Yj<m?Nx1=num2str(x1);y(k+1)=bin2dec(x1);clear x1k=k+1c.蝶形运算程序流程图五、程序代码与实验结果a.FFT程序：%%clear all;close all;clc;%输入数据%A=input('输入x(n)序列','s');A=str2num(A);% A=[1,2,-1,4]; %测试数据%%%%校验序列,%n=length(A);m=log2(n);if (fix(m)~=m)disp('输入序列长度错误，请重新输入!');A=input('输入x(n)序列','s');A=str2num(A);elsedisp('输入正确，请运行下一步')end%%%码位倒置%for k=0:n-1for j=1:m %取M位的二进制数%x1(j)=bitget(k,j); %倒取出二进制数%endx1=num2str(x1); %将数字序列转化为字符串%y(k+1)=bin2dec(x1); %二进制序列转化为十进制数%clear x1endfor k=1:nB(k)=A(y(k)+1); %时间抽取序列%endclear A%%%计算%for L=1:m %分解为M级进行运算%LE=2^L; %第L级群间隔为2^L%LE1=2^(L-1); %第L级中共有2^(L-1)个Wn乘数，进行运算蝶运算的两数序号相隔LE1%W=1;W1=exp(-1i*pi/LE1);for R=1:LE1 %针对第R个Wn系数进行一轮蝶运算，共进行LE1次%for P=R:LE:n %每个蝶的大小为LE% Q=P+LE1;T=B(Q)*W;B(Q)=B(P)-T;B(P)=B(P)+T;endW=W*W1;endendB %输出X(k)%%%验证结果：例6-4b.IFFT程序：%%clear all;close all;clc;%输入数据%A=input('输入X(k)序列','s');A=str2num(A);% A=[6,2+2i,-6,2-2i]; %测试数据%%%%校验序列,%n=length(A);m=log2(n);if (fix(m)~=m)disp('输入序列长度错误，请重新输入!');A=input('输入x(n)序列','s');A=str2num(A);elsedisp('输入正确，请运行下一步')end%%%码位倒置%for k=0:n-1for j=1:m %取M位的二进制数%x1(j)=bitget(k,j); %倒取出二进制数%endx1=num2str(x1); %将数字序列转化为字符串%y(k+1)=bin2dec(x1); %二进制序列转化为十进制数%clear x1endfor k=1:nB(k)=A(y(k)+1); %时间抽取序列%endclear A%%%计算%for L=1:m %分解为M级进行运算%LE=2^L; %第L级群间隔为2^L%LE1=2^(L-1); %第L级中共有2^(L-1)个Wn乘数，进行运算蝶运算的两数序号相隔LE1%W=1;W1=exp(-1i*pi/LE1);for R=1:LE1 %针对第R个Wn系数进行一轮蝶运算，共进行LE1次%for P=R:LE:n %每个蝶的大小为LE%Q=P+LE1;T=B(Q)*W;B(Q)=B(P)-T;B(P)=B(P)+T;endW=W*W1;endendB=conj(B); %取共轭%B=B/n %输出x(n)%验证结果：六、实验心得与结论本次实验借助于Matlab软件，我避开了用C平台进行复杂的复数运算，在一定程度上简化了程序，并添加了简单的检错代码，码位倒置我通过查阅资料，使用了一些函数，涉及到十-二进制转换，数字-文本转换，二-文本转换，相对较复杂，蝶运算我参考了书上了流程图，做些许改动就能直接实现。
信号分析与处理实验
n ωj �
e� n� h � � �
0�n 1� N
ωj
e� H
波滤想理的望希所近逼来 � n� h 列序应响样抽位单长限有用是理原本基的法计设数函窗
� n� ω 数函窗用以所 �的果因非是且 �列序长限无是往往 � n� d h 于由 。 � n� d h 应响样抽位单的器
�应响样抽位单的器波滤字数 RIF 的计设际实到得�理处权加行进并�断截 � n� d h 将
�s
式公
换变在以所�量变参的要紧关无个一了成就 T 期周样采�后标指域字数出给为因 。响影有没果结计设对 s T 变改�答 �么什为�响影无有果结计设对值取的 s T 的中
1� 1�
z � 1 sT z �1 2
�s
式公换变�中程过器波滤字数计设法换变性线双用
题考思�六�
。滑平更形波使�号信用有留保�号信扰干除滤是用作的器波滤�答
5x 4x的 �� x n x� 4 x为 �4 �的 ��7n 为 期 8n是 周 n�� � �因 x 结的确正是 61=N �期 � n� � 7果 x1 � � � n� 7 x 使果结�61 是期周的 � n� 5 x � � n� 56 �n �周
� n� 5 xj � � n� 4 x � � n� 8x .8
s
�
减衰小最带阻� p � 减衰大最带通� s � 率频止截限下带阻� p � 率频止截带通�有数参 。器波滤计设来换变性线双用采能才 �时偿补被能或的许容是真失性线非种这当有只 �答 �性特和数参的器波滤录记�四�
� 2 �1 � N � �
k � 1�
� z � neve H
N � k � 1�
性称对足满须必还 � n� h 则�位相性线求要果如
信号分析与处理实验报告
信号分析与处理实验报告一、实验目的1.了解信号分析与处理的基本概念和方法；2.掌握信号分析与处理的基本实验操作；3.熟悉使用MATLAB进行信号分析与处理。
二、实验原理信号分析与处理是指利用数学和计算机技术对信号进行分析和处理的过程。
信号分析的目的是了解信号的特性和规律，通过对信号的频域、时域和幅频特性等进行分析，获取信号的频率、幅度、相位等信息。
信号处理的目的是对信号进行数据处理，提取信号的有效信息，优化信号的质量。
信号分析和处理的基本方法包括时域分析、频域分析和滤波处理。
时域分析主要是对信号的时变过程进行分析，常用的方法有波形分析和自相关分析。
频域分析是将信号转换到频率域进行分析，常用的方法有傅里叶级数和离散傅里叶变换。
滤波处理是根据信号的特性选择适当的滤波器对信号进行滤波，常用的滤波器有低通滤波器、高通滤波器、带通滤波器和带阻滤波器等。
三、实验内容1.信号的时域分析将给定的信号进行波形分析，绘制信号的时域波形图；进行自相关分析，计算信号的自相关函数。
2.信号的频域分析使用傅里叶级数将信号转换到频域，绘制信号的频域图谱；使用离散傅里叶变换将信号转换到频域，绘制信号的频域图谱。
3.滤波处理选择合适的滤波器对信号进行滤波处理，观察滤波前后的信号波形和频谱。
四、实验步骤与数据1.时域分析选择一个信号进行时域分析，记录信号的波形和自相关函数。
2.频域分析选择一个信号进行傅里叶级数分析，记录信号的频谱；选择一个信号进行离散傅里叶变换分析，记录信号的频谱。
3.滤波处理选择一个信号，设计适当的滤波器对信号进行滤波处理，记录滤波前后的信号波形和频谱。
五、实验结果分析根据实验数据绘制的图像进行分析，对比不同信号在时域和频域上的特点。
观察滤波前后信号波形和频谱的变化，分析滤波效果的好坏。
分析不同滤波器对信号的影响，总结滤波处理的原理和方法。
六、实验总结通过本次实验，我们了解了信号分析与处理的基本概念和方法，掌握了信号分析与处理的基本实验操作，熟悉了使用MATLAB进行信号分析与处理。
武汉工程大学数字信号处理实验二 时域离散系统及系统响应
实验二时域离散系统及系统响应一、实验目的1、掌握求解离散时间系统冲激响应和阶跃响应的方法；2、进一步理解卷积定理，掌握应用线性卷积求解离散时间系统响应的基本方法；3、掌握离散系统的响应特点。
二、实验内容1、请分别用impz 和dstep函数求解下面离散时间系统的冲激响应和阶跃响应。
（1）系统的差分方程为：)ynnxy-=(n-+y+n)2.0866((648.0()1.0)a=[1,-0.8,0.64];b=[0.866,0,0];n=20;hn=impz(b,a,n); %冲激响应gn=dstep(b,a,n); %阶跃响应subplot(2,1,1),stem(hn,'filled'); %显示冲激响应曲线title('系统的单位冲激响应');ylabel('h(n)');xlabel('n');axis([0,n,1.1*min(hn),1.1*max(hn)]);subplot(2,1,2),stem(gn,'filled'); %显示阶跃响应曲线title('系统的单位阶跃响应');ylabel('g(n)');xlabel('n');axis([0,n,1.1*min(gn),1.1*max(gn)]);2468101214161820-0.4-0.200.20.40.60.8系统的单位冲激响应h (n )n246810121416182011.21.41.6系统的单位阶跃响应g (n )n（2）系统的系统函数为：21115.01)(---+--=z z z z Ha=[1,-1,1]; b=[1,-0.5,0]; n=20;hn=impz(b,a,n); %冲激响应 gn=dstep(b,a,n); %阶跃响应subplot(2,1,1),stem(hn,'filled'); %显示冲激响应曲线 title('系统的单位冲激响应');ylabel('h(n)');xlabel('n');axis([0,n,1.1*min(hn),1.1*max(hn)]);subplot(2,1,2),stem(gn,'filled'); %显示阶跃响应曲线 title('系统的单位阶跃响应');ylabel('g(n)');xlabel('n');axis([0,n,1.1*min(gn),1.1*max(gn)]);2468101214161820-1-0.500.51系统的单位冲激响应h (n )n2468101214161820-0.500.511.5系统的单位阶跃响应g (n )n2、运行例题2.3，理解卷积过程和程序中每一句的意义。
信号分析与处理实验报告
实验一图像信号频谱分析及滤波一：实验原理FFT不是一种新的变化，而是DFT的快速算法。
快速傅里叶变换能减少运算量的根本原因在于它不断地把长序列的离散傅里叶变换变为短序列的离散傅里叶变换，在利用的对称性和周期性使DFT运算中的有些项加以合并，达到减少运算工作量的效果。
为了消除或减弱噪声，提取有用信号，必须进行滤波，能实现滤波功能的系统成为滤波器。
按信号可分为模拟滤波器和数字滤波器两大类。
数字滤波器的关键是如何根据给定的技术指标来得到可以实现的系统函数。
从模拟到数字的转换方法很多，常用的有双线性变换法和冲击响应不变法，本实验主要采用双线性变换法。
双线性变换法是一种由s平面到z平面的映射过程，其变换式定义为：数字域频率与模拟频率之间的关系是非线性关系。
双线性变换的频率标度的非线性失真是可以通过预畸变的方法去补偿的。
变换公式有Ωp=2/T*tan(wp/2)Ωs=2/T*tan(ws/2)二：实验内容1．图像信号的采集和显示选择一副不同彩色图片，利用Windows下的画图工具，设置成200*200像素格式。
然后在Matlab软件平台下，利用相关函数读取数据和显示图像。
要求显示出原始灰度图像、加入噪声信号后的灰度图像、滤波后的灰度图像。
2．图像信号的频谱分析要求分析和画出原始灰度图像、加入噪声信号后灰度图像、滤波后灰度图像信号的频谱特性。
3.数字滤波器设计给出数字低通滤波器性能指标:通带截止频率fp＝10000 Hz，阻带截止频率fs＝15000 Hz，阻带最小衰减Rs＝50 dB，通带最大衰减Rp＝3 dB，采样频率40000Hz。
三：实验程序clear allx=imread('D:\lan.jpg');%原始彩色图像的数据读取x1=rgb2gray(x);%彩色图像值转化为灰度图像值[M,N]=size(x1);%数据x1的长度,用来求矩阵的大小x2=im2double(x1);%unit8转化为double型x3=numel(x2);%计算x2长度figure(1);subplot(1,3,1);imshow(x2);title('原始灰度图')z1=reshape(x2,1,x3);%将二维数据转化成一维数据g=fft(z1);%对图像进行二维傅里叶变换mag=fftshift(abs(g));%fftshift是针对频域的，将FFT的DC分量移到频谱中心K=40000;Fs=40000;dt=1/Fs;n=0:K-1;f1=18000;z=0.1*sin(2*pi*f1*n*dt);x4=z1+z;%加入正弦噪声f=n*Fs/K;y=fft(x4,K);z2=reshape(x4,M,N);%将一维图转换为二维图subplot(1,3,2);imshow(z2);title('加入噪声后')g1=fft(x4);mag1=fftshift(abs(g1));%设计滤波器ws=0.75*pi;wp=0.5*pi;fs=10000;wp1=2*fs*tan(wp/2);ws1=2*fs*tan(ws/2);rs=50;rp=3;% [n,wn]=buttord(wp/pi,ws/pi,rp,rs);% [bz,az]=butter(n,wn);[n,wn]=buttord(wp1,ws1,rp,rs,'s');[z,p,k]=buttap(n);[b,a]=zp2tf(z,p,k);[B,A]=lp2lp(b,a,wn);[bz,az]=bilinear(B,A,fs);[h,w]=freqz(bz,az,128,fs);L=numel(z2);z3=reshape(z2,1,L);x6=filter(bz,az,double(z3));x7=reshape(x6,M,N);subplot(1,3,3);imshow(x7);g2=fft(x6);mag2=fftshift(abs(g2));title('滤波后')%建立频谱图figure(2);subplot(1,3,1);plot(mag);title('原始Magnitude')subplot(1,3,2);plot(mag1);title('加噪声Magnitude')subplot(1,3,3);plot(mag2);title('滤波后Magnitude')figure(3);subplot(1,2,1)plot(w,abs(h));xlabel('f');ylabel('h');title('滤波器幅谱');subplot(1,2,2);plot(w,angle(h));title('滤波器相谱');四：实验结果与分析图一图二分析：由图二可以知道加入噪声后的幅值谱和原始图的幅值谱明显多了两条幅值线，而这两条幅值线就是我们对原始灰度图加入的正弦噪声，而相应的图一中的加噪声后的图与原始图相比，出现了明显的变化。
大学信号分析实验报告
一、实验目的1. 理解信号分析的基本概念和原理；2. 掌握信号的时域和频域分析方法；3. 熟悉MATLAB在信号分析中的应用；4. 培养实验操作能力和数据分析能力。
二、实验原理信号分析是研究信号特性的科学，主要包括信号的时域分析和频域分析。
时域分析关注信号随时间的变化规律，频域分析关注信号中不同频率分量的分布情况。
1. 时域分析：通过对信号进行采样、时域卷积、微分、积分等操作，分析信号的时域特性。
2. 频域分析：通过对信号进行傅里叶变换、频域卷积、滤波等操作，分析信号的频域特性。
三、实验内容1. 信号采集与处理（1）采集一段语音信号，利用MATLAB的录音功能将模拟信号转换为数字信号。
（2）对采集到的信号进行采样，选择合适的采样频率，确保满足奈奎斯特采样定理。
（3）绘制语音信号的时域波形图，观察信号的基本特性。
2. 信号频谱分析（1）对采集到的信号进行傅里叶变换，得到信号的频谱。
（2）绘制信号的频谱图，分析信号的频域特性。
3. 信号滤波（1）设计一个低通滤波器，滤除信号中的高频噪声。
（2）将滤波后的信号与原始信号进行对比，分析滤波效果。
4. 信号调制与解调（1）对原始信号进行幅度调制，产生已调信号。
（2）对已调信号进行解调，恢复原始信号。
（3）分析调制与解调过程中的信号变化。
四、实验步骤1. 采集语音信号，将模拟信号转换为数字信号。
2. 对采集到的信号进行采样，确保满足奈奎斯特采样定理。
3. 绘制语音信号的时域波形图，观察信号的基本特性。
4. 对信号进行傅里叶变换，得到信号的频谱。
5. 绘制信号的频谱图，分析信号的频域特性。
6. 设计低通滤波器，滤除信号中的高频噪声。
7. 对滤波后的信号与原始信号进行对比，分析滤波效果。
8. 对原始信号进行幅度调制，产生已调信号。
9. 对已调信号进行解调，恢复原始信号。
10. 分析调制与解调过程中的信号变化。
五、实验结果与分析1. 时域分析通过观察语音信号的时域波形图，可以看出信号的基本特性，如信号的幅度、频率等。
信号分析与处理(第3版)-第3章part1(时域分析)
14
五、离散信号的描述-序列的表示方法
• 集合表示法：
{x(n)}={……, 0，1，2，3, 4，3，2，1，0，……}
n=0
n值规定为自左向右逐一递增
• 公式表示法： x(n) 4 n , n 3
x(n)
• 图形表示法：
4
3
2 1
5 4 3 2 1 0 1 2 3 4 5 n
15
1、单位脉冲序列
奈奎斯特（Nyquist）频率： s 2m
10
2、由抽样信号恢复原连续信号
• 取主频带 X () ：
• 时域卷积定理： X () X s ()H ()
xs (t) x(nTs ) (t nTs ) n
h(t )
c
Sa( ct )
x(t) xs (t) * h(t)
n
c
x(nTs
• 频谱发生了周期延拓，即将原连续信号的频 谱X()分别延拓到以±s， ±2s ……为中心的
频谱，其中s为采样角频率
• 频谱的幅度乘上了因子1/Ts，其中Ts为采样周 期
9
二、时域采样定理
对于频谱受限的信号x(t)，如果其最高频率分量为 ωm，为了保留原信号的全部信息，或能无失真地恢 复原信号，在通过采样得到离散信号时，其采样频 率应满足ω s ≥ 2ωm
• 预习内容：
• 离散信号的频域分析
• 实验1：信号的采样与恢复
34
•即
y(n) {1,1,4,23,32,13,34,21,5,20} 32
7、两序列相关运算
• 序列的相关运算被定义为
xy (n) x(m) y(n m) m
• 可以用卷积符号“＊”来表示相关运算
xy (n) x(n) * y(n)
信号分析与处理实验报告格式(1)概要
北京工商大学信号与系统实验姓名刘冰璇学号 1404010415班级信息142教师陈秀新2015. .实验一基本信号的产生一实验目的学习使用MATLBA产生信号，绘制信号波形，实现信号的基本运算，为信号分析和系统设计奠定基础。
二实验内容1 利用MATLBA产生下列连续信号并作图：(1) x(t)=-2u(t-1), -1<t<5t= -1: 0.02:5; p=t-1; x=-2.*(p>=0); plot(t,x)(2) x(t)=e-0.1t sin(2/3t), 0<t<30t = 0: 0.1: 30; x = exp(-0.1*t); y=sin((2/3).*t); A=x.*y; plot(t,A);(3) x(t)=cos(100t)+cos(3000t), -0.1<t<0.1t=-0.1:0.001:0.1; x=cos(100.*t)+cos(3000.*t); plot(t,x);(4) x(t)=cos(0.1πt)cos(0.8πt),0<t<200>> t=0:0.01:200;x=cos(0.1*pi*t).*cos(0.8*pi*t);>> plot(t,x);xlabel('x=cos(0.1*pi*t).*cos(0.8*pi*t)');axis([0,30,-1,1]);2、利用MA TLBA产生下列信号并作图(1)15k15,0551][≤<-⎩⎨⎧≤-=，其他，kxk= -15: 15; x=[zeros(1,10),ones(1,11),zeros(1,10)]; stem(k,x)(2) x[k]=(0.9)k [sin(0.25k)+cos(0.25k)],-20<k20k=-20:20; x=(0.9.^k) y=sin(0.25*pi.*k)+cos(0.25*pi.*k); A=x.*y; stem(k,A);？？没有图像3、已知序列：x[k]=[1,2,0,-1,3,2;k=-2,-1,0,1,2,3],h[k]=[1,-1,1;k=0,1,2]。
数字信号处理实验报告 (3)
武汉工程大学实验报告实验课程数字信号处理一、实验目的(1)加深对离散傅里叶变换(DFT)基本概念的理解。
(2)了解有限长序列傅里叶变换(DFT)与周期序列傅里叶级数(DFS)、离散时间傅里叶变换(DTFT)的联系。
(3)掌握用MA TLAB语言进行离散傅里叶变换和逆变换的方法。
二、实验内容1.有限长序列的傅里叶变换(DFT)和逆变换(IDFT)2.有限长序列DFT与周期序列DFS的联系3.有限长序列DFT与离散时间傅里叶变换DTFT的联系三、实验环境MA TLAB7.0四丶：实验内容、原理描述及实验结果1.离散时间信号的表示离散时间信号定义为一时间函数，它只在某些离散的瞬时给出函数值，而在其他处无定义。
因此，它是时间上不连续按一定先后次序排列的一组数的集合，故称为时间序列，简称序列，通常表示为{x（n）} -∞<n<+∞(1)单位抽样序列用Matlab编写的实验程序n0=0;n1=-5;n2=5;n=[n1:n2];nc=length(n);x=zeros(1,nc);for i=1:ncif n(i)==n0x(i)=1;endendstem(n,x)xlabel('n');ylabel('x(n)');title('单位抽样序列');grid（2）单位阶跃序列用Matlab编写编写的生成单位阶跃序列的函数n0=0;n1=-5;n2=5;n=[n1:n2];x=[(n-n0)>=0];stem(n,x)xlabel('n');ylabel('x(n)');title('单位阶跃序列');grid图形如下（3）指数序列程序代码如下：n=[0:20];x=(0.78).^n;stem(n,x)xlabel('n');ylabel('x(n)');title('指数序列'); grid图形如下：（4）正余弦序列用matlab编写正弦序列x(n)=5sin(0.1πn+π/3)函数的程序。
信号分析报告与处理实验报告材料(基于某MATLAB)
实用文档武汉工程大学电气信息学院2、四、思考：1、为什么图二中t=0处曲线是间断的，如何使其成为连续的曲线？因为axis函数对纵坐标的的上边界限定过小，使图形在边界处不能完整的显示。
2.3.四、思考：1、代数运算符号*和.*的区别是？*是矩阵相乘，是矩阵A行元素与B的列元素相乘的和.*是数组相乘，表示数组A和数组B中的对应元素相乘实验容实验三连续时间信号的卷积一、实验容1、已知两连续时间信号如下图所示，绘制信号f1(t)、f2(t)及卷积结果f(t)的波形；设时间变化步长dt分别取为0.5、0.1、0.01，当dt取多少时，程序的计算结果就是连续时间卷积的较好近似？2.实验容2.实验容实验五连续时间信号的频域分析一、实验容1、如图5.4所示的奇谐周期方波信号，周期为T1=1，幅度为A=1，将该方波信号展开成三角形式Fourier级数并分别采用频域矩形窗和Hanning窗加权，绘制两种窗函数加权后的方波合成图像。
时间围取为-2~2，步长值取为0.01。
2、将图5.5中的锯齿波展开为三角形式Fourier级数，按(2)式求出Fourier级数的系数，并在频域分别采用矩形窗、Hanning窗和三角窗加权，观察其Gibbs效应及其消除情况。
时间围取为-2~2，步长值取为0.01。
3、选做：编程计算连续时间周期信号的三角形式傅里叶级数展开的系数二、实验方法与步骤1、将方波信号展开成三角形式Fourier级数并分别采用频域矩形窗和Hanning窗加权方波展开的三角式傅立叶级数为：()()tkkt xLk1,5,3,1sin4ωπ⋅∑=∞=采用频域矩形窗加权，则展开式变为：()()()[]tkkt x Kk112sin124ωπ+⋅+∑==采用Hanning窗加权，则展开式变为：()()()()[]tkKkkt x Kk112sin122cos5.05.0124ωππ+⋅⎥⎦⎤⎢⎣⎡++⋅+∑==a0=2/T*int(f,t,0,T); %求函数f对t从0到T的定积分a0=simplify(a0) %得出结果syms kfa=t*cos(k*w*t);fb=t*sin(k*w*t);ak=2/T*int(fa,t,0,T); %求函数fa对t从0到T的定积分bk=2/T*int(fb,t,0,T); %求函数fb对t从0到T的定积分ak=simplify(ak)bk=simplify(bk)三、实验数据与结果分析1.2.3.实验容根据绘制的幅频特性曲线，系统具有低通滤波特性2.根据绘制的幅频特性曲线，系统具有带通滤波特性。
信号分析与处理实验报告
信号分析与处理实验报告
班级_________________________
学生姓名_________________________
学号_________________________
所在专业_________________________
成绩_________________________
上海大学
二0 0 年月日
图1-2 芯片参数设置界面
4. 利用数字公式编程生成正弦波、噪声或三角波等数字信号，可以选择其中一种信号，
图3-1 滤波器的种类
下图是用带通滤波器消除信号钢管无损探伤信号中由于传感器晃动带来的低频干扰，以及由于电磁噪声等带来的高频干扰的例子。
用滤波器消除信号中的干扰
图3-3 滤波器的作用实验
下面是该实验的装配图和信号流图，图中线上的数字为连接软件芯片的软件总线数
图3-4 滤波器的作用实验装配图。
实验二离散时间系统的时域分析实验
数字信号处理——实验二武汉工程大学电气信息学院通信工程红烧大白兔一、实验目的1、在时域中仿真离散时间系统，进而理解离散时间系统对输入信号或延时信号进行简单运算处理，生成具有所需特性的输出信号的方法。
2、仿真并理解线性与非线性、时变与时不变等离散时间系统。
3、掌握线性时不变系统的冲激响应的计算并用计算机仿真实现。
4、仿真并理解线性时不变系统的级联、验证线性时不变系统的稳定特性。
二、实验设备计算机，MATLAB语言环境三、实验根底理论1、系统的线性线性性质表现为系统满足线性叠加原理：假设某一输入是由N个信号的加权和组成的，输出就是由系统对这N个信号中每一个的响应的相应加权和组成的。
设x1〔n〕和〔n〕分别作为系统的输入序列，其输出分别用y1(n)和y2(n)表示，即Y1(n)=T[x1(n)],y2(n)=T[x2(n)]假设满足T[a1x1(n)+a2x2(n)]=a1y1(n)+a2y2(n)x2那么那么该系统服从线性叠加原理，或者称为该系统为线性系统。
2、系统的时不变特性假设系统的变换关系不随时间变化而变化，或者说系统的输出随输入的移位而相应移位但形状不变，那么称该系统为时不变系统。
对于时不变系统，假设y(n)=T[x(n)]那么T[x(n-m)]=y(n-m)3、系统的因果性系统的因果性既系统的可实现性。
如果系统n时刻的输出取决于n时刻及n时刻以前的输入，而和以后的输入无关，那么该系统是可实现的，是因果系统。
系统具有因果性的充分必要条件是h(n)=0,n<04、系统的稳定性稳定系统是指有界输入产生有界输出〔BIBO)的系统。
如果对于输入序列x(n),存在一个不变的正有限值M，对于所有n值满足|x(n)|≤M<∞那么称该输入序列是有界的。
稳定性要求对于每个有界输入存在一个不变的正有限值K，对于所有n值，输出序列y(n)满足|y(n)|≤K<∞系统稳定的充分必要条件是系统的单位取样响应绝对可和，用公式表示为|h(n)|n5、系统的冲激响应设系统输入x(n)=δ(n)，系统输出y(n)的初始状态为零，这时系统输出用即h(n)=T[δ(n)]那么称h(n)为系统的单位脉冲响应。
武汉工程大学matlab实验二离散时间信号的分析实验【范本模板】
武汉工程大学数字信号处理实验报告二专业班级：14级通信03班学生姓名：秦重双学号：1404201114实验时间：2017年5月3日实验地点：4B315指导老师: 杨述斌实验一离散时间信号的分析实验一、实验目的①认识常用的各种信号,理解其数学表达式和波形表示。
②掌握在计算机中生成及绘制数值信号波形的方法。
③掌握序列的简单运算及计算机实现与作用。
④理解离散时间傅里叶变换、Z变换及它们的性质和信号的频域特性。
二、实验设备计算机,MATLAB语言环境。
三、实验基础理论1、序列的相关概念离散时间信号用一个称为样本的数字序列来表示。
一般用｛x［n］｝表示，其中自变量n的取值范围是﹣∞到﹢∞之间的整数。
为了表示方便，序列通常直接用x[n]表示。
离散时间信号可以是一个有限长序列，也可以是一个无限长序列。
有限长（也称为有限时宽)序列仅定义在有限的时间间隔中：﹣∞≤N1 ≤N2 ≤+∝。
有限长序列的长度或时宽为N=N1 -N2+1。
满足x[n+kN］=x[n]（对于所有n）的序列称为周期为N的周期序列,其中N取任意正整数；k取任意整数；2、常见序列常见序列有单位取样值信号、单位阶跃序列、矩形序列、斜变序列、单边指数序列、正弦序列、复指数序列等。
3、序列的基本运算序列的基本运算有加法、乘法、倒置（反转）、移位、尺度变换、卷积等。
4、离散傅里叶变换的相关概念5、Z变换的相关概念四．实验内容与步骤1、知识准备认真复习以上基础理论,理解本实验所用到的实验原理。
2、离散时间信号（序列）的产生利用MATLAB语言编程和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形，以加深对离散信号时域表示的理解。
①单位取样值信号Matlab程序x=0;y=1;stem(x，y）；title（'单位样值’）；axis([—2,2，0,1]）；②单位阶跃序列Matlab程序n0=0;n1=—5;n2=5；n=[n1:n2］;x=［（n—n0）>=0];stem(n，x）；xlabel（'n'）；ylabel（’x(n)’）；title(’单位阶跃序列’）;③指数序列、正弦序列Matlab程序n=［0：10］;x=(1/3)。
《信号分析与处理》实验报告
(此文档为word格式，下载后您可任意编辑修改！)序号：号项目名称：《信号分析与处理》实验报告学生学院：信息工程学院专业班级：学生学号：学生姓名：指导老师：朱铮涛2013年12月25日目录实验一、基本信号的产生和时频域抽样实验 (1)一、实验目的 (1)二、实验内容及所得图表 (1)三、思考题解答 (15)实验二、连续和离散系统分析 (16)一、实验目的 (16)二、实验内容和要求 (16)三、思考题解答 (22)实验三、用FFT实现谱分析实验 (23)一、实验目的 (23)二、实验原理 (23)三、实验内容及实验得到的结果 (23)四、实验结论 (26)五、思考题解答 (26)实验四、IIR数字滤波器设计和应用 (27)一、实验目的 (27)二、实验原理 (27)三、实验内容和结果 (27)四、思考题解答 (33)实验五、FIR数字滤波器设计和应用 (34)一、实验目的 (34)二、FIR数字滤波器的设计基本原理 (34)三、实验内容和实验结果 (37)四、思考题解答 (40)实验一、基本信号的产生和时频域抽样实验一、实验目的1、学习使用matlab产生基本信号波形、实现信号的基本运算；2、熟悉连续信号经理想采样前后的频谱变化关系,加深对时域采样定理的理解;3、加深理解频谱离散化过程中的数学概念和物理概念，掌握频域抽样定理的基本内容。
二、实验内容及所得图表1、用Matlab产生以下序列的样本，并显示其波形：(a):()(0.9)cos(0.2/3),020nx n n nππ=+≤≤(b):)20()5()(---=nununx(c):)*5.0exp()(n nx-=(d):(e):(f):)()sin()(t u tAetx taΩ=-α2 设（a）：求其傅里叶变换；（b）：用频率对进行采样，求出采样所得离散时间信号的傅里 Fs=5000Hz对进行采样，求出采样所得离散时间信号的叶变换；再用频率Fs=1000Hz傅里叶变换；（c）：分别针对（b）中采样所得离散时间信号和，重建出对应的连续时间信号和，并分别与原连续时间信号进行比较；根据抽样定理（即Nyquist 定理）的知识，说明采样频率对信号重建的影响。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
f (t ) f1 (t ) * f 2 (t ) f1 ( ) f 2 (t )d
（二）线性时不变（LTI）系统的单位冲激响应 给定一个连续时间 LTI 系统，在系统的初始条件为零时，用单位冲激信号(t)作用 于系统，此时系统的响应信号称为系统的单位冲激响应（Unit impulse response） ，一般 用 h(t)来表示。需要强调的是，系统的单位冲激响应是在激励信号为 (t)时的零状态响 应（Zero-state response） 。 系统的单位冲激响应是一个非常重要的概念，如果已知一个系统的单位冲激响 应，那么，该系统对任意输入信号的响应信号都可以求得。 （三）卷积的意义 对于 LTI 系统，根据系统的线性和时不变性以及信号可以分解成单位冲激函数可 得，任意 LTI 系统可以完全由它的单位冲激响应 h(t)来确定，系统的输入信号 x(t)和输 出信号 y(t)之间的关系可以用卷积运算来描述，即：
(t )dt 1 (t ) 0,
t0
1.1
1, u (t ) 0,
t0 t0
1.2
下面给出产生单位冲激信号和单位阶跃信号的两个函数，供参考。 产生单位冲激信号的程序为：
function x=delta(t1,t2,t0) dt=0.01; t=t1:dt:t2; n=length(t); %信号的时间间隔 %信号时间样本点向量 %时间样本点向量长度
的说明检查程序的方法。 用上述向量对连续信号进行表示后，就可以用 plot 命令绘制信号的时域波形。命 令如下：
plot(t,x) title(‘x(t)=Sa(t)’) xlabel(‘t’)
axis([-10,10,-0.2,1.2])
绘制的信号波形如图一所示，当把 t 改为：t =-10:0.5:10;则可得到图二。因为 plot 命令将点与点之间用直线连接，当点与点之间距离很小时，绘出的图形就成了光滑的 曲线。但图二在 t=0 时，曲线是间断的。
信号与系统实验指导书
实验二
一、实验目的：
连续时间信号的时域分析
1、掌握用 MATLAB 对连续信号进行基本运算和时域变换的方法； 2、掌握用 MATLAB 表示周期信号的方法。 二、实验原理： （一）连续时间信号的基本运算和时域变换 1、加法：x1(t)+x2(t) 信号的加法运算为对应位置处量值的相加，在 MATLAB 中可用运算符“+”实 现，但要求参与运算的两信号向量的长度必须相等。如果长度不等或者长度相等但采 样位置不同，则不能直接应用该运算符，此时需要先给定参数使序列具有相同的位置 向量和长度。下面给出 sigadd 函数实现任意两信号的加法运算。
x=zeros(1,n); x(1,(t0-t1)/dt+1)=1/dt; stairs(t,x);
%各样本点函数值赋值为零 %在时间 t=-t0 处，样本点赋值为 1/dt
产生单位阶跃信号的程序为：
function x=ut(t) x=(t>0); %t>0 时 x 为 1，否则为 0
在调用该函数表示信号时，需要先定义向量 t，如 t=-1:0.01:3。 对于其他常用信号，可以直接调用 MATLAB 中的内部函数进行定义，例如：正 弦信号：sin( )，余弦信号：cos( )，指数信号：exp( )，符号函数：sign( )；周期方波： square( )，周期锯齿波：sawtooth( )。复指数信号是时间 t 的复函数，需要用模和相角或 实部和虚部来表示复指数信号随时间变化的规律，对应的函数分别为 abs( )，angle( )， real( )，imag( )。 各函数的参数及定义方法可参考 MATLAB 的帮助文件。 三、实验内容： 1、参考示例程序，绘制信号 e 2t cos 3t[u (t ) u (t 3)] 的图形，t 取-1 到 4，步长 值设为 0.01。 2、产生一个指数为 [0.1 ( / 4) * i ] * t 的复指数函数，绘出函数的实部、虚部、幅 度和相位的波形，t 取 0 到 20，步长值设为 0.1。 四、思考： 1、为什么图二中 t=0 处曲线是间断的，如何使其成为连续的曲线？
实验一
一、实验目的：
连续时间信号的时域表示
1、掌握连续时间信号的 MATLAB 表示方法； 2、掌握用 MATLAB 描绘二维图像的方法。 二、实验原理： 信号是消息的载体，是消息的一种表现形式。信号可以是多种多样的，通常表现 为随时间变化的某些物理量，一般用 x(t)或 x(n)来表示。信号按照自变量的取值是否连 续可分为连续时间信号和离散时间信号。 连续时间信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了 有若干不连续点以外，信号都有确定的值与之对应。严格来说，MATLAB 并不能处理 连续信号，而是用等时间间隔点的样值来近似地表示连续信号。当取样时间间隔足够 小时，这些离散的样值就能较好地近似连续信号。在 MATLAB 中通常用向量来表示 连续时间信号，向量需要与时间变量相对应。 对于连续时间信号 x(t)，可用 x、t 两个行向量来表示。其中向量 t 是形如 t = t1:p:t2 的 MATLAB 命令定义的时间范围向量，t1 为信号起始时间，t2 为终止时间，p 为时间 间隔。向量 x 为连续信号 x(t)在向量 t 所定义的时间点上的样值。如产生连续信号
其中 x1 和 x2 为参与加法运算的两信号，t1 和 t2 分别为 x1 和 x2 的时间向量。 2、乘法：x1(t)·x2(t) 序列的乘法运算为对应位置处量值的相乘，在 MATLAB 中由数组运算符“.*”实 现，也受到“+”运算符同样的限制。 3、时移：y(t) = x(t - t0) 其中，t0 为位移量，当 t0>0 时，y(t)为 x(t)右移 t0 时刻之后的结果，当 t0<0 时， y(t)为 x(t)左移|t0|时刻之后的结果。 在 MATLAB 中，时移运算与数学上习惯表达方法完全相同。例：
function [y,t] = sigadd(x1,t1,x2,t2) t = min(min(t1),min(t2)):max(max(t1),max(t2)); y1 = zeros(1,length(t)); y2 = y1; y1(find((t>=min(t1))&(t<=max(t1))==1))=x1; y2(find((t>=min(t2))&(t<=max(t2))==1))=x2; y = y1+y2; %结果的时间向量 %初始化 %在公共区间定义 y1 %在公共区间定义 y2
x(t ) Sa (t ) sin(t ) 可用如下命令实现： t
t =-10:1.5:10; x=sin(t)./ t;
在命令窗口（Command Window）中可得到程序执行的结果即 x、t 的具体值。注 意：在 MATLAB 程序调试过程中，有时程序执行不出结果或虽然出结果但存在一些 问题，MATLAB 都会在 Command 窗口中给出错误说明，掌握利用 Command 窗口中
y (t ) x( )h(t )d
由于系统的单位冲激响应是零状态响应，故按照上式求得的系统响应也是零状态 响应。它是描述连续时间系统输入输出关系的一个重要表达式。 （四）函数说明 利用 MATLAB 的内部函数 conv( 其语法为：y = )可以很容易地完成两个信号的卷积积分运算。
plot(t1,x1) 注：函数左移时，t0<0，即 t-|t0|；函数右移时，t0>0。
4、反折：y(t) = x(-t) 在 MATLAB 中有多种方法可以实现信号的反折运算。 （1）修改绘图函数 plot(t,x)中的时间变量 t，即用-t 代替原来的 t。 （2）直接利用原信号与其反折信号的数学关系式来实现。这种方法最符合信号反 折运算的实际意义。 （3）使用 MATLAB 内部函数 fliplr( )来实现信号的反褶运算。其用法为： y=
Байду номын сангаас
若信号的自变量的范围和 t 的范围相同，则不能用上述方法，如将 x = exp(-0.5*t) 进行左移得到 x1 = exp(-0.5*(t+2))后，还需要对 x1 的时间变量重新定义。由于函数的 平移可看作是函数时间向量的平移，而对应的样值不变，当函数左移时，所有时间序 号都减小|t0|个单位，反之，则增加 t0 个单位。因此可用如下方式实现： t1=t+t0; x1=x;
clear; t = -5:0.01:5; x = exp(-0.5*t).*ut(t); x1 = exp(-0.5*(t+2)).*ut(t+2); subplot(211)
信号与系统实验指导书
plot(t,x) title ('原信号 x(t)') subplot (212) plot (t,x1) title (' x(t)左移 2') xlabel (' t (sec)')
该图添加横、纵坐标说明和标题；有时在一个程序中需要将几个图形绘制在一个窗 口，利用 subplot(m,n,k)函数可以将当前窗口分成 m 行 n 列个子窗口，并在第 k 个子窗 口绘图，窗口的排列顺序为从左至右，从上至下分别为 1,2,…m*n。 以上为几个常用绘图函数的基本用法，有关各函数的其他参数可参考 MATLAB 的帮助文件。 单位冲激信号(t) 和单位阶跃信号 u(t)是信号处理中两个最基本的信号，它们的定 义如下：
fliplr(x)，其中 x 为原信号，而 y 则为 x 的时域反折。需要说明的是，函数 fliplr()对信 号作时域反折，仅仅将信号中各个元素的次序作了一个反转，这种反转处理是独立于 时间变量 t 的。因此，还需要对时间变量 t 进行反折，即 t= -fliplr(x)。 5、展缩：y(t) = x(at) 其中 a 为任意常数。根据 a 的不同取值，尺度变换对信号 x(t)具有不同的影响。当 a > 1 时，y(t) = x(at)，y(t)是将 x(t)在时间轴上压缩得到；当 0 < a < 1 时，y(t) = x(at)， y(t)是将 x(t)在时间轴上扩展得到。在 MATLAB 中，按照数学上的常规方法即能实 现。 （二）周期信号 周期信号是一类非常重要的信号。给定一个信号 x(t)，如果满足 x(t) = x(t+kT)，则 该信号叫做周期信号。其中，k 为任意整数，T 为常数，通常称为信号的基本周期或最