exm040835_3.m

matlab mfiles实例分析

%“零阶”法计算卷积
t2=3;t4=11;T=0.01;
tu=0:T:t2;N2=length(tu);
th=0:T:t4;N4=length(th);
u=exp(-tu);h=th.*exp(-th/2);
tx=0:T:(t2+t4);Nx=length(tx);
yt3=T*conv(u,h);
%把计算结果与理论值比较
t=tx;yyt1=eval(vectorize(char(yt1)));%例4.8.3.5-2第（1）步的计算结果
[dy,kd]=max(abs(yyt1(1:N2)-yt3(1:N2)));
dy13(1)=dy/abs(yyt1(kd));
[dy,kd]=max(abs(yyt1(N2+1:N4)-yt3(N2+1:N4)));
dy13(2)=dy/abs(yyt1(N2+kd));
[dy,kd]=max(abs(yyt1(N4+1:Nx)-yt3(N4+1:Nx)));
dy13(3)=dy/abs(yyt1(N4+kd));   

disp('与理论结果的相对误差')
disp([blanks(4),'[0,3]段   [3,11]段  [11,14]段']),disp(dy13)
plot(t,yyt1,':b',tx,yt3,'r')
set(gca,'Xtick',[0,3,11,14]),grid