samp2_16.m

来自「数字信号处理的Matlab实现的随书光盘」· M 代码 · 共 16 行

%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')
%绘出积分后信号图