信号与系统_实验二
电路、信号与系统(2)实验指导书
描述线性时不变离散系统的差分方程为
编写求解上述方程的通用程序。
[建模]
将方程变形可得(用MATLAB语言表示)
a(1)*y(n)= b(1)*u(n)+…+ b(nb)*u(n-nb+1)- a(2)*y(n-1)-…- a(na)*y(n-na+1)
令us== [u(n),…, u(n-nb+1)]; ys=[y(n-1),…, y(n-na+1)]
x(n)={2,1,-1,3,1,4,3,7}(其中加下划线的元素为第0个采样点)在MATLAB中表示为:
n=[-3,-2,-1,0,1,2,3,4]; x=[2,1,-1,3,1,4,3,7];
当不需要采样位置信息或这个信息是多余的时候,可以只用x向量来表示。
(一)离散信号的MATLAB表述
[问题]
实验一连续时间信号与系统分析
一、实验目的
1、了解连续时间信号的特点;
2、掌握连续时间信号的MATLAB描述;
3、掌握连续LTI系统单位冲激响应的求解方法;
4、掌握连续LTI系统的零状态响应的求解方法。
二、实验内容
严格说来,只有用符号推理的方法才能分析连续系统,用数值方法是不能表示连续信号的,因为它给出的是各个样点的数据。只有当样本点取得很密时才可看成连续信号。所谓很密,是相对于信号变化的快慢而言的。以下均假定相对于采样点密度而言,信号变化足够慢。
elseif lu<lh nh=0; nu=lh-lu;
else nu=0; nh=0;
end
dt=0.1;
lt=lmax;
u=[zeros(1, lt), uls, zeros(1, nu), zeros(1, lt)];
信号与系统实验实验报告
信号与系统实验实验报告一、实验目的本次信号与系统实验的主要目的是通过实际操作和观察,深入理解信号与系统的基本概念、原理和分析方法。
具体而言,包括以下几个方面:1、掌握常见信号的产生和表示方法,如正弦信号、方波信号、脉冲信号等。
2、熟悉线性时不变系统的特性,如叠加性、时不变性等,并通过实验进行验证。
3、学会使用基本的信号处理工具和仪器,如示波器、信号发生器等,进行信号的观测和分析。
4、理解卷积运算在信号处理中的作用,并通过实验计算和观察卷积结果。
二、实验设备1、信号发生器:用于产生各种类型的信号,如正弦波、方波、脉冲等。
2、示波器:用于观测输入和输出信号的波形、幅度、频率等参数。
3、计算机及相关软件:用于进行数据处理和分析。
三、实验原理1、信号的分类信号可以分为连续时间信号和离散时间信号。
连续时间信号在时间上是连续的,其数学表示通常为函数形式;离散时间信号在时间上是离散的,通常用序列来表示。
常见的信号类型包括正弦信号、方波信号、脉冲信号等。
2、线性时不变系统线性时不变系统具有叠加性和时不变性。
叠加性意味着多个输入信号的线性组合产生的输出等于各个输入单独作用产生的输出的线性组合;时不变性表示系统的特性不随时间变化,即输入信号的时移对应输出信号的相同时移。
3、卷积运算卷积是信号处理中一种重要的运算,用于描述线性时不变系统对输入信号的作用。
对于两个信号 f(t) 和 g(t),它们的卷积定义为:\(f g)(t) =\int_{\infty}^{\infty} f(\tau) g(t \tau) d\tau \在离散时间情况下,卷积运算为:\(f g)n =\sum_{m =\infty}^{\infty} fm gn m \四、实验内容及步骤实验一:常见信号的产生与观测1、连接信号发生器和示波器。
2、设置信号发生器分别产生正弦波、方波和脉冲信号,调整频率、幅度和占空比等参数。
3、在示波器上观察并记录不同信号的波形、频率和幅度。
实验二 连续信号与系统的频域分析
(2)绘出f(t)的时域波形及频谱图。
f(t) 1
(2)电路的系统函数为 H(jω)
1 j 1 j
(b)用MATLAB求系统的单位冲激响应。
(c) 当输入为 f (t) sint sin3t t
求系统的稳态响应。
程序清单: syms w t; Hw=(1-j*w)/(1+j*w); ht=ifourier(Hw,t); ft=sin(t)+sin(3*t); Fw=fourier(ft); Yw=Fw*Hw; yt=ifourier(Yw,t);
函数fourier()——傅立叶正变换 函数ifourier()——傅立叶逆变换
3、连续时间系统的响应
已知某电路的系统传递函数为 H(jω)=1/(0.08(jω)2+0.4jω+1)
用MATLAB绘制系统的幅频特性曲线和相频特性曲线,并分
析该系统的频率特性。 系统频率特性
H(j)=|H(j)|ej() |H(j)|——系统的幅频特性
for i=1:9
a(i)=subs(an,n,i); % 计算系数a1~a9,存于数组a中 end
a0=double(a0);a=double(a); %转换成数值型
stem(0,a0,i,a); %绘f(t)的频谱图
2、非周期信号的分析 (1)已知某一连续时间信号为
f t e2 t
试绘出它的时域波形及相应的频谱图。
山科大信号与系统实验二-LTI系统的响应
实验二 LTI 系统的响应一、 实验目的1. 熟悉连续时间系统的单位冲激响应、阶跃响应的意义及求解方法2. 熟悉连续(离散)时间系统在任意信号激励下响应的求解方法3. 熟悉应用MATLAB 实现求解系统响应的方法二、 实验原理1.连续时间系统在MATLAB 中有专门用于求解连续系统冲激响应和阶跃响应, 并绘制其时域波形的函数impulse( ) 和step( )。
如果系统输入为f (t ),冲激响应为h(t),系统的零状态响应为y (t ),则有:()()()y t h t f t =*。
若已知系统的输入信号及初始状态,我们便可以用微分方程的经典时域求解方法,求出系统的响应。
但是对于高阶系统,手工计算这一问题的过程非常困难和繁琐。
在MATLAB 中,应用lsim( )函数很容易就能对上述微分方程所描述的系统的响应进行仿真,求出系统在任意激励信号作用下的响应。
lsim( )函数不仅能够求出连续系统在指定的任意时间范围内系统响应的数值解,而且还能同时绘制出系统响应的时域波形图。
说明:(1)当系统有初始状态时,若使用lsim( )函数求系统的全响应,就要使用系统的状态空间描述法,即首先要根据系统给定的方式,写出描述系统的状态方程和输出方程。
假如系统原来给定的是微分方程或系统函数,则可用相变量法或对角线变量等方法写出系统的状态方程和输出方程。
其转换原理如前面实验四所述。
(2)显然利用lsim( )函数不仅可以分析单输入单输出系统,还可以分析复杂的多输入多输出系统。
例题1: 若某连续系统的输入为e (t ),输出为r (t ),系统的微分方程为:''()5'()6()3'()2()y t y t y t f t f t ++=+① 求该系统的单位冲激响应h (t )及其单位阶跃响应g (t )。
a=[1 5 6];b=[3 2];subplot(2,1,1),impulse(b,a,0:0.01:5); subplot(2,1,2),step(b,a,0:0.01:5);-10123Im pulse ResponseTim e (sec)A m p l i t u d e0.20.40.60.8Step ResponseTim e (sec)A m p l i t u d e② 若2()()tf t e t ε-= 求出系统的零状态响应y(t ) a=[1 5 6];b=[3 2]; t=0:0.01:5; f=exp(-2*t); lsim(b,a,f,t);Linear Sim ulation R esultsTim e (sec)A m p l i t u d e-0.20.20.40.60.811.2例题2 已知一个过阻尼二阶系统的状态方程和输出方程分别为:010'()()()232x t X t f t ⎡⎤⎡⎤=+⎢⎥⎢⎥--⎣⎦⎣⎦, r (t )=[0 1]X (t ) 。
信号实验报告 2
信号与系统实验报告实验一信号与系统的时域分析一、实验目的1、熟悉和掌握常用的用于信号与系统时域仿真分析的MA TLAB函数;2、学会用MA TLAB进行信号基本运算的方法;3、掌握连续时间和离散时间信号的MA TLAB产生,掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB编程。
二、实验内容Q1-1:修改程序Program1_1,将dt改为0.2,再执行该程序,保存图形,看看所得图形的效果如何?dt = 0.01时的程序clear, % Clear all variablesclose all, % Close all figure windowsdt = 0.01; % Specify the step of time variablet = -2:dt:2; % Specify the interval of timex = sin(2*pi*t); % Generate the signalplot(t,x) % Open a figure window and draw the plot of x(t)title('Sinusoidal signal x(t)')xlabel('Time t (sec)')dt = 0.2时的程序clear, % Clear all variablesclose all, % Close all figure windowsdt = 0.2; % Specify the step of time variablet = -2:dt:2; % Specify the interval of timex = sin(2*pi*t); % Generate the signalplot(t,x) % Open a figure window and draw the plot of x(t)title('Sinusoidal signal x(t)')xlabel('Time t (sec)')dt = 0.01时的信号波形dt = 0.2时的信号波形这两幅图形有什么区别,哪一幅图形看起来与实际信号波形更像?答:dt = 0.01的图形比dt = 0.2的图形光滑,dt = 0.01看起来与实际信号波形更像。
信号与系统实验报告
信号与系统实验报告目录1. 内容概要 (2)1.1 研究背景 (3)1.2 研究目的 (4)1.3 研究意义 (4)2. 实验原理 (5)2.1 信号与系统基本概念 (7)2.2 信号的分类与表示 (8)2.3 系统的分类与表示 (9)2.4 信号与系统的运算法则 (11)3. 实验内容及步骤 (12)3.1 实验一 (13)3.1.1 实验目的 (14)3.1.2 实验仪器和设备 (15)3.1.4 实验数据记录与分析 (16)3.2 实验二 (16)3.2.1 实验目的 (17)3.2.2 实验仪器和设备 (18)3.2.3 实验步骤 (19)3.2.4 实验数据记录与分析 (19)3.3 实验三 (20)3.3.1 实验目的 (21)3.3.2 实验仪器和设备 (22)3.3.3 实验步骤 (23)3.3.4 实验数据记录与分析 (24)3.4 实验四 (26)3.4.1 实验目的 (27)3.4.2 实验仪器和设备 (27)3.4.4 实验数据记录与分析 (29)4. 结果与讨论 (29)4.1 实验结果汇总 (31)4.2 结果分析与讨论 (32)4.3 结果与理论知识的对比与验证 (33)1. 内容概要本实验报告旨在总结和回顾在信号与系统课程中所进行的实验内容,通过实践操作加深对理论知识的理解和应用能力。
实验涵盖了信号分析、信号处理方法以及系统响应等多个方面。
实验一:信号的基本特性与运算。
学生掌握了信号的表示方法,包括连续时间信号和离散时间信号,以及信号的基本运算规则,如加法、减法、乘法和除法。
实验二:信号的时间域分析。
在本实验中,学生学习了信号的波形变换、信号的卷积以及信号的频谱分析等基本概念和方法,利用MATLAB工具进行了实际的信号处理。
实验三:系统的时域分析。
学生了解了线性时不变系统的动态响应特性,包括零状态响应、阶跃响应以及脉冲响应,并学会了利用MATLAB进行系统响应的计算和分析。
信号与系统实验报告
信号与系统实验报告一、信号的时域基本运算1.连续时间信号的时域基本运算两实验之一实验分析:输出信号值就等于两输入信号相加(乘)。
由于b=2,故平移量为2时,实际是右移1,符合平移性质。
两实验之二心得体会:时域中的基本运算具有连续性,当输入信号为连续时,输出信号也为连续。
平移,伸缩变化都会导致输出结果相对应的平移伸缩。
2.离散时间信号的时域基本运算两实验之一实验分析:输出信号的值是对应输入信号在每个n值所对应的运算值,当进行拉伸变化后,n值数量不会变,但范围会拉伸所输入的拉伸系数。
两实验之二心得体会:离散时间信号可以看做对连续时间信号的采样,而得到的输出信号值,也可以看成是连续信号所得之后的采样值。
二、连续信号卷积与系统的时域分析1.连续信号卷积积分两实验之一实验分析:当两相互卷积函数为冲激函数时,所卷积得到的也是一个冲激函数,且该函数的冲激t值为函数x,函数y冲激t值之和。
两实验之二心得体会:连续卷积函数每个t值所对应的卷积和可以看成其中一个在k值取得的函数与另外一个函数相乘得到的一个分量函数,并一直移动k值直至最后,最后累和出来的最终函数便是所得到的卷积函数。
3.RC电路时域积分两实验之一实验分析:全响应结果正好等于零状态响应与零输入响应之和。
两实验之二心得体会:具体学习了零状态,零输入,全响应过程的状态及变化,与之前所学的电路知识联系在一起了。
三、离散信号卷积与系统的时域分析1.离散信号卷积求和两实验之一实验分析:输出结果的n值是输入结果的k号与另一个n-k的累和两实验之二心得体会:直观地观察到卷积和的产生,可以看成连续卷积的采样形式,从这个方面去想,更能深入地理解卷积以及采样的知识。
2.离散差分方程求解两实验之一实验分析:其零状态响应序列为0 0 4 5 7.5,零输入响应序列为2 4 5 5.5 5.75,全状态响应序列为2 4 9 10.5 13.25,即全状态=零输入+零状态。
两实验之二心得体会:求差分方程时,可以根据全状态响应是由零输入输入以及零状态相加所得,分开来求,同时也加深了自己对差分方程的求解问题的理解。
实验2 信号卷积实验
实验二 信号卷积实验一、实验目的1. 理解卷积的概念及物理意义;2. 通过实验的方法加深对卷积运算的图解方法及结果的理解。
二、实验设备1. 信号与系统实验箱 1台2. 双踪示波器1台3. 铆孔连接线 若干二、实验原理说明卷积积分的物理意义是将信号分解为冲激信号之和,借助系统的冲激响应,求解系统对任意激励信号的零状态响应。
设系统的激励信号为)t (x ,冲激响应为)t (h ,则系统的零状态响应为)(*)()(t h t x t y =()()x t h t d ττ∞-∞=-⎰。
对于任意两个信号)t (f 1和)t (f 2,两者做卷积运算定义为:()()()12f t f t f t d ττ∞-∞=-⎰=)t (f 1*)t (f 2=)t (f 2*)t (f 1。
1. 两个矩形脉冲信号的卷积过程两信号)t (x 与)t (h 都为矩形脉冲信号,如图10-1所示。
下面由图解的方法(图10-1)给出两个信号的卷积过程和结果,以便与实验结果进行比较。
0≤<∞-t210≤≤t 1≤≤t 41≤≤t ∞<≤t 2124τ(b)(a)(c)(d)(e)(f)(g)(h)(i)2卷积结果图10-1 两矩形脉冲的卷积积分的运算过程与结果2. 矩形脉冲信号与锯齿波信号的卷积信号)t(f1为矩形脉冲信号,)t(f2为锯齿波信号,如图10-2所示。
根据卷积积分的运算方法得到)t(f1和)t(f2的卷积积分结果)t(f,如图10-2(c)所示。
(a)(b)(c)图10-2 矩形脉冲信号与锯齿脉冲信号的卷积积分的结果3. 本实验进行的卷积运算的实现方法在本实验装置中采用了DSP数字信号处理芯片,因此在处理模拟信号的卷积积分运算时,是先通过A/D转换器把模拟信号转换为数字信号,利用所编写的相应程序控制DSP芯片实现数字信号的卷积运算,再把运算结果通过D/A转换为模拟信号输出。
结果与模拟信号的直接运算结果是一致的。
信号与系统分析实验报告
信号与系统分析实验报告信号与系统分析实验报告引言:信号与系统分析是电子工程领域中的重要课程之一,通过实验可以更好地理解信号与系统的基本概念和原理。
本实验报告将对信号与系统分析实验进行详细的描述和分析。
实验一:信号的采集与重构在这个实验中,我们学习了信号的采集与重构。
首先,我们使用示波器采集了一个正弦信号,并通过数学方法计算出了信号的频率和幅值。
然后,我们使用数字信号处理器对采集到的信号进行重构,并与原始信号进行比较。
实验结果表明,重构后的信号与原始信号非常接近,证明了信号的采集与重构的有效性。
实验二:线性系统的时域响应本实验旨在研究线性系统的时域响应。
我们使用了一个线性系统,通过输入不同的信号,观察输出信号的变化。
实验结果显示,线性系统对于不同的输入信号有不同的响应,但都遵循线性叠加的原则。
通过分析输出信号与输入信号的关系,我们可以得出线性系统的传递函数,并进一步研究系统的稳定性和频率响应。
实验三:频域特性分析在这个实验中,我们研究了信号的频域特性。
通过使用傅里叶变换,我们将时域信号转换为频域信号,并观察信号的频谱。
实验结果显示,不同频率的信号在频域上有不同的分布特性。
我们还学习了滤波器的设计和应用,通过设计一个低通滤波器,我们成功地去除了高频噪声,并得到了干净的信号。
实验四:系统辨识本实验旨在研究系统的辨识方法。
我们使用了一组输入信号和对应的输出信号,通过数学建模的方法,推导出了系统的传递函数。
实验结果表明,通过系统辨识可以准确地描述系统的特性,并为系统的控制和优化提供了基础。
结论:通过本次实验,我们深入学习了信号与系统分析的基本概念和原理。
实验结果证明了信号的采集与重构的有效性,线性系统的时域响应的线性叠加原则,信号的频域特性和滤波器的设计方法,以及系统辨识的重要性。
这些知识和技能对于我们理解和应用信号与系统分析具有重要的意义。
通过实验的实际操作和分析,我们对信号与系统的理论有了更深入的理解,为我们今后的学习和研究打下了坚实的基础。
2010信号与系统实验2
实验二:连续和离散系统的频域分析一:实验目的1:学习傅里叶正变换和逆变换,理解频谱图形的物理含义2:了解连续和离散时间系统的单位脉冲响应3:掌握连续时间系统的频率特性二:实验原理1. 傅里叶正变换和逆变换公式 正变换:()()j t F f t e dt ωω∞--∞=⎰逆变换:1()()2j t f t F e d ωωωπ∞-∞=⎰2. 频域分析t j tj e d d e t e ωωωπωωωπ⎰⎰∞∞-∞∞-E =E =)(21)(21)(将激励信号分解为无穷多个正弦分量的和。
⎰∞∞-H E =ωωωπωd e t r tj zs )()(21)(,R(ω)为)(t r zs 傅里叶变换;πωωd )(E 各频率分量的复数振幅 激励单位冲激响应时的零状态响应→ )(t δ)(t h单位阶跃响应时的零状态响应激励→)(t u )(t g3 各函数说明:(1)impulse 冲激响应函数:[Y ,X,T]=impulse(num,den);)1()2()1()1()2()1()()()(11++++++++==--n a s a s a m b s b s b s A s B s H n n m m num 分子多项式系数; num=[b(1) b(2) … b(n+1)]; den 分母多项式系数; den=[a(1) a(2) … a(n+1)];Y ,X,T 分别表示输出响应,中间状态变量和时间变量; 如:352)(2+++=s s s s H ,等价于)(2)()(3)(5)(t e t e t r t r t r +=++ 定义den=[1 5 3];num=[1 2]; [Y ,X,T]=impulse(num,den);(2)step 阶跃响应函数:[Y,X,T]=step(num,den);num 分子多项式;den 分母多项式 Y ,X,T 分别表示输出响应,中间状态变量和时间变量;如:352)(2+++=s s s s H ,den=[1 5 3];num=[1 2];[Y ,X,T]= step (num,den);(3)impz 数字滤波器的冲激响应 [h,t] = impz(b,a,n) b 分子多项式系数;a 分母多项式系数;n 采样样本h 离散系统冲激响应;t 冲激时间,其中t=[0:n-1]', n=length(t)时间样本数(4)freqs 频域响应 [h,w] = freqs(b,a,f) b,a 定义同上,f 频率点个数 h 频域响应,w 频域变量)1()2()1()1()2()1()()()(11++++++++==--m a s a s a n b s b s b s A s B s H m m n n三.实验内容1 周期信号傅里叶级数 已知连续时间信号()()2/8cos 3/4coscos )(321ππ++++=t A t A wt A t x ,其中321,,A A A 取值如下:(X 为学号的后两位)]10,1[,5.02321∈⎪⎩⎪⎨⎧===X X A X A X A ]20,11[,55321∈⎪⎩⎪⎨⎧+==-=X X A XA X A ⎪⎩⎪⎨⎧=-=-=X A X A X A 32151020,>X 要求画出信号的时域波形和频域波形(幅度谱和相位谱)。
信号与系统的实验报告
信号与系统的实验报告信号与系统的实验报告引言:信号与系统是电子工程、通信工程等领域中的重要基础学科,它研究的是信号的传输、处理和变换过程,以及系统对信号的响应和特性。
在本次实验中,我们将通过实际操作和数据分析,深入了解信号与系统的相关概念和实际应用。
实验一:信号的采集与重构在这个实验中,我们使用了示波器和函数发生器来采集和重构信号。
首先,我们通过函数发生器产生了一个正弦信号,并将其连接到示波器上进行观测。
通过调整函数发生器的频率和幅度,我们可以观察到信号的不同特性,比如频率、振幅和相位等。
然后,我们将示波器上的信号通过数据采集卡进行采集,并使用计算机软件对采集到的数据进行处理和重构。
通过对比原始信号和重构信号,我们可以验证信号的采集和重构过程是否准确。
实验二:信号的时域分析在这个实验中,我们使用了示波器和频谱分析仪来对信号进行时域分析。
首先,我们通过函数发生器产生了一个方波信号,并将其连接到示波器上进行观测。
通过调整函数发生器的频率和占空比,我们可以观察到方波信号的周期和占空比等特性。
然后,我们使用频谱分析仪对方波信号进行频谱分析,得到信号的频谱图。
通过分析频谱图,我们可以了解信号的频率成分和能量分布情况,进而对信号的特性进行深入研究。
实验三:系统的时域响应在这个实验中,我们使用了函数发生器、示波器和滤波器来研究系统的时域响应。
首先,我们通过函数发生器产生了一个正弦信号,并将其连接到滤波器上进行输入。
然后,我们通过示波器观测滤波器的输出信号,并记录下其时域波形。
通过改变滤波器的参数,比如截止频率和增益等,我们可以观察到系统对信号的响应和滤波效果。
通过对比输入信号和输出信号的波形,我们可以分析系统的时域特性和频率响应。
实验四:系统的频域响应在这个实验中,我们使用了函数发生器、示波器和频谱分析仪来研究系统的频域响应。
首先,我们通过函数发生器产生了一个正弦信号,并将其连接到系统中进行输入。
然后,我们通过示波器观测系统的输出信号,并记录下其时域波形。
《信号与系统》实验报告
信号与系统实验报告班级:姓名:信息与通信工程学院实验一 系统的卷积响应实验性质:提高性 实验级别:必做 开课单位:信息与通信工程学院 学 时:2一、实验目的:深刻理解卷积运算,利用离散卷积实现连续卷积运算;深刻理解信号与系统的关系,学习MATLAB 语言实现信号通过系统的仿真方法。
二、实验设备: 计算机,MATLAB 软件 三、实验原理: 1、 离散卷积和: 调用函数:conv ()∑∞-∞=-==i i k f i f f f conv S )()(1)2,1(为离散卷积和,其中,f1(k), f2 (k) 为离散序列,K=…-2, -1, 0 , 1, 2, …。
但是,conv 函数只给出纵轴的序列值的大小,而不能给出卷积的X 轴序号。
为得到该值,进行以下分析:对任意输入:设)(1k f 非零区间n1~n2,长度L1=n2-n1+1;)(2k f 非零区间m1~m2,长度L2=m2-m1+1。
则:)(*)()(21k f k f k s =非零区间从n1+m1开始,长度为L=L1+L2-1,所以S (K )的非零区间为:n1+m1~ n1+m1+L-1。
2、 连续卷积和离散卷积的关系:计算机本身不能直接处理连续信号,只能由离散信号进行近似: 设一系统(LTI )输入为)(t P ∆,输出为)(t h ∆,如图所示。
)t)()(t h t P ∆∆→)()(lim )(lim )(0t h t h t P t =→=∆→∆∆→∆δ若输入为f(t):∆∆-∆=≈∑∞-∞=∆∆)()()()(k t P k f t f t f k得输出:∆∆-∆=∑∞-∞=∆∆)()()(k t hk f t y k当0→∆时:⎰∑∞∞-∞-∞=∆→∆∆→∆-=∆∆-∆==ττδτd t f k t P k f t f t f k )()()()(lim)(lim )(0⎰∑∞∞-∞-∞=∆→∆∆→∆-=∆∆-∆==τττd t h f k t hk f t y t y k )()()()(lim)(lim )(0所以:∆∆-∆=-==∑⎰→∆)()(lim)()()(*)()(212121k t f k fd t f f t f t f t s τττ如果只求离散点上的f 值)(n f ∆])[()()()()(2121∑∑∞-∞=∞-∞=∆-∆∆=∆∆-∆∆=∆k k k n f k f k n f k fn f所以,可以用离散卷积和CONV ()求连续卷积,只需∆足够小以及在卷积和的基础上乘以∆。
信号与系统的实验报告(2)
信号与系统实验报告——连续时间系统的复频域分析班级:05911101学号:**********姓名:***实验五连续时间系统的复频域分析——1120111487 信息工程(实验班)蒋志科一、实验目的①掌握拉普拉斯变换及其反变换的定义,并掌握MA TLAB 实现方法 ②学习和掌握连续时间系统系统函数的定义及其复频域分析方法③掌握系统零极点的定义,加深理解系统零极点分布与系统特性的关系。
二、实验原理与方法 1、拉普拉斯变换连续时间信号x(t)的拉普拉斯变换定义为:X s =x (t )e −st dt +∞−∞拉普拉斯反变换为:x t =12πj X (s )e st ds σ+j ∞σ−j ∞在MA TLAB 中可以采用符号数学工具箱中的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 。
2、连续时间系统的系统函数连续时间系统的系统函数是系统单位冲激响应的拉氏变换H s =ℎ(t )e −st dt +∞−∞此外,连续时间系统的系统函数还可以由系统输入和输出信号的拉氏变换之比得到H s =Y(s)/X(s) 单位冲激响应h(t)反映了系统的固有性质,而H(s)从复频域反映了系统的固有性质。
对于H(s)描述的连续时间系统,其系统函数s 的有理函数H s =b M s M +b M−1s M−1+⋯+b 0a n s n +a n −1s M−1+⋯+a 03、连续时间系统的零极点分析系统的零点指使式H s 的分子多项式为零的点,极点指使分母多项式为零的点,零点使系统的值为零,极点使系统函数的值无穷大。
信号与系统实验报告
《信号与系统》实验报告湖南工业大学电气与信息工程学院实验一用同时分析法观测50Hz非正弦周期信号的分解与合成一、实验目的1、用同时分析法观测50Hz非正弦周期信号的频谱,并与傅立叶级数各项的频率与系数作比较。
2、观测基波和其谐波的合成。
二、实验设备1、信号与系统实验箱:TKSS -A型或TKSS -B 型TKSS -C 型;2、双踪示波器三、实验原理1、 一个非正弦周期函数可以用一系列频率成整数倍的正弦函数来表示,其中与非正弦具有相同频率的成分称为基波或一次谐波,其他成分则根据其频率为基波频率的2、3、4、…、n 等倍数分别称为二次、三次、四次、…、n 次谐波,其幅度将随着谐波次数的增加而减小,直至无穷小。
2、 不同频率的谐波可以合成一个非正弦周期波,反过来,一个非正弦周期波也可以分解为无限个不同频率的谐波成分,3、 一个非正弦周期函数可以用傅立叶级数来表示,级数各项系数之间的关系可用一个频谱来表示,不同的非正弦周期函数具有不同的频谱图,各种不同波形及其傅氏级数表达式见表2-1,方波频谱图如图2-1表示Um1351/91/51/71/3790ωωωωωω图1-1 方波频谱图表2-1 各种不同波形的傅立叶级数表达式UmtTU 2τ方波UmTU 2τ正弦整流全波UmTU 2τ三角波Um0T2τ正弦整流半波t tUm0tT U 2τ矩形波U1、方波 ())7sin 715sin 513sin 31(sin 4 ++++=t t t t u t u m ωωωωπ 2、三角波())5sin 2513sin 91(sin 82++-=t t t u t u mωωωπ3、半波())4cos 1512cos 31sin 421(2 +--+=t t t u t u m ωωωππ 4、全波 ())6cos 3514cos 1512cos 3121(4 +---=t t t u t u m ωωωπ5、 矩形波())3cos 3sin 312cos 2sin 21cos (sin 2 ++++=t T t T t T U T U t u m m ωτπωτπωτππτ实验装置的结构如图1-2所示DC20f f f f f f 3456图1-2信号分解于合成实验装置结构框图图中LPF 为低通滤波器,可分解出非正弦周期函数的直流分量。
信号与系统期末试卷及答案
读书破万卷下笔如有神实验二利用DFT分析离散信号频谱一、实验目的应用离散傅里叶变换(DFT),分析离散信号的频谱。
深刻理解DFT分析离散信号频谱的原理,掌握改善分析过程中产生的误差的方法。
二、实验原理根据信号傅里叶变换建立的时域与频域之间的对应关系,可以得到有限长序列的离散傅里叶变换(DFT)与四种确定信号傅里叶变换之间的关系(见教材),实现由DFT分析其频谱。
三、实验内容?3的频谱;1.利用FFT分析信号x(310),nn?,1,...,n)?cos(8(1)、确定DFT计算的参数;N=32;n=0:N-1;x=cos(3*pi/8*n);X=fft(x,N);subplot(2,1,1);stem(n,abs(fftshift(X)));ylabel('Magnitude');xlabel('Frequency (rad)');title('朱艺星杨婕婕'); subplot(2,1,2);stem(n,angle(fftshift(X)));ylabel('Phase');xlabel('Frequency(rad)');读书破万卷下笔如有神进行理论值与计算值比较,讨论信号频谱分析过程中误差原因及改善2)(方法。
在频谱分析过程中由于取样频率过低或者由于信号的截取长度不当将会答:产生误差。
可以适当提高取样率,增加样点数,可能会产生混频现象,取样频率过低,来减少混叠对频谱分析所造成的误差。
对于连续周期信号,其时域取样必须kfo,即(其中K≥2*N+1N为最高谐波分量)其取样点数满足时域取样定理:2fm+fo。
≥≥2Nfo+fo;fs截取信号长度不当,会产生功率泄露,对周期序列进行频谱分析时,为避免泄露应做到:截取的长度应取一个基本周期或基本周期的整数倍,若待分析的周期信号事先不知道其确切的周期,则可截取较长时间长度的样点进行分析,以减少功率泄露误差。
信号与系统 matlab实验报告
信号与系统 matlab实验报告《信号与系统 Matlab实验报告》摘要:本实验报告通过使用 Matlab 软件进行信号与系统实验,探讨了信号与系统在数字领域的应用。
实验结果表明,Matlab 软件具有强大的信号处理和系统分析功能,能够有效地进行信号与系统的模拟和分析。
引言:信号与系统是电子工程领域中的重要基础课程,它研究了信号的产生、传输和处理,以及系统对信号的响应和影响。
在数字领域,信号与系统的理论和方法也得到了广泛的应用。
Matlab 软件作为一种强大的数学计算工具,为信号与系统的模拟和分析提供了便利和高效的途径。
实验一:信号的生成与显示在本实验中,我们首先使用 Matlab 软件生成了几种常见的信号,包括正弦信号、方波信号和三角波信号。
通过调整信号的频率、幅度和相位等参数,我们观察了信号的变化,并利用 Matlab 的绘图功能将信号图形显示出来。
实验结果表明,Matlab 软件能够方便地生成各种类型的信号,并能够直观地显示信号的波形和特性。
实验二:信号的采样与重构在本实验中,我们使用 Matlab 软件对信号进行了采样和重构。
我们首先对一个连续信号进行了离散采样,然后利用 Matlab 的插值函数对采样信号进行了重构。
实验结果表明,采样和重构过程中存在信号失真和频率混叠等问题,但通过适当的采样和重构方法,我们能够有效地还原原始信号。
实验三:系统的响应与分析在本实验中,我们使用 Matlab 软件对系统的响应进行了分析。
我们构建了几种常见的系统模型,包括线性时不变系统和滤波器系统,然后利用 Matlab 的系统分析工具对系统的频率响应、相位响应和单位脉冲响应等进行了分析。
实验结果表明,Matlab 软件能够有效地进行系统的模拟和分析,为系统设计和优化提供了有力的支持。
结论:通过本实验,我们深入了解了信号与系统在数字领域的应用,并掌握了使用 Matlab 软件进行信号与系统模拟和分析的方法。
信号与系统实验一、二
chapter1实验内容:1、画出以下连续时间信号的波形1-0)f(t)=cos(2πt)代码如下:pi=3.14159;t=0:0.01:8;fa=cos(2*pi*t);plot(t,fa);1-1)f (t)=sin(2πt)代码如下:pi=3.14159;t=0:0.01:8;fa=sin(2*t*pi); plot(t,fa);2-0)f (t)=Sa(t/π) 代码如下:pi=3.14159;t=0:0.01:8;fa=sinc(t/pi); plot(t,fa);3-0)f (t)=2[u(t 3)- u(t 5)] 代码如下:t=-1:0.01:10;ft=2*((t>=3)-(t>=5)); plot(t,ft);axis([-1,10,0,3]);4-1)f (t)=e t 代码如下:t=0:0.01:10; ft=exp(t);plot(t,ft);4-2)f (t)=e-t u(t) 代码如下:t=0:0.01:10;f1=(t>=0);f2=exp(-t);plot(t,f1.*f2);5-0)f(t)=2e j(π/4)t,画出实部、虚部、模和相角的波形代码如下:t=0:0.01:10;ft=2*exp(j*(pi/4)*t);h=real(ft); %实部g=imag(ft); %虚部r=abs(ft); %模a=angle(ft); %相角subplot(2,2,1),plot(t,h),title('实部') subplot(2,2,3),plot(t,g),title('虚部') subplot(2,2,2),plot(t,r),title('模')subplot(2,2,4),plot(t,a),title('相角')7)f (t) = u(t)代码如下:t=-1:0.01:5ft=(t>=0);plot(t,ft);axis([-1,5,0,1.5]);8)f (t) =δ(t)代码如下:t=-1:0.01:5;ft=(t>=0)-(t>=0.1); plot(t,ft);axis([-1,1,0,1.1]);9)f9为周期矩形信号,其幅度从-1 到1,占空比为75% 代码如下:pi=3.14159;t=-10:0.01/pi:10;ft=square(t,75);plot(t,ft);2、信号本身运算画出f1(t)为宽度是4,高为1,斜度为0.5 的三角脉冲,然后画出f1(-t),f1(2t),f1(2-2t)的波形以及f1(t)的微分和积分波形。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
实 验 报 告姓名:周松江 学号:201201050435 时间:2013年11月11日实验名称:LTI 系统的响应一、 实验目的1. 熟悉连续时间系统的单位冲激响应、阶跃响应的意义及求解方法2. 熟悉连续(离散)时间系统在任意信号激励下响应的求解方法3. 熟悉应用MATLAB 实现求解系统响应的方法二、 实验原理1.连续时间系统对于连续的LTI 系统,当系统输入为f (t ),输出为y (t ),则输入与输出之间满足如下的线性常系数微分方程:()()0()()nmi j i j i j a yt b f t ===∑∑,当系统输入为单位冲激信号δ(t )时产生的零状态响应称为系统的单位冲激响应,用h(t)表示。
若输入为单位阶跃信号ε(t )时,系统产生的零状态响应则称为系统的单位阶跃响应,记为g(t),如下图所示。
系统的单位冲激响应h (t )包含了系统的固有特性,它是由系统本身的结构及参数所决定的,与系统的输入无关。
我们只要知道了系统的冲激响应,即可求得系统在不同激励下产生的响应。
因此,求解系统的冲激响应h(t )对我们进行连续系统的分析具有非常重要的意义。
在MATLAB 中有专门用于求解连续系统冲激响应和阶跃响应, 并绘制其时域波形的函数impulse( ) 和step( )。
如果系统输入为f (t ),冲激响应为h(t),系统的零状态响应为y (t ),则有:()()()y t h t f t =*。
若已知系统的输入信号及初始状态,我们便可以用微分方程的经典时域求解方法,求出系统的响应。
但是对于高阶系统,手工计算这一问题的过程非常困难和繁琐。
在MATLAB 中,应用lsim( )函数很容易就能对上述微分方程所描述的系统的响应进行仿真,求出系统在任意激励信号作用下的响应。
lsim( )函数不仅能够求出连续系统在指定的任意时间范围内系统响应的数值解,而且还能同时绘制出系统响应的时域波形图。
以上各函数的调用格式如下: ⑴ impulse( ) 函数函数impulse( )将绘制出由向量a 和b 所表示的连续系统在指定时间范围内的单位冲激响应h (t )的时域波形图,并能求出指定时间范围内冲激响应的数值解。
impulse(b,a) 以默认方式绘出由向量a 和b 所定义的连续系统的冲激响应的时域波形。
impulse(b,a ,t0) 绘出由向量a 和b 所定义的连续系统在0 ~ t0时间范围内冲激响应的时域波形。
impulse(b,a,t1:p:t2) 绘出由向量a 和b 所定义的连续系统在t1 ~ t2时间范围内,并且以时间间隔p 均匀取样的冲激响应的时域波形。
y=impulse(b,a,t1:p:t2) 只求出由向量a 和b 所定义的连续系统在t1 ~ t2时间范围内,并且以时间间隔p 均匀取样的冲激响应的数值解,但不绘出其相应波形。
⑵ step( ) 函数函数step( )将绘制出由向量a 和b 所表示的连续系统的阶跃响应,在指定的时间范围内的波形图,并且求出数值解。
和impulse( )函数一样,step( )也有如下四种调用格式:step( b,a) step(b,a,t0)step(b,a,t1:p:t2) y=step(b,a,t1:p:t2)上述调用格式的功能和impulse( )函数完全相同,所不同只是所绘制(求解)的是系统的阶跃响应g (t ),而不是冲激响应h (t )。
⑶ lsim( )函数根据系统有无初始状态,lsim( )函数有如下两种调用格式:①系统无初态时,调用lsim( )函数可求出系统的零状态响应,其格式如下:lsim(b,a,x,t) 绘出由向量a 和b 所定义的连续系统在输入为x 和t 所定义的信号时,系统零状态响应的时域仿真波形,且时间范围与输入信号相同。
其中x 和t 是表示输入信号的行向量,t 为表示输入信号时间范围的向量,x 则是输入信号对应于向量t 所定义的时间点上的取样值。
y=lsim(b,a,x,t) 与前面的impulse 和step 函数类似,该调用格式并不绘制出系统的零状态响应曲线,而只是求出与向量t 定义的时间范围相一致的系统零状态响应的数值解。
②系统有初始状态时,调用lsim( )函数可求出系统的全响应,格式如下:lsim(A,B,C,D,e,t,X0) 绘出由系数矩阵A,B,C,D 所定义的连续时间系统在输入为e 和t 所定义的信号时,系统输出函数的全响应的时域仿真波形。
t 为表示输入信号时间范围的向量,e 则是输入信号e(t)对应于向量t 所定义的时间点上的取样值,X0表示系统状态变量X=[x1,x2,…..xn]'在t=0时刻的初值。
[Y,X]= lsim(A,B,C,D,e,t,X0) 不绘出全响应波形,而只是求出与向量t 定义的时间范围相一致的系统输出向量Y 的全响应以及状态变量X 的数值解。
显然,函数lsim( )对系统响应进行仿真的效果取决于向量t 的时间间隔的密集程度,t 的取样时间间隔越小则响应曲线越光滑,仿真效果也越好。
说明:(1)当系统有初始状态时,若使用lsim( )函数求系统的全响应,就要使用系统的状态空间描述法,即首先要根据系统给定的方式,写出描述系统的状态方程和输出方程。
假如系统原来给定的是微分方程或系统函数,则可用相变量法或对角线变量等方法写出系统的状态方程和输出方程。
其转换原理如前面实验四所述。
(2)显然利用lsim( )函数不仅可以分析单输入单输出系统,还可以分析复杂的多输入多输出系统。
例题1: 若某连续系统的输入为e (t ),输出为r (t ),系统的微分方程为:''()5'()6()3'()2()y t y t y t f t f t ++=+①求该系统的单位冲激响应h (t )及其单位阶跃响应g (t )。
②若2()()tf t e t ε-= 求出系统的零状态响应y(t )分析: ① 求冲激响应及阶跃响应的MATLAB 程序:a=[1 5 6];b=[3 2];subplot(2,1,1), impulse(b,a,4) subplot(2,1,2), step(b,a,4)运行结果如右:② 求零状态响应的MATLAB 程序:a=[1 5 6];b=[3 2];p1=0.01; %定义取样时间间隔为0.01 t1=0:p1:5; %定义时间范围 x1=exp(-2*t1); %定义输入信号lsim(b,a,x1,t1), %对取样间隔为0.01时系统响应进行仿真hold on; %保持图形窗口以便能在同一窗口中绘制多条曲线 p2=0.5; %定义取样间隔为0.5 t2=0:p2:5; %定义时间范围 x2=exp(-2*t2); %定义输入信号lsim(b,a,x2,t2), hold off %对取样间隔为0.5时系统响应进行仿真并解除保持运行结果如下:例题2 已知一个过阻尼二阶系统的状态方程和输出方程分别为:010'()()()232x t X t f t ⎡⎤⎡⎤=+⎢⎥⎢⎥--⎣⎦⎣⎦, r (t )=[0 1]X (t ) 。
若系统初始状态为X (0)=[4 -5]T , 求系统在4()3()tf t e t ε-=作用下的全响应。
求全响应程序如下:A=[0 1 ; -2 -3] ;B=[0 2]';C=[0 1];D=[0];X0=[4 -5]'; %定义系统初始状态 t=0: 0.01:10;E =[3*exp(-4*t).*ones(size(t))]'; %定义系统激励信号[r , x]=lsim(A,B,C,D,E,t,X0); %求出系统全响应的数值解 plot(t,r) %绘制系统全响应波形 运行结果如右。
2.离散时间系统LTI 离散系统中,其输入和输出的关系由差分方程描述:00()()n mi ji j a y k i bf k j ==+=+∑∑ (前向差分方程)()()nmi ji j a y k i bf k n j ==-=-+∑∑ (后向差分方程)当系统的输入为单位序列δ(k )时产生的零状态响应称为系统的单位函数响应,用h (k )表示。
当输入为 ε(k )时产生的零状态响应称为系统的单位阶跃应,记为:g (k ),如下图所示。
如果系统输入为e (k ),冲激响应为h (k ),系统的零状态响应为y(k ),则有:()()()y k h k f k =*。
与连续系统的单位冲激响应h (t )相类似,离散系统的单位函数响应h (k )也包含了系统的固有特性,与输入序列无关。
我们只要知道了系统的单位函数响应,即可求得系统在不同激励信号作用下产生的响应。
因此,求解系统的单位函数响应h (k )对我们进行离散系统的分析也同样具有非常重要的意义。
MATLAB 中为用户提供了专门用于求解离散系统单位函数响应, 并绘制其时域波形的函数impz( )。
同样也提供了求离散系统响应的专用函数filter( ),该函数能求出由差分方程所描述的离散系统在指定时间范围内的输入序列作用时,产生的响应序列的数值解。
当系统初值不为零时,可以使用dlsim( )函数求出离散系统的全响应,其调用方法与前面连续系统的lsim( )函数相似。
另外,求解离散系统阶跃响应可以通过如下两种方法实现:一种是直接调用专用函数dstep( ),其调用方法与求解连续系统阶跃响应的专用函数step( )的调用方法相似;另一种方法是利用求解离散系统零状态响应的专用函数filter( ),只要将其中的激励信号看成是单位阶跃信号ε(k )即可。
函数的调用格式分别如下: ⑴ impz( )函数impz(b,a) 以默认方式绘出由向量a 和b 所定义的离散系统单位函数响应的时域波形。
impz(b,a,n) 绘出由向量a 和b 所定义的离散系统在0 ~ n (n 必须为整数)的离散时间范围内单位函数响应的时域波形。
impz(b,a,n1:n2) 绘出由向量a 和b 所定义的离散系统在n1 ~ n2 (n1、n2必须为整数)的离散时间范围内单位函数响应的时域波形。
y=impz(b,a,n1:n2) 求出由向量a 和b 所定义的离散系统在n1 ~ n2 (n1、n2必须为整数)的离散时间范围内单位函数响应的数值解,但不绘出波形。
⑵ filter( ) 函数filter(b,a,x) 其中a 和b 与前面相同,x 是包含输入序列非零样值点的的行向量。
此命令将求出系统在与x 的取样时间点相同的输出序列样值。
例题:已知描述离散系统的差分方程为:()0.25(1)0.5(2)()(1)y k y k y k f k f k --+-=+-,且已知系统输入序列为12()()()kf k k ε=,① 求出系统的单位函数响应h (k )在-3 ~10离散时间范围内响应波形。