%Samp2_16
dt=0.01;            %信號的采樣間隔
t=0:dt:4*pi;        %時間序列
y1=sin(t);          %原始正弦函數值序列
subplot(3,1,1),plot(t,y1),title('原始信號')
%繪出原始正弦信號
y2=diff(y1)/dt;   %對信號進行微分
subplot(3,1,2),plot(t,[0,y2]) ,title('微分后信號')
%繪出微分后的信號，因為微分后信號比原信號少一個元素，用零補。
for i=1:length(y2)
y3(i)=sum(y2(1:i))*dt;
end
%將原始信號微分后的y2進行積分
subplot(3,1,3),plot(t,[0 y3]) ,title('積分后信號'),xlabel('時間/s')
%繪出積分后信號圖