mscohere
震级的相干性
Sintaxis
Descripcion
cxy= mscohere (x，y）cxy， de las señales de entrada，xey．
如果
xey两个两个矢量，一个是纵向的，一个是纵向的。Si una de las señales es una矩阵通过la otra向量，la longitude del vector通过ser = número de filas de la matrix。La función扩展了el向量，通过发展unmatriz de estimaciones de coherencia de magnitude cuadrada column - por column。
如果
xey子矩阵conel mismo número de filas, pero conente número de columns，mscoheredevelop una matriz de coherencia múltiple。洛杉矶一列圆柱米德cxy中国大陆estimación del grado de correlación中心今天将通过señales de entrada y la señal de salida米．Para obtener más información，咨询震级的相干性．如果
xey子矩阵del mismo tamaño，mscohere歌剧院的圆柱;Cxy (:，n) = mscohere(x(:，n)，y(:，n))．相干矩阵的参数计算múltiple, añada“再分配”A la lista de argumenttos。
mscohere (＿＿＿）从理论上讲，我们可以从estimación的角度出发，我们可以从实际的角度出发，我们可以从实际的角度出发。
包括
Estimación de la coherencia de dos secencia
计算代表了la estimación de la coherencia中心dos secuencias de ruido de color。
generere una señal que consista en rudo blanco gaussiano。
R = randn(16384,1);
Para crear la premiersecuencia, utiltunfilter, paso banda en la señal。Diseñe unfiltro de 16。按顺序排列的频率为0,2π到0,4π / 1。特别是una atenuación de la banda de parada de 60 dB。过滤la señal原创。
Dx = designfilt(“bandpassiir”，“FilterOrder”, 16岁,…“StopbandFrequency1”, 0.2,“StopbandFrequency2”, 0.4,…“StopbandAttenuation”、60);X = filter(dx,r);
Para crear la secunda secencia, diseñe unfilter de 16。按顺序排列的频率为0,6π到0,8π /m。特别是una ondulación de banda de paso de 0,1 dB。过滤la señal原创。
Dy = designfilt(“bandstopiir”，“FilterOrder”, 16岁,…“PassbandFrequency1”, 0.6,“PassbandFrequency2”, 0.8,…“PassbandRipple”, 0.1);Y = filter(dy,r);
Estime la coherencia de la magnitude cuadrada dexey．汉明的公共设施，512个城市。具体来说，500个太阳能建筑将被分配到2048个DFT建筑中。
[cxy,fc] = mscohere(x,y,hamming(512)，500,2048);
用superponga表示función de cohercia(相干性)。
[qx,f] = freqz(dx);Qy = freqz(dy);情节(fc /π,cxy)在情节(f /π,abs(季度)、f /π,abs (qy))从
Coherencia múltiple y Coherencia ordinaria
generere una señal运河管理局，x．Genere otra señal;y， filtrando en paso bajo los DOS运河通过sumándolos。特别适用于FIR 30。orden con una frecucia de corte de 0,3π y diseñado con una ventana矩形。
H = fir1(30,0.3,rectwin(31));X = randn(16384,2);Y = sum(filter(h,1,x)，2);
微积分estimación de la coherencia múltiple dexey．Aplique una ventana de Hann de 1024 muestras a las señales。具体来说，512个单元的solapapamiito是由1024个单元的DFT组成的。代表la estimación。
重叠= 512;NFFT = 1024;mscohere (x, y,损害(nfft) noverlap, nfft,“再分配”）
比较了estimación de la coherencia与la respuesta，甚至是freucia del filter。Las caídas de coherencia se对应于los ceros de la respuesta en freucia。
[H,f] = freqz(H);持有在yyaxis正确的情节(f /π,20 * log10 (abs (H)))从
计算表示为estimación de la coherencia de magnitude cuadrada ordinaria dexey．La estimación no lega a 1 en ninguno de los canales。
图mscohere (x, y,损害(nfft) noverlap, nfft)
系统MIMO的相干性
generere dos señales multianal, cadada una de ellas muestreada 1千赫持续2秒。La primera señal, La de entrada, est
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* *f*[1 3])* t -(π/ 2π/ 4));Oupt = Oupt +randn(size(Oupt));
esttime el grado de correlación entre todas as señales de entrada y cado de los canales de salida。海明经度为100度。mscohere发展una función de coherencia para cada canal de salida。Las funcciones de coherencia alcanzan los máximos en Las frerecuencias comppartidas por la entrada by la salida。
[Cxy f] = mscohere (inpt oupt,汉明(100),[],[],fs,“再分配”）;为k = 1:尺寸(oupt, 2)次要情节(大小(oupt, 2), 1, k)情节(f, Cxy (:, k))标题([“输出”int2str (k)'，所有输入'])结束
Cambie las señales de entrada y salida y culculla función de cohercoherence múltiple。用海明的工具。No hay correlación entre la entrada y la salida a 480 Hz。可怜的tanto，没有干草picos en la tercera función de correlación。
[Cxy f] = mscohere (oupt inpt,汉明(100),[],[],fs,“再分配”）;为k = 1:尺寸(inpt, 2)次要情节(大小(inpt, 2), 1, k)情节(f, Cxy (:, k))标题([“输入”int2str (k)'，所有输出'])结束
Repita el cálculo, utilitando la functionalidad de representación gráfica demscohere．
clf mscohere (oupt inpt,汉明(100),[],[],fs,“再分配”）
微积分función第二阶段的一般一致性señal y los dos primeros canales de la primera señal。从不同的数据中获取不同的数据función一致性múltiple。
[Cxy f] = mscohere (oupt inpt(:,(1 2)),汉明(100),[],[],fs);情节(f, Cxy)
研究结果表明:不同的光谱计算方法和不同的光谱计算方法(ángulo)和不同的光谱计算方法(máxima相干)。
Pxy =运行cpsd (oupt inpt(:,(1 2)),汉明(100),[],[],fs);[~，mxx] = max(cxx);为K = 1:2 fprintf('相位滞后%d = %5.2f*pi\n', k,角(Pxy (mxx (k), k)) / pi)结束
相位滞后1 = -0.51*pi相位滞后2 = -0.22*pi
修改ungráfica de coherencia de magnitude cuadrada
generere dos señales sinusoidales必须在1 kHz的频率下持续1秒。在250赫兹以下的频率上加正弦信号。Una de las señales tiene undesfase deπ/3弧度相对于a弧度。Inserte ambas señales在统一变量分析的基础上。
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));
UtilicemscoherePara calculation y representation la coherencia de magnitude cuadrada de las señales。
mscohere(嗯,联合国 ,[],[],[], fs)
修改el título de la gráfica, la etiqueta del ejexY los límites del ejey．
标题(“平方一致性”)包含(“f (Hz)”) ylim([0 .1])
UtilicegcaPara otener:联合国身份证件。Cambie la ubicación de las marcas de verificación。消除“礼仪”y．
Ax = gca;斧子。XTick = 0:250:500; ax.YTick = 0:0.25:1; ax.YLabel.String = [];
我很高兴孩子们Del identiidador para cambiar el color由el ancho de la línea代表。
ln = ax.Children;ln。Color = [0.8 0 0]; ln.LineWidth = 1.5;
tamambisaman puede utilitar集y得到Para modiar las proedades de la línea。
集(get (gca),“孩子”)，“颜色”，[0 0.4 0]，“线型”，“——”，“线宽”, 1)
争论的焦点
论据
Más acerca de
Referencias
[1] Gómez González, A.， J. Rodríguez, X. Sagartzazu, A. Schumacher, I. Isasa。非平稳信号下噪声和振动传输路径的时域多相干分析。2010国际噪声与振动工程学术会议论文集，ISMA2010-USD2010．3927 - 3941页。
好，史蒂文。现代光谱估计。Englewood Cliffs, NJ: Prentice-Hall, 1988。
bbb Rabiner, Lawrence R.和Bernard Gold。数字信号处理理论与应用“，”Englewood Cliffs, NJ: Prentice-Hall, 1975。
斯托伊卡、彼得和伦道夫·摩西。信号的频谱分析。上马鞍河，新泽西州:普伦蒂斯厅，2005。
[5] Welch, Peter D.“快速傅立叶变换在功率谱估计中的应用:一种基于时间平均的方法，在短的，修改的周期图上。”IEEE®音频与电声学学报。AU-15卷，1967年，第70-73页。
Capacidades ampliadas
历史版本
R2006a的引入
