啁啾信号chirp（扫频余弦信号）
Chirp信号
Generate swept-frequency cosine (chirp) signal
春天景象
Chirp信号是⼀个典型的⾮平稳信号，在通信、声纳、雷达等领域具有⼴泛的应⽤，好看的爱情小说
Chirp译名：啁啾（读⾳：“周纠”），是通信技术有关编码脉冲技术中的⼀种术语，是指对脉冲进⾏编码时，其载频在脉冲持续时间内线性地增加，当将脉冲变到⾳频地，会发出⼀种声⾳，听起来像鸟叫的啁啾声，故名“啁啾”。
后来就将脉冲传输时中⼼波长发⽣偏移的现象叫做“啁啾”。例如在光纤通信中由于激光⼆极管本⾝不稳定⽽使传输单个脉冲时中⼼波长瞬时偏移的现象，也叫“啁啾”。
Chirp信号的表达式如下：
式中f0称作起始频率。
u0为调频率
对相位进⾏求导，得到⾓频率以及频率随时间的线性变化关系 f = f0+u0*t
matlab⽣成chirp信号
chirp函数
y = chirp(t,f0,t1,f1,’method’,phi,'shape')
根据指定的⽅法在时间t上产⽣余弦扫频信号，f0为初始时刻的瞬时频率，f1为t1时刻的瞬时频率，f0和f1单位都为Hz。如果未指定，f0默认为e-6（对数扫频⽅法）或0（其他扫频⽅法），t1为1，f1为100Hz。
扫频⽅法有linear线性扫频、quadratic⼆次扫频、logarithmic对数扫频；
phi允许指定⼀个初始相位（以°为单位），默认为0，如果想忽略此参数，直接设置后⾯的参数，可以指定为0或[]；
销售激励口号
shape指定⼆次扫频⽅法的抛物线的形状，凹还是凸，值为concave或convex，如果此信号被忽略，则根据f0和f1的相对⼤⼩决定是凹还是凸。
蔡和蔡线性扫频
t = 0:0.001:1; % 1 cs @ 1kHz sample rate
y = chirp(t,0,1,50); % Start @ DC 0Hz, % cross 50Hz at t=1 c
plot(t, y)
figure
% spectrogram(x,window,noverlap,nfft,fs)
spectrogram(y,256,250,1024,1000,'yaxis') %% nfft=1024, fs=1000Hz, window = 256, numoverlap = 250
%specgram(y,1024,1000,256,250) % nfft=1024, fs=1000Hz, window = 256, numoverlap = 250
短时傅⾥叶变换
specgram(y, nfft, fs, window, numoverlap) 或 spectrogram(y, window, numoverlap, nfft, fs, 'yaxis') nfft是fft的长度，越长的话，频域分辨率越⾼。
fs就是采样率；
window是指窗的长度，⼀般和nfft相同即可；
姨侄numoverlap是指nfft减去步长，越⼤越好，但运算量越⼤。⼀般取nfft的3/4效果就⽐较好了。
⼆次扫频
t = 0:0.001:1; % 1 cs @ 1kHz sample rate
y = chirp(t,0,1,50, 'quadratic', 0, 'convex'); % Start @ DC 0Hz, % cross 50Hz at t=1 c
%y = chirp(t,0,1,50, 'quadratic', 0, 'concave'); % Start @ DC 0Hz, % cross 50Hz at t=1 c
怎么提高自信
plot(t, y)
figure
% spectrogram(x,window,noverlap,nfft,fs)
spectrogram(y,256,250,1024,1000,'yaxis') %% nfft=1024, fs=1000Hz, window = 256, numoverlap = 250
对数扫频
t = 0:0.001:1; % 1 cs @ 1kHz sample rate
%y = chirp(t,0,1,50, 'quadratic', 0, 'convex'); % Start @ DC 0Hz, % cross 50Hz at t=1 c
y = chirp(t,1e-6,1,50, 'logarithmic'); % Start @ DC 1e-6Hz, % cross 50Hz at t=1 c
plot(t, y)
figure品茗杯
xxx按摩
% spectrogram(x,window,noverlap,nfft,fs)
spectrogram(y,256,250,1024,1000,'yaxis') %% nfft=1024, fs=1000Hz, window = 256, numoverlap = 250
