Chirp信號(hào)

Generate swept-frequency cosine (chirp) signal

Chirp信號(hào)是一個(gè)典型的非平穩(wěn)信號(hào)，在通信、聲納、雷達(dá)等領(lǐng)域具有廣泛的應(yīng)用，

Chirp譯名:啁啾(讀音：“周糾”)，是通信技術(shù)有關(guān)編碼脈沖技術(shù)中的一種術(shù)語，是指對(duì)脈沖進(jìn)行編碼時(shí)，其載頻在脈沖持續(xù)時(shí)間內(nèi)線性地增加，當(dāng)將脈沖變到音頻地，會(huì)發(fā)出一種聲音，聽起來像鳥叫的啁啾聲，故名“啁啾”。

后來就將脈沖傳輸時(shí)中心波長發(fā)生偏移的現(xiàn)象叫做“啁啾”。例如在光纖通信中由于激光二極管本身不穩(wěn)定而使傳輸單個(gè)脈沖時(shí)中心波長瞬時(shí)偏移的現(xiàn)象，也叫“啁啾”。

Chirp信號(hào)的表達(dá)式如下：

式中f0稱作起始頻率。

u0為調(diào)頻率

對(duì)相位進(jìn)行求導(dǎo)，得到角頻率以及頻率隨時(shí)間的線性變化關(guān)系 f = f0+u0*t

matlab生成chirp信號(hào)

chirp函數(shù)

y = chirp(t,f0,t1,f1,’method’,phi,'shape')

根據(jù)指定的方法在時(shí)間t上產(chǎn)生余弦掃頻信號(hào)，f0為初始時(shí)刻的瞬時(shí)頻率，f1為t1時(shí)刻的瞬時(shí)頻率，f0和f1單位都為Hz。如果未指定，f0默認(rèn)為e-6(對(duì)數(shù)掃頻方法)或0(其他掃頻方法)，t1為1，f1為100Hz。

掃頻方法有l(wèi)inear線性掃頻、quadratic二次掃頻、logarithmic對(duì)數(shù)掃頻；

phi允許指定一個(gè)初始相位(以°為單位)，默認(rèn)為0，如果想忽略此參數(shù)，直接設(shè)置后面的參數(shù)，可以指定為0或[]；

shape指定二次掃頻方法的拋物線的形狀，凹還是凸，值為concave或convex，如果此信號(hào)被忽略，則根據(jù)f0和f1的相對(duì)大小決定是凹還是凸。

線性掃頻

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

短時(shí)傅里葉變換

specgram(y,?nfft, fs, window,?numoverlap) 或 spectrogram(y, window, numoverlap, nfft, fs,?'yaxis')

nfft是fft的長度，越長的話，頻域分辨率越高。

fs就是采樣率；

window是指窗的長度，一般和nfft相同即可；

numoverlap是指nfft減去步長，越大越好，但運(yùn)算量越大。一般取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

對(duì)數(shù)掃頻

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

總結(jié)

以上是生活随笔為你收集整理的matlab啁啾信号,啁啾信号chirp（扫频余弦信号）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。