q6014.m

来自「电路与系统专业的matlab程序。有很多经典的实例。」· M 代码 · 共 16 行

% 本例演示卷积的时域积分法：
% 已知系统冲激响应h(t)=(exp(-t/T))/T ，求输入u(t)=exp(-t)下的输出响应。
syms T t tao;
ut=exp(-t);							%定义系统输入
ht=exp(-t/T)/T;						%定义系统冲激响应
uh_tao=subs(ut,t,tao)*subs(ht,t,t-tao);	%运用变量替换指令形成被积函数
yt=int(uh_tao,tao,0,t);					%实施卷积
yt=simple(yt)                           %化简
yt1=subs(yt,T,'2')
ezplot(yt1,[-2,2])
% 本例演示通过变换和反变换求取卷积。系统冲激响应、输入同上例，求输出。
syms s;
yt=ilaplace(laplace(ut,t,s)*laplace(ht,t,s),s,t);
yt2=simple(yt)