exa1_31.m

新版Matlab 7辅助信号处理技术与应用

% exa1-31_chirp.m for example 1-31
% chirp function

clear all;
fs=1000; %sample frequency
t=0:1/fs: 1;
f0=10; %10Hz at t=0
t1=1;
f1=100; %cross 100Hz at t=1sec
y=chirp(t,f0,t1,f1);%generate a linear chirp
y1=chirp(t,f0,t1,f1,'q');%generate quadratic chirp
subplot(211);plot(t,y)
title('线性调频信号')
subplot(212); plot(t,y1)
title('二次调频信号')