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

实验一 时域离散信号与系统变换域分析

一、实验目的

1．了解时域离散信号的产生及基本运算实现。 2．掌握离散时间傅里叶变换实现及系统分析方法。 3. 熟悉离散时间傅里叶变换性质。 4. 掌握系统Z 域分析方法。

5. 培养学生运用软件分析、处理数字信号的能力。 二、实验内容 1. 序列的基本运算

1.1 产生余弦信号)04.0cos()(n n x π=及带噪信号)(

2.0)04.0cos()(n w n n y +=π 0<=n<=50（噪声采用randn 函数）

1.2 已知12)(1-=n n x 51≤≤n ，22)(2-=n n x 62≤≤n ,求两个序列的和、乘积、序列x1的移位序列（右移2位）,序列x2的翻褶序列，画出原序列及运算结果图。

2. 序列的傅里叶变换

2.1 已知序列)()5.0()(n u n x n =。试求它的傅里叶变换，并且画出其幅度、相角、实部和虚部的波形，并分析其含有的频率分量主要位于高频区还是低频区。

2.2 令||1000)(t a e t x -=，求其傅立叶变换)(Ωj X a 。分别用kHz f s 1=和kHz f s 5=对其进行采样，求出离散时间傅立叶变换)(ωj e X ，写出程序，并画出相应频谱，分析结果的不同及原因。

3. 序列的傅里叶变换性质分析

3.1 已知序列n j e n x )9.0()(3/π=，100≤≤n ，求其傅里叶变换，并讨论其傅里叶变换的周期性和对称性。

3.2 已知序列n n x )9.0()(-=，55≤≤-n ，求其傅里叶变换，并讨论其傅里叶变换的周期性和对称性。

为了方便，考虑在两个周期，例如[ππ2,2-]中2M+1个均匀频率点上计算FT ，并且观察其周期性和对称性。为此给出function 文件如下，求解FT 变换：

function [X,w]=ft1(x,n,k) w=(pi/abs(max(k)/2))*k

X=x*(exp(-j*pi/abs(max(k)/2))).^(n'*k)

3.3 编写程序验证序列傅里叶变换频移性质，时域卷积定理（时域卷积后的频域特性）。（所需信号自行选择）

4. 时域差分方程的求解

4.1求解差分方程y(n)＋a1y(n-1)＋a2y(n-2)=b0x(n)＋b1x(n-1)的零状态响应和全响应。已知X(n)为单位取样序列，y(-1)=1,y(-2)=2,a1=0.5,a2=0.06,b0=2,b1＝3。

5. 离散系统的Z 域分析

5.1 利用系统函数)(z H 分析系统的稳定性。假设系统函数如下式：

5147

.13418.217.198.33)

3)(9()(2

34-++--+=

z z z z z z z H ，试判断系统是否稳定。 5.2 已知线性时不变系统的系统函数2

11

12.08.013.01.0)(-----+=z

z z z H ，编写程序求其单位取样响应，频率响应及系统零极点，并画出相应图形。

6. 创新训练拓展内容

6.1 利用Matlab 自带的录音功能，或利用Goldwave 等音频编辑软件，对语音或其他音频信号进行采集并保存为*.wav 文件。

要求：（1）采用不同的采样频率（2000Hz ，4000Hz ，8000Hz ，16000Hz 等）。

（2）对采集得到的信号进行播放。

（3）分析在不同采样频率下得到的信号有何不同。

6.2 设定一个连续时间信号，进行抽样和恢复，要求分析不同采样频率对恢复结果的影响，给出实验程序及各关键步骤图形结果。

6.3 设计内容：

设计一个离散系统，给定系统函数或差分方程，设定激励及初始条件。要求： （1）绘制系统函数零极点图，判断稳定性； （2）求单位序列响应h （n ）；

（3）求系统零输入响应及零状态响应，要求零状态响应采样三种方法求解（卷积的方法、迭代解法、差分方程求解函数方法），激励自定；

（4）分析系统频响特性，画出频响函数幅频曲线和相频曲线。 三、试验要求

第一部分：验证试验内容

根据给定的试验内容，部分试验给出了参考程序段，见下面各段程序。请基于Matlab 环境进行验证试验。

第二部分：编程试验内容

对于给定的试验内容中，没有参考程序段的部分，进行编程，并给出试验结果，进行相应的分析。

第三部分：创新训练拓展内容

此部分内容，要求根据个人能力，进行选作。

1.序列的基本运算

%1.单位取样序列 x(n)=delta(n-n0) 要求n1<=n0<=n2

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

n=[n1:n2]; x=[(n-n0)==0]; == 是逻辑判断

%2.单位阶跃序列 x(n)=u(n-n0) 要求n1<=n0<=n2

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

n=[n1:n2]; x=[(n-n0)>=0];

%3.信号加 y(n)=x1(n)+x2(n)

%find函数：找出非零元素的索引号

%x1:第一个序列的值，n1:序列x1的索引号

%x2:第二个序列的值，n2:序列x2的索引号

function[y,n]=sigadd(x1,n1,x2,n2)

n=min(min(n1),min(n2)):max(max(n1),max(n2));

y1=zeros(1,length(n)); y2=y1;

y1(find((n>=min(n1))&(n<=max(n1))==1))=x1;

y2(find((n>=min(n2))&(n<=max(n2))==1))=x2;

y=y1+y2;

%4.信号乘 y(n)=x1(n)*x2(n)

function[y,n]=sigmult(x1,n1,x2,n2)

n=min(min(n1),min(n2)):max(max(n1),max(n2));

y1=zeros(1,length(n)); y2=y1;

y1(find((n>=min(n1))&(n<=max(n1))==1))=x1;

y2(find((n>=min(n2))&(n<=max(n2))==1))=x2;

y=y1.*y2;

%5.移位 y(n)=x(n-n0)

function[y,n]=sigshift(x,m,n0)

n=m+n0; y=x;

%6.翻褶 y(n)=x(-n)

function[y,n]=sigfold(x,n)

y=fliplr(x); n=-fliplr(n);

2.序列的傅里叶变换

%7. 求序列)

n

x n

的傅里叶变换

u

(n

(

)5.0(

)

w=[0:1:500]*pi/500

X=exp(j*w)./(exp(j*w)-0.5*ones(1,501)) magX=abs(X)

angX=angle(X)

realX=real(X)

imagX=imag(X)

subplot(2,2,1)

plot(w/pi,magX)

grid

xlabel('frequency in pi units')

title('Magnitude Part')

ylabel('Magnitude')

subplot(2,2,3)

plot(w/pi,angX)

grid

xlabel('frequency in pi units')

title('Angle Part')

ylabel('Radians')

subplot(2,2,2)

plot(w/pi,realX)

grid

xlabel('frequency in pi units')

title('Real Part')

ylabel('Real')

subplot(2,2,4)

plot(w/pi,imagX)

grid

xlabel('frequency in pi units')

title('Imaginary Part')

ylabel('Imaginary')

程序执行结果：

0.5

1

0.5

11.5

2

frequency in pi units

Magnitude Part

M a g n i t u d e

0.5

1

-0.8

-0.6-0.4-0.2

0frequency in pi units

Angle Part

R a d i a n s

0.5

1

0.5

11.5

2frequency in pi units Real Part

R e a l

0.5

1

-0.8

-0.6-0.4-0.2

0frequency in pi units

Imaginary Part

I m a g i n a r y

%8 令||1000)(t a e t x -=，绘制其傅立叶变换)(Ωj X a 。用不同频率对其进行采样，分别画出

)(ωj e X 。

Dt=0.00005; %步长为0.00005s t=-0.005:Dt:0.005;

xa=exp(-1000*abs(t)); %取时间从-0.005s 到0.005s 这段模拟信号 Wmax=2*pi*2000; %信号最高频率为2π*2000 K=500; %频域正半轴取500个点进行计算 k=0:1:K;

W=k*Wmax/K; %K

W k max *=Ω 求模拟角频率

Xa=xa*exp(-j*t'*W)*Dt; %计算连续时间傅立叶变换（利用矩阵运算实现） Xa=real(Xa); %取实部

W=[-fliplr(W),W(2:501)]; %将角频率范围扩展为从-到+ Xa=[fliplr(Xa),Xa(2:501)]; subplot(2,2,1);

plot(t*1000,xa); %画出模拟信号，横坐标为时间（毫秒），纵坐标为幅度 xlabel('time(millisecond)');ylabel('xa(t)'); title('anolog signal'); subplot(2,2,2);

plot(W/(2*pi*1000),Xa*1000); %画出连续时间傅立叶变换 xlabel('frequency(kHZ)'); %横坐标为频率（kHz ） ylabel('xa(jw)'); %纵坐标为幅度 title('FT');

%下面为采样频率5kHz 时的程序

T=0.0002; %采样间隔为s f s

0002.01

=

n=-25:1:25;

x=exp(-1000*abs(n*T)); %离散时间信号 K=500;k=0:1:K;w=pi*k/K; %w 为数字频率

X=x*exp(-j*n'*w); %计算离散时间傅立叶变换（序列的傅立叶变换） X=real(X);

w=[-fliplr(w),w(2:K+1)]; X=[fliplr(X),X(2:K+1)]; subplot(2,2,3);

stem(n*T*1000,x); %画出采样信号（离散时间信号） xlabel('time(millisecond)'); ylabel('x1(n)');

title('discrete signal'); subplot(2,2,4);

plot(w/pi,X); %画出离散时间傅立叶变换 xlabel('frequency(radian)'); %横坐标为弧度 ylabel('x1(jw)');title('DTFT');

3. 序列的傅里叶变换性质分析

%9 已知序列n j e n x )9.0()(3/π=，100≤≤n ，求其傅里叶变换，并讨论其傅里叶变换的周期性和对称性。 n=0:10

x=(0.9*exp(j*pi/3)).^n k=-200:200

[X,w]=ft1(x,n,k) magX=abs(X) angX=angle(X) subplot(2,1,1) plot(w/pi,magX) grid

xlabel('frequency in pi units') ylabel('/X/')

title('Magnitude Part') subplot(2,1,2)

plot(w/pi,angX/pi) grid

xlabel('frequency in pi units') ylabel('Radians/pi') title('Angle Part')

-2.5

-2

-1.5

-1

-0.500.51

1.5

2

0246

frequency in pi units

/X /

-2.5

-2-1.5-1

-0.500.51 1.52

-0.50

0.5

frequency in pi units

R a d i a n s /p i

Angle Part

由图可见，序列n j e n x )9.0()(3/π=的傅里叶变换对w 是周期的，但不是共轭对称的。 %10、已知序列n n x )9.0()(-=，55≤≤-n ，求其傅里叶变换，并讨论其傅里叶变换的周期性和对称性。

n=-5:5

x=(-0.9).^n k=-200:200

[X,w]=ft1(x,n,k) magX=abs(X) angX=angle(X) subplot(2,1,1) plot(w/pi,magX) grid

xlabel('frequency in pi units') ylabel('/X/')

title('Magnitude Part') subplot(2,1,2)

plot(w/pi,angX/pi) grid

xlabel('frequency in pi units') ylabel('Radians/pi') title('Angle Part')

-2.5

-2-1.5-1

-0.500.51 1.52

0510

frequency in pi units

/X /

-2.5

-2

-1.5

-1

-0.500.51

1.5

2

-1-0.500.5

1frequency in pi units

R a d i a n s /p i

Angle Part

由图可见，序列n n x )9.0()(-=的傅里叶变换对w 是周期的，是共轭对称的。

4. 时域差分方程的求解

采用filter 函数实现线性常系数差分方程的递推求解，函数调用格式如下：

● yn=filter(B,A,xn) 计算输入信号xn 的零状态响应yn

● yn=filter(B,A,xn,xi) 计算输入信号xn 的全响应yn ，xi 为等效初始条件的输入序列

● xi=filtic(B,A,ys,xs) 由初始条件计算xi 的函数

4.1求解差分方程y(n)＋a1y(n-1)＋a2y(n-2)=b0x(n)＋b1x(n-1)的零状态响应和全响应。已知X(n)为单位取样序列，y(-1)=1,y(-2)=2,a1=0.5,a2=0.06,b0=2,b1＝3。

程序：

xn=[1 zeros(1,20)] B=[2,3] A=[1,0.5,0.06] ys=[1,2]

xi=filtic(B,A,ys) yn1=filter(B,A,xn)

yn2=filter(B,A,xn,xi) subplot(2,1,1) n1=0:length(yn1)-1 stem(n1,yn1,'.') axis([0,21,-3,3]) subplot(2,1,2) n2=0:length(yn2)-1 stem(n2,yn2,'.')

4.2 结果图形

上图为零状态响应、下图为全响应。

5. 离散系统的Z 域分析

%11 利用系统函数)(z H 分析系统的稳定性。假设系统函数如下式：

2

4

6

8

10

12

14

16

18

20

2

4

6

8

10

12

14

16

18

20

5147

.13418.217.198.33)

3)(9()(234-++--+=

z z z z z z z H ，试判断系统是否稳定。

解：

%调用roots 函数求极点， 并判断系统的稳定性

A=［3， －3.98， 1.17， 2.3418， －1.5147］；

%H(z)的分母多项式系数

p=roots(A) %求H(z)的极点

pm=abs(p)； %求H(z)的极点的模

if max(pm)<1 disp(′系统因果稳定′)， else ， disp(′系统不因果稳定′)， end

程序运行结果如下：

极点： －0.7486 0.6996-0.7129i 0.6996+0.7129i 0.6760

pm = 0.7486 0.9988 0.9988 0.6760

由极点分布判断系统因果稳定。

四、编程练习题

1. 已知序列)()5.0()(n u n x n -=。试求它的傅里叶变换，并且画出其幅度、相角、实部和虚部的波形，并分析其含有的频率分量主要位于高频区还是低频区。

2. 下面四个二阶网络的系统函数具有一样的极点分布：

1）211

29425.06.113.01)(---+--=z z z z H 2）2

1139425.06.118.01)(---+--=z z z z H 3）21139425.06.118.01)(---+--=z z z z H 4）2

12

149425.06.118.06.11)(----+-+-=z z z z z H

请分析研究零点分布对于单位脉冲响应的影响。 要求：

（1） 分别画出各系统的零、 极点分布图；

（2） 分别求出各系统的单位脉冲响应，并画出其波形；

（3） 分析零点分布对于单位脉冲响应的影响。

信号与系统实验

《信号与系统及MATLAB实现》实验指导书

前言 长期以来，《信号与系统》课程一直采用单一理论教学方式，同学们依靠做习题来巩固和理解教学内容，虽然手工演算训练了计算能力和思维方法，但是由于本课程数学公式推导较多，概念抽象，常需画各种波形，作题时难免花费很多时间，现在，我们给同学们介绍一种国际上公认的优秀科技应用软件MA TLAB，借助它我们可以在电脑上轻松地完成许多习题的演算和波形的绘制。 MA TLAB的功能非常强大，我们此处仅用到它的一部分，在后续课程中我们还会用到它，在未来地科学研究和工程设计中有可能继续用它，所以有兴趣的同学，可以对MA TLAB 再多了解一些。 MA TLAB究竟有那些特点呢？ 1．高效的数值计算和符号计算功能，使我们从繁杂的数学运算分析中解脱出来； 2．完备的图形处理功能，实现计算结果和编程的可视化； 3．友好的用户界面及接近数学表达式的自然化语言，易于学习和掌握； 4．功能丰富的应用工具箱，为我们提供了大量方便实用的处理工具； MA TLAB的这些特点，深受大家欢迎，由于个人电脑地普及，目前许多学校已将它做为本科生必须掌握的一种软件。正是基于这些背景，我们编写了这本《信号与系统及MA TLAB实现》指导书，内容包括信号的MA TLAB表示、基本运算、系统的时域分析、频域分析、S域分析、状态变量分析等。通过这些练习，同学们在学习《信号与系统》的同时，掌握MA TLAB的基本应用，学会应用MA TLAB的数值计算和符号计算功能，摆脱烦琐的数学运算，从而更注重于信号与系统的基本分析方法和应用的理解与思考，将课程的重点、难点及部分习题用MA TLAB进行形象、直观的可视化计算机模拟与仿真实现，加深对信号与系统的基本原理、方法及应用的理解，为学习后续课程打好基础。另外同学们在进行实验时，最好事先预习一些MA TLAB的有关知识，以便更好地完成实验，同时实验中也可利用MA TLAB的help命令了解具体语句以及指令的使用方法。

离散系统的时域及变换域分析

实验1 离散系统的时域及变换域分析 一、实验目的： 1.加深对离散系统的差分方程、单位抽样响应和卷积分析方法的理解。 2.加深对离散系统的频率响应分析和零、极点分布的概念理解。 二、实验原理： 1.时域 离散系统 其输入、输出关系可用以下差分方程描述： ∑∑==-=-M m m N k n m n x b k n y a )()( 输入信号分解为冲激信号， ∑∞ -∞ =-= m m n m x n x )()()(δ 系统单位抽样序列h （n ）， 则系统响应为如下的卷积计算式： ∑∞ -∞ =-= *=m m n h m x n h n x n y )()()()()( 当0 0≠a N k a k ,...2,1,0==时，h(n)是有限长度的（n ：[0，M]），称系统为FIR 系统；反之，称系统为IIR 系统。 在MATLAB 中，可以用函数y=filter(b,a,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积。 2.变换域 离散系统的时域方程为 ∑∑==-=-M m m N k n m n x b k n y a )()(

其变换域分析方法如下： X(z)H(z) Y(z) )()()()()(=?-= *=∑∞ -∞ =m m n h m x n h n x n y 系统函数为 N N M M z a z a a z b z b b z X z Y z H ----++++++= =......)()()(110110 分解因式 ∏∏∑∑=-=-=-=---== N k k M m m N k k k M m m m z d z c K z a z b z H 1 1 11 ) 1() 1()( ， 其中 m c 和 k d 称为零、极点。 在MATLAB 中，可以用函数[z ，p ，K]=tf2zp （num ，den ）求得有理分式形式的系统函数的零、极点，用函数zplane （z ，p ）绘出零、极点分布图；也可以用函数zplane （num ，den ）直接绘出有理分式形式的系统函数的零、极点分布图。使用h=freqz(num,den,w)函数可求系统的频率响应，w 是频率的计算点，如w=0:pi/255:pi, h 是复数，abs(h)为幅度响应，angle(h)为相位响应。另外，在MATLAB 中，可以用函数 [r ，p ，k]=residuez （num ，den ）完成部分分式展开计算；可以用函数sos=zp2sos （z ，p ，K ）完成将高阶系统分解为2阶系统的串联。 三 、实验内容 1.时域 （1.）编制程序求解下列系统的单位抽样响应，并绘出其图形。 )1()()2(125.0)1(75.0)(--=-+-+n x n x n y n y n y 解 用MATLAB 计算程序如下： N=15; n=0:N-1; b=[1,-1]; a=[1,0.75,0.125]; x=[n==0]; y=filter(b,a,x); subplot(3,2,1); stem(n,y,'.'); axis([0,N,-1,2]); ylabel('y(n)');

信号与系统分析实验报告连续系统的时域分析

北京联合大学 实验报告 课程（项目）名称： 信号与系统分析 学 院： 自动化学院 专 业： 信息处理与智能技术 班 级： 0910030204 学 号：2009100302440 姓 名： 韩禹辉 成 绩： 2011年 5 月 21 日 实验二 连续系统的时域分析 冲激响应与阶跃响应实验 一、实验目的 1.观察典型二阶电路的阶跃响应与冲激响应的波形和相关参数，并研究参数变化对响应状态的影响. 2.掌握系统阶跃响应与冲激响应的观测方法. 3.理解系统阶跃响应与冲激响应的关系. 二、实验设备 PC 机一台，TD-SAS 系列教学实验系统一套. 三、实验原理 本实验是观察典型的二阶系统的阶跃响应和冲激响应的三种不同状态.二阶系统的微分方程通式为： 2()2()()()n n y t ay t y t f t ωω'''++= 其特征根为： 1,2a λ=-对于不同的a 和n ω值，特征根四种不同的情况，如表2-1-1所示，分别对应两个不等实根、两个相等实根、共轭复根和共轭虚根.相应的冲激响应和阶跃响应波形如图2-1-1所示. 表2-1-1 二阶系统的冲激响应和阶跃响应

图2-1-1二阶系统的冲激响应和阶跃响应 本实验电路采用由运放组成的典型二阶电路，如图2-1-2所示，它与RLC 串联电路构成二阶系统完成如图2-1-3所示的功能.实验中通过调节器Rp 便可以使系统处于不同的状态. 图2-1-2 由运放构成的二阶电路 图2-1-3 RLC 二阶电路 通过电路图可以得到该系统的微分方程为： 从公式可以得到：

由上式得到系统响应的三种状态： （1）当n a ω>时，即Rp>4K Ω时，系统有两个不等实根，处于过阻尼状态； （2）当n a ω=时，即Rp=4K Ω时，系统有两个相等实根处于临界阻尼状态； （3）当n a ω<时，即Rp<4K Ω时，系统有一对共轭复根，处于欠阻尼状态. 四、实验步骤 本实验在阶跃与冲激响应单元完成. 1．阶跃响应观察 （1）使信号发生器输出幅值2V 、频率为1Hz 、占空比为50%的脉冲信号，其中每个高电平作为一次阶跃输入.将脉冲信号接入IN 端. （2）用示波器同时测量IN 和OUT 两端，记录当电位器Rp 值分别为1.5K 、4K 和8K 时OUT 端的波形. 使用万用表测量电位器阻值时，先关闭实验箱电源开关，将短路块N 断开，这样电位器就从电路中断开，并且测量时应当注意表笔的正负端应和测量点的正负端一致.然后再打开实验箱电源开关，测量完后将短路块闭合，使电位器重新接入电路. （3）分别保存Rp 值在上述取值时的阶跃响应波形，并加以比较看是否满足图2-1-1（b ）所述.

实验6离散时间系统的z域分析

实验6 离散时间系统的z 域分析 一、实验目的 1.掌握z 变换及其反变换的定义，并掌握MATLAB 实现方法。 2.学习和掌握离散时间系统系统函数的定义及z 域分析方法。 3.掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理 1. Z 变换 序列x(n)的z 变换定义为 ()()n n X z x n z +∞ -=-∞ = ∑ Z 反变换定义为 1 1 ()()2n r x n X z z dz j π-= ? 在MATLAB 中，可以采用符号数学工具箱的ztrans 函数和iztrans 函数计算z 变换和z 反变换： Z=ztrans(F) 求符号表达式F 的z 变换。 F=ilaplace(Z) 求符号表达式Z 的z 反变换。 2.离散时间系统的系统函数 离散时间系统的系统函数H(z)定义为单位抽样响应h(n)的z 变换 ()()n n H z h n z +∞ -=-∞ = ∑ 此外，连续时间系统的系统函数还可以由系统输入和输出信号的z 变换之比得到 ()()/()H z Y z X z =

由上式描述的离散时间系统的系统函数可以表示为 101101()M M N N b b z b z H z a a z a z ----+++= +++…… 3.离散时间系统的零极点分析 离散时间系统的零点和极点分别指使系统函数分子多项式和分母多项式为零的点。在MATLAB 中可以通过函数roots 来求系统函数分子多项式和分母多项式的根，从而得到系统的零极点。 此外，还可以利用MATLAB 的zplane 函数来求解和绘制离散系统的零极点分布图，zplane 函数调用格式为： zplane(b,a) b,a 为系统函数的分子、分母多项式的系数向量（行向量）。 zplane(z,p) z,p 为零极点序列（列向量）。 系统函数是描述系统的重要物理量，研究系统函数的零极点分布不仅可以了解系统单位抽样响应的变化，还可以了解系统的频率特性响应以及判断系统的稳定性： ①系统函数的极点位置决定了系统单位抽样响应h(n)的波形，系统函数零点位置只影响冲激响应的幅度和相位，不影响波形。 ②系统的频率响应取决于系统的零极点，根据系统的零极点分布情况，可以通过向量分析系统的频率响应。 ③因果的离散时间系统稳定的充要条件是H(z)的全部极点都位于单位圆内。 三、实验内容 (1)已知因果离散时间系统的系统函数分别为： ①23221()0.50.0050.3 z z H z z z z ++=--+

信号与系统分析实验信号的频谱分析

实验三信号的频谱分析 1方波信号的分解与合成实验 1实验目的 1. 了解方波的傅立叶级数展开和频谱特性。 2. 掌握方波信号在时域上进行分解与合成的方法。 3. 掌握方波谐波分量的幅值和相位对信号合成的影响。 2 实验设备 PC机一台，TD－SAS系列教学实验系统一套。 3 实验原理及内容 1. 信号的傅立叶级数展开与频谱分析 信号的时域特性和频域特性是对信号的两种不同的描述方式。对于一个时域的周期信号f(t)，只要满足狄利克莱条件，就可以将其展开成傅立叶级数： 如果将式中同频率项合并，可以写成如下形式： 从式中可以看出，信号f(t)是由直流分量和许多余弦（或正弦）分量组成。其中第一项A0/2是常数项，它是周期信号中所包含的直流分量；式中第二项A1cos(Ωt+φ1)称为基波，它的角频率与原周期信号相同，A1是基波振幅，φ1是基波初相角；式中第三项A2cos(Ωt+φ2)称为二次谐波，它的频率是基波的二倍，A2是基波振幅，φ2是基波初相角。依此类推，还有三次、四次等高次谐波分量。 2. 方波信号的频谱 将方波信号展开成傅立叶级数为： n=1，3，5…

此公式说明，方波信号中只含有一、三、五等奇次谐波分量，并且其各奇次谐波分量的幅值逐渐减小，初相角为零。图3-1-1为一个周期方波信号的组成情况，由图可见，当它包含的分量越多时，波形越接近于原来的方波信号，还可以看出频率较低的谐波分量振幅较大，它们组成方波的主体，而频率较高的谐波分量振幅较小，它们主要影响波形的细节。 （a）基波（b）基波＋三次谐波 （c）基波＋三次谐波＋五次谐波 （d）基波＋三次谐波＋五次谐波＋七次谐波 （e）基波＋三次谐波＋五次谐波＋七次谐波＋九次谐波 图3-1-1方波的合成 3. 方波信号的分解 方波信号的分解的基本工作原理是采用多个带通滤波器，把它们的中心频率分别调到被测信号的各个频率分量上，当被测信号同时加到多路滤波器上，中心频率与信号所包含的某次谐波分量频率一致的滤波器便有输出。在被测信号发生的实际时间内可以同时测得信号所包含的各频率分量。本实验便是采用此方法，实验中共有5路滤波器，分别对应方波的一、 三、五、七、九次分量。 4. 信号的合成 本实验将分解出的1路基波分量和4路谐波分量通过一个加法器，合成为原输入的方波信号，信号合成电路图如图3-1-2所示。 图3-1-2

离散信号与系统时域分析

目录 第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) 绘制结果图形并对仿真结果进行详细的分析。

实验一 时域离散信号与系统变换域分析(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 函数。 序列傅里叶变换的性质：

离散系统的变换域分析

实验2 离散系统的变换域分析 一、实验目的 1、熟悉对离散系统的频率响应分析方法； 2、加深对零、极点分布的概念理解。 二、实验原理 离散系统的时域方程为 其变换域分析方法如下： 频域: 系统的频率响应为: Z域: 系统的转移函数为: 分解因式: ， 其中 和 称为零、极点。 三、预习要求 1. 在MATLAB中，熟悉函数tf2zp、zplane、freqz、residuez、 zp2sos的使用，其中：[z，p，K]=tf2zp（num，den）求得有 理分式形式的系统转移函数的零、极点；zplane（z，p）绘 制零、极点分布图；h=freqz(num,den,w)求系统的单位频率 响应；[r，p，k]=residuez（num，den）完成部分分式展开 计算；sos=zp2sos（z，p，K）完成将高阶系统分解为2阶系 统的串联。 2. 阅读文中的范例，学习频率分析法在MATLAB中的实现； 3. 编程实现系统参数输入，绘出幅度频率响应和相位响应曲线

和零、极点分布图。 四、实验内容 求系统 的零、极点和幅度频率响应和相位响应。 五、范例 求下列直接型系统函数的零、极点，并将它转换成二阶节形式 解:用MATLAB计算程序如下： num=[1 -0.1 -0.3 -0.3 -0.2]; den=[1 0.1 0.2 0.2 0.5]; [z,p,k]=tf2zp(num,den); m=abs(p); disp('零点');disp(z); disp('极点');disp(p); disp('增益系数');disp(k); sos=zp2sos(z,p,k); disp('二阶节');disp(real(sos)); zplane(num,den) 输入到“num”和“den”的分别为分子和分母多项式的系数。计算求得零、极点增益系数和二阶节的系数： 零点: 0.9615 -0.5730 -0.1443 + 0.5850i -0.1443 - 0.5850i 极点: 0.5276+0.6997i 0.5276-0.6997i -0.5776+0.5635i -0.5776-0.5635i 增益系数: 1 二阶节: 1.0000 -0.3885 -0.5509 1.0000 1.15520 0.6511

信号与系统MATLAB实验

2016-2017学年第一学期 信号与系统实验报告 班级： 姓名： 学号： 成绩： 指导教师：

实验一常见信号的MATLAB 表示及运算 一．实验目的 1．熟悉常见信号的意义、特性及波形 2．学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法 二．实验原理 信号一般是随时间而变化的某些物理量。按照自变量的取值是否连续，信号分为连续时间信号和离散时间信号，一般用()f t 和()f k 来表示。若对信号进行时域分析，就需要绘制其波形，如果信号比较复杂，则手工绘制波形就变得很困难，且难以精确。MATLAB 强大的图形处理功能及符号运算功能，为实现信号的可视化及其时域分析提供了强有力的工具。 根据MATLAB 的数值计算功能和符号运算功能，在MATLAB 中，信号有两种表示方法，一种是用向量来表示，另一种则是用符号运算的方法。在采用适当的MATLAB 语句表示出信号后，就可以利用MATLAB 中的绘图命令绘制出直观的信号波形了。下面分别介绍连续时间信号和离散时间信号的MATLAB 表示及其波形绘制方法。 1.连续时间信号 所谓连续时间信号，是指其自变量的取值是连续的，并且除了若干不连续的点外，对于一切自变量的取值，信号都有确定的值与之对应。从严格意义上讲，MATLAB 并不能处理连续信号。在MATLAB 中，是用连续信号在等时间间隔点上的样值来近似表示的，当取样时间间隔足够小时，这些离散的样值就能较好地近似出连续信号。在MATLAB 中连续信号可用向量或符号运算功能来表示。 ⑴ 向量表示法 对于连续时间信号()f t ，可以用两个行向量f 和t 来表示，其中向量t 是用形如12::t t p t 的命令定义的时间范围向量，其中，1t 为信号起始时间，2t 为终止时间，p 为时间间隔。向量f 为连续信号()f t 在向量t 所定义的时间点上的样值。 说明：plot 是常用的绘制连续信号波形的函数。 严格说来，MATLAB 不能表示连续信号，所以，在用plot()命令绘制波形时，要对自变量t 进行取值，MATLAB 会分别计算对应点上的函数值，然后将各个数据点通过折线连接起来绘制图形，从而形成连续的曲线。因此，绘制的只是近似波形，而且，其精度取决于t 的取样间隔。t 的取样间隔越小，即点与点之间的距离越小，则近似程度越好，曲线越光滑。例如：图1-1是在取样间隔为p=0.5时绘制的波形，而图1-2是在取样间隔p=0.1时绘制的波形，两相对照，可以看出图1-2要比图1-1光滑得多。

数字信号处理实验3-离散系统的变换域分析

实验3 离散系统的变换域分析 一、实验目的： 加深对离散系统的频率响应分析和零、极点分布的概念理解。 二、实验原理： 离散系统的时域方程为 ∑∑==-=-M m m N k n m n x b k n y a 00)()( 其变换域分析方法如下： X(z)H(z)Y(z) )()()()()(=?-= *=∑∞-∞=m m n h m x n h n x n y 系统函数为 N N M M z a z a a z b z b b z X z Y z H ----++++++==......)()()(110110 分解因式 ∏∏∑∑=-=-=-=---==k k M m m k k k M m m m z d z c K z a z b z H 111100)1() 1()( ， 其中 m c 和 k d 称为零、极点。 在MATLAB 中，可以用函数[z ，p ，K]=tf2zp （num ，den ）求得有理分式形式的系统函数的零、极点，用函数zplane （z ，p ）绘出零、极点分布图；也可以用函数zplane （num ，den ）直接绘出有理分式形式的系统函数的零、极点分布图。使h=freqz(num,den,w)函数可求系统的频率响应，w 是频率的计算点，如w=0:pi/255:pi, h 是复数，abs(h)为幅度响应，angle(h)为相位响应。另外，在MATLAB 中，可以用函数 [r ，p ，k]=residuez （num ，den ）完成部分分式展开计算；可以用函数sos=zp2sos （z ，p ，K ）完成将高阶系统分解为2阶系统的串联。 （在实验报告中对这几种函数的使用方法及参数含义做出说明，这一部分手写） 三、实验内容 例1 求下列直接型系统函数的零、极点，并将它转换成二阶节形式 解 用MATLAB 计算程序如下： num=[1 -0.1 -0.3 -0.3 -0.2]; den=[1 0.1 0.2 0.2 0.5]; [z,p,k]=tf2zp(num,den); disp('零点');disp(z); disp('极点');disp(p); disp('增益系数');disp(k); sos=zp2sos(z,p,k);

信号与系统实验实验报告

信号与系统实验实验报 告 文档编制序号：[KKIDT-LLE0828-LLETD298-POI08]

实验五连续系统分析一、实验目的 深刻理解连续时间系统的系统函数在分析连续系统的时域特性、频域特性及稳定性中的重要作用及意义，掌握根据系统函数的零极点设计简单的滤波器的方法。掌握利用MATLAB分析连续系统的时域响应、频响特性和零极点的基本方法。 二、实验原理 MATLAB提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频域响应等分析函数。 三、实验内容 1.已知描述连续系统的微分方程为，输入，初始状态 ，计算该系统的响应，并与理论结果比较，列出系统响应分析的步骤。 实验代码： a=[1 10]; b=[2]; [A B C D]=tf2ss(b,a); sys=ss(A,B,C,D); t=0: :5; xt=t>0; sta=[1]; y=lsim(sys,xt,t,sta); subplot(3,1,1); plot(t,y); xlabel('t'); title('系统完全响应 y(t)'); subplot(3,1,2); plot(t,y,'-b'); hold on yt=4/5*exp(-10*t)+1/5; plot(t,yt,' : r'); legend('数值计算','理论计算'); hold off xlabel('t'); subplot(3, 1 ,3); k=y'-yt; plot(t,k); k(1) title('误差');

实验结果： 结果分析： 理论值 y(t)=0. 8*exp(-10t)*u(t)+ 程序运行出的结果与理论预期结果相差较大误差随时间增大而变小，初始值相差最大，而后两曲线基本吻合，表明该算法的系统响应在终值附近有很高的契合度，而在初值附近有较大的误差。 2.已知连续时间系统的系统函数为，求输入分别为，， 时，系统地输出，并与理论结果比较。 a=[1,3,2,0]; b=[4,1]; sys=tf(b,a); t=0: :5; x1=t>0; x2=(sin(t)).*(t>0); x3=(exp(-t)).*(t>0); y1=lsim(sys,x1,t); y2=lsim(sys,x2,t); y3=lsim(sys,x3,t); subplot(3,1,1); plot(t,y1); xlabel('t'); title('X(t)=u(t)'); subplot(3,1,2); plot(t,y2); xlabel('t'); title('X(t)=sint*u(t)'); subplot(3, 1 ,3); plot(t,y3); xlabel('t'); title('X(t)=exp(-t)u(t)'); 实验结果： 结果分析： a=[1,3,2,0]; b=[4,1]; sys=tf(b,a); t=0: :5; x1=t>0; x2=(sin(t)).*(t>0); x3=(exp(-t)).*(t>0); y1=lsim(sys,x1,t); y2=lsim(sys,x2,t); y3=lsim(sys,x3,t); subplot(3,1,1); plot(t,y1,'-b');

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

实验报告 实验二 信号、系统及系统响应，离散系统的时域分析 一、实验目的 （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 ωωω=-式中

信号与系统分析实验

2019年信号与系统分析实验报告 班级： 通信工程170x 班 姓名：苏斌斌 学号：2017xxx 成绩： ================================================================ 实验1: MATLAB 中信号的表示、运算与分析 实验任务：1. 对信号进行数学表示、波形绘制及信号参数观察分析； 2．对信号进行翻转、时移与尺度变换，微分与积分等基本运算； 3．对周期信号进行分解与合成分析； 4．对信号进行频域分析。 实验目的：能够运用合适的数学方法和工具软件，完成信号在时域和频域的分析与 运算，取得有效结论，加深对信号分析基本概念和基本原理的理解，掌握信号分析的基本方法。 实验内容： 1. 完成预习PPT 中MA TLAB 相关程序的学习和练习。 2. 已知x (t )波形如下图所示，使用利用MATLAB 画出的2个由x (t )位移、尺度、翻转结合产生的波形。 t x (t )

3. 已知x= [1,2,4,6,8]，n=0,1,2,3,4 ，h= [1,1,1,1,1]，n=-1,0,1,2,3 ，使用MATALAB 画出x，h以及y=x*h的序列。（使用stem语句） 序列h波形 序列y=x*h波形 4. 通过MATALAB观察课本P66例3.2.4周期方波傅里叶级数展开效果。

5. 周期三角波信号傅里叶级数展开式如下所示，分析其频谱并通过MATLAB 验证。 由对称性知 0n b =，可得：002200 2 4)2 ()cos d /2 0T T n A n n a f t n t t T A n πω-?-= =?=?? ，为奇数 ， 其中002/T ωπ=，所以02241()cos , 1,3,5,2n A A f t n t n n ωπ?? = -= ??? ∑ 可得f (t )幅度谱为：22 4, 1,3,5, 02 n A n n A A n π?=??=??=?? 利用MATLAB 分析幅度为1v ，频率为1KHz 周期三角波的幅度谱和谐波合成波形； 实验任务：1. 使用MA TLAB 求解连续LTI 系统的单位冲激响应和零状态响应；

数字信号处理matlab实验4 离散系统的变换域分析

实验4离散系统的变换域分析 实验目的：加深对离散系统的频率响应分析和零、极点分布的概念理解。 实验原理：离散系统的时域方程为 ∑∑==-=-M k k N k k k n x p k n y d 00) ()(其变换域分析方法如下： 频域)()()(][][][][][ΩΩ=Ω?-= *=∑∞ -∞=H X Y m n h m x n h n x n y m 系统的频率响应为Ω -Ω-Ω-Ω-++++++=ΩΩ=ΩjN N j jM M j e d e d d e p e p p D p H ......)()()(1010Z 域)()()(][][][][][z H z X z Y m n h m x n h n x n y m =?-=*=∑∞ -∞=系统的转移函数为N N M M z d z d d z p z p p z D z p z H ----++++++==......)()()(110110分解因式∏-∏-=∑∑==-=-=-=-N i i M i i N i i k M i i k z z K z d z p z H 11110 0)1()1()(λξ，其中i ξ和i λ称为零、极点。 在MATLAB 中，可以用函数[z，p，K]=tf2zp（num，den）求得有理分式形式的系统转移函数的零、极点，用函数zplane（z，p）绘出零、极点分布图；也可以用函数zplane（num，den）直接绘出有理分式形式的系统转移函数的零、极点分布图。使h=freqz(num,den,w)函数可求系统的频率响应，w 是频率的计算点，如w=0:pi/255:pi,h 是复数，abs(h)为幅度响应，angle(h)为相位响应。另外，在MATLAB 中，可以用函数[r，p，k]=residuez（num，den）完成部分分式展开计算；可以用函数sos=zp2sos（z，p，K）完成将高阶系统分解为2阶系统的串联。例1 求下列直接型系统函数的零、极点，并将它转换成二阶节形式 解用MATLAB 计算程序如下： num=[1-0.1-0.3-0.3-0.2]; den=[10.10.20.20.5];

信号与系统实验报告

实验一二阶有源电路滤波系统设计以及频率响应特性分析一、实验目的 通过定性观察不同输入信号下的电路输出，理解输入信号的频谱分布与系统频率响应的相对关系对信号的影响；通过二阶系统频率响应特性的定量测试，考察系统参数对典型二阶系统滤波特性的影响。 二、实验设备 信号与系统实验电路模板、信号发生器、双综示波器、直流电源等。 三、实验内容与要求 1. 已知一有源二阶系统如下图： ①理论分析其作为微分电路、积分电路或比例网络的条件： ②为验证等效微分电路，选择一组元件参数，如R1＝1k、R2＝1k、C1＝0.1uF、C2＝ 0. 01uF，τ1= 0.0001s,τ2 = 0.00001s，用信号发生器产生频率为100Hz，占空比为1:1 的周期方波信号，实验观察记录输入、输出波形并做出适当解释； ③类比②，适当选择元件参数以及输入信号来验证等效积分电路和比例电路。要求画出 所选元件参数下的电路频率特性、记录输入、输出波形并做出合理解释。 2、二阶系统频率响应特性分析 ①从四种二阶有源滤波参考电路任选1-2种，自选参数，分析计算系统的幅频响应特性；

②将正弦信号输入二阶系统，改变输入信号频率，用示波器观察并记录输出信号的波形、幅值，并与理论计算比较。说明二阶谐振系统的阻尼系数a 或品质因数Q 对系统频率特性以及系统稳定性的影响。 四．理论分析 1.首先，在实际电路系统中并不存在理想的微分、积分或比例电路；只能是针对一定频率范围的信号，特定的电路系统呈现出相应的微分、积分或比例特性。以上电路的传递函数为 ) /1(1 )/1()(2211ττ++- =s C s R s s H ，其中τ1＝R 1C 1, τ2＝R 2C 2。可以看出，这是 一个一阶高通与一个一阶低通环节的串连，其中τ1>τ2。幅频特性、相频特性为： 2 2222121)/1(1 )/1()(τωτωωω++= C R j H ()()()2190ωτωτω?arctg arctg --?=。当输入信号的角频率ω<<1/τ1时， ∣H(j ω)∣≒C 1R 2ω、φ（ω）≒90，即等效于微分电路；当输入信号的角频率ω>>1/τ2时，∣H(j ω)∣≒1/(R 1C 2ω)、φ（ω）≒-90，即等效于积分电路；当输入信号的角频率1/τ2>>ω>>1/τ1时，∣H(j ω)∣≒R 2/R 1、φ（ω）≒0，即等效比例电路。 2.低通滤波电路如下： a) 低通滤波器

信号与系统实验报告

学生实验报告 （理工类） 课程名称：信号与系统实验专业班级：电子信息（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 加深对零，极点分布的概念的理解。 实验内容： 程序： b=[0.036 0.143 0.214 0.143 0.036]; a=[1 -1.035 0.826 0.260 0.040]; [z,p,k]=tf2zp(b,a); disp('零点');disp(z); disp('极点');disp(p); disp('增益系数');disp(k); sos=zp2sos(z,p,k); disp('二阶节');disp(real(sos)); zplane(b,a); k=256; b=[0.036 0.143 0.214 0.143 0.036]; a=[1 -1.035 0.826 0.260 0.040]; w=0:pi/k:pi; h=freqz(b,a,w); subplot(2,2,1); plot(w/pi,real(h));grid title('实部') xlabel('\omega/\pi'); subplot(2,2,2); plot(w/pi,imag(h));grid title('虚部') xlabel('\omega/\pi'); subplot(2,2,3); plot(w/pi,abs(h));grid title('幅度谱') xlabel('\omega/\pi');ylabel('幅值') subplot(2,2,4) plot(w/pi,angle(h));grid title('相位谱') xlabel('\omega/\pi');ylabel('相位') 实验结果： 零点 -0.9861 + 0.1661i -0.9861 - 0.1661i -1.0000 -1.0000

信号与系统实验指导书

信号与系统实验指导书 赵欣、王鹏 信息与电气工程学院 2006.6.26

前言 “信号与系统”是无线电技术、自动控制、生物医学电子工程、信号图象处理、空间技术等专业的一门重要的专业基础课，也是国内各院校相应专业的主干课程。 当前，科学技术的发展趋势既高度综合又高度分化，这要求高等院校培养的大学生，既要有坚实的理论基础，又要有严格的工程技术训练，不断提高实验研究能力、分析计算能力、总结归纳能力和解决各种实际问题的能力。21世纪要求培养“创造型、开发型、应用型”人才，即要求培养智力高、能力强、素质好的人才。 由于该课程核心的基本概念、基本理论和分析方法都非常重要，而且系统性、理论性很强，为此在学习本课程时，开设必要的实验，对学生加深理解、深入掌握基本理论和分析方法，培养学生分析问题和解决问题的能力，以及使抽象的概念和理论形象化、具体化，对增强学习的兴趣有极大的好处，做好本课程的实验，是学好本课程的重要教学辅助环节。 在做完每个实验后，请务必写出详细的实验报告，包括实验方法、实验过程与结果、心得和体会等。

目录 实验一无源和有源滤波器 (1) 实验二方波信号的分解 (6) 实验三用同时分析法观测方波信号的频谱 (8) 实验四二阶网络状态轨迹的显示 (10) 实验五二阶网络函数的模拟 (14) 实验六抽样定理 (18) 附录 (22)

实验一无源和有源滤波器 一、实验目的 1、了解RC无源和有源滤波器的种类、基本结构及其特性。 2、对比研究无源和有源滤波器的滤波特性。 3、学会列写无源和有源滤波器网络函数的方法。 二、基本原理 1、滤波器是对输入信号的频率具有选择性的一个二端口网络，它允许某些频率(通常是某个频带范围)的信号通过，而其它频率的信号受到衰减或抑制，这些网络可以是由RLC元件或RC元件构成的无源滤波器，也可以是由RC元件和有源器件构成的有源滤波器。 2、根据幅频特性所表示的通过或阻止信号频率范围的不同，滤波器可分为低通滤波器(LPF)、高通滤波器(HPF)、带通滤波器(BPF)和带阻滤波器(BEF)四种。我们把能够通过的信号频率范围定义为通带，把阻止通过或衰减的信号频率范围定义为阻带。而通带与阻带的分界点的频率f，称为截止频率或称转折频率。图1-1中的A up为通带的电压放大倍数，f0为中心频率，f CL和f CH分别为低端和高端截止频率。 A A up f C f f C f f CL f CH f f CL f CH f 图1-1 各种滤波器的理想幅频特性 四种滤波器的实验线路如图1-2所示：