c=square(pi*t1,50);
t2=0:0.01:1-0.01; y1=t2; y2=1-t2+0.01; y=[y1,y2]; subplot(221); stem(t,a); xlabel('t'); ylabel('x(t)/h(t)'); title('输入方波'); subplot(222); stem(n,b); xlabel('t'); ylabel('y(t)');
title('卷积结果（实际）'); subplot(223); plot(t1,c,'r'); xlabel('t'); ylabel('x(t)/h(t)'); title('输入方波'); subplot(224);
plot(t1,y,'r'); axis([0,4,0,1]); xlabel('t'); ylabel('y(t)');
title('卷积结果（理论）');
验证e指数函数卷积的程序如下： close all; clc; clear; t=0:0.1:10; a=exp(-2*t).*(t>=0); b=[ones(1,100)]; c=conv(a,b)*0.1;
m=(length(c)-1)*0.1; n=0:0.1:m; t1=0:0.1:10;
y1=0.5*(1-exp(-2*t1)); t2=10:0.1:20;
y2=-0.5*(exp(-20)-exp(-2*(t2-10))); subplot(211); stem(t,a); xlabel('t'); ylabel('h(t)'); title('系统'); axis([0,10,0,1]); subplot(212); stem(n,c); hold on;
plot(t1,y1,'r',t2,y2,'r'); xlabel('t'); ylabel('y(t)'); title('卷积结果'); axis([0,20,0,1]);
实验二 用FFT作谱分析
下图给出了主程序框图，供参考。本实验提供FFT子程序和通用绘图子程序。
调用信号产生子程序产生实验信号 读入长度N 开始
