实验一时域离散信号与系统变换域分析(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.52
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.20frequency 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.51frequency 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 结果图形
上图为零状态响应、下图为全响应。
2
4
6
8
10
12
14
16
18
20
-3-2-101230
2
4
6
8
10
12
14
16
18
20
-2-10123
5. 离散系统的Z 域分析
%11 利用系统函数)(z H 分析系统的稳定性。假设系统函数如下式:
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) 分析零点分布对于单位脉冲响应的影响。
北京交通大学信号与系统时域分析
【研讨题目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.加深对离散系统的频率响应分析和零、极点分布的概念理解。 二、实验原理: 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)');
§8-1 引 言 一、 变换域分析的目的 变换域分析的目的,在于将原来的求解问题简化。 对于连续时间系统,通过L.T.,可以将原来求解微分方程的问题转变为求解代数方程的问题; 对于离散时间系统,通过Z 变换(Z.T.),可以将原来求解差分方程的问题转变为求解代数方程的问题。 二、 Z 变换的发展史 十八世纪,DeMoivre 提出生成函数,并应用于概率论; 十九世纪Laplace 、二十世纪Seal 对其进行了进一步深入研究; 二十世纪六十年代起,由于计算机技术和控制技术的飞速发展,抽样控制理论的应用,离散信号处理和数字信号处理得到了广泛应用。作为离散时间系统分析的重要工具,Z.T.得到了很大的发展,其用途甚至超过了L.T. 三、 离散时间序列的频域分析方法 离散时间系统和离散时间序列也可以通过正交分解的方法,在频域进行分析。离散系统也有频率响应(对各种频率的离散正弦信号的响应)。傅利叶变换的离散形式——离散傅利叶变换(DFT )——在离散时间系统分析中同样占用很重要的地位,而DFT 的快速算法——FFT ——的提出使得DFT 在各种信号处理场合得到的广泛的应用。 除了DFT 以外,其信号分析方法,如沃尔什变换等,在离散信号处理中同样得到的很广泛的应用。 §8-2 Z 变换及其性质 一、 Z 变换的定义 Z 变换的定义可以从纯数学的角度进行,也可以通过信号分解的 角度提出。后者更加容易理解。本课程中,通过连续时间系统的F.T.,导出Z.T.。 离散时间信号f(k)可以看成是连续时间信号通过抽样而得到的冲激序列: )(k f ——> ∑+∞ -∞ =-= k kT t k f t f ) ()()(δδ 对其)(t f δ 进行F.T.: () ∑ ∑ ∑? ∑??∑ ? ∞ +-∞ =-∞+-∞ =-∞+-∞=∞ +∞ --∞ +-∞=∞+∞ --∞+∞--∞ +-∞=+∞ ∞ --= = ? ? ? ? ??-=-=??? ? ????-==k kT j k kT j k t j k t j t j k t j e k f e k f dt e kT t k f dt e kT t k f dt e kT t k f dt e t f j F ω ωωωωωδδδδω)()()()()()()()()()( 根据Dirichlet 条件,只有在信号满足绝对可积条件——这里可以
目录 第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) 绘制结果图形并对仿真结果进行详细的分析。
实验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.了解时域离散信号的产生及基本运算实现。 2.掌握离散时间傅里叶变换实现及系统分析方法。 3. 熟悉离散时间傅里叶变换性质。 4. 掌握系统Z 域分析方法。 5. 培养学生运用软件分析、处理数字信号的能力。 二、实验设备 1、计算机 2、Matlab7.0以上版本 三、实验内容 1、对于给定的时域离散信号会进行频谱分析,即序列的傅里叶变换及其性质分析。 2、对于离散系统会进行频域分析及Z 域分析。包括频谱特性、零极点画图、稳定性分析。 3、对于差分方程会用程序求解,包括求单位冲击序列响应,零输入响应、零状态响应、全响应,求其系统函数,及其分析。 4、信号时域采样及其频谱分析,序列恢复。 5、扩展部分主要是关于语音信号的读取及其播放。 四、实验原理 1、序列的产生及运算 在Matlab 中自带了cos 、sin 、exp (指数)等函数,利用这些函数可以产生实验所需序列。 序列的运算包括序列的加法、乘法,序列)(n x 的移位)(0n n x -,翻褶)(n x -等。序列的加法或乘法指同序号的序列值逐项对应相加或相乘,但Matlab 中“+”“.*”运算是对序列的值直接进行加或乘,不考虑两序列的序号是否相同,因此编程时考虑其序号的对应。 2、序列的傅里叶变换及其性质 序列的傅里叶变换定义:)(|)(|)()(ω?ωωω j j n n j j e e X e n x e X ==∑∞-∞=-,其幅度特性为|)(|ωj e X , 在Matlab 中采用abs 函数;相位特性为)(ω?,在Matlab 中采用angle 函数。 序列傅里叶变换的性质:
实验一信号与系统的时域分析 一、实验目的 1、熟悉与掌握常用的用于信号与系统时域仿真分析的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),根
一.典型连续信号和离散信号的时域波形。 1.单边指数信号)()(t u Ae t y t α=; 2.单位冲激信号)()(0t t t y +=δ; 3.单位阶跃信号)()(0t t u t y +=; 4.矩形脉冲信号)]()([)(21t t u t t u A t y +-+?=; 5.正弦信号)()sin()(t u t A t y ω?=; 6.单位序列)()(0n n n y +=δ; 7.单位阶跃序列)()(0n n u n y +=; 8.单位矩形序列)()()(21n n u n n u n y +-+=; 9.指数序列)()(n u a A n y n ?=; 10.正弦序列)()sin()(n u n A n y ω?=。
单边指数信号 function zhishu(A,a,t1,t2,dt) t1=0 t2=10 A=1 A=-0.4 dt=0.01 t=t1:dt:t2; y=A*exp(a*t); plot(t,y) axis([t1,t2,0,1.2]) xlabel('t') ylabel('y(t)') title(' 单边指数信号') 单位冲激信号 function chongji(t1,t2,t0) dt=0.01; t1=10; t2=-5; t=t1:dt:t2; n=length(t); x=zeros(1,n); x(1,(-t0-t1)/dt+1)=1/dt; stairs(t,x); axis([t1,t2,0,1.2/dt]) xlabel('t') ylabel('y(t)') title('单位冲激信号')
实验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);
实验一 时域离散信号的产生与基本运算 一、实验目的 1、了解常用的时域离散信号及其特点。 2、掌握MATLAB 产生常用时域离散信号的方法。 3、掌握时域离散信号简单的基本运算方法。 二、实验内容 1、自己设定参数,分别表示并绘制单位抽样序列、单位阶跃序列、正弦序列、 实指数序列、随机序列。 2、自己设定参数,分别表示并绘制信号移位、信号相加、信号相乘、信号翻转、 信号和、信号积、信号能量。 3、已知信号 (1) 描绘)(n x 序列的波形。 (2) 用延迟的单位脉冲序列及其加权和表示)(n x 序列。 (3) 描绘以下序列的波形:)2()(),2(2)(),2(2)(321n x n x n x n x n x n x -=+=-= 三、实现步骤 1、自己设定参数,分别表示并绘制单位抽样序列、单位阶跃序列、正弦序列、 实指数序列、随机序列。 (1)单位抽样序列 程序: x=zeros(1,10);
x(2)=1; stem(x,'filled') axis([0,10,-0.2,1]); title('μ¥??3é?ùDòáD'); -0.20 0.2 0.4 0.6 0.8 图 1 (2)单位阶跃序列 程序: N=10; u=ones(1,N); stem(u,'filled') axis([-10,10,0,1]); title('μ¥???×??DòáD');
00.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 单位阶跃序列 图 2 (3)正弦序列 程序: x=-20:1:20; y=sin(0.2*pi.*x+0.5*pi); stem(x,y,'filled'); axis([-20,20,-2,2]); title('?y?òDòáD');
实验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
实验1用M A T L A B产生时域离散信号 一、.实验目的: 1、了解常用时域离散信号及其特点 2、掌握用MATLAB产生时域离散信号的方法 二、实验内容及步骤 1、阅读并上机验证实验原理部分的例题程序,理解每一条语句的含义。 改变例题中的有关参数(如信号的频率、周期、幅度、显示时间的取值范围、采样点数等),观察对信号波形的影响。 2、编写程序,产生以下离散序列: n1=-3;n2=4;n0=0; n=n1:n2; x=[n==n0]; stem(n,x,'filled'); axis([n1,n2,0,*max(x)]); xlabel('时间(n)');ylabel('幅度x(n)'); title('单位脉冲序列'); (2)n1=-5;n2=5;n0=0; n=n1:n2; x=[n>=n0]; stem(n,x,'filled') axis([n1,n2,0,*max(x)]); xlabel('时间(n)');ylabel('幅度x(n)'); title('单位阶跃序列'); n1=20;a=;w=*pi; n=0:n1; x=exp((a+j*w)*n); subplot(2,2,1);plot(n,real(x)); title('复指数信号的实部'); subplot(2,2,3);stem(n,real(x),'filled'); title('复指数序列的实部'); subplot(2,2,2);plot(n,imag(x)); title('复指数信号的虚部'); subplot(2,2,4);stem(n,imag(x),'filled'); title('复指数序列的虚部');
实验一---时域离散信号与系统变换域分析(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 函数。
实验一时域离散信号与系统变换域分析(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)
实验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];
信号与系统 课程实验报告 实验四 实验名称 信号的时域分析 及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('脉冲响应') 结果:
FFT对连续信号和时域离散信号进行谱分析
————————————————————————————————作者:————————————————————————————————日期:
一、实验目的与要求 学习用FFT对连续信号和时域离散信号进行谱分析的方法,了解可能出现的分析误差及其原因,以便正确应用FFT。 二、实验原理 用FFT对信号作频分析是学习数字信号处理的重要内容,经常需要进行分析的信号是模拟信号的时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D和分析误差。频谱分辨率直接和FFT的变换区间N有关,因为FFT能够实现的频率分辨率是2π/N,因此要求2π/N 小于等于D。可以根据此式选择FFT的变换区间N。误差主要来自于用FFT作频谱分析时,得到的是离散谱,而信号(周期信号除外)是连续谱,只有当N较大时,离散谱的包络才能逼近连续谱,因此N要适当选择大一些。 三、实验步骤及内容 (1)对以下序列进行FFT分析: x1(n)=R4(n) n+1 0≤n≤3 x2(n)={ 8-n 4≤n≤7 0 其它n 4-n 0≤n≤3 X3(n)={ n-3 4≤n≤7 0 其它n 选择FFT的变换区间N为8和16两种情况进行频谱分析,分别打印出幅频特性曲线,并进行讨论、分析与比较 xn1=[1 1 1 1]; Xk18=fft(xn1,8); yn11=abs(Xk18); n11=0:length(yn11)-1; Xk116=fft(xn1,16); yn12=abs(Xk116); n12=0:length(yn12)-1; n=0:3; x21=n+1; x31=4-n; n=4:7; x22=8-n; x32=n-3; xn2=[x21,x22]; Xk28=fft(xn2,8); yn21=abs(Xk28); n21=0:length(yn21)-1; Xk216=fft(xn2,16); yn22=abs(Xk216); n22=0:length(yn22)-1; xn3=[x31,x32]; Xk38=fft(xn3,8);
实验一 连续时间信号与系统的时域分析及MATLAB 实现 实验目的: 1. 掌握MATLAB 相关函数的调用格式及作用; 2. 会用MATLAB 生成和实现连续时间信号,并对连续系统进行分析。 实验内容: 1. 先在MATLAB 帮助菜单下查看以下函数的调用格式及作用; Plot,ezplot,sym,subplot,stepfun,diff,int,conv,impulse,step,lsim,roots. 2.绘出下列信号的时域波形。 (1)()sin t f t t =;(2)()5sin 4f t t π??= ??? ; (3)()52t f t e -=;(4)()42j t f t e π=(实部、虚部、相角及模) 3.完成下列信号相应的运算和变换。 (1)已知信号()()()11222f t t u t u t ? ?=+?+--?? ?? ???,绘制时域波形并用MA TLAB 求()()()()2,,2,f t f t f t f t +--及其波形。 (2)已知信号()()()()144f t t u t u t =-+?--????和()2sin(2)f t t π=,用MATLAB 计算()()12f t f t +;()1d f t dt ,()1t f d ττ-∞?并绘制波形。; (3)计算信号()()()12f t u t u t =--和()()()24f t u t u t =--的卷积。 4.已知描述某连续连续系统的微分方程为: ()()()()()''''282y t y t y t f t f t ++=+ 试用MATLAB 绘出该系统的冲激响应和阶跃响应,并求出系统在()()2t f t e u t -=的激励下的零状态响应。
离散系统变换域分析 实验目的 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
实验一 时域离散信号与系统变换域分析 一、实验目的 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 函数。 序列傅里叶变换的性质: