信号与系统实验报告-实验3--周期信号的频谱分析
- 格式:doc
- 大小:177.00 KB
- 文档页数:13
实验三信号的频谱分析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. 方波信号的分解方波信号的分解的基本工作原理是采用多个带通滤波器,把它们的中心频率分别调到被测信号的各个频率分量上,当被测信号同时加到多路滤波器上,中心频率与信号所包含的某次谐波分量频率一致的滤波器便有输出。
实验报告一、实验目的和要求谱分析即求信号的频谱。
本实验采用DFT/FFT技术对周期性信号进行谱分析。
通过实验,了解用X(k)近似地表示频谱X(ejω)带来的栅栏效应、混叠现象和频谱泄漏,了解如何正确地选择参数(抽样间隔T、抽样点数N)。
二、实验内容和步骤2-1 选用最简单的周期信号:单频正弦信号、频率f=50赫兹,进行谱分析。
2-2 谱分析参数可以从下表中任选一组(也可自定)。
对各组参数时的序列,计算:一个正弦周期是否对应整数个抽样间隔?观察区间是否对应整数个正弦周期?2-3 对以上几个正弦序列,依次进行以下过程。
2-3-1 观察并记录一个正弦序列的图形(时域)、频谱(幅度谱、频谱实部、频谱虚部)形状、幅度谱的第一个峰的坐标(U,V)。
2-3-2 分析抽样间隔T、截断长度N(抽样个数)对谱分析结果的影响;2-3-3 思考X(k)与X(e jω)的关系;2-3-4 讨论用X(k)近似表示X(ejω)时的栅栏效应、混叠现象、频谱泄漏。
三、主要仪器设备MATLAB编程。
四、操作方法和实验步骤(参见“二、实验内容和步骤”)五、实验数据记录和处理clc;clf;clear;%清除缓存%第一组数据的MATLAB程序(之后几组只需要将参数改变即可) T=0.000625;length=32;n=0:length-1;t=0:0.0001:31;%原序列和采样序列xn=sin(2*pi*50*n*T);xt=sin(2*pi*50*t);%画第一幅图(原序列和采样序列)figure(1);subplot(2,1,1);plot(t,xt);xlabel('t');ylabel('xt');axis([0,0.2,-1.1,1.1]);title('原序列时域');subplot(2,1,2);stem(n,xn ,'filled');xlabel('n');ylabel('xn');axis([0,length,-1.1,1.1]);title('采样后序列时域');%画第二幅图(采样序列实部、虚部、模和相角)figure(2);subplot(2,2,1);stem(n,real(xn) ,'filled');xlabel('n');ylabel('real(xn)');axis([0,length,-1.1,1.1]);title('采样序列的实部');subplot(2,2,2);stem(n,imag(xn) ,'filled');xlabel('n');ylabel('imag(xn)');axis([0,length,-1.1,1.1]);title('采样序列的虚部');subplot(2,2,3);stem(n,abs(xn) ,'filled');xlabel('n');ylabel('abs(xn)');axis([0,length,-1.1,1.1]);title('采样序列的模');subplot(2,2,4);stem(n,angle(xn) ,'filled');xlabel('n');ylabel('angle(xn)');axis([0,length,-(pi+0.5),pi+0.5]);title('采样序列的相角');%计算DFTDFT=fft(xn,length);%画第三幅图(DFT的幅度、实部和虚部)figure(3);subplot(3,1,1);stem(n,abs(DFT) ,'filled');xlabel('k');%DFT后的频域变量为kylabel('abs(DFT)');title('DFT 幅度谱');subplot(3,1,2);stem(n,real(DFT) ,'filled');xlabel('k');ylabel('real(DFT)');title('DFT的实部');subplot(3,1,3);stem(n,imag(DFT) ,'filled');xlabel('k');ylabel('imag(DFT)');title('DFT的虚部');六、实验结果与分析实验结果:第一组数据:实验名称:DFT/FFT的应用之一 确定性信号谱分析姓名:张清学号:3110103952 P.4第二组数据:第三组数据:第四组数据:第五组数据:第六组数据:6-1 实验前预习有关概念,并根据上列参数来推测相应频谱的形状、谱峰所在频率(U)和谱峰的数值(V)、混叠现象和频谱泄漏的有无。
信号与系统实验报告一、实验目的(1) 理解周期信号的傅里叶分解,掌握傅里叶系数的计算方法;(2)深刻理解和掌握非周期信号的傅里叶变换及其计算方法;(3) 熟悉傅里叶变换的性质,并能应用其性质实现信号的幅度调制;(4) 理解连续时间系统的频域分析原理和方法,掌握连续系统的频率响应求解方法,并画出相应的幅频、相频响应曲线。
二、实验原理、原理图及电路图(1) 周期信号的傅里叶分解设有连续时间周期信号()f t ,它的周期为T ,角频率22fT,且满足狄里赫利条件,则该周期信号可以展开成傅里叶级数,即可表示为一系列不同频率的正弦或复指数信号之和。
傅里叶级数有三角形式和指数形式两种。
1)三角形式的傅里叶级数:01212011()cos()cos(2)sin()sin(2)2cos()sin()2n n n n a f t a t a t b t b t a a n t b n t 式中系数n a ,n b 称为傅里叶系数,可由下式求得:222222()cos(),()sin()T T T T nna f t n t dtb f t n t dtTT2)指数形式的傅里叶级数:()jn tn nf t F e式中系数n F 称为傅里叶复系数,可由下式求得:221()T jn tT nF f t edtT周期信号的傅里叶分解用Matlab进行计算时,本质上是对信号进行数值积分运算。
Matlab中进行数值积分运算的函数有quad函数和int函数。
其中int函数主要用于符号运算,而quad函数(包括quad8,quadl)可以直接对信号进行积分运算。
因此利用Matlab进行周期信号的傅里叶分解可以直接对信号进行运算,也可以采用符号运算方法。
quadl函数(quad系)的调用形式为:y=quadl(‘func’,a,b)或y=quadl(@myfun,a,b)。
其中func是一个字符串,表示被积函数的.m文件名(函数名);a、b分别表示定积分的下限和上限。
信号与系统实验报告-实验3--周期信号的频谱分析信号与系统实验报告实验三周期信号的频谱分析实验三周期信号的频谱分析实验目的:1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法;2、观察截短傅里叶级数而产生的“Gibbs现象”,了解其特点以及产生的原因;3、掌握各种典型的连续时间非周期信号的频谱特征。
实验内容:(1)Q3-1 编写程序Q3_1,绘制下面的信号的波形图:其中,0 = 0.5π,要求将一个图形窗口分割成四个子图,分别绘制cos(0t)、cos(30t)、cos(50t) 和x(t) 的波形图,给图形加title,网格线和x坐标标签,并且程序能够接受从键盘输入的和式中的项数。
程序如下:clear,%Clear all variablesclose all,%Close all figure windowsdt = 0.00001; %Specify the step of time variable t = -2:dt:4; %Specify the interval of timew0=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;endsubplot(221)plot(t,x1)%Plot x1axis([-2 4 -2 2]);grid on,title('signal cos(w0.*t)')subplot(222)plot(t,x2)%Plot x2axis([-2 4 -2 2]); grid on,title('signal cos(3*w0.*t))')subplot(223)plot(t,x3)%Plot x3axis([-2 4 -2 2])grid on,title('signal cos(5*w0.*t))')subplot(224)plot(t,x)%Plot xtaxis([-2 4 -2 2])grid on,title('signal xt')(2)给程序3_1增加适当的语句,并以Q3_2存盘,使之能够计算例题1中的周期方波信号的傅里叶级数的系数,并绘制出信号的幅度谱和相位谱的谱线图。
实验一 抽样定理与信号恢复一、实验目的1. 观察离散信号频谱,了解其频谱特点;2. 验证抽样定理并恢复原信号。
二、实验原理1. 离散信号不仅可从离散信号源获得,而且也可从连续信号抽样获得。
抽样信号 Fs (t )=F (t )·S (t )。
其中F (t )为连续信号(例如三角波),S (t )是周期为Ts 的矩形窄脉冲。
Ts 又称抽样间隔,Fs=1Ts 称抽样频率,Fs (t )为抽样信号波形。
F (t )、S (t )、Fs (t )波形如图1-1。
t-4T S -T S 0T S 4T S8T S 12T S tt02/1τ1τ2/31τ2/1τ1τ2/31τ2/1τ-(a)(b)(c)图1-1 连续信号抽样过程将连续信号用周期性矩形脉冲抽样而得到抽样信号,可通过抽样器来实现,实验原理电路如图1-2所示。
2. 连续周期信号经周期矩形脉冲抽样后,抽样信号的频谱()∑∞∞--∙=m s s m m SaTsA j )(22s F ωωπδτωτω 它包含了原信号频谱以及重复周期为fs (f s =πω2s 、幅度按ST A τSa (2τωs m )规律变化的原信号频谱,即抽样信号的频谱是原信号频谱的周期性延拓。
因此,抽样信号占有的频带比原信号频带宽得多。
以三角波被矩形脉冲抽样为例。
三角波的频谱 F (j ω)=∑∞-∞=-K k k sa E )2()2(12τπωδππ抽样信号的频谱Fs (j ω)=式中 取三角波的有效带宽为31ω18f f s =作图,其抽样信号频谱如图1-3所示。
图1-2 信号抽样实验原理图)(2(212s m k s m k k Sa m Sa TS EA ωωωδπτωτπ--∙∙∑∞-∞=-∞=111112ττπω==f 或(b) 抽样信号频谙图1-3 抽样信号频谱图如果离散信号是由周期连续信号抽样而得,则其频谱的测量与周期连续信号方法相同,但应注意频谱的周期性延拓。
信号频谱分析实验报告信号频谱分析实验报告引言:信号频谱分析是一种重要的信号处理技术,可以帮助我们了解信号的频率成分和能量分布情况。
通过对信号频谱的分析,我们可以更好地理解信号的特性,并在实际应用中进行优化和改进。
本实验旨在通过实际操作,探究信号频谱分析的原理和方法。
实验设备和步骤:实验中我们使用了信号发生器、示波器和频谱分析仪作为主要设备。
首先,我们将信号发生器连接到示波器,通过调节信号发生器的频率和幅度,产生不同特性的信号。
然后,将示波器的输出信号连接到频谱分析仪上,通过频谱分析仪对信号进行频谱分析。
在实验过程中,我们记录了不同信号频谱的变化情况,并进行了数据的整理和分析。
实验结果:在实验中,我们产生了多种不同频率和幅度的信号,并对其进行了频谱分析。
通过观察频谱图,我们可以清晰地看到不同频率成分的能量分布情况。
实验结果表明,信号的频谱在不同频率范围内具有不同的能量分布,且能量峰值对应着信号的主要频率成分。
此外,我们还观察到信号的幅度对频谱的形态有着重要影响,幅度较大的信号在频谱图上表现出更强的峰值。
讨论与分析:通过对实验结果的分析,我们可以得出以下几点结论:1. 信号频谱分析可以帮助我们了解信号的频率成分和能量分布情况,从而更好地理解信号的特性。
2. 不同频率成分的能量分布情况在频谱图上呈现为峰值,峰值对应着信号的主要频率成分。
3. 信号的幅度对频谱的形态有着重要影响,幅度较大的信号在频谱图上表现出更强的峰值。
4. 通过对信号频谱的分析,我们可以优化和改进信号的特性,以满足实际应用的需求。
实验的局限性和改进方向:在本实验中,我们只使用了简单的信号发生器和示波器进行频谱分析,实验结果可能受到设备本身的限制。
为了更准确地分析信号的频谱,可以考虑使用更高精度的频谱分析仪和信号源。
此外,我们在实验中只观察了信号频谱的静态特性,对于动态信号的分析还需要进一步研究。
结论:通过本次实验,我们深入了解了信号频谱分析的原理和方法,并通过实际操作获得了实验结果。
一、实验目的1. 理解信号频谱分析的基本原理和重要性。
2. 掌握使用MATLAB进行信号频谱分析的方法和步骤。
3. 通过实验验证不同信号类型(如连续信号、离散信号)的频谱特性。
4. 学习如何利用频谱分析进行信号处理和滤波。
二、实验原理信号频谱分析是将信号从时域转换到频域的一种方法,它可以帮助我们了解信号的频率成分、幅度分布和相位特性。
常见的频谱分析方法包括傅里叶变换(FT)、快速傅里叶变换(FFT)等。
傅里叶变换将信号分解为不同频率的正弦波和余弦波的线性组合,从而揭示信号的频率成分。
FFT是一种高效的傅里叶变换算法,广泛应用于信号处理领域。
三、实验仪器与软件1. 仪器:信号发生器、示波器、计算机2. 软件:MATLAB四、实验步骤1. 信号生成:使用信号发生器生成不同的信号,如正弦波、方波、三角波等。
2. 信号采集:使用示波器采集信号的时域波形,并将数据导入MATLAB进行后续处理。
3. 频谱分析:- 使用MATLAB的FFT函数对采集到的信号进行傅里叶变换。
- 绘制信号的频谱图,观察信号的频率成分、幅度分布和相位特性。
4. 滤波:- 根据实验需求,设计合适的滤波器(如低通、高通、带通等)。
- 对信号进行滤波处理,观察滤波效果。
5. 结果分析:- 分析不同信号类型的频谱特性,如正弦波、方波、三角波等。
- 分析滤波器对信号的影响,如信号失真、噪声抑制等。
五、实验结果与分析1. 正弦波频谱分析:- 正弦波的频谱只有一个频率成分,即其本身频率。
- 频谱图上,该频率处的幅度为最大值,其余频率处的幅度为零。
2. 方波频谱分析:- 方波的频谱包含多个频率成分,包括基波及其整数倍谐波。
- 频谱图上,基波频率处的幅度最大,谐波频率处的幅度逐渐减小。
3. 三角波频谱分析:- 三角波的频谱包含基波及其整数倍谐波。
- 频谱图上,基波频率处的幅度最大,谐波频率处的幅度逐渐减小,且衰减速度比方波慢。
4. 滤波效果分析:- 滤波器可以有效抑制不需要的频率成分,保留需要的频率成分。
信号与系统实验报告一、信号的时域基本运算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,即全状态=零输入+零状态。
两实验之二心得体会:求差分方程时,可以根据全状态响应是由零输入输入以及零状态相加所得,分开来求,同时也加深了自己对差分方程的求解问题的理解。
一,实验目的四,心得体会了解信号频谱和信号频域,掌握其特性。
一,实验原理实验主要分为四个部分,分别分析了连续和离散信号的周期、非周期情况下特性。
1.连续周期信号的频谱分析首先手算出信号的傅里叶级数,得出信号波形,然后通过代码画出信号波形图。
2.连续非周期信号的频谱分析先由非周期信号的时域信号得到它的频谱X(w),再通过MATLAB求出其傅里叶变换并绘出图形。
X=fourier(x)x=ifourier(x)①符号运算法syms t②数值积分法quad(fun,a,b)③数值近似法3.离散周期信号的频谱分析X=fft(x)4.离散非周期信号的频谱分析可以化为两个相乘的矩阵,从而由MATLAB实现。
三,实验内容(1)已知x(t)是如图周期矩形脉冲信号。
1).计算该信号的傅里叶级数。
2).利用MATLAB绘出由前N次谐波合成的信号波形,观察随着N的变化合成信号波形的变化规律。
3).利用MATLAB绘出周期矩形脉冲信号的频谱,观察参数T和τ变化时对频谱波形的影响。
思考下列问题:①什么是吉伯斯现象?产生吉伯斯现象的原因是什么?②以周期矩形脉冲信号为例,说明周期信号的频谱有什么特点。
③周期矩形脉冲信号参数τ/T的变化,其频谱结构(如频谱包络形状、过零点、频谱间隔等)如何变化?(2)已知x(t)是如图所示矩形脉冲信号。
1).求该信号的傅里叶变幻。
2). 利用MATLAB绘出周期矩形脉冲信号的频谱,观察参数T和τ变化时对频谱波形的影响。
3). 让矩形脉冲宽度始终等于一,改变矩形脉冲宽度,观察矩形脉冲信号时域波形和频谱随矩形脉冲宽度的变化趋势。
①比较矩形脉冲信号和周期矩形脉冲信号的频谱,两者之间有何异同。
②让矩形脉冲的面积始终等于一,改变矩形脉冲的宽度,观察矩形脉冲信号时域波形和频谱波形随矩形脉冲宽度的变化趋势。
(1)已知x(t)是如图所示的周期矩形脉冲信号①,计算该信号的傅里叶级数答:由图中x(t)波形可知信号为通过计算,可以知道所以x(t)的傅里叶级数为。
学院:电子工程学院班级:姓名:学号:信号与系统测试实验总结当前,科学技术都向两极化发展,既向微观发展又向宏观发展。
各学科之间既高度综合又高度分化。
这就要求了我们当代的大学生既要有坚实的理论基础,又还必须具备极强的动手能力和解决各种实际问题的能力。
而这个学期所开展的信号与系统测试的实验课程给我提供了一个很好的机会和平台。
本学期的第八周到十二周期间,我们有幸做了四次信号与系统测试实验。
这四次实验分别为:信号的分类与观察、非正弦周期信号的频谱分析、信号的抽样与恢复(PAM)和模拟滤波器实验。
通过四次印象深刻的实验,不仅在理论上加深了我的理论概念知识,更是通过实践锻炼我们的动手能力,学会使用示波器、信号发生器、频谱仪、信号与系统试验箱等实验仪器。
由于第一次做实验,所以对于实验室里面的很多仪器都感到很新奇,给我留下的印象也很是深刻。
我们目前学的信号与系统基本都是一些数学理论内容,实在是抽象的紧,缺乏和实际的联系。
而这门课程的关键就是在于补充这一方面,让我们的理论和实际得到一定印证。
第一次实验中,实验过程比较简单,稍微复杂的是在于函数图像的绘制上。
而实验之后的理论计算则是让我费了一番功夫,这也让我体会到了理论和实际结合的重要性。
第二个实验是非正弦信号的频谱分析。
在这次试验中,我们接触到了频谱仪这个很重要的工具。
这一次实验中,实验的操作很简单,但实验的原理倒是颇为复杂,这也导致我们进行理论计算时十分的麻烦。
第三次实验做了信号的抽样与恢复。
这是一个很有意义的实验,它向我们展示了现代通信技术的基础,也正是它才使得信息可以有效地传递。
这次实验,我们主要通过矩形脉冲对正弦信号进行抽样,再把它还原回来,最后用还原的图形与原图形对比,分析实验并总结。
试验中,抽样后的波形不稳定,很难根据示波器上的图形进行图形描绘,老师便告诉了我们一个办法,即用手机把图形拍下来再进行绘制,这一环节中,老师如果不醒那么早,让我们自己去思考寻找解决办法将是更好的一个考验。
信号与系统分析实验报告信号与系统分析实验报告引言:信号与系统分析是电子工程领域中的重要课程之一,通过实验可以更好地理解信号与系统的基本概念和原理。
本实验报告将对信号与系统分析实验进行详细的描述和分析。
实验一:信号的采集与重构在这个实验中,我们学习了信号的采集与重构。
首先,我们使用示波器采集了一个正弦信号,并通过数学方法计算出了信号的频率和幅值。
然后,我们使用数字信号处理器对采集到的信号进行重构,并与原始信号进行比较。
实验结果表明,重构后的信号与原始信号非常接近,证明了信号的采集与重构的有效性。
实验二:线性系统的时域响应本实验旨在研究线性系统的时域响应。
我们使用了一个线性系统,通过输入不同的信号,观察输出信号的变化。
实验结果显示,线性系统对于不同的输入信号有不同的响应,但都遵循线性叠加的原则。
通过分析输出信号与输入信号的关系,我们可以得出线性系统的传递函数,并进一步研究系统的稳定性和频率响应。
实验三:频域特性分析在这个实验中,我们研究了信号的频域特性。
通过使用傅里叶变换,我们将时域信号转换为频域信号,并观察信号的频谱。
实验结果显示,不同频率的信号在频域上有不同的分布特性。
我们还学习了滤波器的设计和应用,通过设计一个低通滤波器,我们成功地去除了高频噪声,并得到了干净的信号。
实验四:系统辨识本实验旨在研究系统的辨识方法。
我们使用了一组输入信号和对应的输出信号,通过数学建模的方法,推导出了系统的传递函数。
实验结果表明,通过系统辨识可以准确地描述系统的特性,并为系统的控制和优化提供了基础。
结论:通过本次实验,我们深入学习了信号与系统分析的基本概念和原理。
实验结果证明了信号的采集与重构的有效性,线性系统的时域响应的线性叠加原则,信号的频域特性和滤波器的设计方法,以及系统辨识的重要性。
这些知识和技能对于我们理解和应用信号与系统分析具有重要的意义。
通过实验的实际操作和分析,我们对信号与系统的理论有了更深入的理解,为我们今后的学习和研究打下了坚实的基础。
实验一典型连续时间信号和离散时间信号一、实验目的掌握利用Matlab画图函数和符号函数显示典型连续时间信号波形、典型时间离散信号、连续时间信号在时域中的自变量变换。
二、实验内容1、典型连续信号的波形表示(单边指数信号、复指数信号、抽样信号、单位阶跃信号、单位冲击信号)1)画出教材P28习题1-1(3) ()[(63)(63)]t=----的波形图。
f t e u t u t2)画出复指数信号()()j t f t e σω+=当0.4, 8σω==(0<t<10)时的实部和虚部的波形图。
t=0:0.01:10;f1='exp(0.4*t)*cos(8*t)';f2='exp(0.4*t)*sin(8*t)';figure(1)ezplot(f1,t);grid on;figure(2)ezplot(f2,t);grid on;3)画出教材P16图1-18,即抽样信号Sa(t)的波形(-20<t<20)。
t=-10:0.01:10;f='sin(t)/t';ezplot(f,t);grid on;4)用符号函数sign画出单位阶跃信号u(t-3)的波形(0<t<10)。
t=0:0.01:10;f='(sign(t-3)+1)/2';ezplot(f,t);grid on;5)单位冲击信号可看作是宽度为∆,幅度为1/∆的矩形脉冲,即t=t 1处的冲击信号为11111 ()()0 t t t x t t t otherδ∆⎧<<+∆⎪=-=∆⎨⎪⎩画出0.2∆=, t 1=1的单位冲击信号。
t=0:0.01:2;f='5*(u(t-1)-u(t-1.2))';ezplot(f,t);grid on;axis([0 2 -1 6]);2、典型离散信号的表示(单位样值序列、单位阶跃序列、实指数序列、正弦序列、复指数序列)编写函数产生下列序列:1)单位脉冲序列,起点n0,终点n f,在n s处有一单位脉冲。
信号与系统实验报告-实验3--周期信号的频谱分析信号与系统实验报告实验三周期信号的频谱分析实验三周期信号的频谱分析实验目的:1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法;2、观察截短傅里叶级数而产生的“Gibbs现象”,了解其特点以及产生的原因;3、掌握各种典型的连续时间非周期信号的频谱特征。
实验内容:(1)Q3-1 编写程序Q3_1,绘制下面的信号的波形图:其中,0 = 0.5π,要求将一个图形窗口分割成四个子图,分别绘制cos(0t)、cos(30t)、cos(50t) 和x(t) 的波形图,给图形加title,网格线和x坐标标签,并且程序能够接受从键盘输入的和式中的项数。
程序如下:clear,%Clear all variablesclose all,%Close all figure windowsdt = 0.00001; %Specify the step of time variable t = -2:dt:4; %Specify the interval of timew0=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;endsubplot(221)plot(t,x1)%Plot x1axis([-2 4 -2 2]);grid on,title('signal cos(w0.*t)')subplot(222)plot(t,x2)%Plot x2axis([-2 4 -2 2]); grid on,title('signal cos(3*w0.*t))')subplot(223)plot(t,x3)%Plot x3axis([-2 4 -2 2])grid on,title('signal cos(5*w0.*t))')subplot(224)plot(t,x)%Plot xtaxis([-2 4 -2 2])grid on,title('signal xt')(2)给程序3_1增加适当的语句,并以Q3_2存盘,使之能够计算例题1中的周期方波信号的傅里叶级数的系数,并绘制出信号的幅度谱和相位谱的谱线图。
程序如下:% Program3_1 clear, close allT = 2;dt = 0.00001;t = -2:dt:2;x1 = ut(t) - ut(t-1-dt);x = 0;for m = -1:1x = x + ut(t-m*T) - ut(t-1-m*T-dt);endw0 = 2*pi/T;N = 10;L = 2*N+1;for k = -N: N;ak(N+1+k) = (1/T)*x1*exp(-j*k*w0*t')*dt; endphi = angle(ak);subplot(211)'k = -10:10;stem (k,abs(ak),'k');axis([-10,10,0,0.6]);grid on;title('fudupu');subplot(212);k = -10:10stem(k,angle(ak),'k');axis([-10,10,-2,2]);grid on;titie('xiangweipu');xlabel('Frequency index x');(3)反复执行程序Program3_2,每次执行该程序时,输入不同的N值,并观察所合成的周期方波信号。
通过观察,你了解的吉伯斯现象的特点是:程序如下:clear,close allT = 2;dt = 0.00001;t = -2:dt:2;x1 = ut(t)-ut(t-1-dt);x = 0; for m = -1:1x = x + ut(t-m*T) - ut(t-1-m*T-dt);endw0 = 2*pi/T;N = input('Type in the number of the harmonic components N = :');L = 2*N+1;for k = -N:1:N;ak(N+1+k) = (1/T)*x1*exp(-j*k*w0*t')*dt; endphi = angle(ak);y=0;for q = 1:L;y = y+ak(q)*exp(j*(-(L-1)/2+q-1)*2*pi*t/T); end;subplot(221),plot(t,x),title('The original signal x(t)'),axis([-2,2,-0.2,1.2]),subplot(223),plot(t,y),title('The synthesis signal y(t)'),axis([-2,2,-0.2,1.2]),xlabel('Time t'),subplot(222)k=-N:N;stem(k,abs(ak),'k.'),title('The amplitude |ak| of x(t)'),axis([-N,N,-0.1,0.6])subplot(224)stem(k,phi,'r.'),title('The phase phi(k) of x(t)'),axis([-N,N,-2,2]),xlabel('Index k')N=1N=3通过观察我们了解到:如果一个周期信号在一个周期有内断点存在,那么,引入的误差将除了产生纹波之外,还将在断点处产生幅度大约为9%的过冲(Overshot),这种现象被称为吉伯斯现象(Gibbs phenomenon)。
即信号在不连续点附近存在一个幅度大约为9%的过冲,且所选谐波次数越多,过冲点越向不连续点靠近。
(4)计算如图的傅里叶级数的系数程序如下:clc,clear,close allT=2;dt=0.00001;t=-3:dt:3;x=(t+1).*(u(t+1)-u(t))-(t-1).*(u(t)-u(t-1));x1=0; for m=-2:2x1=x1+(t+1-m*T).*(u(t+1-m*T)-u(t-m*T))-(t-1-m*T).*(u(t-m *T)-u(t-1-m*T));endw0=2*pi/T;N=10;L=2*N+1;for k=-N:N;ak(N+1+k)=(1/T)*x*exp(-j*k*w0*t')*dt;endphi=angle(ak);plot(t,x1);axis([-4 4 0 1.2]);grid on;title('The signal x1(t)'); xlabel('Time t (sec)');ylabel('signal x1(t)');(5)仿照程序3_1,编写程序Q3_5,以计算x2(t) 的傅里叶级数的系数(不绘图)。
程序如下:clc,clear,close allT=2;dt=0.00001;t=-3:dt:3;x=ut(t+0.2)-ut(t-0.2-dt);x2=0;for m=-1:1x2=x2+ut(t+0.2-m*T)-ut(t-0.2-m*T)-ut(t-0.2-m*t-dt); endw0=2*pi/T;N=10;L=2*N+1for k=-N:N;ak(N+1+k)=(1/T)*x*exp(-j*k*w0*t')*dt;endphi=angle(ak);plot(t,x2);axis([-2.5 2.5 0 1.2]);grid on;title('The signal x2(t)');xlabel('Time t (sec)');ylabel('signal x2(t)');(6)仿照程序3_2,编写程序Q3_6,计算并绘制出原始信号x1(t) 的波形图,用有限项级数合成的y1(t) 的波形图,以及x1(t) 的幅度频谱和相位频谱的谱线图。
程序如下:clc,clear,close allT=2;dt=0.00001;t=-3:dt:3;x=(t+1).*(ut(t+1)-ut(t))-(t-1).*(ut(t)-ut(t-1));x1=0;for m=-2:2x1=x1+(t+1-m*T).*(ut(t+1-m*T)-ut(t-m*T))-(t-1-m*T).*(ut( t-m*t)-ut(t-1-m*t));endw0=2*pi/T;N=10;L=2*N+1;for k=-N:N;ak(N+1+k)=(1/T)*x*exp(-j*k*w0*t')*dt; endphi=angle(ak);y=0;for q=1:L;y=y+ak(q)*exp(j*(q-1-N)*w0*t);end;subplot(221)plot(t,x)%plot xaxis([-3 3 -0.2 1.2]);grid on;title('The original signal x(t)'); subplot(223)plot(t,y)%Plot yaxis([-3 3 -0.2 1.2]);grid on;title('The synthesis signal y(t)'); subplot(222);xlabel('Time i (sec)');subplot(222);实验心得:通过这次实验,了解了连续时间周期信号的傅里叶级数的物理意义,观察了截短傅里叶级数而产生的“Gibbs现象”,了解其特点以及产生的原因,了解掌握了各种典型的连续时间非周期信号的频谱特征。
从开始的不了解,到后来通过看书,上网查找资料做出这个实验,我学到了很多东西,虽然花费了不少时间,但是却获得了编程经验。