信号与系统的MATLAB仿真

1 信号与系统的MATLAB 仿真 一、信号生成与运算的实现 1.1 实现)3(sin)()(ttttStfa )(sin)sin()sin(sin)()(tctttttttStfa m 11.m t=-3*pi:0.01*pi:3*pi; % 定义时间范围向量 t f=sinc(t/pi); % 计算 Sa(t)函数 plot(t,f); % 绘制 Sa(t)的波形 运行结果： 1.2 实现)10()sin()(sin)(ttttctf m 12.m t=-10:0.01:10; % 定义时间范围向量 t f=sinc(t); % 计算 sinc(t)函数 plot(t,f); % 绘制 sinc(t)的波形 运行结果： 1.3 信号相加：tttf20cos18cos)( m 13.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)的波形 运行结果： 2 1.4 信号的调制：tttf50cos)4sin22()( m 14.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 信号相乘：)20cos()(sin)(ttctf m 15.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)'); % 加注波形标题 运行结果： 3 二、系统时域的仿真分析 2.1 实现卷积)(*)(thtf，其中：)2()()()],1()([2)(ttthtttf 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) su bplot(3,1,1),stairs(nf,f); % 绘制 f(t)的波形 title('f(t)');ax is([0 3 0 2.1]); su bplot(3,1,2),stairs(nh,h); % 绘制 h(t)的波形 title('h(t)');ax is([0 3 0 1.1]); su bplot(3,1,3),plot(k,y ); % 绘制 y (t)=f(t)*h(t)的波形 title('y (t)=f(t)*h(t)');ax is([0 3 0 2.1]); 子程序 sconv.m % 此函数用于计算连续信号的卷积 y (t)=f(t)*h(t) fu nction [y ,k]=sconv (f,h,nf,nh,p) % y :卷积积分 y (t)对应的非零样值向量 % k:y (t)对应的时间向量 % f:f(t)对应的...