1 MHzで10 ms間サンプリングした2の線形チャプを生成します。

S/N比が40 dBとなるホワaaplトガウスノaaplズを付加します。

nSamp = 10000;Fs = 1000e3;信噪比= 40;t = (0:nSamp-1)'/Fs;X1 = chirp(t,150e3,t(end)，350e3);x1 = x1+randn(size(x1))*std(x1)/db2mag(SNR);X2 = chirp(t,200e3,t(end)，300e3);x2 = x2+randn(size(x2))*std(x2)/db2mag(SNR);

2 .。信号を200サンプルのセグメントに分割し，各セグメントにハミングウィンドウを適用します。隣接するセグメント間の80サンプルのオーバーラップ,および1024サンプルのDFT長を指定します。

xspectrogram (x1, x2,汉明(200),80年,1024年,Fs,“桠溪”）

2番目のチャープを変更し,測定中に周波数が50千赫から350 kHzに増加するようにします。形状係数 $$\beta ＝5$$で500サンプルのカザウィンドウを使用して，セグメントにウィンドウを適用します。450サンプルのオ，バ，ラップおよび256のDFT長を指定します。クロススペクトログラムを計算してプロットします。

X2 = chirp(t,50e3,t(end)，350e3);x2 = x2+randn(size(x2))*std(x2)/db2mag(SNR);xspectrogram (x1, x2,凯瑟(500 5),450256年,Fs,“桠溪”）

両方の場合にいて，関数は2。

44100 Hzでサンプリングされた2の音声信号を含むファルを読み込みます。

2の信号をプロットします。2番目の信号を垂直方向にオフセットして、両方が表示されるようにします。

负载(“voice.mat”）要收听，输入soundsc(transform,fs)，pause(2)，soundsc(reform,fs)T = (0:length(reform)-1)/fs;情节(t,变换,t,改革+ 0.3)传说(”“变换函数”，“改革司法”的）

2の信号のクロススペクトログラムを計算します。信号を1000サンプルのセグメントに分割し，ハミングウィンドウを適用します。隣接するセグメント間に800個のサンプルのオ，バ，ラップを指定します。4 kHzまでの周波数のみを含めます。

Nwin = 1000;NVLP = 800;Fint = 0:4000;[s,f,t] = xspectrogram(transform,reform,hamming(nwin)，nvlp,fint,fs);网格(t,f,20*log10(s))视图(2)轴紧

クロススペクトログラムは，信号がより多くの周波数成分を共通しても時間間隔を強調表示します。“form”という音節が特に顕著になります。

それぞれ1 kHzで2秒間サンプリングされた2の2次チャプを生成します。両方のチャ，プの初期周波数は100hzで，測定の半ばには200hzに増加します。2番目のチャ，プの最初のチャ，プに対する位相差は23°です。

Fs = 1e3;T = 0:1/fs:2;Y1 = chirp(t,100, 1200，“二次”, 0);Y2 = chirp(t,100, 1200，“二次”、23);

チャ，プの複素クロススペクトログラムを計算し，その間の位相シフトを抽出します。信号を128サンプルのセグメントに分割します。隣接するセグメント間に120個のサンプルのオバラップを指定します。形状係数がβ= 18のカイザーウィンドウを使用して各セグメントをウィンドウ処理し,128サンプルのDFT長を指定します。xspectrogramのプロット機能を使用して，クロススペクトログラムを表示します。

[~ f, xt, c] = xspectrogram (y1, y2,皇帝(128年,18),120128年,fs);xspectrogram (y1、y2、皇帝(128年,18),120128年,fs,“桠溪”）

クロススペクトログラムの最大エネルギの時間—周波数リッジを抽出して表示します。

[tfr，~，lridge] = tfridge(c,f);持有在情节(xt,总和生育率,“k”，“线宽”, 2)从

位相シフトは，リッジに沿った時変クロススペクトルの虚数部と実数部の比です。位相シフトを計算して，度単位で表します。平均値を表示します。

PSHFT =角(c(lridge))*180/pi;意思是(pshft)

Ans = -23.0000

3 kHzでそれぞれ1秒間サンプリングされた2の信号を生成します。最初の信号は2次チャープで,測定中に周波数が300 Hzから1300 Hzに増加します。チャプはホワトガウスノズに組み込まれます。2番目の信号もホワイト ノイズに組み込まれますが、これは正弦関数的に変化する周波数成分をもつチャープです。

Fs = 3000;T = 0:1/fs:1-1/fs;X1 = chirp(t,300,t(end)，1300，“二次”) + randn(大小(t)) / 100;x2 = exp (2 j *π* 100 * cos(2 *π* 2 * t)) + randn(大小(t)) / 100;

2 .の信号のクロススペクトログラムを計算し，プロットします。隣接するセグメント間で255サンプルがオーバーラップする256サンプルのセグメントに信号を分割します。形状係数がβ= 30のカ。既定のDFTの点の数を使用します。クロススペクトログラムの中央を周波数ゼロに揃えます。

Nwin = 256;xspectrogram (x1, x2, kaiser (nwin 30), nwin-1, [], fs,“中心”，“桠溪”）

パワ，スペクトル密度の代わりにパワ，スペクトルを計算します。-40 dBより小さい値をゼロに設定します。プロットの中央をナ@ @キスト周波数に揃えます。

xspectrogram (x1, x2, kaiser (nwin 30), nwin-1, [], fs,.．.“权力”，“MinThreshold”, -40,“桠溪”)标题(二次啁啾与复啁啾的交叉谱图）

しきい値処理により，共通周波数の領域がさらに強調表示されます。

2 .。

各シケンスを4096サンプルの長さに指定します。

N = 4096;

最初のシーケンスを作成するために,ホワイトガウスノイズに組み込まれる凸二次チャープを生成し,バンドパスフィルターを適用します。

rx = chirp(0:N-1,0.1/2,N,0.8/2，“二次”[],“凸”) + randn (N - 1) / 100;Dx = designfilt(“bandpassiir”，“FilterOrder”, 16岁,.．.“StopbandFrequency1”, 0.2,“StopbandFrequency2”, 0.4,.．.“StopbandAttenuation”、60);X = filter(dx,rx);

2番目のシーケンスを作成するために,ホワイトガウスノイズに組み込まれる線形チャープを生成し,バンドストップフィルターを適用します。

ry = chirp(0:N-1,0.9/ 2,n,0.1/2)'+randn(N,1)/100;Dy = designfilt(“bandstopiir”，“FilterOrder”, 16岁,.．.“PassbandFrequency1”, 0.6,“PassbandFrequency2”, 0.8,.．.“PassbandRipple”1);Y = filter(dy,ry);

2 .のシ，ケンスをプロットします。2番目のシーケンスを垂直方向にオフセットして、両方が表示されるようにします。

情节([x y + 2])

xとyのクロススペクトログラムを計算し，プロットします。512サンプルのハミングウィンドウを使用します。隣り合ったセグメント間のオ，バ，ラップを500サンプル，dft点を2048に指定します。

xspectrogram (x, y,汉明(512),500年,2048年,“桠溪”）

-50 dBより小さいクロススペクトログラム値をゼロに設定します。

xspectrogram (x, y,汉明(512),500年,2048年,“MinThreshold”, -50,“桠溪”）

スペクトログラムに，フィルタ，によって強調または抑制された周波数領域が示されます。