Chirp信号
Generate swept-frequency cosine (chirp) signal
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 secs @ 1kHz sample rate y = chirp(t,0,1,50); % Start @ DC 0Hz, % cross 50Hz at t=1 sec 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 secs @ 1kHz sample rate y = chirp(t,0,1,50, 'quadratic', 0, 'convex'); % Start @ DC 0Hz, % cross 50Hz at t=1 sec %y = chirp(t,0,1,50, 'quadratic', 0, 'concave'); % Start @ DC 0Hz, % cross 50Hz at t=1 sec 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 secs @ 1kHz sample rate %y = chirp(t,0,1,50, 'quadratic', 0, 'convex'); % Start @ DC 0Hz, % cross 50Hz at t=1 sec y = chirp(t,1e-6,1,50, 'logarithmic'); % Start @ DC 1e-6Hz, % cross 50Hz at t=1 sec 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
今天的文章chirp信号的频谱图_单频余弦信号是什么分享到此就结束了,感谢您的阅读,如果确实帮到您,您可以动动手指转发给其他人。
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。
如需转载请保留出处:https://bianchenghao.cn/45688.html