生成两个有色噪声信号和情节交叉功率谱密度。指定一个长度- 1024 500点FFT和三角窗没有重叠。

r = randn (16384 1);0.2 hx = fir1(30日,rectwin (31);x =过滤器(hx 1 r);hy = (10) /√(10);y =过滤器(已有1 r);运行cpsd (x, y, triang(500), 250年,1024年)

生成两个双通道正弦信号采样1 kHz 1秒。第一个信号通道的频率200赫兹和300赫兹。第二个信号通道的频率300赫兹和400赫兹。这两个信号嵌入在unit-variance高斯白噪声。

fs = 1 e3;t = (0:1 / fs: 1 - 1 / fs) ';q = 2 * sin(2 *π* (200 300)。* t);q = q + randn(大小(q));r = 2 * sin(2 *π* (300 400)。* t);r = r + randn(大小(r));

计算交叉两个信号的功率谱密度。使用256 -样本Bartlett窗口把信号分成段和窗口的段。指定128个样本之间的重叠相邻段和DFT 2048点。使用内置的功能运行cpsd绘制结果。

运行cpsd (q, r,巴特利特(256),128年,2048年,fs)

默认情况下,运行cpsd作品column-by-column矩阵输入相同的大小。每个通道的峰值频率的正弦曲线。

重复计算,但现在追加“再分配”的参数列表。

运行cpsd (q, r,巴特利特(256),128年,2048年,fs,“再分配”)

当被称为“再分配”选项,运行cpsd返回一个三维数组,其中包含交叉功率谱密度估计输入列的所有组合。第二通道的估计问和第一个频道r显示了一个增强的峰值在300赫兹的频率。

生成两个100 Hz正弦信号在1千赫采样296 ms。一个正弦曲线落后于其他2.5毫秒,相当于的相位滞后π/ 2。两个信号嵌入在高斯白噪声的方差1/4²。

Fs = 1000;t = 0:1 / Fs: 0.296;x = cos(2 *π* t * 100) + 0.25 * randn(大小(t));τ= 1/400;y = cos(2 *π* 100 * (t-tau)) + 0.25 * randn(大小(t));

计算和情节交叉功率谱密度的大小。使用默认设置运行cpsd。大小山峰在频率信号之间存在显著的一致性。

运行cpsd (x, y, [] [], [], Fs)

情节平方相干函数和交叉谱的阶段。的纵坐标high-coherence频率对应于正弦信号之间的相位滞后。

[Cxy F] = mscohere (x, y, []、[] [], Fs);[Pxy F] =运行cpsd (x, y, []、[] [], Fs);次要情节(2,1,1)情节(F, Cxy)标题(“平方一致性”次要情节(2,1,2)情节(F,角(Pxy))在情节(F, 2 *π* 100 *τ*的(大小(F)),“——”)举行从包含(“赫兹”)ylabel (“\θ(f)”)标题(“交叉谱阶段”)

生成两个 $$N$$样本指数序列, $$x_{\mathrm{一个}}={\mathrm{一个}}^n$$和 $$x_b=b^n$$, $$n\ge 0$$。指定 $$\mathrm{一个}=0。8$$, $$b=0。9$$,和一个小 $$N$$——的尺寸效应。

N = 10;n = 0: n - 1;一个= 0.8;b = 0.9;xa = a。^ n;xb = b。^ n;

计算的交叉功率谱密度和情节完整的归一化频率间隔序列, $$(-\; -\; -\; -\; -\; -\pi ,\pi ]$$。指定一个矩形窗口的长度 $$N$$和段之间没有重叠。

w = -π:1/1000:π;风= rectwin (N);小说= 0;[pxx f] =运行cpsd (xa, xb,风、诺w);

跨两个序列的功率谱的解析表达式 $$N$$:

这个表达式转换成一个十字架功率谱密度除以 $$2\pi N$$。比较结果。的涟漪运行cpsd结果是一个窗口的结果。

nfac = 2 *π* N;X = 1. / (1 a * exp (1 j * w));Y = 1。/ (1 b * exp (1 j * w));R = x * Y / nfac;semilogy (f /π,abs (pxx))在semilogy (w /π,abs (R))从传奇(“运行cpsd”,“分析”)

重复的计算 $$N=25$$。曲线同意六位数 $$N$$小100。

N = 25;n = 0: n - 1;xa = a。^ n;xb = b。^ n;风= rectwin (N);[pxx f] =运行cpsd (xa, xb,风、诺w);R = x * Y /(2 *π* N);semilogy (f /π,abs (pxx))在semilogy (w /π,abs (R))从传奇(“运行cpsd”,“分析”)

使用交叉功率谱密度来确定一个高度的基调。

产生的声音信号,当你拨一个号码或符号在数字电话的正弦曲线频率来自两个不同的组。每一对音调低组的包含一个频率(697 Hz, 770 Hz, 852 Hz, 941 Hz)和一个频率高集团(1209 Hz, 1336 Hz, 1477 Hz)。

产生信号对应的所有符号。样品在4 kHz每个音调半秒钟。准备一个引用表。

fs = 4 e3;t = 0:1 / fs: 0.5 - 1 / f;nms = [' 1 ';' 2 ';“3”;“4”;“5”;“6”;“7”;“8”;“9”;‘*’;' 0 ';“#”];版本= (697 770 852 941);贺南洪= (1209 1336 1477);v =长度(版本);h =长度(小时);为k = 1: v为l = 1: h idx = h * (k - 1) + l;语气=总和(罪(2 *π*[版本(k),贺南洪(l)]。* t)) ';音调(:,idx) =基调;结束结束

情节的韦尔奇周期图和注释组件频率信号。使用200 -样本汉明窗将信号划分为互不重叠的段段和窗口。

[pxx f] = pwelch(音调,汉明(200),0,[],fs);为k = 1: v为l = 1: h idx = h * (k - 1) + l;ax =次要情节(v、h、idx);情节(f, 10 * log10 (pxx (:, idx))) ylim(-80[0])标题(nms (idx)) tx =[版本(k),贺南洪(l)];斧子。XTick = tx;斧子。XTickLabel = int2str (tx);结束结束

发送一个信号由拨号数字8通过噪声信道。接收到的信号是如此腐败,不能被检查。

mys =总和(罪(2 *π*[版本(3),贺南洪(2)]。* t))的+ 5 * randn(大小(t));%,类型soundsc(岬,fs)

计算交叉破坏信号的功率谱密度和参考信号。窗口信号用一个样品512 - Kaiser窗形状系数β= 5。每个光谱图的大小。

[pxy f] =运行cpsd(岬,音调,凯瑟(512 5),100年,[],fs);为k = 1: v为l = 1: h idx = h * (k - 1) + l;ax =次要情节(v、h、idx);情节(f, 10 * log10 (abs (pxy (:, idx)))) ylim(-80[0])标题(nms (idx)) tx =[版本(k),贺南洪(l)];斧子。XTick = tx;斧子。XTickLabel = int2str (tx);结束结束

损坏的数字信号的频谱最高的山峰和均方根值最高。

(~,loc) = max (rms (abs (pxy)));数字= nms (loc)

数字= ' 8 '