2011信与系统实验指导

1、信号与系统实验指导实验一熟悉Matlab软件环境一、实验目的：熟悉Matlab函数创建、基本计算及绘图。二、试验内容及要求：编制Matlab程序，实现以下内容1 熟悉变量、矩阵及创建矩阵、向量等基本概念；创建基本矩阵函数：n=3、m=2ones(n,m)，zeros(n,m)，rand(n,m)，randn(n,m)，eye(n)2 创建函数：求10个数值的平均值、总和。3 绘制图形：y=sint、y=exp(-0.1*t)*sin(t+1)三、实验步骤：打开电脑，双击桌面Matlab程序图标，打开文件菜单，点击新建M文件，在M文件编辑器中编写程序，编写完毕后按F5键运行调试。编写函数时注意以

2、下几点：.扩展名为.m;.函数必须以关键字“functio开头；.函数文件的第一行为函数说明语句：function返回参数1，返回参数2，=函数名(传入参数1，).函数文件保存的文件名应与用户定义的函数名一致，并且函数名不能与Matlab自带的函数名一致。四、试验报告：试验目的、内容、步骤、程序、图形(可打印)等。附录%试验一,内容1A=ones(3,2)B=zeros(3,2)C=rand(3,2)D=randn(3,2)E=eye(3)%试验一,内容2x=input('x=')y,z=qiuhe(x);%函数见qiuhe.mfprintf('N个数总和为')

3、,yfprintf('N个数平均为'),z%试验一,内容3t=0:0.1:30;y1=sin(t);y2=exp(-0.1.*t).*sin(t+1);subplot(2,1,1),plot(t,y1),grid;title('正弦信号');subplot(2,1,2),plot(t,y2),grid;title('衰减正弦信号');functiony,z=qiuhe(x)m,n=size(x);y=sum(x);z=sum(x)/n;functionr=polyadd(p,q)lp=length(p);lq=length(q);k=lp-lq;

4、ifk>=0r=p+zeros(1,k),q;elser=zeros(1,-k),p+q;end16试验二连续系统的时域分析一、试验目的：熟悉掌握利用Matlab产生基本连续时间信号及时域分析方法。二、试验内容及要求：编制Matlab程序，实现以下内容1利用信号处理工具箱提供函数生成单位冲激函数、单位阶跃函数、三角波、方波。2已知描述系统的微分方程和初始状态如下，求其零输入响应。(1) y''(t)+5y'(t)+6y(t)=f(t);y(O_)=l,y'(O_)=-l;(2) y'、'(t)+4y1'(t)+5y'(t)+

5、2y(t)=f(t),y(0_)=0,y'(0_)=-1,y、(0)=-1；d2ydy3某LTI系统的微分方程为+2+8y二x(t)，求冲激响应；dt2dt若输入为x(t)=cos(0.1t)，求其零状态响应。三、实验步骤：打开电脑，双击桌面Matlab程序图标，打开文件菜单，点击新建M文件，在M文件编辑器中编写程序，编写完毕后按F5键运行调试。四、实验原理：对任何系统依据产生系统响应的原因可以将完全响应分解为零输入响应和零状态响应两个分量，系统完全响应=零输入响应+零状态响应，即r(t)=r.(t)+r(t)，其中r.(t)为零输入响应；r(t)为零状态响应。zizszizsdndt

6、ndn1r(t)+Cr(t)+zi1dtn1zi零输入响应：没有外加激励信号的作用，只由起始状态(t=0)产生的响应，满足：+C£r(t)+Cr(t)二0n1dtzinzi及起始状态r(k)(0)(k=0,l,2,.,n-l)的解；它是齐次解中的一部分：rzi(t)=乙Azikeakt,由于没有外界激励，系统状态不会发生变化，即:r(k)(O_)=r(k)(O+)k=1零状态响应：不考虑起始时刻(t=0)系统储能的作用，仅由系统的外加激励信号产生响应，满足：dndn-1dCr(t)+Cr(t)+Cr(t)+Cr(t)0dtnzs1dtn-1zsn-1dtzsnzsdmdm-1d=Ee

7、(t)+Ee(t)+Ee(t)+Ee(t)0dtm1dtm-1m-1dtm及起始状态r(k)(0)=0(k=0,l,.,n-1)的解；它是由特解和一部分齐次解构成：r(t)=yAeakt+B(t)zszskk=1零状态响应也可由卷积的方法得出：r(t)=e(t)*h(t)zs五、试验报告：试验报告应包含以下内容试验目的、内容、步骤、原理、程序、图形(可打印)等。附录%试验二,内容1%产生单位冲激信号clearall;t0=0;t1=1;dt=0.05;t=0:0.05:5;st=length(t);n1=floor(t1-t0)/dt);x1=zeros(1,st);x1(n1)=1/dt;s

8、ubplot(4,1,1),stairs(t-0.95,x1),grid;axis(-1,4,-1,22)%产生单位阶跃信号x2=zeros(1,n1-1),ones(1,st-n1+1);subplot(4,1,2),stairs(t-0.95,x2),grid;axis(-1,4,-0.1,1.1)%产生周期的三角波clearall;t=-4*pi:0.01:4*pi;a=sawtooth(t,0.5);subplot(4,1,3),plot(t,a),axis(-4*pi,4*pi,-1.1,1.1),grid;%产生周期的方波b=square(t,50);subplot(4,1,4),

9、plot(t,b),axis(-4*pi,4*pi,-1.1,1.1),grid;%实验二，内容2a=input('输入分母系数向量a=al,a2,.=');n=length(a)-1;y0=input('输入初始条件向量yO=yO,dyO,.=');p=roots(a);v=rot90(vander(p);c=v/y0;dt=0.2;tf=8;t=0:dt:tf;y=zeros(l,length(t);fork=l:ny=y+c(k)*exp(p(k)*t);end;plot(t,y),grid;%实验二，内容3a=input('输入多项式分母系数向量

10、a=a1,a2,.=');b=input(输入多项式分子系数向量b=b1,b2,.=');u=input('输入信号u(t)=');dt=0.2;tf=8;t=0:dt:tf;r,p,k=residue(b,a);h=zeros(1,length(t);fori=1:length(a)-1;h=h+r(i)*exp(p(i)*t);endsubplot(2,1,1),plot(t,h),grid;title('冲激响应')；y=conv(u,h);subplot(2,1,2),plot(t,y(1:length(t),grid;title(

11、9;零状态响应');试验三连续系统的频域分析一、试验目的：分析傅里叶级数以及利用傅氏变换研究信号频谱。二、试验内容及要求：编制Matlab程序，实现以下内容1傅氏级数:将周期为T=2兀,幅度为1的对称方波分解为多次正弦波之和，用Matlab演示谐波合成情况。(取到7次谐波)2利用Matlab分析矩形脉冲信号的频谱；求f(t)=u(t+0.5)-u(t-0.5)的傅立叶变换三、实验步骤:打开电脑，双击桌面Matlab程序图标，打开文件菜单，点击新建M文件，在M文件编辑器中编写程序，编写完毕后按F5键运行调试。四、实验原理:傅式级数所选项数越多，相加后波形越逼进原信号，二者方均误差越小；当

12、信号是脉冲信号时，其高频分量主要影响脉冲的跳变沿，低频分量主要影响脉冲的顶部，即信号波形变化越剧烈，包含的高频成分越丰富；变化越缓慢，包含的低频成分越丰富；当信号中任一频谱分量的幅度或相位发生相对变化时，输出波形一般要发生失真。傅立叶变换对:F(®)=f(t)eTotdt傅立叶正变换;gf(t)二1貯fGoL®td®2兀g傅立叶逆变换，简写fG)oFeo)傅立叶变换的物理意义在于信号可分解为无穷多个幅度为无穷小的连续指数信号d“I之和,整个频域-8T8由于幅度无穷小，所以频谱不能再用幅度表示，而改用密度函数表示。五、试验报告要求：试验报告应包含以下内容试验目的、内

13、容、步骤、原理、程序、图形(可打印)等附录%实验三，内容1functionh=my301(a);t=0:0.01:2*pi;nn=a;y1=0;fori=1:nn;xb=i*2-1;subplot(2,1,1),y=sin(xb*t)/xb;plot(t,y),grid,holdon;title('对称方波傅式级数');subplot(2,1,2)y1=y1+y;endplot(t,y1),gridh=0;%实验三，内容2clear,tf=1;N=256;dt=1/N;t=-N+1:N*dt;f=zeros(1,N/2),ones(1,N),zeros(1,N/2);wf=50

14、;nf=64;w1=linspace(0,wf,nf);dw=wf/(nf-1);f1=f*exp(-j*t'*w1)*dt;w=-fliplr(w1),w1(2:nf);f2=-fliplr(f1),f1(2:nf);subplot(l,2,l),plot(t,f,'linewidth',1.5),axis(-l,l,-0.1,l.l),grid;title('时域方波');subplot(l,2,2),plot(w,abs(f2),'linewidth',1.5),grid;title('方波频谱')；试验四连续系统的

15、复频域分析一、试验目的1观察系统的零极点分布对系统稳定性的影响以及系统的串并联。2练习建立Simulink模型。二、试验内容及要求：1编制Matlab程序，实现以下内容：系统的串、并联及反馈：两个单入单出系统W=10(0鬆+»、w=一-一,A5s2+2s+1Bs(s+1)绘制出构成串联、并联及反馈系统后的零极点分布图，并判断稳定性。2试用Simulink建立上题两系统串联、并联的模型；设输入信号为阶跃信号，用示波器观察输出。三、实验步骤：打开电脑，双击桌面Matlab程序图标，打开文件菜单，点击新建M文件，在M文件编辑器中编写程序，编写完毕后按F5键运行调试。点击Simulink图标

16、，运行仿真工具箱，在对话框中拖拽各模块建立系统模型，完成之后按F5键运行调试四、实验原理：根据系统函数极点的分布，可判别系统的稳定性(零点分布只影响原函数幅度和相位)。系统极点为一阶极点时，极点位于s平面左侧，h(t)收敛，系统稳定；极点位于s平面右侧，h(t)发散，系统不稳定；极点位于s平面虚轴上，h(t)等幅振荡，系统处于临界稳定状态。系统极点为二阶极点时，极点位于s平面左侧，h(t)收敛，系统稳定；极点位于s平面右侧或虚轴上，h(t)发散，系统不稳定。LTIS互联的系统函数：LTI系统的并联：H(s)二H(s)+h2(s)LTI系统的级联：H(s)二H(s)-H2(s)H(s)H(s)L

17、TI系统的反馈连接：(s)+h(s)H(s)12建立Simulink仿真框图参见下图例双击模块可更改参数)五、试验报告要求：试验报告应包含以下内容试验目的、内容、步骤、原理、程序、图形(可打印)等附录%实验四，内容1a%计算传递函数fA=5,10;gA=5,2,1;fB=4;gB=1,1,0;%串联fh1=conv(fA,fB);gh1=conv(gA,gB);disp('串联后传递函数')printsys(fh1,gh1,'s')%并联fh2=polyadd(conv(fA,fB),conv(gA,gB);gh2=conv(gA,gB);disp('并

18、联后传递函数')printsys(fh2,gh2,'s')%反馈后传递函数fh3=conv(fA,gB);gh3=polyadd(conv(fA,fB),conv(gA,gB);disp('反馈后传递函数')printsys(fh3,gh3,'s')%实验四，内容1b%画零极点图num=input('输入分子系数sum=al,a2,.=');den=input('输入分母系数den=b1,b2,.=');z,p=tf2zp(num,den);zplane(z,p),grid试验五信号的调制与解调一、试验目的

19、：1、了解用Matlab实现信号调制与解调的方法。2、利用仿真工具箱实现调制与解调系统的框图。二、试验内容：1编制Matlab程序，实现信号的调制：将一低频信号g（t）=cos5t进行调制，载波频率叫=100,并绘制频谱图。g（t）f（t）=g（t）cos®0t相乘COS®0t2利用Simulink建立内容1中所示框图的调制与解调模型（同步解调），并用示波器观察各点的时域波形，绘制频谱图。三、实验步骤：打开电脑，双击桌面Matlab程序图标，打开文件菜单，点击新建M文件，在M文件编辑器中编写程序，编写完毕后按F5键运行调试。点击Simulink图标，运行仿真工具箱，在对话框

20、中拖拽各模块建立系统模型，完成之后按F5键运行调试四、实验原理：在通信系统中，信号从发射端传输到接收端，为实现信号的传输，往往要进行调制和解调。调制是将信号的频谱搬移到任何所需的较高频段上的过程。调制有很多种比如调幅分为振幅调制（AM）、抑制载波振幅调制（AM-SC）、单边带调制（SSB）、残留边带调制（VSB）等，按控制载波的其他参数也有调频（FM）与调相（PM）等。1、调制：我们这里介绍抑制载波振幅调制(AM-SC),将调制信号g(t)与载波信号cos叫t(叫)相乘，如图所示，得到已调信号f(t)=g(t)cos叫t。g(t)根据傅氏变换性质知时域相乘，频域f(t)=g(t)cos

21、4;0t卷积可求出已调信号f(t)的频谱COS®0tF(o)=1lG(o，这200样就将调制信号g(t)的频谱G(®)搬移到了比较高的一个频段上，搬移过程见下图。解调是将已调信号恢复成原始信号的过程，与调制相反，如下图所示(同步解调)f(t)=g(t)cosw0t这种解调方式要求接收段的本地载波与发送端的载波同频同相，所以成为同步解调，已调信号经乘法器输出g0(t)g(t)=g(t)cos2(tot)=1g(t)l1+cos(2®tU=.!g(t)+g(t)cos(2®t)0020220对应频谱为G0（）二2G（）+4G（2°0）+4F（+2&

22、#176;0）,g0（t）通过截止频率®m<®c<2®0-®m的低通滤波器就可以恢复出原调制信号g（t）,完成解调，频谱图如下：“G(w)1/20滤波器1/4m2®0卞二:增益为2J1plGn(TO)0TOto1/4严2toto0m五、试验报告要求：试验报告应包含以下内容试验目的、内容、步骤、原理、程序、图形（可打印）等附录%实验五，内容1%信号调制cleart=linspace(0,2*pi,1001);u=0.5,0.5;u1=cos(5*t);u2=cos(100*t);u31=u(1)*cos(95*t);u32=u(2)*

23、cos(105*t);u3=u31+u32;subplot(3,1,1),plot(t,u1),grid;title('调制');subplot(3,1,2),plot(t,u2),gridsubplot(3,1,3),plot(t,u3),grid试验六调幅信号通过带通滤波器一、实验目的：1掌握调幅信号通过带通滤波器的原理，以及其分析过程;2练习建立Simulink模型。二、试验内容：1编制Matlab程序，实现以下内容调幅信号通过带通滤波器：已知带通滤波器系统函数为：H(s)=虬(S+1)2+1002激励电压U(t)=(l+cost)coslOOt,求(1) 带通滤波器的频率响应；(2) 输出的稳态响应。2试用Simulink建立上题的模型；(观察滤波器前后的波形)讨论结果。三、实验步骤：打开电脑，双击桌面Matlab程序图标，打开文件菜单，点击新建M文件，在M文件编辑器中编写程序，编写完毕后按F5键运行调试。点击Simulink图标，运行仿真工具箱，在对话框中拖拽各模块建立系统模型，完成