赌场送彩金35|临海赌场看场子
查看: 723|回复: 1|关注: 0
打印 上一主题 下一主题

[未答复] matlab相干函数(功率谱方法)及窗函数相关问题

[复制链接]

?#29575;?/font>

5 麦片

财富积分


050


1

主题

3

帖子

0

最佳答案
跳转到指定楼层
1#
发表于 2019-4-11 12:44:49 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
本帖最后由 smilemat 于 2019-4-11 12:46 编辑

大神们,初来乍到,最近在做描述函数识别,遇到几个问题:
首先,原理利用的是 ,即利用输入输出信号的互谱、自谱进行计算两者的线性相关性。
我做了以下尝试:
1、使用mscohere函数,得到

图1

类似噪声很大的图像,不知道是否正确,这个相关系数应该如何看?是均值还是最高值的曲线呢?
2、想使用题目中的公式,用matlab带的频谱函数计算,与mscohere对比。
  1. fs=100;
  2. L1=1024;
  3. window=hamming(L1);
  4. noverlap=L1/2;
  5. [pcc_jz,w]=pwelch(K(2,:),window,noverlap);%Pcc
  6. [puu_jz,w]=pwelch(K(4,:),window,noverlap,w);%Puu 纠结的地方开始
  7. [pcu_jz,w]=cpsd(K(2,:),K(4,:),window,noverlap,w);%pcu
  8. rou=power(abs(pcu_jz),2)./pcc_jz./puu_jz;
  9. plot(w/pi,rou)
  10. figure
  11. mscohere(K(2,:),K(4,:),hamming(1024),514,1024)
复制代码
发现在%纠结的地方,这几行起,如果使用上述代码,plot出来的与图1有倍数的差别(看纵坐标的相关性),如图

图2

但是如果代码中
  1. [pcc_jz,w]=pwelch(K(2,:),window,noverlap);%Pcc
  2. [puu_jz,w]=pwelch(K(4,:),window,noverlap,w);%Puu 纠结的地方开始
  3. [pcu_jz,w]=cpsd(K(2,:),K(4,:),window,noverlap,w);%pcu
复制代码
变成
  1. [pcc_jz,w]=pwelch(K(2,:),window,noverlap);%Pcc
  2. puu_jz=pwelch(K(4,:),window,noverlap);%Puu  差别是,这一次不按照之前得出的w,去固定的后面函数的w
  3. pcu_jz=cpsd(K(2,:),K(4,:),window,noverlap);%pcu
复制代码


图3

不知道是为啥?本人是没有学过信号相关的课程,不知道属于什么学科内容呢?或者只是matlab设置参数中的内容?请指教!

3、最基本的问题,就是:

  1. rou=power(abs(pcu_jz),2)./pcc_jz./puu_jz;
复制代码
这一行代码是否符合题目中的那个公式呢?

理论上,附件里的K数据,希望计算出来的是接近1的数,现在算起来只有0.5,不知道是这个代码写的不对,或者就是数据本身相关性不好呢?

(在网上看到说,计算功率谱定义窗函数?#38498;?#20250;有能量损失??)

?#34892;桓行唬?/div>

附件(即代码中的K):

jz.mat (609.4 KB, 下载次数: 2)

?#29575;?/font>

5 麦片

财富积分


050


1

主题

3

帖子

0

最佳答案
2#
 楼主| 发表于 2019-4-18 12:21:24 | 只看该作者
给自己顶一下
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表
赌场送彩金35 南特足球俱乐部缩写 塞维利亚三字代码 江苏时时彩走势图 湖南十分快乐开奖结果 圣斗士星矢极乐净土篇 古罗马建筑 阿曼达·贝尼斯 加拿大卑斯快乐8走势图 瓦伦西亚战责 那不勒斯的黎明鞠婧祎