HYPERLINK"http://www.wordwendang.com"http://www.wordwendang.com更多相关文档免费下载请登录:HYPERLINK"http://www.wordwendang.com"http://www.wordwendang.com-中文word文档库信号与系统的MATLAB仿真一、信号生成与运算的实现本文档由【HYPERLINK"http://www.wordwendang.com"中文word文档库】HYPERLINK"http://www.wordwendang.com"www.wordwendang.com提供，转载分发敬请保留本信息；HYPERLINK"http://www.wordwendang.com"中文word文档库免费提供海量范文、教育、学习、政策、报告和经济类word文档。<ahref='http://www.wordwendang.com/'>word文档</a>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)ezplot(f,[0pi]);%绘制f(t)的波形运行结果：1.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)');%加注波形标题运行结果：二、系统时域的仿真分析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);%绘制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:取样时间间隔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实现卷积，其中：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(n)的值[y,k]=sconv(f,h,