题主的问题出在f和Y的数据长度不一致。
可以通过修改t和N值，就可以正常出图。即下列两处
T=1e-2;
N=400; 增加
t=linspace(0,2,N/2);%t=0:0.001:2;
。。。
subplot(212);
%N=200;
修改后运行，可以得到如下结果
