北京理工大学信号与系统实验报告5-连续时间系统的复频域分析
实验5连续时间系统的复频域分析
（综合型实验）
一、实验目的
1）掌握拉普拉斯变换及其反变换的定义并掌握MATLAB 实现方法。
2）学习和掌握连续时间系统函数的定义及复频域分析方法。
3）掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换
连续时间信号x(t)的拉普拉斯变换定义为
(s)(t)e st X x dt
+∞
--∞
=
?
（1）
拉普拉斯反变换为1
(t)(s)e 2j st j x X ds
j σσπ+∞
-∞
=?
（2）
MATLAB 中相应函数如下：
(F)
L laplace = 符号表达式F 拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。
(F,t)L laplace =用t 替换结果中的变量s 。
()
F ilaplace L =以s 为变量的符号表达式L 的拉氏反变换，返回时间变量为t 的结果表达式。
(,)
F ilaplace L x =用x 替换结果中的变量t 。
的连续时间系统，其系统函数为s 的有理函数
110
110
...(s)...M M M M N N N N b s b s b H a s a s a ----+++=
+++ （7）
3.连续时间系统的零极点分析
系统的零点指使式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s 平面上，零点用O 表示，极点用?表示，这样得到的图形为零极点分布图。可以通过利用MATLAB 中的求多项式根的roots 函数来实现对（7）分子分母根的求解，调用格式如下：
r=roots(c)，c 为多项式的系数向量，返回值r 为多项式的根向量。
求取零极点以及绘制系统函数的零极点分布图可以采用pzmap 函数，调用格式如下： pzmap(sys)绘出由系统模型sys 描述的系统的零极点分布图。
[p,z]=pzmap(sys)这种调用方式返回极点与零点，不绘出零极点分布图。
还有两个专用函数tf2zp 和zp2tf 可实现系统的传递函数模型和零极点增益模型的转换。调用格
式如下：
[z,p,k]=tf2zp(b,a)
[b,a]=tf2zp(z,p,k)
研究系统函数的零极点分布不仅可以了解系统冲激响应的形式，还可以了解系统的频率特性以及判断系统的稳定性。
1）零极点分布与冲激响应的关系
系统的极点位置决定着系统冲激响应h(t)的波形，冲激响应的幅值是由系统函数的零点和极点共同确定的，系统的零点位置只影响冲激响应的幅度和相位，不影响波形。
2）零极点分布与系统频率响应的关系
系统函数的零极点分布不仅决定了系统函数H(s)，也决定了系统的频率响应()
H ，根据系统的零极点分布情况，可以由几何矢量法分析系统的频率响应。
3）零极点分布与系统稳定性的关系
稳定性是系统的固有性质，与激励信号无关，由
于系统函数(s)
H包含了系统的所固有的性质，因而可以根据系统函数的零极点分布判断系统的
稳定性。因果系统稳定的充要条件是(s)
H的全部极点位于s的左半平面。
三．实验内容
（1）已知系统的冲激响应(t)u(t)u(t 2)h =--，输入信号
(t)u(t)
x =，试采用复频域的方法求解系统的响应，
编写MATLAB 程序实现。 代码：
%DFTfifth_2_1.m syms t
h=heaviside(t)-heaviside(t-2); x=heaviside(t); H=laplace(h); X=laplace(x); Y=H*X; y=ilaplace(Y)
>> DFTfifth_2_1 y =
t - heaviside(t - 2)*(t - 2)
所以系统的响应为y(t)=t-(t-2)*u(t-2)
（2）已知因果连续时间系统的系统函数分别如下： 1）3
21
(s)221
H s s s =+++
2）5
4321
(s)23332
H s
s s s s =+-+++
试采用MATLAB 绘出其零极点分布图，求解系统的冲激响应h(t)和频率响应()H ，并判断系统是否稳定。 1） >> b=[1]; >> a=[1 2 2 1]; >> sys=tf(b,a); >> [p,z]=pzmap(sys) p =
-1.0000 -0.5000 + 0.8660i -0.5000 - 0.8660i z =
Empty matrix:
0-by-1
>> pzmap(sys)
所有极点都位于s 平面的左半平面，所以系统是
稳定的。 >> syms s
Real Axis (seconds -1)
I m a g i n a r y A x i s (s e c o n d s -1)
>> Hs=1/(s^3+2*s^2+2*s+1); >> h=ilaplace(Hs) h =
exp(-t) - exp(-t/2)*(cos((3^(1/2)*t)/2) - (3^(1/2)*sin((3^(1/2)*t)/2))/3)
所以系统的冲激响应为2
(t)[e t t t
h e -
-=-
绘制时域和频域的曲线：
b=[1]; a=[1 2 2 1]; sys=tf(b,a); subplot(311); impulse(sys); xlabel('t'); title('h(t)'); subplot(312); [H,w]=freqs(b,a);
plot(w,abs(H)); xlabel('w'); ylabel('Magnitude'); title('abs(H)'); subplot(313); plot(w,angle(H)); xlabel('w'); ylabel('phase'); title('phase(H)');
t (seconds)
A m p l i t u d e
w
M a g n i t u d e
abs(H)
w
p h a s e
phase(H)
2）
>> b=[1 0 1];
>> a=[1 2 -3 3 3 2]; >> sys=tf(b,a)
sys =
s^2 + 1 ------------------------------------- s^5 + 2 s^4 - 3 s^3 + 3 s^2 + 3 s + 2
Continuous-time transfer function.
>> [p,z]=pzmap(sys)
p =
-3.1704 0.9669 + 0.9540i 0.9669 - 0.9540i -0.3817 + 0.4430i -0.3817 - 0.4430i z =
0 + 1.0000i 0 - 1.0000i
>> pzmap(sys)
由于s 平面有半平面有极点，所以是不稳定系
统。
绘制冲激响应和频域响应的图形 方法同上一题 图形如下：
（3）已知连续时间系统函数的极点位置分别如下所示（设系统无零点）：
Real Axis (seconds -1
)
I m a g i n a r y A x i s (s e c o n d s -1)
w M a g n i t u d e
abs(H)
w
p h a s e
phase(H)
28
h(t)
t (seconds)
A m p l i t u d e
分别绘制以下六种不同情况下，系统函数的零极点分布图，并绘制相应冲激响应的时域波形，观察并分析系统函数极点位置对冲激响应时域特性的影响。
1）p=0
>> b=[1];
>> a=[1 0]; >> sys=tf(b,a)
sys =
1
-
s Continuous-time transfer function. >> pzmap(sys)
1
(s)(t)u(t) H h
s
=?=
>> syms t
>> h=heaviside(t); >> ezplot(h,[-5 5]) >> title('h(t)')
Real Axis (seconds-1)
I
m
a
g
i
n
a
r
y
A
x
i
s
(
s
e
c
o
n
d
s
-
1
)
2） p=-2 >> b=[1]; >> a=[1 2]; >> sys=tf(b,a) sys =
1 ----- s +
2 Continuous-time transfer function.
>> pzmap(b,a)
21
(s)(t)e *(t)2
t H h u s -=
?=+
>> syms t >>
h=exp(-2*t)*heavisi
t
h(t)
P ole-Zero Map
Real Axis (seconds -1)
I m a g i n a r y A x i s (s e c o n d s -1)
>> ezplot(h)
3)p=2 >> b=[1]; >> a=[1 -2]; >> sys=tf(b,a) sys =
1 ----- s - 2
Continuous-time transfer function.
>> pzmap(b,a)
21
(s)(t)*(t)2
t H h e u s =
?=-
>> syms t >>
h=exp(2*t)*heavisid
t
P ole-Zero Map
Real Axis (seconds -1)
I m a g i n a r y A x i s (s e c o n d s -1)
>> ezplot(h)
4) 1
22,2p
j p j
==-
>> b=[1]; >> a=[1 0 4]; >> sys=tf(b,a) sys =
1 ------- s^
2 + 4
Continuous-time transfer function.
>> pzmap(b,a)
2
11
(s)(t)sin(2t)*(t)42
H h u s =
?=+
>> syms t >>
h=(1/2)*sin(2*t)*he
t
P ole-Zero Map
Real Axis (seconds -1
)
I m a g i n a r y A x i s (s e c o n d s -1)
>> ezplot(h)
5) 1
214,14p
j p j
=-+=--
>> b=[1]; >> a=[1 2 17]; >> sys=tf(b,a)
sys =
1 -------------- s^
2 + 2 s + 17
Continuous-time transfer function. >> pzmap(b,a)
2
11(s)(t)sin(2)(t)2172
t
H h e t u s s -=
?=++
>> syms t >>
h=(1/2)*exp(-t)*sin
t
Real Axis (seconds -1
)
I m a g i n a r y A x i s (s e c o n d s -1)
>> ezplot(h)
6) 1
214,14p
j p j
=+=-
>> b=[1]; >> a=[1 -2 17]; >> sys=tf(b,a) sys =
1 -------------- s^
2 - 2 s + 17
Continuous-time
transfer function.
>> pzmap(b,a)
2
11(s)(t)sin(2)(t)2172
t
H h e t u s s =
?=-+
>> syms t
t
Real Axis (seconds -1
)
I m a g i n a r y A x i s (s e c o n d s -1)
>>
h=exp(t)*sin(2*t)*h eaviside(t)/2; >> ezplot(h)
极点在左半平面时呈衰减趋势，在左半平面坐标轴上时呈指数衰减，在非坐标轴位置上时成衰减振荡；在右半平面时成增加趋势，在右半平面坐标轴上时呈增加趋势，在非坐标轴上时呈增幅振荡；在纵轴上时，在非原点时呈等幅振荡，在原点时为单位阶跃响应。
（4）已知三个连续时间系统的系统函数，极点相同，零点不同，试用MATLAB 分别绘制系统的零极点分布图及相应冲激响应的时域波形，观察并分析系统函数零点位置对冲激响应时域特性的影响。 1）2
1
(s)217
H s
s =++
>> b=[1]; >> a=[1 2 17];
>> sys=tf(b,a)
t
(sin(2 t) exp(t) heaviside(t))/2
sys =
1
-------------- s^2 + 2 s + 17
Continuous-time transfer function.
>> pzmap(b,a)
2
11(s)(t)sin(4t)u(t)2174
t
H h e s s -=
?=++
>> syms t >>
h=(1/4)*exp(-t)*sin (4*t)*heaviside(t); >> ezplot(h) >> axis tight
2）2
8(s)217
s H s
s +=++
>> b=[1 8]; >> a=[1 2 17];
Real Axis (seconds -1)
I m a g i n a r y A x i s (s e c o n d s -1
)
t
(sin(4 t) exp(-t) heaviside(t))/4
>> sys=tf(b,a) sys =
s + 8
-------------- s^2 + 2 s + 17
Continuous-time transfer function.
>> pzmap(b,a)
>> syms s
>> H=(s+8)/(s^2+2*s+17); >> h=ilaplace(H) h =
exp(-t)*(cos(4*t) + (7*sin(4*t))/4)
2
87(s)(t)[cos(4t)sin(4t)]u(t)2174
t
s H h e s s -+=
?=+++
>> clear
Real Axis (seconds -1)
I m a g i n a r y A x i s (s e c o n d s -1)
>> syms t
>>h=exp(-t)*(cos(4*t)+(7*sin(4*t))/4)*heaviside(t); >> ezplot(h) >> axis tight
3）2
8(s)217
s H s
s -=++
>> b=[1 -8]; >> a=[1 2 17]; >> sys=tf(b,a) sys =
s - 8
--------------
s^2 + 2 s + 17
Continuous-time transfer function.
>> pzmap(b,a)
t
exp(-t) heaviside(t) (cos(4 t) + (7 sin(4 t))/4)
Real Axis (seconds -1)
I m a g i n a r y A x i s (s e c o n d s -1)
实验五 信号与系统的复频域分析 王靖 08通信 12号 实验目的 （1）掌握利用MA TLAB 进行连续时间信号与系统的复频域分析。 （2）掌握利用MA TLAB 进行离散系统的复频域分析。 实验环境 安装MATLAB7.0以上版本的计算机 实验内容 1. 利用help 命令了解以下命令的基本用法 residue ，roots ，pzmap ，cart2pol ，residuez ，tf2zp ，zplane 2. 部分分式展开的MATLAB 实现 用部分分式展开法求X(s)的反变换。 2321 ()452s X s s s s +=+++ 步骤一：建立新的m 文件，保存并命名为program1.m 。 步骤二：输入以下命令，理解每条命令的含义。 %program1，部分分式展开法求反变换 [10 1];[1452];[,,](,) n u m d en r p k resid u e n u m d en === 步骤三：保存程序并运行，记录得到的结果。 如右图所示 步骤四：由得到的结果可以直接获得X(s)展开表示式 25 4 2 ()21(1)X s s s s =-++++： 步骤五：由此可得到X(s)反变换的原函数，记录。 X(t)=(5exp(-2*t)-4exp(-t)+2texp(-t)) 思考：将其转换成极坐标形式，应该如何使用cart2pol 命令？离散系统的部分分式展开，如何使用命 令residuez ，得到的结果如何利用？ 将笛卡尔坐标转化为极坐标用 [angle,mag]=cart2pol(real(r),imag(r)) [r,p,k] = residuez(nun,,den)
北京理工大学信号与系统实验报告5-连续时间系统的复频域分析
实验5连续时间系统的复频域分析 （综合型实验） 一、实验目的 1）掌握拉普拉斯变换及其反变换的定义并掌握MATLAB 实现方法。 2）学习和掌握连续时间系统函数的定义及复频域分析方法。 3）掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换 连续时间信号x(t)的拉普拉斯变换定义为 (s)(t)e st X x dt +∞ --∞ = ? （1） 拉普拉斯反变换为1 (t)(s)e 2j st j x X ds j σσπ+∞ -∞ =? （2） MATLAB 中相应函数如下： (F) L laplace = 符号表达式F 拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。 (F,t)L laplace =用t 替换结果中的变量s 。 () F ilaplace L =以s 为变量的符号表达式L 的拉氏反变换，返回时间变量为t 的结果表达式。 (,) F ilaplace L x =用x 替换结果中的变量t 。
的连续时间系统，其系统函数为s 的有理函数 110 110 ...(s)...M M M M N N N N b s b s b H a s a s a ----+++= +++ （7） 3.连续时间系统的零极点分析 系统的零点指使式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s 平面上，零点用O 表示，极点用?表示，这样得到的图形为零极点分布图。可以通过利用MATLAB 中的求多项式根的roots 函数来实现对（7）分子分母根的求解，调用格式如下： r=roots(c)，c 为多项式的系数向量，返回值r 为多项式的根向量。 求取零极点以及绘制系统函数的零极点分布图可以采用pzmap 函数，调用格式如下： pzmap(sys)绘出由系统模型sys 描述的系统的零极点分布图。 [p,z]=pzmap(sys)这种调用方式返回极点与零点，不绘出零极点分布图。 还有两个专用函数tf2zp 和zp2tf 可实现系统的传递函数模型和零极点增益模型的转换。调用格
电子工程系 信号与系统课程实验报告 2011-----2012学年第一学期 专业: 电子信息工程技术班级: 学号 : 姓名: 指导教师: 实常用连续时间信号的实现
一、实验目的 （1）了解连续时间信号的特点； （2）掌握连续时间信号表示的向量法和符号法； （3）熟悉MATLAB Plot函数等的应用。 二、实验原理 1、信号的定义 信号是随时间变化的物理量。信号的本质是时间的函数。 2、信号的描述 1）时域法 时域法是将信号表示成时间的函数f(t)来对信号进行描述的方法。信号的时间特性指的是信号的波形出现的先后，持续时间的长短，随时间变化的快慢和大小，周期的长短等。 2）频域（变换域）法 频域法是通过正交变换，将信号表示成其他变量的函数来对信号进行描述的方法。一般常用的是傅立叶变换。信号的频域特性包括频带的宽窄、频谱的分布等。 信号的频域特性与时域特性之间有着密切的关系。 3、信号的分类 按照特性的不同，信号有着不同的分类方法。 （1）确定性信号：可以用一个确定的时间函数来表示的信号。 随机信号：不可以用一个确定的时间函数来表示，只能用统计特性加以描述的信号。 （2）连续信号：除若干不连续的时间点外，每个时间点在t上都有对应的数值信号。离散信号：只在某些不连续的点上有数值，其他时间点上信号没有定义的信号。 （3）周期信号：存在T，使得等式f(t+T)=f(t)对于任意时间t都成立的信号。非周期信号：不存在使得等式f(t+T)=f(t)对于任意时间t都成立的信号。 绝对的周期信号是不存在的，一般只要在很长时间内慢走周期性就可以了。 （4）能量信号：总能量有限的信号。 功率信号:平均功率有限切非零的信号。 （5）奇信号：满足等式f(t)=--f(--t)的信号。偶信号：满足等式f(t)=f(--t)的信号。 三、涉及的MATLAB函数 1、plot函数 功能：在X轴和Y轴方向都按线性比例绘制二维图形。 调用格式： Plot（x,y）:绘出相x对y的函数线性图。 Plot（x1,y1,x2,y2,…..）:会出多组x对y的线性曲线图。 2、ezplot函数 功能：绘制符号函数在一定范围内的二维图形。简易绘制函数曲线。 调用格式： Ezplot (fun):在[-2π，2π]区间内绘制函数。 Ezplot (fun，［min，max］)：在［min,max］区间内绘函数。 Ezplot (funx,funy):定义同一曲面的函数，默认的区间是[0, 2π]。】 3、sym函数 功能：定义信号为符号的变量。 调用格式：sym(fun):fun为所要定义的表达式。 4、subplot函数
1.周期信号的频谱 周期信号在满足一定条件时，可以分解为无数三角信号或指数之和。这就是周期信号的傅里叶级数展开。在三角形式傅里叶级数中，各谐波分量的形式为()1cos n n A n t ω?+；在指数形式傅里叶级数中，分量的形式必定为1j n t n F e ω 与1-j -n t n F e ω 成对出现。为了把周期信号所具有的各 次谐波分量以及各谐波分量的特征（如模、相角等）形象地表示出来，通常直接画出各次谐波的组成情况，因而它属于信号的频域描述。 以周期矩形脉冲信号为lifenxi 周期信号频谱的特点。周期矩形信号在一个周期（-T/2,T/2）内的时域表达式为 ,2 0,>2 ()A t T t f t ττ ≤?=?? （2-6） 其傅里叶复数系数为 12 n n A F Sa T ωττ?? = ??? （2-7） 由于傅里叶复系数为实数，因而各谐波分量的相位为零（n F 为正）或为π±（n F 为负），因此不需要分别画出幅度频谱n F 与相位频谱n φ。可以直接画出傅里叶系数n F 的分布图。 如图2.4.1所示。该图显示了周期性矩形脉冲信号()T f t 频谱的一些性质，实际上那个也是周期性信号频谱的普遍特性： ① 离散状频谱。即谱线只画出现在1ω的整数倍频率上，两条谱线的间隔为1ω（等于2π/t ）。 ② 谱线宽度的包络线按采样函数()1/2a S n ωτ的规律变化。如图2.4.2所示。但1ω 为 2π τ 时，即( )2m π ωτ =（m=1，2，……）时，包络线经过零点。在两相邻 零点之间，包络线有极值点，极值的大小分别为-0.212()2A T τ，
实验5 连续时间系统的复频域分析 一、实验目的 1.掌握拉普拉斯变换及其反变换的定义，并掌握MATLAB 实现方法。 2.学习和掌握连续时间系统系统函数的定义及复频域分析方法。 3.掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换 连续时间信号)(t x 的拉普拉斯变换定义为 )1.....(..........)()(dt e t x s X st ? +∞ ∞ --= 拉普拉斯反变换定义为 )2....(..........)(21)(ds e s X j t x j j st ?∞ +∞ -=σσπ 在MATLAB 中，可以采用符号数学工具箱的laplace 函数和ilaplace 函数进行拉氏变换和反拉氏变换。 L=laplace(F)符号表达式F 的拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。 L=laplace(F,t)用t 替换结果中的变量s 。 F=ilaplace(L)以s 为变量的符号表达式L 的拉氏反变换，返回时间变量为t 的结果表达式。 F=ilaplace(L,x)用x 替换结果中的变量t 。 除了上述ilaplace 函数，还可以采用部分分式法，求解拉普拉斯逆变换，具体原理如下： 当 X (s )为有理分式时，它可以表示为两个多项式之比： )3.(..........)()()(0 110 11a s a s a b s b s b s D s N s X N N N N M M M M +?+++?++==---- 式（3）可以用部分分式法展成一下形式 )4.....(.............)(2211N N p s r p s r p s r s X -++-+-= 通过查常用拉普拉斯变换对，可以由式（1-2）求得拉普拉斯逆变换。 利用 MATLAB 的residue 函数可以将 X (s )展成式（1-2）所示的部分分式展开式，该 函数的调用格式为：[r,p,k] = residue(b,a) 其中b 、a 为分子和分母多项式系数向量，r 、p 、k 分别为上述展开式中的部分分式系数、极点和直项多项式系数。 2.连续时间系统的系统函数
实验六信号与系统复频域分析 一、实验目的 1.学会用MATLAB进行部分分式展开； 2.学会用MATLAB分析LTI系统的特性； 3.学会用MATLAB进行Laplace正、反变换。 4.学会用MATLAB画离散系统零极点图； 5.学会用MATLAB分析离散系统的频率特性； 二、实验原理及内容 1．用MATLAB进行部分分式展开 用MATLAB函数residue可以得到复杂有理分式F(s)的部分分式展开式，其调用格式为 其中，num,den分别为F(s)的分子和分母多项式的系数向量，r为部分分式的系数，p为极点，k为F(s)中整式部分的系数，若F(s)为有理真分式，则k为零。 例6-1 用部分分式展开法求F(s)的反变换 解:其MATLAB程序为 format rat; num=[1,2]; den=[1,4,3,0]; [r,p]=residue(num,den) 程序中format rat是将结果数据以分数形式显示
F(s)可展开为 210.536()13 F s s s s --=++++ 所以，F(s)的反变换为 3211()()326t t f t e e u t --??=--???? 2．用MATLAB 分析LTI 系统的特性 系统函数H （s ）通常是一个有理分式，其分子和分母均为多项式。计算H （s ）的零极点可以应用MATLAB 中的roots 函数，求出分子和分母多项式的根，然后用plot 命令画图。 在MATLAB 中还有一种更简便的方法画系统函数H （s ）的零极点分布图，即用pzmap 函数画图。其调用格式为 pzmap(sys) sys 表示LTI 系统的模型，要借助tf 函数获得，其调用格式为 sys=tf(b,a) 式中，b 和a 分别为系统函数H （s ）的分子和分母多项式的系数向量。 如果已知系统函数H （s ），求系统的单位冲激响应h(t)和频 率响应H ω（j ）可以用以前介绍过的impulse 和freqs 函数。 例6-2 已知系统函数为 321221 s s s +++H(s)= 试画出其零极点分布图，求系统的单位冲激响应h(t)和频率响应H ω（j ），并判断系统是否稳定。 解：其MATLAB 程序如下： num=[1];
实验一：信号的时域分析 一、实验目的 1.观察常用信号的波形特点及产生方法 2.学会使用示波器对常用波形参数的测量 二、实验仪器 1.信号与系统试验箱一台（型号ZH5004） 2.40MHz双踪示波器一台 3.DDS信号源一台 三、实验原理 对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(t)信号、钟形信号、脉冲信号等。 1、信号：指数信号可表示为f(t)=Ke at。对于不同的a取值，其波形表现为不同的形式，如下图所示： 图1―1 指数信号 2、信号：其表达式为f（t）=Ksin（ωt＋θ），其信号的参数：振幅K、角频率ω、与初始相位θ。其波形如下图所示：
图1－2 正弦信号 3、指数衰减正弦信号：其表达式为其波形如下图： 图1－3 指数衰减正弦信号 4、Sa(t)信号：其表达式为：。Sa(t)是一个偶函数，t= ±π,±2π,…,±nπ时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：
图1－4 Sa（t）信号 5、钟形信号（高斯函数）：其表达式为：其信号如下图所示： 图1－5 钟形信号 6、脉冲信号：其表达式为f(t)=u(t)-u(t-T)，其中u(t)为单位阶跃函数。其信号如下图所示： 7、方波信号：信号为周期为T，前T/2期间信号为正电平信号，后T/2期间信号为负电平信号，其信号如下图所示 U(t)
周期信号的时域及其频域分析 姓名：张敏靓学号：1007433014 一、实验目的 1.掌握Multisim软件的应用及用虚拟仪器对周期信号的频谱测量 2.掌握选频电平表的使用，对信号发生器输出信号（方波、矩形波、 三角波等）频谱的测量 二、实验原理 周期信号的傅里叶级数分析法，可以把周期信号表示为三角傅里叶级数或指数傅里叶级数，其中周期信号满足。 1. 周期信号表示为三角傅里叶级数 2. 周期信号表示为指数傅里叶级数 其中， 周期矩形信号的频谱
三、实验内容 1.在Multisim上实现周期信号的时域、频域测量及分析 （1）绘制测量电路 （2）周期信号时域、频域（幅度频谱）的仿真测量 虚拟信号发生器分别设置如下参数： 周期方波信号：周期T=100μs，脉冲宽度τ=50μs，脉冲幅度 V P=5V； 周期矩形信号：周期T=100μs，脉冲宽度τ=20μs，脉冲幅度 V P=5V； 周期三角波信号：周期T=200μs，脉冲幅度V P=5V； 采用虚拟示波器及虚拟频谱仪分别测量上述信号的时域、频域波形并保存测试波形及数据。
2.周期信号时域、频域（幅度频谱）的测量 信号发生器、示波器、选频电平表的连线如上图所示。信号发生器的输出信号分别为周期分别信号、周期矩形信号、周期三角波信号，参数设置同仿真测量。采用示波器及选频电平表对信号发生器的输出信号分别测量，并将测量数据记录下表中。
四、实验总结 1.在周期矩形信号的实验中，信号频率减小，频谱减小；信号占空 比减小，频谱减小；幅度值减小，频谱减小。 2.未安装Origin绘图软件，Excel绘图未能达到理想效果。
实验八系统的复频域 分析
一、实验目的 1、掌握系统的复频域分析方法。 2、掌握测试系统的频率响应的方法。 二、预习内容 1、系统频响的方法。（见第四章波特图的介绍） 三、实验原理 1. N 阶系统系统的传递函数 用微分方程描述的N 阶系统为： 根据零状态响应（起始状态为零），则对其进行拉氏变换有： 则系统传递函数可表达为： 用差分方程描述的N 阶系统为： 根据零状态响应（起始状态为零），则对其进行拉氏变换有： 则系统传递函数可表达为： 2.根据系统传递函数的零极点图分析系统 零点：传递函数分子多项式的根。 极点：传递函数分母多项式的根。 根据零极点图的不同分布分析系统。 3.涉及到的Matlab 函数 (1)freqz 函数：实验六中出现过，可用来求单位圆上的有理z 变换的值。调用格式：同实验六 (2)zplane 函数：得到有理z 变换的零极点图。 调用格式：zplane(num,den)
其中，num和 den是按z ?1 的升幂排列的、z 变换分子分母多项式系数的行向量。 (3)roots 函数：求多项式的根。 调用格式：r=roots(c), c 为多项式系数向量；r 为根向量。 四、实验内容 1.系统零极点的求解 (1)求解系统和的零极点，验 证下面程序的运行结果，根据系统零极点图分析系统性质。 b=[1,0,-1]; a=[1,2,3,2]; zr=roots(b); pr=roots(a); plot(real(zr),imag(zr),'go',real(pr),imag(pr),'mx','markersize',12,'linewidth',2); grid; legend('零点','极点'); figure; zplane(b,a); (2)参考上述程序，绘制系统和 的零极点图，并分析系统性质。与用zplane 函数直接绘制系统零极点图（注:圆心的圆圈并非系统的零点）做比较。
实验二 连续时间信号的频域分析 一、实验目的 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法； 2、观察截短傅里叶级数而产生的“Gibbs 现象”，了解其特点以及产生的原因； 3、掌握连续时间傅里叶变换的分析方法及其物理意义； 4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质； 5、学习掌握利用Matlab 语言编写计算CTFS 、CTFT 和DTFT 的仿真程序，并能利用这些程序对一些典型信号进行频谱分析，验证CTFT 、DTFT 的若干重要性质。 基本要求：掌握并深刻理傅里叶变换的物理意义，掌握信号的傅里叶变换的计算方法，掌握利用Matlab 编程完成相关的傅里叶变换的计算。 二、原理说明 1、连续时间周期信号的傅里叶级数CTFS 分析 任何一个周期为T 1的正弦周期信号，只要满足狄利克利条件，就可以展开成傅里叶级数。 三角傅里叶级数为： ∑∞ =++=1 000)]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 。也可以反过来理解三角傅里叶级数：用无限多个正弦谐波分量可以合成一个任意的非正弦周期信号。 指数形式的傅里叶级数为：
实验四：连续系统的复频域分析 一、实验目的： 1、掌握连续与离散时间系统的正反复频域与Z域变换 2、掌握利用MATLAB进行零极点分析，进一步了解零极点对整个系统的影响 3、掌握simulink环境下系统建模与仿真以及系统求解。 二、实验内容： 1、已知某连续系统的系统函数为： （1）利用[r, p, k]=residue(num, den),求H（s）的极零点以及多项式系数; （2）画出系统的零极点分布图，判断系统得稳定性。 （3）求h(t)，判断系统得稳定性。 2、已知某离散系统的系统函数为：， （1）利用[r, p, k]=residuez(num, den)求H（z）的极零点以及多项式系数; （2）画出零极点分布图，判断系统得稳定性。 （3）求单位函数响应用impz(b, a)，判断系统是否稳定； 3、已知线性时不变微分方程 在Simulink环境下搭建起系统的仿真模型，并查看仿真结果曲线。（1）写出传递函数H（s），绘出系统模拟框图； （2）当f(t)分别为，，的零状态响应；且当与课本P81的结果进行比较（3）方程的初值为, ,求全响应； 4、已知某信号，n(t)为正态噪声干扰且服从N(0,0.22)分布，对此信号进行采样，采样间隔为0.001s，之后对此信号进行Botterworth低通滤波，从信号中过滤10HZ的输出信号，试对系统进行建模与仿真。 三、实验数据处理与结果分析： 第一题：题1_1：
>> num=[2,5]; den=[1,1,3,2]; [r,p,k]=residue(num,den) r = -0.5750 - 0.7979i -0.5750 + 0.7979i 1.1499 p =-0.1424 + 1.6661i -0.1424 - 1.6661i -0.7152 k =[]
学生实验报告 （理工类） 课程名称：信号与系统实验专业班级：电子信息（1）班学生学号：1005101058 学生姓名：严生生 所属院部：信息技术学院指导教师：杨婧 20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制
实验报告书写要求 实验报告原则上要求学生手写，要求书写工整。若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 （1）细致观察，及时、准确、如实记录。 （2）准确说明，层次清晰。 （3）尽量采用专用术语来说明事物。 （4）外文、符号、公式要准确，应使用统一规定的名词和符号。 （5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细，一律用红色笔批改。实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。 实验报告装订要求 实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。
实验项目名称：常用连续信号的表示实验学时： 1 同组学生姓名：实验地点： B402 实验日期：实验成绩： 批改教师：杨婧批改时间： 一、实验目的和要求 熟悉MATLAB软件，利用MATLAB软件，绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机，装有MATLAB7.0软件。 三、实验过程 1，绘制正弦信号f（t）=Asin（ωt+ψ），其中A=1，ω=2π, ψ=π/6; 2,绘制指数信号f（t）=Ae^at，其中A=1，a=-0.4； 3，绘制矩形脉冲信号，脉冲宽度为2； 4，绘制三角波脉冲信号，脉冲宽度为4；斜度为0.5； 5，对上题三角波脉冲信号进行尺度变换，分别得出f（2t），f（2-2t）； 6，绘制抽样函数Sa(t),t取值在-3π到+3π之间； 7，绘制周期矩形脉冲信号，参数自定； 8，绘制周期三角脉冲信号，参数自定； 1，打开MATLAB界面，建立新文件。 2，根据实验要求，编写程序。
信号与系统 实验报告 实验三周期信号的频谱分析 实验报告评分：_______ 实验三周期信号的频谱分析 实验目的： 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法； 2、观察截短傅里叶级数而产生的“Gibbs现象”，了解其特点以及产生的原因；
3、掌握各种典型的连续时间非周期信号的频谱特征。 实验内容： （1）Q3-1 编写程序Q3_1，绘制下面的信号的波形图： 其中，0 = 0.5π，要求将一个图形窗口分割成四个子图，分别绘制cos( 0t)、cos(3 0t)、cos(5 0t)和x(t) 的波形图，给图形加title，网格线和x坐标标签，并且程序能够接受从键盘输入的和式中的项数。 程序如下: clear,%Clear all variables close all,%Close all figure windows dt = 0.00001; %Specify the step of time variable t = -2:dt:4; %Specify the interval of time w0=0.5*pi; x1=cos(w0.*t); x2=cos(3*w0.*t); x3=cos(5*w0.*t); N=input('Type in the number of the harmonic components N='); x=0; for q=1:N; x=x+(sin(q*(pi/2)).*cos(q*w0*t))/q; end subplot(221) plot(t,x1)%Plot x1 axis([-2 4 -2 2]); grid on, title('signal cos(w0.*t)') subplot(222) plot(t,x2)%Plot x2 axis([-2 4 -2 2]); grid on, title('signal cos(3*w0.*t))') subplot(223) plot(t,x3)%Plot x3 axis([-2 4 -2 2])
实验三 常见信号的MATLAB 表示及运算 一、实验目的 1．熟悉常见信号的意义、特性及波形 2．学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法 二、实验原理 根据MATLAB 的数值计算功能和符号运算功能，在MA TLAB 中，信号有两种表示方法，一种是用向量来表示，另一种则是用符号运算的方法。在采用适当的MA TLAB 语句表示出信号后，就可以利用MA TLAB 中的绘图命令绘制出直观的信号波形了。 1.连续时间信号 从严格意义上讲，MATLAB 并不能处理连续信号。在MATLAB 中，是用连续信号在等时间间隔点上的样值来近似表示的，当取样时间间隔足够小时，这些离散的样值就能较好地近似出连续信号。在MATLAB 中连续信号可用向量或符号运算功能来表示。 ⑴ 向量表示法 对于连续时间信号()f t ，可以用两个行向量f 和t 来表示，其中向量t 是用形如12::t t p t =的命令定义的时间范围向量，其中，1t 为信号起始时间，2t 为终止时间，p 为时间间隔。向量f 为连续信号()f t 在向量t 所定义的时间点上的样值。 ⑵ 符号运算表示法 如果一个信号或函数可以用符号表达式来表示，那么我们就可以用前面介绍的符号函数专用绘图命令ezplot()等函数来绘出信号的波形。 ⑶ 常见信号的MATLAB 表示 单位阶跃信号 单位阶跃信号的定义为：10 ()0 t u t t >?=? 0); %定义函数体，即函数所执行指令
信号与系统 实验报告 实验五连续系统的复频域分析 实验五连续系统的复频域分析 一、实验目的 1. 深刻理解拉普拉斯变换、逆变换的定义，掌握用MATLAB实现拉普拉斯变换、逆变换的方法。 2会求几种基本信号的拉氏变换。 3 掌握用MATLAB绘制连续系统零、极点的方法。 4 求解系统函数H(s)。 二
1已知连续时间信号f（t）=sin（t）u（t）、求出该信号的拉普拉斯变换，并用MATLAB 绘制拉普拉斯变换的曲面图。 syms t; ft=sin(t)*heaviside(t); Fs=Laplace(ft); a=-0.5:0.08:0.5; b=-2:0.08:2; [a,b]=meshgrid(a,b); c=a+i*b; d=ones(size(a)); c=c.*c; c=c+d; c=1./c; c=abs(c); mesh(a,b,c); surf(a,b,c) axis([-0.5,0.5,-2,2,0,10]) colormap(hsv
) 2求[(1-e^(-at))]/t的拉氏变换。 syms t s a f1=(1-exp(-a*t))/t; F=laplace(f1,t,s) F = log(s+a)-log(s) 3求F（s）=-log（s）+ log（s+a）的拉氏逆变换syms t s a F =log(s+a)-log(s); f1=ilaplace(F,s,t) f1 = (1-exp(-a*t))/t
4已知某连续系统的系统函数为： H（s）=(s^2+3s+2)/(8s^4+2s^3+3s^2+5)试用MATLAB求出该系统的零极点，画出零极点分布图。 b=[1 3 2]; a=[8 2 3 0 5]; zs=roots(b); ps=roots(a); hold on plot(real(zs),imag(zs),'o'); plot(real(ps),imag(ps),'x'); grid axis([-2.5,1,-1,1]) 5已知H(s)=（s+1）/（s^2+s+1）,绘制阶跃响应图形，冲激响应图形，频率激响应图形。 syms t s H=(s+1)/(s^2+s+1); f1=ilaplace(H,s,t); f2=heaviside(t);
计算机与信息工程学院设计性实验报告 一、实验目的 1.掌握用matlab 分析系统时间响应的方法 2.掌握用matlab 分析系统频率响应的方法 3.掌握系统零、极点分布与系统稳定性关系 二、实验原理 1.系统函数H(s) 系统函数:系统零状态响应的拉氏变换与激励的拉氏变换之比. H(s)=R(s)/E(s) 在matlab 中可采用多种方法描述系统,本文采用传递函数(系统函数)描述法. 在matlab 中, 传递函数描述法是通过传递函数分子和分母关于s 降幂排列的多项式系数来表示的.例如,某系统传递函数如下 )1(8 .03.11 )(2+++=s s s s H 则可用如下二个向量num 和den 来表示: num=[1,1]；den=[1,1.3,0.8] 2.用matlab 分析系统时间响应 1)脉冲响应 y=impulse(num,den,T) T:为等间隔的时间向量,指明要计算响应的时间点. 2)阶跃响应 y=setp(num,den,T) T 同上. 3)对任意输入的响应 y=lsim(num,den,U,T) U:任意输入信号. T 同上. 3.用matlab 分析系统频率响应特性 频响特性: 系统在正弦激励下稳态响应随信号频率变化的特性. ()()() ()j s j H j H s H j e φωω ωω=== |H(j ω)|:幅频响应特性. ?(ω):相频响应特性(或相移特性).
Matlab 求系统频响特性函数freqs 的调用格式: h=freqs(num,den,ω) ω:为等间隔的角频率向量,指明要计算响应的频率点. 4.系统零、极点分布与系统稳定性关系 系统函数H(s)集中表现了系统的性能,研究H(s)在S 平面中极点分布的位置,可很方面地判断系统稳定性. 1) 稳定系统: H(s)全部极点落于S 左半平面(不包括虚轴),则可以满足 0)]([lim =∞ →t h t 系统是稳定的. 2)不稳定系统: H(s)极点落于S 右半平面,或在虚轴上具有二阶以上极点,则在足够长时间后,h(t)仍继续增长, 系统是不稳定的. 3)临界稳定系统: H(s)极点落于S 平面虚轴上,且只有一阶,则在足够长时间后,h(t)趋于一个非零数值或形成一个等幅振荡. 系统函数H(s)的零、极点可用matlab 的多项式求根函数roots()求得. 极点:p=roots(den) 零点:z=roots(num) 根据p 和z 用plot()命令即可画出系统零、极点分布图,进而分析判断系统稳定性. 三、实验内容 设()(1)(2) s H s s p s p = -- 设①p1=-2,p2=-30; ②p1=-2,p2=3 1. 针对极点参数①②, 画出系统零、极点分布图, 判断该系统稳定性. 2. 针对极点参数①②,绘出系统的脉冲响应曲线,并观察t →∞时, 脉冲响应变化趋势. 3. 针对极点参数①, 绘出系统的频响曲线. 四、实验要求 1．预习实验原理； 2．对实验内容编写程序(M 文件),上机运行； 3．绘出实验内容的各相应曲线或图。 五、实验设备 1.装MATLAB 软件的计算机 1台
理工大学信号与系统实验报告连续时间系统的 复频域分析 Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】
实验5连续时间系统的复频域分析 （综合型实验） 一、实验目的 1）掌握拉普拉斯变换及其反变换的定义并掌握MATLAB 实现方法。 2）学习和掌握连续时间系统函数的定义及复频域分析方法。 3）掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换 连续时间信号x(t)的拉普拉斯变换定义为(s)(t)e st X x dt +∞ --∞ =? （1） 拉普拉斯反变换为1 (t)(s)e 2j st j x X ds j σσπ+∞ - ∞ = ? （2） MATLAB 中相应函数如下： (F)L laplace = 符号表达式F 拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。 (F,t)L laplace =用t 替换结果中的变量s 。 ()F ilaplace L =以s 为变量的符号表达式L 的拉氏反变换，返回时间变量 为t 的结果表达式。 (,)F ilaplace L x =用x 替换结果中的变量t 。 拉氏变换还可采用部分分式法，当(s)X 为有理分式时，它可以表示为两个多项式之比： 110 1 10 ...(s)(s)(s)...M M M M N N N N b s b s b N X D a s a s a ----+++==+++ （3）
上式可以采用部分分式法展成以下形式 1212(s)...N N r r r X s p s p s p = +++--- （4） 再通过查找常用拉氏变换对易得反变换。 利用residue 函数可将X(s)展成（4）式形式，调用格式为： [r,p,k]residue(b,a)=其中b 、a 为分子和分母多项式系数向量，r 、p 、k 分 别为上述展开式中的部分分式系数、极点和直项多项式系数。 2.连续时间系统的系统函数 连续时间系统的系统函数是指系统单位冲激响应的拉氏变换 (s)(t)e st H h dt +∞ --∞ = ? （5） 连续时间系统的系统函数还可以由系统输入与输出信号的拉氏变换之比得到。 (s)(s)/X(s)H Y = （6） 单位冲激响应(t)h 反映了系统的固有性质，而(s)H 从复频域反映了系统的固有性质。由（6）描述的连续时间系统，其系统函数为s 的有理函数 110 1 10 ...(s)...M M M M N N N N b s b s b H a s a s a ----+++=+++ （7） 3.连续时间系统的零极点分析 系统的零点指使式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s 平面上，零点用O 表示，极点用?表示，这样得到的图形为零极点分布图。可以通过利用MATLAB 中的求多项式根的roots 函数来实现对（7）分子分母根的求解，调用格式如下：
实验一:信号得时域分析 一、实验目得 1.观察常用信号得波形特点及产生方法 2.学会使用示波器对常用波形参数得测量 二、实验仪器 1.信号与系统试验箱一台（型号ＺH５0０4) 2.40MHｚ双踪示波器一台 3.ＤDS信号源一台 三、实验原理 对于一个系统特性得研究,其中重要得一个方面就是研究它得输入输出关系，即在一特定得输入信号下,系统对应得输出响应信号.因而对信号得研究就是对系统研究得出发点,就是对系统特性观察得基本手段与方法.在本实验中,将对常用信号与特性进行分析、研究。 信号可以表示为一个或多个变量得函数，在这里仅对一维信号进行研究,自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(ｔ)信号、钟形信号、脉冲信号等。 1、信号:指数信号可表示为ｆ(t）＝Ke at。对于不同得ａ取值，其波形表现为不同得形式,如下图所示: 图１―1 指数信号 2、信号:其表达式为f(t）=Ｋsin(ωt+θ),其信号得参数:振幅K、角频率ω、与初始相位θ。其波形如下图所示:
图1-2 正弦信号 3、指数衰减正弦信号:其表达式为其波形如下图: 图１-３指数衰减正弦信号 ４、Sａ(t）信号:其表达式为:。Sa（t）就是一个偶函数，t= ±π,±２π,…,±ｎπ时,函数值为零。该函数在很多应用场合具有独特得运用。其信号如下图所示: 图1－4 Ｓa(t)信号 5、钟形信号(高斯函数):其表达式为：其信号如下图所示：
图1-5 钟形信号 6、脉冲信号:其表达式为ｆ(t)=u(t)-ｕ（t—T),其中u（t)为单位阶跃函数。其信号如下图所示: f(ｔ） ? ……??…… 0 t 图1-6脉冲信号 7、方波信号:信号为周期为T,前T/2期间信号为正电平信号,后T/２期间信号为负电平信号，其信号如下图所示 U(t) ………… ?0?t 图1－７方波信号 四、实验内容及主要步骤 下列实验中信号产生器得工作模式为11 1、指数信号观察 通过信号选择键1，设置Ａ组输出为指数信号(此时信号输出指示灯为0００000）。用示波器测量“信号Ａ组”得输出信号。 输出波形为：
a=-0.5:0.08:0.5; b=-2:0.08:2; [a,b]=meshgrid(a,b); d=ones(size(a)); c=a+i*b; c=c.*c; c=c+d; c=1./c c=abs(c); surf(a,b,c); axis=([-0.5,0.5,-2,2,0.15]); title('单边正弦信号拉氏变换图'); colormap(hsv); 2. a=0:0.5:5; b=-20:0.1:20; [a,b]=meshgrid(a,b); c=a+i*b; c=(1-exp(-2*c)./c); c=abs(c) mesh(a,b,c); sufr(a,b,c); view(-10,20); axis([-0.5:-20,20,0.2]); title('拉氏变换S域像函数'); w=-20:0.1:20; Fw=(2*sin(w).*exp(i*w)); plot(w,abs(Fw)) title=('傅里叶变换'); xlabel('频率w'); 3. a=[1,2,-3,2,1]; b=[1,4]; sjdt(a,b); a=[1,5,16,30]; b=[5,20,25,0]; sjdt(a,b); 4. a=[8,2,3,15]; b=[1,3,2]; [p,q]=sjdt(a,b); 5.
a=[1,0]; b=[1]; impulse(b,a) a=[1,2]; b=[1]; impulse(b,a) a=[1,-2]; b=[1]; impulse(b,a) a=[1,1,16,25]; b=[1]; impulse(b,a) a=[1,0,16]; b=[1]; impulse(b,a) a=[1,-1,16,25]; b=[1]; impulse(b,a) 6. q=[0,0]; p=[-100,100]; f1=0; f2=1000; k=0.1; splxy(f1,f2,k,p,q) q=[0,0]; f1=0; f2=1000; k=0.1; p=[-500,-1000]; splxy(f1,f2,k,p,q) q=[0,0]; f1=0; f2=1000; k=0.1; p=[-2000,-4000]; splxy(f1,f2,k,p,q) 7. 8
实验一连续时间信号的时域基本运算 一、实验目的 （1）掌握连续时间信号时域运算的基本方法； （2）掌握相关函数的调用格式及作用； （3）掌握连续信号的基本运算。 二、实验原理 信号的基本运算包括信号的相加（减）和相乘(除).信号的时域变换包括信号的平移、翻转、倒相以及尺度变换。 （1）加减： f(t)=f1(t)±f2(t)（2）乘: f(t)=f1(t)×f2(t) （3）延时或平移：f(t)→(t-t0) t0>0时右移；t0<0时左移 （4）翻转： f(t)→f(-t) （5）尺度变换：f(t)→ f(at) |a|>1时尺度缩小；|a|<1时尺度放大；a<0时，尺度翻转。 （6）标量乘法：f(t)→af(t) （7）倒相： f(t)→-f(t) （8）微分： f(t)→df(t)/dt （9）积分： f(t)→∫t -∞f(t)d(t) 三、涉及的MATLAB函数及其运算 1、stepfun函数 功能：产生一个阶跃信号。 调用格式： Stepfun(t,t 0)其中，t是时间区间，在该区间内阶跃信号一定会产生；t 是信号 发生从0到1跳跃的时刻。 2、diff函数 调用格式： diff (f) : 求函数f对预设独立变数的一次微分值。 diff (f, ’t’) : 求函数f对独立变数t的一次微分值。 3、int函数 调用格式: Int（f）: 函数f对预设独立变数的积分值。 Int(f,’t’): 函数f对独立变数t的积分值。 4、heaviside函数 Heaviside（t）:产生没有移位的阶跃信号。 Heaviside（t-k）:产生向右平移K单位的阶跃信号。 四、实验内容与方法1、验证性实验 （1）移位 实现连续信号的移位，即f(t-t0),或者f(t+t0);常数t0>0。 MATLAB程序 clear all t=0:0.0001:2 y=sin(2*pi*(t)); y1=sin(2*pi*(t-0.2)); plot(t,y,'-',t,y1,'--'); ylabel('f(t)');xlabel('t');title('信号的移位') 运行结果：