信号和系统的MATLAB仿真.doc

信号与系统的MATLAB仿真 m11.m t=-3*pi:0.01*pi:3*pi; % 定义时间范围向量t f=sinc(t/pi); % 计算Sa(t)函数 plot(t,f); % 绘制Sa(t)的波形 运行结果： 1.2 实现 m12.m t=-10:0.01:10; % 定义时间范围向量t f=sinc(t); % 计算sinc(t)函数 plot(t,f); % 绘制sinc(t)的波形 运行结果： 1.3 信号相加： m13.m syms t; % 定义符号变量t f=cos(18*pi*t)+cos(20*pi*t); % 计算符号函数f(t)=cos(18*pi*t)+cos(20*pi*t) ezplot(f,[0 pi]); % 绘制f(t)的波形 运行结果： 1.4 信号的调制： m14.m syms t; % 定义符号变量t f=(2+2*sin(4*pi*t))*cos(50*pi*t) % 计算符号函数f(t)=(2+2*sin(4*pi*t))*cos(50*pi*t) ezplot(f,[0 pi]); % 绘制f(t)的波形 运行结果： 1.5 信号相乘： m15.m t=-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.m p=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([0 3 0 2.1]); subplot(3,1,2),stairs(nh,h); % 绘制h(t)的波形 title(h(t));axis([0 3 0 1.1]); subplot(3,1,3),plot(k,y); % 绘制y(t)=f(t)*h(t)的波形 title(y(t)=f(t)*h(t));axis([0 3 0 2.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)