资源描述：
《信号与系统(第四版 燕庆明)(含习题解答)信号与系统的MATLAB仿真.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、信号与系统的MATLAB仿真一、信号生成与运算的实现1.1实现m11.mt=-3*pi:0.01*pi:3*pi;%定义时间范围向量tf=sinc(t/pi);%计算Sa(t)函数plot(t,f);%绘制Sa(t)的波形运行结果：1.2实现m12.mt=-10:0.01:10;%定义时间范围向量tf=sinc(t);%计算sinc(t)函数plot(t,f);%绘制sinc(t)的波形运行结果：1.3信号相加：m13.msymst;%定义符号变量tf=cos(18*pi*t)+cos(20*pi*t);%计算符号函数f(t)=cos(18*pi*t)+cos(20*pi*t)ezpl
2、ot(f,[0pi]);%绘制f(t)的波形运行结果：271.4信号的调制：m14.msymst;%定义符号变量tf=(2+2*sin(4*pi*t))*cos(50*pi*t)%计算符号函数f(t)=(2+2*sin(4*pi*t))*cos(50*pi*t)ezplot(f,[0pi]);%绘制f(t)的波形运行结果：1.5信号相乘：m15.mt=-5:0.01:5;%定义时间范围向量f=sinc(t).*cos(20*pi*t);%计算函数f(t)=sinc(t)*cos(20*pi*t)plot(t,f);%绘制f(t)的波形title('sinc(t)*cos(20*pi*t
3、)');%加注波形标题运行结果：27二、系统时域的仿真分析2.1实现卷积，其中：m21.mp=0.01;%取样时间间隔nf=0:p:1;%f(t)对应的时间向量f=2*((nf>=0)-(nf>=1));%序列f(n)的值nh=0:p:2;%h(t)对应的时间向量h=(nh>=0)-(nh>=2);%序列h(n)的值[y,k]=sconv(f,h,nf,nh,p);%计算y(t)=f(t)*h(t)subplot(3,1,1),stairs(nf,f);%绘制f(t)的波形title('f(t)');axis([0302.1]);subplot(3,1,2),stairs(nh,h);
4、%绘制h(t)的波形title('h(t)');axis([0301.1]);subplot(3,1,3),plot(k,y);%绘制y(t)=f(t)*h(t)的波形title('y(t)=f(t)*h(t)');axis([0302.1]);子程序sconv.m%此函数用于计算连续信号的卷积y(t)=f(t)*h(t)function[y,k]=sconv(f,h,nf,nh,p)%y:卷积积分y(t)对应的非零样值向量%k:y(t)对应的时间向量%f:f(t)对应的非零样值向量%nf:f(t)对应的时间向量%h:h(t)对应的非零样值向量%nh:h(t)对应的时间向量%p:取样时
5、间间隔y=conv(f,h);%计算序列f(n)与h(n)的卷积和y(n)y=y*p;%y(n)变成y(t)left=nf(1)+nh(1)%计算序列y(n)非零样值的起点位置right=length(nf)+length(nh)-2%计算序列y(n)非零样值的终点位置k=p*(left:right);%确定卷积和y(n)非零样值的时间向量运行结果：272.2实现卷积，其中：m22.mp=0.01;%取样时间间隔nf=0:p:2;%f(t)对应的时间向量f=2*((nf>=0)-(nf>=2));%序列f(n)的值nh=0:p:4;%h(t)对应的时间向量h=exp(-nh);%序列h
6、(n)的值[y,k]=sconv(f,h,nf,nh,p);%计算y(t)=f(t)*h(t)subplot(3,1,1),stairs(nf,f);%绘制f(t)的波形title('f(t)');axis([0602.1]);subplot(3,1,2),plot(nh,h);%绘制h(t)的波形title('h(t)');axis([0601.1]);subplot(3,1,3),plot(k,y);%绘制y(t)=f(t)*h(t)的波形title('y(t)=f(t)*h(t)');axis([0602.1]);运行结果：272.3设方程，试求零状态响应m23.m：yzs=ds
7、olve('D2y+5*Dy+6*y=2*exp(-t)','y(0)=0,Dy(0)=0')ezplot(yzs,[08]);运行结果：yzs=exp(-t)+exp(-3*t)-2*exp(-2*t)即：2.4已知二阶系统方程对下列情况分别求，并画出其波形。a.b.c.d.m24.m:R=input('电阻R=');%以交互方式输入电阻R的值L=input('电感L=');%以交互方式输入电阻L的值C=input('电容C=');%以交互方
