mscohere
平方的一致性
语法
描述
cxy= mscohere (x,y)cxy输入信号的x和y。
如果
x和y都是向量,它们必须具有相同的长度。如果一个信号是一个矩阵,另一个是一个向量,向量的长度必须等于矩阵的行数。功能扩展的向量,并返回一个矩阵column-by-column平方相干估计。
如果
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 /样品。指定一个阻带衰减的60分贝。过滤器原始信号。
dx = designfilt (“bandpassiir”,“FilterOrder”,16岁,…“StopbandFrequency1”,0.2,“StopbandFrequency2”,0.4,…“StopbandAttenuation”、60);x =过滤器(dx, r);
创建第二个序列,设计一个16阶滤波器停止之间的归一化频率0.6π和0.8πrad /样品。指定一个通带波纹为0.1 dB。过滤器原始信号。
dy = designfilt (“bandstopiir”,“FilterOrder”,16岁,…“PassbandFrequency1”,0.6,“PassbandFrequency2”,0.8,…“PassbandRipple”,0.1);y =过滤器(dy, r);
估计的平方连贯性x和y。使用512 - sample汉明窗。指定500个样本之间的重叠相邻段和DFT 2048点。
[cxy, fc] = mscohere (x, y,汉明(512),500年,2048年);
情节的相干函数和覆盖过滤器的频率响应。
[qx、f] = freqz (dx);qy = freqz (dy);情节(fc /π,cxy)在情节(f /π,abs(季度)、f /π,abs (qy))从
多个相干和普通的连贯性
生成一个随机的双通道信号,x。生成另一个信号,y通过低通滤波两个渠道和添加在一起。指定一个30阶数字滤波器截止频率为0.3π和设计使用一个矩形窗口。
0.3 h = fir1(30日,rectwin (31);x = randn (16384 2);y =总和(过滤器(h, 1, x), 2);
计算multiple-coherence估计的x和y。窗口信号以1024 -样本损害窗口。指定512个样本之间的重叠相邻段和DFT 1024点。情节的估计。
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)
相干分布式天线系统
生成两个多路信号,每1千赫采样2秒。第一信号输入,包括三个正弦曲线频率为120 Hz, 360 Hz, 480 Hz。第二信号输出,由两个正弦信号的频率120赫兹和360赫兹。一个正弦曲线滞后π/ 2的第一个信号。另一个正弦信号的滞后π/ 4。信号都是嵌入在高斯白噪声。
fs = 1000;f = 120;t = (0:1 / fs: 2 - 1 / fs) ';inpt =罪(2 *π* f * [1 3 4]。* t);inpt = inpt + randn(大小(inpt));oupt =罪(2 *π* f * 3 [1]。* t -(π/ 2π/ 4));oupt = oupt + randn(大小(oupt));
估计之间的相关程度,所有的输入信号和每一个输出通道。利用汉明窗长度100窗口数据。mscohere返回一个相干函数为每一个输出通道。频率的相干函数达到极大值的输入和输出。
[Cxy f] = mscohere (inpt oupt,汉明(100),[],[],fs,“再分配”);为k = 1:尺寸(oupt, 2)次要情节(大小(oupt, 2), 1, k)情节(f, Cxy (:, k))标题([“输出”int2str (k)”,所有输入”])结束
开关的输入和输出信号,计算出多个相干函数。使用相同的汉明窗。没有输入和输出之间的相关性在480赫兹。因此没有第三相关函数峰值。
[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流(的相位滞后% d = % 5.2 f *π\ n”,k,角(Pxy (mxx (k), k)) / pi)结束
相位滞后1 = -0.51 *π相位滞后2 = -0.22 *π
修改平方连贯情节
生成两个正弦信号采样1秒1 kHz。每一个正弦信号的频率为250赫兹。落后于其他阶段的信号之一π/ 3弧度。嵌入单元的信号在高斯白噪声方差。
fs = 1000;f = 250;t = 0:1 / fs: 1 - 1 / f;嗯=罪(2 *π* f * t) +兰德(大小(t));联合国=罪(2 *π* f * t-pi / 3) +兰德(大小(t));
使用mscohere计算和情节的平方相干信号。
mscohere(联合国,[]、[][],fs)
修改图的标题的标签x设在和的极限y设在。
标题(“平方一致性”)包含(“f (Hz)”1.1)ylim ([0])
使用gca获得处理当前轴。改变刻度线的位置。删除的标签y设在。
甘氨胆酸ax =;斧子。XTick = 0:250:500; ax.YTick = 0:0.25:1; ax.YLabel.String = [];
调用孩子们财产处理改变颜色和画线的宽度。
ln = ax.Children;ln。Color = [0.8 0 0]; ln.LineWidth = 1.5;
另外,使用集和得到修改属性。
集(get (gca),“孩子”),“颜色”(0.4 0 0),“线型”,“——”,“线宽”,1)
输入参数
输出参数
更多关于
引用
[1]戈麦斯冈萨雷斯。,J. Rodríguez, X. Sagartzazu, A. Schumacher, and I. Isasa. “Multiple Coherence Method in Time Domain for the Analysis of the Transmission Paths of Noise and Vibrations with Non-Stationary Signals.”《2010国际会议的噪音和振动工程、isma2010 -美元- 2010。3927 - 3941页。
[2]凯,史蒂文。现代谱估计。恩格尔伍德悬崖,台北:普伦蒂斯·霍尔出版社,1988年。
[3]拉宾,劳伦斯·R。,伯纳德·金。数字信号处理的理论和应用。恩格尔伍德悬崖,台北:普伦蒂斯·霍尔出版社,1975年。
[4]斯托伊卡,那和伦道夫摩西的。光谱分析的信号。上台北:Prentice Hall出版社,2005年。
[5]韦尔奇,彼得·d·“使用快速傅里叶变换功率谱的估计:基于时间平均在短方法,改进的周期图。”IEEE®交易音频和电声学。AU-15卷,1967年,页70 - 73。
扩展功能
版本历史
之前介绍过的R2006a
