mscohere
平方的一致性
语法
描述
cxy= mscohere (x，y）cxy的输入信号，x而且y．
如果
x而且y都是向量，它们的长度必须相同。如果其中一个信号是一个矩阵，另一个是一个向量，那么向量的长度必须等于矩阵中的行数。该函数展开矢量并返回逐列的平方相干估计矩阵。
如果
x而且y矩阵的行数相同但列数不同吗mscohere返回一个多重相干矩阵。的米的第Th列cxy包含了所有输入信号与的相关程度的估计值米输出信号。看到平方的一致性获取更多信息。如果
x而且y那么矩阵的大小是否相等呢mscohere操作列:Cxy (:，n) = mscohere(x(:，n)，y(:，n))．若要获得多重相干矩阵，请添加“再分配”到参数列表。
mscohere (＿＿＿）在没有输出参数的情况下，在当前图形窗口中绘制平方相干估计。
例子
两个序列的相干性估计
计算并绘制两个彩色噪声序列之间的相干估计。
生成一个由高斯白噪声组成的信号。
R = randn(16384,1);
为了创建第一个序列，对信号进行带通滤波。设计一个16阶滤波器，通过归一化频率在0.2π和0.4π rad/sample之间。指定阻带衰减为60 dB。对原始信号进行滤波。
Dx = designfilt(“bandpassiir”，“FilterOrder”, 16岁,.．.“StopbandFrequency1”, 0.2,“StopbandFrequency2”, 0.4,.．.“StopbandAttenuation”、60);X = filter(dx,r);
为了创建第二个序列，设计一个16阶滤波器，停止0.6π和0.8π rad/sample之间的归一化频率。指定0.1 dB的通带纹波。对原始信号进行滤波。
Dy = designfilt(“bandstopiir”，“FilterOrder”, 16岁,.．.“PassbandFrequency1”, 0.6,“PassbandFrequency2”, 0.8,.．.“PassbandRipple”, 0.1);Y = filter(dy,r);
估计的平方相干性x而且y．使用512个样本的汉明窗口。指定相邻段之间的重叠样本500个，DFT点2048个。
[cxy,fc] = mscohere(x,y,hamming(512)，500,2048);
绘制相干函数，叠加滤波器的频率响应。
[qx,f] = freqz(dx);Qy = freqz(dy);情节(fc /π,cxy)在情节(f /π,abs(季度)、f /π,abs (qy))从
多重相干与普通相干
生成一个随机的双通道信号，x．生成另一个信号，y，通过低通滤波两个通道并将它们相加。指定一个30阶FIR滤波器，其截止频率为0.3π，并使用矩形窗口设计。
H = fir1(30,0.3,rectwin(31));X = randn(16384,2);Y = sum(filter(h,1,x)，2);
计算的多相干估计x而且y．用1024样本Hann窗口对信号进行窗口化。指定相邻段之间重叠的512个样本和1024个DFT点。绘制估算图。
Noverlap = 512;NFFT = 1024;mscohere (x, y,损害(nfft) noverlap, nfft,“再分配”）
将相干估计与滤波器的频率响应进行比较。相干性的下降对应于频率响应的零点。
[H,f] = freqz(H);持有在yyaxis正确的情节(f /π,20 * log10 (abs (H)))从
计算和绘制的普通平方相干估计x而且y．任何通道的估计值都不会达到1。
图mscohere (x, y,损害(nfft) noverlap, nfft)
MIMO系统的相干性
生成两个多通道信号，每个信号以1千赫采样2秒。第一个信号，即输入信号，由三个频率分别为120hz、360hz和480hz的正弦波组成。第二个信号，即输出信号，由两个频率分别为120hz和360hz的正弦波组成。其中一个正弦信号滞后于第一个信号π/2。另一个正弦波有π/4的滞后。两个信号都嵌入高斯白噪声中。
Fs = 1000;F = 120;T = (0:1/fs:2-1/fs)';Inpt = sin(2*pi*f*[1 3 4].*t);Inpt = Inpt +randn(size(Inpt));Oupt = sin(2*pi*f*[1 3]。* t -(π/ 2π/ 4));Oupt = Oupt +randn(size(Oupt));
估计所有输入信号和每个输出通道之间的相关程度。使用长度为100的汉明窗口来显示数据。mscohere为每个输出通道返回一个相干函数。相干函数在输入和输出共享的频率处达到最大值。
[Cxy f] = mscohere (inpt oupt,汉明(100),[],[],fs,“再分配”）;为k = 1:尺寸(oupt, 2)次要情节(大小(oupt, 2), 1, k)情节(f, Cxy (:, k))标题([“输出”int2str (k)'，所有输入']）结束
切换输入输出信号，计算多重相干函数。使用相同的汉明窗口。在480hz时，输入输出无相关性。因此在第三个相关函数中没有峰值。
[Cxy f] = mscohere (oupt inpt,汉明(100),[],[],fs,“再分配”）;为k = 1:尺寸(inpt, 2)次要情节(大小(inpt, 2), 1, k)情节(f, Cxy (:, k))标题([“输入”int2str (k)'，所有输出']）结束
的绘图功能重复计算mscohere．
clf mscohere (oupt inpt,汉明(100),[],[],fs,“再分配”）
计算所述第二信号和所述第一信号的前两个信道的普通相干函数。非峰值值不同于多重相干函数。
[Cxy f] = mscohere (oupt inpt(:,(1 2)),汉明(100),[],[],fs);情节(f, Cxy)
通过计算交叉谱在最大相干点处的夹角来求相位差。
Pxy =运行cpsd (oupt inpt(:,(1 2)),汉明(100),[],[],fs);[~，mxx] = max(Cxy);为K = 1:2 fprintf('相位滞后%d = %5.2f*pi\n', k,角(Pxy (mxx (k), k)) / pi)结束
相位滞后1 = -0.51*pi相位滞后2 = -0.22*pi
修改平方相干图
生成两个正弦信号采样1秒，每个在1千赫。每个正弦波的频率为250hz。其中一个信号在相位上比另一个滞后π/ 3弧度。将两个信号都嵌入单位方差的高斯白噪声中。
Fs = 1000;F = 250;T = 0:1/fs:1-1/fs;Um = sin(2*pi*f*t)+rand(size(t));UN = sin(2*pi*f*t-pi/3)+rand(size(t));
使用mscohere计算和绘制信号的平方相干度。
mscohere(嗯,联合国 ,[],[],[], fs)
修改标题的情节，标签的x的极限y设在。
标题(“平方一致性”)包含(“f (Hz)”) ylim([0 1.1])
使用gca获取当前轴的句柄。更改标记的位置。的标签取下y设在。
Ax = gca;斧子。XTick = 0:250:500; ax.YTick = 0:0.25:1; ax.YLabel.String = [];
调用孩子们属性更改绘制的线条的颜色和宽度。
ln = ax.儿童;ln。Color = [0.8 0 0]; ln.LineWidth = 1.5;
另外,使用集而且得到修改行属性。
集(get (gca),“孩子”)，“颜色”，[0 0.4 0]，“线型”，“——”，“线宽”，1)
输入参数
输出参数
更多关于
参考文献
[1] Gómez González, A.， J. Rodríguez, X. Sagartzazu, A. Schumacher，和I. Isasa。非平稳信号噪声和振动传播路径分析的时域多重相干法2010年噪声与振动工程国际会议论文集，ISMA2010-USD2010．3927 - 3941页。
[2]凯，史蒂文M。现代谱估计。恩格尔伍德悬崖，新泽西州:Prentice-Hall, 1988。
Rabiner, Lawrence R.和Bernard Gold。数字信号处理理论与应用“，”恩格尔伍德悬崖，新泽西州:Prentice-Hall, 1975。
Stoica, Petre和Randolph Moses。信号的光谱分析。上马鞍河，新泽西州:Prentice Hall, 2005。
[5] Welch, Peter D.“使用快速傅里叶变换估计功率谱:一种基于时间平均的短修正周期图的方法。”IEEE®音频与电声学汇刊。卷AU-15, 1967，第70-73页。
扩展功能
您也可以从以下列表中选择一个网站:
