exm05835_2.m

精通matlab

function yt1=exm05835_2(flag)
%exm05835_2.m
%flag取 1 时，只进行符号计算
%flag取 2 时，只进行SIMULINK仿真
%flag取 3 时，把SIMULINK仿真结果与解析结果比较

switch flag
case 1			%（1）符号卷积法得解析结果
   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氏反变换得卷积的理论解
case 2		   %（2）SIMULINK卷积法
   exm05835_2_2
   yt1=[];
case 3			%（3）计算结果比较
   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]=sim('exm05835_2_2');
   yyt1=eval(vectorize(char(yt1)));	%理论解的数值序列
   [dy,kd]=max(abs(yyt1-yt2(:,2)));
   dy12=dy/abs(yyt1(kd))
end