资源描述:
《用matlab进行信号与系统的时频域分析》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、用Matlab进行信号与系统的时频域分析专业学生姓名班级学号指导教师指导单位曰期一、实验目的1、了解matlab中相关函数的调用,以及其参数设定方式。学会调试代码运行过程中出现的错误。2、熟悉掌握时域、频域的相关知识,加深对书本知识的理解。二、实验任务1、完成实验内容的三个小题,分析解决调试代码过程中出现的问题。2、对第三个小题做拓展,完成思考题并做小结。3、认真完成本次实验小结,总结概括实验心得体会。三、主要仪器设备硬件:微型计算机软件:Matlab四、实验内容(一)连续时间信号的卷积上机题1.已知两个信号七0)=冰一1)-卽-2)和x2(r)=£(z)
2、-冲-1),试分别画出A⑴,x2⑴和卷积yO)=七⑺*x2⑺的波形。图形:代码:T=0.01;tl=l;t2=2;t3=0;t4=l;t3=0;t4=l;t=0:T:(t2+t4);t=0:T:(t2-K4);xl=ones(size(t)).*((t>t1)-(t>t2):xl=ones(size(t)).*((t>tl)-(t>t2)):x2=ones(size(t)).*((t>t3)-(t>t4));y=conv(xl,x2)*T;subplot(3,I,1),plot(1,xl);subplot(3,1,1),plot(t,xl);ylabelC
3、xl(t)');subplot(3,1,2),plot(t,x2);ylabelCx2(t)*);subplot(3,1,3),plot(t,y(l:(t2+t4)/T+l));ylabelCy(t)=xl*x2*);xlabeK*^t/s*);上机题2.已知两个信号=⑺和/2⑴=冲),试用数值计算法求卷积,并分别画出和卷积y(r)=x(0*始)的波形。»t=sym(’V’positive’):???Errorusing==>sym.sym>char2syjnNotavalidsymbolicexpression.Errorin==>sy».symat92
4、S=char2sym(x);»t=symCtrpositive');»xt=exp(-t/2);>>xh垂tao=subs(xt,t,tao)*subs(ht,t,t-tao???Undefinedfunctionorvariable’tao’»symstao;»t=symC'positive');»xh.tao=subs(xt,t,tao)*subs(ht,t,t-tao???Undefinedfunctionorvariable’ht’.»sym5tao;»txs/mC*positive');»xt=exp(-t):»ht=t*exp(-t/2);»x
5、h-tao=subs(xt,t,tao)*subs(ht,t,t-tao»yt=int(xh-tao,tao,0,t)4*exp(-t)+2*t*exp(-1/2*t)-4*exp(-l/2*t)图形:ixure1lilt£du父IntertZ^liIiiU的Da■o代码:»t2=3;t4=ll;»T=0.01;»t=O:T:(t2-K4);»x=exp(-t0.*((t>0)-(t>t2));???x^exp(-t0.*((t>0)-(t>t2));IError:Incompleteormisformedexpressionorstatement.»x=
6、exp(-t).*((t>0)-(t>t2));»h=t.*exp(-t/2)•*((t>0)-(t>t4)〉:»y=conv(x,h)*T;»yt=4*exp(-l)+2*t.exp(-l/2*t)-4*exp(-l/2*t);???Attempttoreferencefieldofnon-structurearray.»yt=4*exp(-l)+2*t.*exp(-l/2*t)-4*exp(-l/2*t);»subplot(3,1,1〉,plot(t,x);»ylabelCx(t)*);»sub(3,1,2),plot(t,h);???Undefine
7、dcommand/function’sub’.»subplot(3,1,2),plot(t,h):»ylabel(’h(t〉’〉:»subplot(3,1,3),plot(t,y(l:(t2+t4)/T+l〉,yt,’--r*);»legend(bynumericaV,•Theoretical'):»ylabelCy=x*h*);»xlabel(>t/s.);(二)信号的频域分析上机题3.求周期矩形脉冲信号的频谱图,己知A=l,r=0.15,7=0.55图形:代码:»a=l;tao=0.1;t=0.5;»nO=t/tao;»n=0:2*n0;»fn_p=a
8、*tao/1*(sin(n*pi*tao/t+eps*(n==0)