多个汤普森多窗口平均功率谱密度估计

3视图(30天)

显示旧的评论

菲尔。 2014年8月2日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/144481-average-of-multiple-thompson-multitaper-power-spectral-density-estimates

编辑: 图像分析 2014年8月3日

答:接受图像分析

我正在做信号处理由摩擦产生的声音。我设法画出令牌在相同的情节,但是我也想阴谋的平均光谱,根据其他PSDs 。当然,我想如果我也能画出不同颜色的光谱平均!任何帮助将不胜感激。

创建当前光谱我有使用功能:

                         Hs = spectrum.mtm;
                        
                         fs = 44100%这是录音的质量信号;
                        
                         x = wavread (“文件名”);
                        
                         y = wavread (“otherfilename”);
                        
                         z = wavread (“lastfilename”);
                        
                         x1 = psd (Hs x,“fs”fs);
                        
                         y2 = psd (Hs, y,“fs”fs);
                        
                         z1 = psd (Hs, z,“fs”fs);
                        
                         持有在;
                        
                         情节(x1);
                        
                         情节(y1);
                        
                         情节(z1);

我想阴谋平均频谱上的光谱 (黑色,如果可能的话)!再次感谢您的任何帮助!

我已经附加在. zip . wav文件文件夹。

菲尔。

1评论
显示隐藏 None

图像分析 2014年8月2日

//www.tianjin-qmedu.com/matlabcentral/answers/13205-tutorial-how-to-format-your-question-with-markup

登录置评。

在回答这个问题。

接受的答案

图像分析 2014年8月2日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/144481-average-of-multiple-thompson-multitaper-power-spectral-density-estimates answer_147212

所以你就不能取平均值和情节吗?我的意思是

                             meanSpectrum = (x1 + y1 + z1) / 3;
                            
                             情节(meanSpectrum“k -”,“线宽”3);
                            
                             传奇(x1的,“日元”,“z1”,“meanSpectrum”);

我似乎真的微不足道,或者遗漏了什么东西? ? ?

16条评论
显示隐藏 15年长的评论

图像分析 2014年8月3日

编辑:图像分析 2014年8月3日

谁告诉你的平均频率的PSD ?不是我。我不会做,因为这是无稽之谈。看看我告诉过你什么。好了,看到这段代码:

                                  Hs = spectrum.mtm;
                                 
                                  fs = 44100%这是录音的质量信号;
                                 
                                  x = wavread (“C: \ \ \ \文件临时用户\ sheppard_1.wav ');
                                 
                                  y = wavread (“C: \ \ \ \文件临时用户\ sheppard_2.wav ');
                                 
                                  z = wavread (“C: \ \ \ \文件临时用户\ sheppard_3.wav ');
                                 
                                  xPsd = psd (Hs x,“fs”fs);
                                 
                                  yPsd = psd (Hs, y,“fs”fs);
                                 
                                  zPsd = psd (Hs, z,“fs”fs);
                                 
                                  持有在;
                                 
                                  情节(xPsd。xPsd.Data频率,的r -,“线宽”2);
                                 
                                  情节(yPsd。yPsd.Data频率,“g -”,“线宽”2);
                                 
                                  情节(zPsd。zPsd.Data频率,“b -”,“线宽”2);
                                 
                                  meanPsd = (xPsd。数据+ yPsd。数据+ zPsd.Data) / 3;
                                 
                                  情节(xPsd。meanPsd频率,“k -”,“线宽”5);
                                 
                                  %全屏放大图。
                                 
                                  集(gcf,“单位”,“归一化”,“OuterPosition”,(0 0 1 1));
                                 
                                  网格在;
                                 
                                  传奇(“x”,“y”,“z”,“的意思是”);

菲尔。 2014年8月3日

example.bmp

显然我没有你告诉我的事情要做,我试过频率为了解决这个问题,但这并没有帮助。这是我,因为我做你建议的时候,我得到一个平坦的线在0轴:

                                  一个= wavread (“sheppard_1”);
                                 
                                  b = wavread (“sheppard_2”);
                                 
                                  c = wavread (“sheppard_3”);
                                 
                                  d = wavread (“sheppard_4”);
                                 
                                  e = wavread (“sheppard_5”);
                                 
                                  f = wavread (“sheppard_6”);
                                 
                                  g = wavread (“sheppard_7”);
                                 
                                  h = wavread (“sheppard_8”);
                                 
                                  我= wavread (“sheppard_9”);
                                 
                                  Hs = spectrum.mtm;
                                 
                                  fs = 44100;
                                 
                                  apsd = psd(海关,“fs”,fs);
                                 
                                  bpsd = psd (h、b“fs”,fs);
                                 
                                  运行cpsd = psd (h, c,“fs”,fs);
                                 
                                  dpsd = psd (Hs d“fs”,fs);
                                 
                                  epsd = psd (Hs, e,“fs”,fs);
                                 
                                  fpsd = psd (f Hs,“fs”,fs);
                                 
                                  gpsd = psd (Hs, g,“fs”,fs);
                                 
                                  hpsd = psd (h, h,“fs”,fs);
                                 
                                  ipsd = psd(海关,我,“fs”,fs);
                                 
                                  meanpsd = (apsd。+ bpsd的数据。数据+运行cpsd。数据+ dpsd。数据+ epsd。数据+ fpsd。数据+ gpsd。数据+ hpsd。数据+ ipsd.data) / 9;
                                 
                                  持有在
                                 
                                  情节(apsd);
                                 
                                  情节(bpsd);
                                 
                                  情节(运行cpsd);
                                 
                                  情节抛物型);
                                 
                                  情节(epsd);
                                 
                                  情节(fpsd);
                                 
                                  情节(gpsd);
                                 
                                  情节(hpsd);
                                 
                                  情节(ipsd);
                                 
                                  情节(apsd。meanpsd频率,“k -”,“线宽”2);
                                 
                                  持有从

我附加的图像结果。

菲尔。 2014年8月3日

我共享的代码用于创建的阴谋。我明白你的情节,但它不是形似谱片像我,这就是我想做的!

好吧,我试着下面,给我我想要的,但是x轴不再是频率(赫兹):

                                  一个= wavread (“sheppard_1”);
                                 
                                  b = wavread (“sheppard_2”);
                                 
                                  c = wavread (“sheppard_3”);
                                 
                                  d = wavread (“sheppard_4”);
                                 
                                  e = wavread (“sheppard_5”);
                                 
                                  f = wavread (“sheppard_6”);
                                 
                                  g = wavread (“sheppard_7”);
                                 
                                  h = wavread (“sheppard_8”);
                                 
                                  我= wavread (“sheppard_9”);
                                 
                                  fs = 44100;
                                 
                                  [pxxa, fa] = pmtm(一个3 (a)长度,fs);
                                 
                                  [pxxb, fb] = pmtm (b, 3,长度(b), fs);
                                 
                                  [pxxc, fc] = pmtm (c 3 (c)长度,fs);
                                 
                                  [pxxd, fd] = pmtm (d、3、长度(d), fs);
                                 
                                  (pxxe、铁)= pmtm (e, 3、长度(e)、fs);
                                 
                                  [pxxf, ff] = pmtm (f, 3,长度(f), fs);
                                 
                                  [pxxg, fg] = pmtm (g, 3、长度(g)、fs);
                                 
                                  [pxxh,跳频]= pmtm (h, 3,长度(h)、fs);
                                 
                                  [pxxi, fi] = pmtm(我3长度(i), fs);
                                 
                                  meanpxx = (pxxa + pxxb + pxxc + pxxd + pxxe + pxxf + pxxg + pxxh + pxxi) / 9;
                                 
                                  持有在
                                 
                                  情节(fa, 10 * log10 (pxxa));
                                 
                                  情节(fb, 10 * log10 (pxxb));
                                 
                                  情节(fc, 10 * log10 (pxxc));
                                 
                                  情节(fd, 10 * log10 (pxxd));
                                 
                                  情节(铁、10 * log10 (pxxe));
                                 
                                  情节(ff, 10 * log10 (pxxf));
                                 
                                  情节(fg, 10 * log10 (pxxg));
                                 
                                  情节(fh, 10 * log10 (pxxh));
                                 
                                  情节(fi 10 * log10 (pxxi));
                                 
                                  情节(fa, 10 * log10 (meanpxx),“k -”,“线宽”2);
                                 
                                  持有从

我附带结果,唯一的问题是规模有点奇怪,它读取作为轴值X 10 ^ 4,这不是我想要的……

图像分析 2014年8月3日

编辑:图像分析 2014年8月3日

这是我的代码。原来你正在策划第一9元素数量的函数,而不是频率,只和PSD是与频率,因此一切都搞砸了。

                                  Hs = spectrum.mtm;
                                 
                                  fs = 44100%这是录音的质量信号;
                                 
                                  w1 = wavread (“C: \ \ \ \文件临时用户\ sheppard_1.wav ');
                                 
                                  w2 = wavread (“C: \ \ \ \文件临时用户\ sheppard_2.wav ');
                                 
                                  w3 = wavread (“C: \ \ \ \文件临时用户\ sheppard_3.wav ');
                                 
                                  w4 = wavread (“C: \ \ \ \文件临时用户\ sheppard_4.wav ');
                                 
                                  w5 = wavread (“C: \ \ \ \文件临时用户\ sheppard_5.wav ');
                                 
                                  将= wavread (“C: \ \ \ \文件临时用户\ sheppard_6.wav ');
                                 
                                  支w7 = wavread (“C: \ \ \ \文件临时用户\ sheppard_7.wav ');
                                 
                                  w8 = wavread (“C: \ \ \ \文件临时用户\ sheppard_8.wav ');
                                 
                                  w9 = wavread (“C: \ \ \ \文件临时用户\ sheppard_9.wav ');
                                 
                                  w1 psd1 = psd(海关,“fs”fs);
                                 
                                  psd2 = psd (Hs w2,“fs”fs);
                                 
                                  psd3 = psd (Hs, w3,“fs”fs);
                                 
                                  w4 psd4 = psd(海关,“fs”fs);
                                 
                                  psd5 = psd (Hs w5,“fs”fs);
                                 
                                  将psd6 = psd (Hs“fs”fs);
                                 
                                  支w7 psd7 = psd(海关,“fs”fs);
                                 
                                  w8 psd8 = psd(海关,“fs”fs);
                                 
                                  w9 psd9 = psd(海关,“fs”fs);
                                 
                                  持有在;
                                 
                                  情节(psd1。psd1.Data频率,“- - -”,“颜色”,兰德(1、3)“线宽”2);
                                 
                                  情节(psd2。psd2.Data频率,“- - -”,“颜色”,兰德(1、3)“线宽”2);
                                 
                                  情节(psd3。psd3.Data频率,“- - -”,“颜色”,兰德(1、3)“线宽”2);
                                 
                                  情节(psd4。psd4.Data频率,“- - -”,“颜色”,兰德(1、3)“线宽”2);
                                 
                                  情节(psd5。psd5.Data频率,“- - -”,“颜色”,兰德(1、3)“线宽”2);
                                 
                                  情节(psd6。psd6.Data频率,“- - -”,“颜色”,兰德(1、3)“线宽”2);
                                 
                                  情节(psd7。psd7.Data频率,“- - -”,“颜色”,兰德(1、3)“线宽”2);
                                 
                                  情节(psd8。psd8.Data频率,“- - -”,“颜色”,兰德(1、3)“线宽”2);
                                 
                                  情节(psd9。psd9.Data频率,“- - -”,“颜色”,兰德(1、3)“线宽”2);
                                 
                                  meanPsd = (psd1。数据+ psd2。数据+ psd3。数据+ psd4。数据+ psd5。数据+ psd6。数据+ psd7。数据+ psd8。数据+ psd9.Data) / 9;
                                 
                                  情节(psd1。meanPsd频率,“k -”,“线宽”5);
                                 
                                  %全屏放大图。
                                 
                                  集(gcf,“单位”,“归一化”,“OuterPosition”,(0 0 1 1));
                                 
                                  网格在;
                                 
                                  传奇(“x”,“y”,“z”,“的意思是”);

菲尔。 2014年8月3日

基于你在做什么我改变了功能:

                                  [pxxa, fa] = pmtm(一个3 (a)长度,fs);
                                 
                                  [pxxb, fb] = pmtm (b, 3,长度(b), fs);
                                 
                                  [pxxc, fc] = pmtm (c 3 (c)长度,fs);
                                 
                                  [pxxd, fd] = pmtm (d、3、长度(d), fs);
                                 
                                  (pxxe、铁)= pmtm (e, 3、长度(e)、fs);
                                 
                                  [pxxf, ff] = pmtm (f, 3,长度(f), fs);
                                 
                                  [pxxg, fg] = pmtm (g, 3、长度(g)、fs);
                                 
                                  [pxxh,跳频]= pmtm (h, 3,长度(h)、fs);
                                 
                                  [pxxi, fi] = pmtm(我3长度(i), fs);

而不是psd的函数!我不确定为什么它工作,但是它,然后我可以跟随你的最初建议的添加,除以总,这给了我一个情节就像我,这正是我想要的情节。得到我想要的轴的方式,我必须手动更改它。我要试着找出一些代码自动化。

谢谢你花时间帮助我!我肯定没有你的建议不能算出来。

谢谢!

登录置评。

类别

信号处理信号处理工具箱光谱分析谱估计

找到更多的在谱估计在帮助中心和文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

多个汤普森多窗口平均功率谱密度估计

1评论
显示隐藏 None

接受的答案

16条评论
显示隐藏 15年长的评论

更多的答案(0)

另请参阅

类别

标签

社区寻宝

多个汤普森多窗口平均功率谱密度估计

1评论 显示隐藏 None

接受的答案

16条评论 显示隐藏 15年长的评论

更多的答案(0)

另请参阅

类别

标签

社区寻宝

1评论
显示隐藏 None

16条评论
显示隐藏 15年长的评论