exm040835_2.m

来自「MATLAB 6.5 从基础讲解其功能」· M 代码 · 共 11 行

syms tao;t=sym('t','positive');		%把t定义为“取正”符号变量
US1=laplace(exp(-t));					%u(t)的L氏变换
HS1=laplace(t*exp(-t/2))				%h(t)的L氏变换
yt1=simple(ilaplace(US1*HS1))			%L氏反变换得卷积的理论解  


t=yt2(:,1);				%exm040835_2_2.mdl示波器所保存的仿真采样时间序列
yyt1=eval(vectorize(char(yt1)));		%理论解的数值序列
[dy,kd]=max(abs(yyt1-yt2(:,2)));
dy12=dy/abs(yyt1(kd))