信号与系统的matlab仿真

《信号与系统的matlab仿真》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、信号与系统的MATLAB仿真信号与系统的MATLAB仿真一、信号生成与运算的实现1.1实现EMBEDEquation.3EMBEDEquation.3m11.mt=-3*pi:0.01*pi:3*pi;%定义时间范围向量tf=sinc(t/pi);%计算Sa(t)函数plot(t,f);%绘制Sa(t)的波形运行结果：1.2实现EMBEDEquation.3m12.mt=-10:0.01:10;%定义时间范围向量tf=sinc(t);%计算sinc(t)函数plot(t,f);%绘制sinc(t)的

2、波形运行结果：1.3信号相加：EMBEDEquation.3m13.msymst;%定义符号变量tf=cos(18*pi*t)+cos(20*pi*t);%计算符号函数f(t)=cos(18*pi*t)+cos(20*pi*t)ezplot(f,[0pi]);%绘制f(t)的波形运行结果：1.4信号的调制：EMBEDEquation.3m14.msymst;%定义符号变量tf=(2+2*sin(4*pi*t))*cos(50*pi*t)%计算符号函数f(t)=(2+2*sin(4*pi*t))*co

3、s(50*pi*t)ezplot(f,[0pi]);%绘制f(t)的波形运行结果：1.5信号相乘：EMBEDEquation.3m15.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)');%加注波形标题运行结果：二、系统时域的仿真分析2.1实现卷积EMBEDEquation.3，其中：EMBEDEqua

4、tion.3m21.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),sta

5、irs(nh,h);%绘制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)对应的非零样值

6、向量%nf:f(t)对应的时间向量%h:h(t)对应的非零样值向量%nh:h(t)对应的时间向量%p:取样时间间隔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)非零样值的时间向量运行结果：2.2实现卷积EMBEDEquation.3

7、，其中：EMBEDEquation.3m22.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(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,

8、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]);运行结果：2.3设方程EMBEDEquation.3，试求零状态响应EMBEDEquation.3m23.m：yzs=dsolve('D2y+5*Dy+6*y=2*exp(-t)','y(0)=0,