mscohere
幅度平方相干
语法
描述
cxy=mscohere(x，y）cxy，输入信号的，x和y．
如果
x和y都是向量，它们的长度必须相同。如果一个信号是矩阵，另一个是向量，那么向量的长度必须等于矩阵中的行数。该函数展开向量并返回一个逐列幅值平方相干估计矩阵。
如果
x和y矩阵的行数相同但列数不同，那么mscohere返回多重一致性矩阵。这个米第列cxy包含所有输入信号与输入信号之间的相关程度的估计值米输出信号。看到平方的一致性了解更多信息。如果
x和y那么矩阵的大小是相等的吗mscohere操作列:cxy（：，n）=mscohere（x（：，n），y（：，n））．要获得多重相干矩阵，请附加“mimo”到参数列表。
mscohere(＿＿＿）在没有输出参数的情况下，在当前图形窗口中绘制幅值平方相干估计。
例子
两个序列的相干估计
计算并绘制两个有色噪声序列之间的相干性估计。
产生由高斯白噪声组成的信号。
r=randn（16384,1）；
要创建第一个序列，带通滤波器的信号。设计一个16阶滤波器，通过0.2 ~ 0.4π rad/sample之间的归一化频率。指定一个60 dB的阻带衰减。过滤原始信号。
dx=设计过滤器(“bandpassiir”，“过滤器订单”,16,...“阻带频率1”, 0.2,“阻带频率2”,0.4,...“止损绷带”、60);x =过滤器(dx, r);
为了创建第二个序列，设计一个16阶滤波器，使归一化频率停止在0.6π和0.8πrad/sample之间。指定0.1 dB的通带纹波。过滤原始信号。
dy=设计过滤器(“bandstopiir”，“过滤器订单”,16,...“通带频率1”,0.6,“通带频率2”,0.8,...“PassbandRipple”, 0.1);y =过滤器(dy, r);
估计的振幅的平方相干性x和y.使用512个样本汉明窗口。指定相邻线段和2048个DFT点之间重叠的500个样本。
[cxy, fc] = mscohere (x, y,汉明(512),500年,2048年);
绘制相干函数并叠加滤波器的频率响应。
[qx，f]=freqz（dx）；qy=freqz（dy）；绘图（fc/pi，cxy）保持在…上情节(f /π,abs(季度)、f /π,abs (qy))关
多重连贯与普通连贯
生成随机双通道信号，x. 产生另一个信号，y，通过对两个通道进行低通滤波并将其相加。指定一个截止频率为0.3π并使用矩形窗口设计的30阶FIR滤波器。
h=fir1（30,0.3，rectwin（31））；x=randn（16384,2）；y=和（滤波器（h，1，x），2）；
计算的多相干估计x和y.使用1024个样本Hann窗口对信号进行窗口处理。指定相邻段和1024个DFT点之间重叠的512个样本。绘制估计值。
noverlap=512；nfft=1024；mscohere（x，y，hann（nfft），noverlap，nfft，“mimo”）
将相干性估计值与滤波器的频率响应进行比较。相干性下降对应于频率响应的零点。
[H，f]=freqz（H）；保持在…上yyaxis正确的图（f/pi，20*log10（abs（H））保持关
计算并绘制地震波的普通震级平方相干估计x和y. 任何信道的估计值均未达到1。
图mscohere (x, y,损害(nfft) noverlap, nfft)
MIMO系统的相干性
产生两个多通道信号，每个采样在1 kHz 2秒。第一个信号，即输入，由三个频率为120 Hz、360 Hz和480 Hz的正弦波组成。第二个信号是输出信号，由两个频率为120 Hz和360 Hz的正弦波组成。其中一个正弦波滞后于第一个信号π/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-[pi/2 pi/4]）；oupt=oupt+randn（size（oupt））；
估计所有输入信号和每个输出通道之间的相关度。使用长度为100的汉明窗口对数据进行窗口处理。mscohere为每个输出通道返回一个相干函数。相干函数在输入和输出共享的频率处达到最大值。
[Cxy f] = mscohere (inpt oupt,汉明(100),[],[],fs,“mimo”);为k = 1:尺寸(oupt, 2)次要情节(大小(oupt, 2), 1, k)情节(f, Cxy (:, k))标题([“输出”int2str（k）'所有输入'])结束
切换输入和输出信号并计算多重相干函数。使用相同的汉明窗口。在480 Hz时，输入和输出之间没有相关性。因此，第三个相关函数中没有峰值。
[Cxy，f]=mscohere（oupt，inpt，hamming（100），[]，[]，fs，“mimo”);为k=1:size（inpt，2）子地块（size（inpt，2），1，k）绘图（f，Cxy（：，k））标题([“输入”int2str（k）'所有输出'])结束
使用的打印功能重复计算mscohere．
clf mscohere (oupt inpt,汉明(100),[],[],fs,“mimo”）
计算第二个信号的普通相干函数和第一个信号的前两个通道。非峰值不同于多重相干函数。
[Cxy，f]=mscohere（oupt，inpt（：，[12]），hamming（100），[]，[]，fs）；plot（f，Cxy）
通过计算最大相干性点处的交叉谱角来求相位差。
Pxy=cpsd（输出，输入（：，[12]），hamming（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
修正幅度平方相干图
生成两个正弦信号，每个信号以1kHz的频率采样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（嗯，联合国，[]，[]，[]，财政司司长）
修改绘图的标题、绘图的标签x-轴，以及y-轴心国。
头衔(“平方一致性”)xlabel(“f (Hz)”)ylim（[01.1]）
使用gca获取当前坐标轴的句柄。更改标记的位置。去掉标签y-轴心国。
ax=gca；ax.XTick=0:250:500；ax.YTick=0:0.25:1；ax.YLabel.String=[]；
打电话给孩子们属性更改打印线的颜色和宽度。
ln=ax.儿童；ln.Color=[0.80]；ln.线宽=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-2010.第3927-3941页。
Kay, Steven M。现代谱估计。恩格尔伍德悬崖，新泽西州：普伦蒂斯大厅，1988年。
[3] 拉宾纳、劳伦斯R.和伯纳德·戈尔德。数字信号处理的理论与应用。恩格尔伍德悬崖，新泽西州：普伦蒂斯大厅，1975年。
[4] 斯托伊卡、彼得和伦道夫·摩西。信号的频谱分析。上鞍河，新泽西州：普伦蒂斯大厅，2005年。
[5] 快速傅里叶变换在功率谱估计中的应用：一种基于短时间平均的方法，修正的周期图IEEE®音频和电声学学报。《非洲联盟》第15卷，1967年，第70-73页。
扩展能力
您还可以从以下列表中选择网站：
