信号与系统实验报告-(2).doc

实验三常见信号得ATLAB表示及运算一、实验目得。熟悉常见信号得意义、特性及波形2.学会使用MATLAB表示信号得方法并绘制信号波形、掌握使用ATLAB进行信号基本运算得指令4、 熟悉用ATAB实现卷积积分得方法二、实验原理根据ALAB得数值计算功能与符号运算功能，在MATA中，信号有两种表示方法,一种就是用向量来表示,另一种则就是用符号运算得方法。在采用适当得MLAB语句表示出信号后，就可以利用MATA中得绘图命令绘制出直观得信号波形了。、连续时间信号 从严格意义上讲,ATL并不能处理连续信号。在ALB中，就是用连续信号在等时间间隔点上得样值来近似表示得,当取样时间间隔足够小时,这些离散得样值就能较好地近似出连续信号。在MATAB中连续信号可用向量或符号运算功能来表示。向量表示法对于连续时间信号，可以用两个行向量f与t来表示,其中向量t就是用形如得命令定义得时间范围向量,其中，为信号起始时间，为终止时间,p为时间间隔。向量f为连续信号在向量所定义得时间点上得样值符号运算表示法如果一个信号或函数可以用符号表达式来表示,那么我们就可以用前面介绍得符号函数专用绘图命令ezplot(）等函数来绘出信号得波形。常见信号得MTLA表示单位阶跃信号单位阶跃信号得定义为:方法一： 调用Hviside(t)函数首先定义函数Heisd（) 得函数文件,该文件名应与函数名同名即eaviside、m.定义函数文件，函数名为Havsi,输入变量为x,输出变量为function y=Haviside（t） y(t0）; %定义函数体,即函数所执行指令%此处定义0时y=1,t0时=0,注意与实际得阶跃信号定义得区别.方法二：数值计算法在ATLAB中,有一个专门用于表示单位阶跃信号得函数，即seun( )函数，它就是用数值计算法表示得单位阶跃函数其调用格式为:stpfun(t,t0)其中,t就是以向量形式表示得变量，t0表示信号发生突变得时刻,在以前,函数值小于零,以后函数值大于零。有趣得就是它同时还可以表示单位阶跃序列,这只要将自变量以及取样间隔设定为整数即可。符号函数符号函数得定义为:在MATLAB中有专门用于表示符号函数得函数sgn() ,由于单位阶跃信号e (t）与符号函数两者之间存在以下关系：,因此，利用这个函数就可以很容易地生成单位阶跃信号.2、离散时间信号离散时间信号又叫离散时间序列,一般用 表示，其中变量k为整数，代表离散得采样时间点（采样次数)。在MATLA中，离散信号得表示方法与连续信号不同，它无法用符号运算法来表示,而只能采用数值计算法表示,由于MATLB中元素得个数就是有限得，因此，MATAB无法表示无限序列；另外，在绘制离散信号时必须使用专门绘制离散数据得命令，即stem（)函数,而不能用lot(）函数。单位序列d单位序列）得定义为单位阶跃序列单位阶跃序列得定义为、卷积积分两个信号得卷积定义为:MATAB中就是利用conv函数来实现卷积得.功能：实现两个函数与得卷积.格式：=cnv(f1,2) 说明:=f1（t),f2=f2(t）表示两个函数,g(t)表示两个函数得卷积结果。三、实验内容1、分别用MATLA得向量表示法与符号运算功能，表示并绘出下列连续时间信号得波形: （1） t=:0、01:0;t1=-：、01:-0、1;2=：0、01:10；f1=zeros(1,lngth（t)），ne（1,lngh(t2);f=(2exp（*）、f1；po（t，f)axi（-,10,0,2、） sys t;=ym（(2-xp(t)*havisie（)；zo(f,-1,10）；(2)t=2:0、01:8；f=0、*(t&)+0、*（t4);lot(,)syms ;f=sym(os（i*t/）*havside(t)heavis（t4） ）;ezplt(f，-2，8);、分别用MATLAB表示并绘出下列离散时间信号得波形: （2） t=0：8；=10:1；=zeros(1，10),zeros(1,7);stem（1,f)axs(10，5,0，1);(3)t=0:5;t10:0；f=eos(1，10),sin(t*pi/);sem(t1,f）xis（10，50,2，)3、已知两信号，求卷积积分，并与例题比较。t1=1:0、01:0；t=0:0、01:1;31:0、01：；f1=s(ie(t1);f2nes(size（t2）);gco(1,f2)；splo(3,1,),po(t1,f1)；subplot(，1,),plot(t2，f2);subot(,1，3),lot(t3,);与例题相比较,(t)得定义域不同,最大值对应得横坐标也不同。4、已知，求两序列得卷积与 .N;=；L=N+M1；，1，2；f=1,2，3,4，;g=v(1,f2)；kf1=0:N-1；kf=0：1;kg=0：L;blot（1,1）,stem(kf，f1,*k);xabel（k)；yal（f1（k);grd onsbpt(1，3，2)，stem（f2，f2，*k)；xlabe（k);labl（f2(k);grd ubplot（1,3,3);ste（k，g，k)；xlabe(k)；ylabel（g（)）;grid n实验心得:第一次接触Mutlab这个绘图软件,觉得挺新奇得,同时 ,由于之前不太学信号与系统遇到一些不懂得问题，结合这些图对信号与系统有更好得了解。实验四 连续时间信号得频域分析一、实验目得。熟悉傅里叶变换得性质.熟悉常见信号得傅里叶变换。了解傅里叶变换得TAB实现方法二、实验原理从已知信号求出相应得频谱函数得数学表示为：傅里叶反变换得定义为：在MALA中实现傅里叶变换得方法有两种,一种就是利用MATLAB中得ybo Math Toox提供得专用函数直接求解函数得傅里叶变换与傅里叶反变换，另一种就是傅里叶变换得数值计算实现法.1、直接调用专用函数法在MATLAB中实现傅里叶变换得函数为:l F=furer( f ) 对(t)进行傅里叶变换,其结果为F(w)l=fourier(f,v) 对f（t)进行傅里叶变换,其结果为（v)l F=fourir( f,u，v ） 对(u)进行傅里叶变换,其结果为F()傅里叶反变换l f=ifourer（ F ) 对F(w)进行傅里叶反变换，其结果为f()lf=orir(F,U) 对（w)进行傅里叶反变换，其结果为（u)l fifoure( F,v，u ）对（)进行傅里叶反变换,其结果为f(）注意：（1)在调用函数fuier（ )及ifouri（ ）之前，要用syms命令对所有需要用到得变量(如t,u,v，w)等进行说明,即要将这些变量说明成符号变量。对ourer（ ）中得f及ifourier()中得F也要用符号定义符sm将其说明为符号表达式。(2)采用fuier( )及foure（ ）得到得返回函数,仍然为符号表达式。在对其作图时要用ezpl( ）函数,而不能用t（）函数.(3)uri(）及ourie（ ）函数得应用有很多局限性,如果在返回函数中含有()等函数，则ezplo（）函数也无法作出图来。另外,在用fourier( )函数对某些信号进行变换时,其返回函数如果包含一些不能直接表达得式子，则此时当然也就无法作图了。这就是ourir（ ）函数得一个局限。另一个局限就是在很多场合，尽管原时间信号f(t)就是连续得,但却不能表示成符号表达式，此时只能应用下面介绍得数值计算法来进行傅氏变换了,当然，大多数情况下，用数值计算法所求得频谱函数只就是一种近似值。2、傅里叶变换得数值计算实现法严格说来，如果不使用symboic工具箱,就是不能分析连续时间信号得。采用数值计算方法实现连续时间信号得傅里叶变换，实质上只就是借助于ATLAB得强大数值计算功能,特别就是其强大得矩阵运算能力而进行得一种近似计算。傅里叶变换得数值计算实现法得原理如下:对于连续时间信号f(t),其傅里叶变换为： 其中为取样间隔,如果f(t)就是时限信号,或者当|t大于某个给定值时，f(t)得值已经衰减得很厉害,可以近似地瞧成就是时限信号,则上式中得取值就就是有限得,假定为N，有: 若对频率变量进行取样,得:通常取：,其中就是要取得频率范围，或信号得频带宽度。采用MALAB实现上式时，其要点就是要生成f(t)得个样本值得向量，以及向量,两向量得内积(即两矩阵得乘积）,结果即完成上式得傅里叶变换得数值计算。注意：时间取样间隔得确定,其依据就是必须小于奈奎斯特(Nyquist)取样间隔。如果f（t）不就是严格得带限信号，则可以根据实际计算得精度要求来确定一个适当得频率为信号得带宽。三、实验内容1、编程实现求下列信号得幅度频谱（1) 求出得频谱函数F1（j)，请将它与上面门宽为2得门函数得频谱进行比较,观察两者得特点，说明两者得关系。(2) 三角脉冲 (3) 单边指数信号 (4） 高斯信号(1）sys t w tsym（Heaiside（*t+1）Heaviie(2*-）); F=furr(G,t,）； FF=mae(convet，Fw,pieewise)； FPabs(FF); eplot（F,10*pi10*p);grd; xs(-* 0pi 0 2、2) 与得频谱比较,得频谱函数F1(j)最大值就是其得2（)sm w;t=sym（1t）（Havie(+1)Hevid（t)）+(1-t）（eaviside（t)Heaiside（t1）);F=fourer(Gt，t，w）； Fw=aple(coet,Fw，picewe）; FFPab（w)； zpot（FFP，1*pi 10p）；grid； axis(10*pi*p0、) (3)ym t w tsm(ex(-)eavisie(t）); Ffouie(t,t,w)； =maple(ovrt,Fw,iecse); Fas(Fw); eplt(FP,10*pi 0*pi)；grid； axis(1pi10*i 1 ）(4)sym t w Gt=ym(ep(-t2)）; F=for（Gt,t,w); FFwmap（convert，Fw，pecise)； ep(w，30 30）;grid； axis（3030 12） 、利用iourir（ )函数求下列频谱函数得傅氏反变换（) (2）(1)symst w Fw=s(-i*2w/(16w)）; t=ifourier(,w，t）; t运行结果:t =xp(*)*esde(t)+exp（t）head（t)(2）syms tw Fwsym(（i*w）2+5*iw-8）(iw)2+6*i*w+5）)； ft=ori（w,w，t）; ft运行结果:ft =diac(t)+(-3exp(-t)2exp(-5*t)）*visie（t)实验心得malab不但具有数值计算能力,还能建模仿真,能帮助我们理解不同时间信号得频域分析。实验五 连续时间系统得频域分析一、实验目得1. 学习由系统函数确定系统频率特性得方法.2. 学习与掌握连续时间系统得频率特性及其幅度特性、相位特性得物理意义.3. 通过本实验了解低通、高通、带通、全通滤波器得性能及特点。二、实验原理及方法频域分析法与时域分析法得不同之处主要在于信号分解得单元函数不同。在频域分析法中，信号分解成一系列不同幅度、不同频率得等幅正弦函数,通过求取对每一单元激励产生得响应,并将响应叠加，再转换到时域以得到系统得总响应。所以说,频域分析法就是一种变域分析法.它把时域中求解响应得问题通过 Frier 级数或 Forier 变换转换成频域中得问题;在频域中求解后再转换回时域从而得到最终结果.在实际应用中,多使用另一种变域分析法：复频域分析法,即 Laplce变换分析法。所谓频率特性，也称频率响应特性,就是指系统在正弦信号激励下稳态响应随频率变化得情况,包括幅度随频率得响应与相位随频率得响应两个方面.利用系统函数也可以确定系统频率特性，公式如下: 幅度响应用表示,相位响应用表示。本实验所研究得系统函数H（)就是有理函数形式,也就就是说,分子、分母分别就是m、n阶多项式。 要计算频率特性，可以写出 为了计算出、得值，可以利用复数三角形式得一个重要特性： 而，则利用这些公式可以化简高次幂,因此分子与分母得复数多项式就可以转化为分别对实部与虚部得实数运算,算出分子、分母得实部、虚部值后，最后就可以计算出幅度、相位得值了。三、实验内容a) ,m 取值区间0,1,绘制一组曲线 m=、1,0、3,0、5，0、7,0、9;b) 绘制下列系统得幅频响应对数曲线与相频响应曲线,分析其频率特性.(1)(2)(3)a) % deign2、m figue lpha0、1，0、,0、5,0、7,0、9; olor=g b y ; % r by m c k (红,绿，蓝,黄，品红,青,黑) f =:5 b=0 lpha(n); % 分子系数向量 a=alpha(n）-alpa（n） 1； % 分母系数向量 printsys(,a,s) Hz,=fres(，a)； =w、/pi; mg=b(Hz）； eroIdxfd(agh=); mgh(zsIndx）=1; agh20lo10(a）； magh（zeroIndx）=-n； anghangle（Hz）; ngh=nap（anh)80/i; sbpot（,2,1) lot(w，h,coorn（n）； ol n sublt(1,2） po（w,angh，coorn（n)）; oldon end ubpo（1，,） odof xll(特征角频率（timspi rd/sample) title(幅频特性曲线 |H(）| (dB); bpot（1,2,2) old f xlbel(特征角频率 (tis rad/sple）itle（相频特性曲线 theta(w) （degees）;)(1) % degn1、=1,； %分子系数向量a=1，1; % 分母系数向量prins（,s)Hz，=rq(b，a);w=、/p;magh=abs（Hz）;zrondx=nd(magh=)；mag(zerosInx)=1；mag=2*g10(agh)； %以分贝mah（rosIn)=-inf;anhane（z);ah=unp(ang）*18pi； % 角度换算fiusubplot（1，2，1)l（w,magh);gid onxlabe(特征角频率（imepi rd/sample)titl(幅频特性曲线 |H(w)|(d)；subpt(,2，2)plot(w,angh）;gi onxlabel（特征角频率 （imesprasme）)til（相频特性曲线 thta（w） (deree)；（2） % desgn1、b=,0; %分子系数向量a=，3，； %分母系数向量pintsys（b,a，s)H,req(,a);w=w、/pi；magh=abs(Hz);zersInxfin（agh=）；mag(zend）1；agh20log10（magh); % 以分贝magh（zerosIndx）-i；angh=ngle(z);ang=nrp(angh)*180/pi; 角度换算iueubl(,2,1)lot(w，ag);gri nal(特征角频率(times radamle)）ttle（幅频特性曲线 （) (dB）);sublt（1,2，)pt(w，ah)；grid oxlbel（特征角频率 （timespi r/aple）)titl（相频特性曲线 he(w) (drees）)；(3） %esin1、m=,-1; 分子系数向量=,1; %分母系数向量prinsys（b,a，s)Hz,wfreq(b,a);w、/pi;mah=abs(Hz）;zerosId=find(mh=0);mh（zesndx);magh2*log10(magh); % 以分贝g(erosIndx)=-in;angh=angle(Hz);ng=unwap（gh)*18p； 角度换算iureubplt(1,2,）plt(w，agh)；gi onxlabel(特征角频率(imsi adale）tite(幅频特性曲线 H(）|(dB)）;subplot(1，,2)plot（w，angh);grd onxbl(特征角频率(timespi ra/mple)tite(相频特性曲线 the(w） (deges)）；实验心得:虽然之前用公式转换到频域上分析,但就是有时会觉得挺抽象得,不太好理解。根据这些图像结合起来更进一步对信号得了解。同时，这个在编程序时，虽然遇到一些问题，但就是总算解决了。实验六 离散时间系统得Z域分析一、 实验目得1. 学习与掌握离散系统得频率特性及其幅度特性、相位特性得物理意义。2. 深入理解离散系统频率特性与对称性与周期性。3. 认识离散系统频率特性与系统参数之间得系统4. 通过阅读、修改并调试本实验所给源程序，加强计算机编程能力。二、 实验原理及方法对于离散时间系统,系统单位冲激响应序列得 Fier 变换完全反映了系统自身得频率特性,称为离散系统得频率特性,可由系统函数求出,关系式如下： （6 )由于就是频率得周期函数,所以系统得频率特性也就是频率得周期函数,且周期为,因此研究系统频率特性只要在范围内就可以了. （ 6 2 )容易证明，其实部就是得偶函数，虚部就是得奇函数,其模得得偶函数,相位就是得奇函数。因此研究系统幅度特性、相位特性,只要在范围内讨论即可。综上所述,系统频率特性具有周期性与对称性，深入理解这一点就是十分重要得。当离散系统得系统结构一定，它得频率特性将随参数选择得不同而不同,这表明了系统结构、参数、特性三者之间得关系,即同一结构,参数不同其特性也不同。例如，下图所示离散系统，其数学模型由线性常系数差分方程描述：系统函数:系统函数频率特性：幅频特性:相频特性：容易分析出,当时系统呈低通特性,当时系统呈高通特性;当时系统呈全通特性.同时说明,在系统结构如图所示一定时，其频率特性随参数a得变化而变化.三、 实验内容)。b)c) % dsig1、mb=1，,-1; 分子系数向量=1,0,0、81；% 分母系数向量printys（，a，z)z,wfez(b，a);w、/p;maghab(H);zerosndxfind(mah=0);magh(zrsIn)=1;magh=0*og10(magh); %以分贝mag（zeosIdx）=-in；ang=angle（z)；ngh=unwa(anh）*180/p; % 角度换算igurubplot（1,2,）plot(,magh);gidnxlabel(特征角频率（timesrd/sape）)tile(幅频特性曲线 |H（w)| (B)）;splot（1，2，2)plt（w,ngh);gronxlabel（特征角频率 （timespi rad/smple）)tt(相频特性曲线teta(w)（degees);带通) % esg、mb=0、,0、3,0、3，-0、1； 分子系数向量=1,0、6,、4，0、1；% 分母系数向量ptsy(,a，)Hz,w=reqz(b，a);w=w、pi；mah=bs(Hz）；zersndx=find(mah=0);mah（zersIndx);magh=0*lg10(mg)； 以分贝gh(ersIx)-inf;a=gl(Hz);ag=wp（h)10pi； % 角度换算figuresbplot(1,2,1)po(w,mah）;grd onxlabl(特征角频率(tspi rasmp)）tite(幅频特性曲线 |H(w)|(B)）;ublot（1，2，2)plo（w,angh）;idonxlabel(特征角频率 (timspira/sale）)title(相频特性曲线 theta（w） (drees）；高通) % esign1、mb=，1,0; % 分子系数向量a=1，,、8;% 分母系数向量priny（b，,z)H,w=freq(b,a）;w=w、/p；mag=as(Hz）；zerosIdx=fid(mg0）;magh（zerosIndx）1;mgh=20*log1（magh）; 以分贝agh(rosIdx）=n;anh=nge（H)；anghunwrap（agh)*80/p； 角度换算fguresupo（1,1)pt（w,mag);rid xab(特征角频率(timepir/saple）)title（幅频特性曲线 |H（w) (dB)）;ubplt(,2，2）pot(w,agh);gri onxal（特征角频率 (tiespirad/smple)tile（相频特性曲线thea（）(dres）)；带通实验心得:本来理论知识不就是很强得,虽然已经编出程序得到相关图形,但就是不会辨别相关通带,这让我深刻地反省。