例2-22 正弦信号的微分运算和积分运算

%2017/8/21日  星期一
%工程北629
%功能：数字信号  微分  积分
%思想：diff函数  cumsum函数或cumtrapz函数
%注意：这里直接使用t作为横坐标，没有使用n

clear all; clc

%输入信号
dt = 0.01;         %采样间隔 
t  = 0: dt: 4*pi;  %时间序列, 4*pi=12.5664

%信号处理
y1 = sin(t);        %采样间隔必须小于pi，这是采样定理要求的，可以从原始式子y = A*sin(2*pi*fc*t)得到。
y2 = diff(y1)/dt;   %对信号进行微分
y3 = cumsum(y2)*dt; %对信号进行积分

%输出信号
subplot(3,1,1),plot(t,y1),     title('原始信号');    grid on
subplot(3,1,2),plot(t,[0, y2]),title('微分后信号');  grid on  %微分后的信号比原信号少一个元素，用0补上
subplot(3,1,3),plot(t,[0, y3]),title('积分后信号');  grid on
xlabel('时间t/s')

%总结：将正弦信号微分得到余弦信号，将余弦信号微分得到正弦信号。
%      改变采样间隔，不会导致信号的周期的改变，仅仅会改变图形的采样的点数。图形大致的形状都没有非常大的改变