samp2_12.m

数字信号处理的Matlab实现的随书光盘

%Samp2_12
t= 0:0.001:1;                % 将1s的时间段以1kHz的采样频率采样
fo=2;f1=10;                  % 频率渐增函数chirp从2Hz增加到10Hz
x=chirp(t,fo,1,f1);          %在时间序列t上产生频率渐增信号--鸟鸣信号
m=1:length(x);               %为原chirp信号的序号序列
subplot(2,1,1),plot((m-1)*0.001,x)  %绘出原始随时间变化chirp信号，m为序号，因此必须减1与采样间隔相乘才得到时间序列。
xlim([0 1.5]),title('原始信号')
[y,n]=sigshift(x,m,500);    %将序号序列为m，值序列为x的信号时移得到序号序列为n，值序列为y的新序列
subplot(2,1,2),plot((n-1)*0.001,y)             %绘制新序列信号
xlim([0 1.5]),title('时移后信号'),xlabel('时间/s')