功率谱(PS)和功率谱密度(PSD)是否一样、谱估计谱图纵轴数值为什么为负
不一样。
目录 1.功率谱和功率谱密度的区别 2.”Spectral estimation——MATLAB”谱图纵轴单位含义 3.功率谱密度的意义 4.谱图纵轴单位为负的解释
1.功率谱和功率谱密度的区别
2.”Spectral estimation——MATLAB”谱图纵轴单位含义
3.功率谱密度的意义
4.谱图纵轴单位为负的解释 )
1.功率谱和功率谱密度的区别[1][2]
功率谱(PS)和功率谱密度(PSD)是不同的。若能量为E,时间为T,频带为F,则功率谱是表示为E/T;而功率谱密度是表示为E/T/F。所以它们的量纲和单位是不同的,表示了不同的物理量。但又由于常把功率谱当作功率谱密度的简称,所以经常容易混淆。
例如:对同一个信号做谱估计(Welch法),分别设定纵轴输出为功率谱’power’和功率谱密度’psd’, 两者得到的谱图近似,同一频率下的功率谱(PS)和功率谱密度(PSD)数值和单位均不一样。前者为(3.65625,-28.1157),单位dB; 后者为(3.65625,-43.5767),单位dB/Hz。所以“功率谱就是功率谱密度”是错误的说法。
总之,功率谱密度(PSD)的单位是带有Hz的,而功率谱(PS)不带,它表示平常意义上的功率(特定频率下的)。
2.”Spectral estimation——MATLAB”谱图纵轴单位含义[3]
使用的MATLAB代码,即MATLAB官网”Spectral estimation——MATLAB”上的代码,虽然名为“谱估计Spectral estimation”(或翻译为功率谱估计),但是在对相关方程做的解释是“功率谱密度(PSD)估计”,如下图:
使用的5个Corresponding Function的描述:”…spectral density estimate…”,其中pwelch( )可以设定纵轴输出为功率谱(PS)或者功率谱密度(PSD),其他4个只能输出功率谱密度(PSD)。
因为大家常把功率谱当作功率谱密度的简称(实际上两者不一样),或者可能在实际应用应用中两者区分不大,所以在”Spectral estimation——MATLAB”中虽然名为功率谱估计(Spectral estimation), 但是实际对象是功率谱密度(PSD)而不是功率谱(PS)。不难解释谱图的纵坐标单位:dB/Hz,dB/Hz就是功率谱密度(PSD)的单位,即Power/frequency(dB/Hz或dB·s)。
3.功率谱密度的意义[4]
在物理学中,信号通常是波的形式,例如电磁波、随机振动或者声波。当波的频谱密度乘以一个适当的系数后将得到每单位频率波携带的功率,这被称为信号的功率谱密度(power spectral density, PSD)或者谱功率分布(spectral power distribution, SPD)。功率谱密度的单位通常用每赫兹的瓦特数(W/Hz)表示,或者使用波长而不是频率,即每纳米的瓦特数(W/nm)来表示。实验中用dB/Hz来表示。
4.谱图纵轴单位为负的解释
为了更好的展示谱图,所以将真实功率进行转换,如25kHz附近的峰值在未转换时无法观察到,在转换后可以观察到。
pxx为真实功率,经过公式10log10(pxx)的转换后,当pxx数值较小时会显示为负数。
同样的代码未转换:
同样的代码转换后:
计算过程:
对于第一张图坐标点(656.25, 9.85621e-06)和第二张图坐标点(656.25, -50.0629)
10log10(9.85621e-06) = -50.0629,与数值吻合。
MATLAB验证代码:
clc;
clear all;
close all;
% Fs=1000; %采样频率
% n=0:1/Fs:3;
% xn=cos(2*pi*n*200)+randn(size(n));
[y,Fs]=audioread('*.wav')
%去直流、线性趋势项
xn = detrend(y)%Welch法
nfft=4096;
M=4096;
window=hann(M);
noverlap=M*0.75;figure(1);
pwelch(xn,window,noverlap,nfft,Fs,'power');
title('power-Welch法');
grid on
grid minorfigure(2);
pwelch(xn,window,noverlap,nfft,Fs,'psd');
title('psd-Welch法');
grid on
grid minor
参考资料:
[1] 功率谱和功率谱密度 – MATLAB中文论坛 (ilovematlab.cn)
[2] (26条消息) 功率谱: 功率谱密度(PSD)、功率谱_Chipei Kung的博客-CSDN博客_功率谱密度
[3] Spectral estimation - MATLAB spectrum - MathWorks 中国
[4] 功率谱密度_百度百科 (baidu.com)
功率谱(PS)和功率谱密度(PSD)是否一样、谱估计谱图纵轴数值为什么为负相关推荐
- 在线PS(PhotoShop),打开PSD文件,图像处理
在线PS(PhotoShop),打开PSD文件,图像处理 打开并编辑PSD,XCF,Sketch(Photoshop,Gimp和Sketch)或任何其他图像文件. 检查PSD和Sketch文件. 将S ...
- ps脚本解析psd文本字号
ps脚本解析psd中文本字号最简单的方法,通过DOM获取:var size = activeDocument.activeLayer.textItem.size; 但是当该文本被拉伸过后,DOM解析出 ...
- 淘宝天猫电商活动大促PS模板合集,附带预览图、源文件
淘宝天猫电商活动大促PS模板,附带预览图.源文件 素材包含源文件PSD.jpg预览图,数量十分齐全 2G多 希望大家能够喜欢 备份(可不限速下载):分秒帧 影音协作 从此无界
- PSD/PS: Matlab估计PSD的方法
Matlab代码: %FFT 获得功率频谱密度估计 N = length(signalData.sampleVolt); xdft = fft(signalData.sampleVolt); xdft ...
- 脑电信号特征提取常用算法(共空间模式CSP、小波变换DWT、功率谱密度PSD、AR模型)
1 共空间模式CSP 原理:共空间模式(CSP)是一种对两分类任务下的空域滤波特征提取算法,能够从多通道的脑机接口数据里面提取出每一类的空间分布成分.公共空间模式算法的基本原理是利用矩阵的对角化,找到 ...
- 大前端求人不如求己系列工具:如何PSD文件自动标注切图和真机预览
前言 在很多时候,UI会给一部分切图,当你使用一个新的技术,需要的不再是png而是webp或者SVG时候,沟通变得异常复杂,沟通成本急剧增加. 很多时候在想,是否有一些工具,可以帮助我们解决这个尴尬的 ...
- Sketch 和 PS中的设计图如何实现“自动切图”?
切图是很多UI设计师的一项日常工作.平时做完设计图,要将设计稿切成便于制作成页面的图片,并标注好尺寸和间距,交付给前端来完成html+css布局的静态页面,有利于交互,形成良好的视觉感. 但有的认为前 ...
- 图片自动切图java源代码_Sketch 和 PS中的设计图如何实现“自动切图”?
切图是很多UI设计师的一项日常工作.平时做完设计图,要将设计稿切成便于制作成页面的图片,并标注好尺寸和间距,交付给前端来完成html+css布局的静态页面,有利于交互,形成良好的视觉感. 但有的认为前 ...
- 图片插值数据_结合PS用这招来增强ArcGIS插值图出图效果
1 概述 在一些科研领域,经常会遇到使用插值的方式进行处理,并生成最终的插值图.插值图在ArcGIS里面非常容易生成,只要具备了采用点数据,通过ArcToolBox里面的插值的工具可以直接生成最终的插 ...
最新文章
- Mysql 小数点后最多保留两位数
- mysql mysqlhotcopy_MySQL备份之mysqlhotcopy与注意事项
- python装饰器-究竟什么是装饰器?python中的装饰器介绍
- 地址总线与数据传送机制
- 论文笔记:Missing Value Imputation for Multi-view UrbanStatistical Data via Spatial Correlation Learning
- matlab padarray函数详解
- ListView下拉刷新、上拉载入更多之封装改进
- Java编辑环境搭建
- windows2008开机占用多少内存_如何提升电脑开机速度?
- java中实体类中实现Serializable接口的作用
- Python的三种基本数据类型
- 机器视觉算法(数据结构)
- 乌班图Linux程序自动启动,通过update-rc.d来管理Ubuntu系统的自动启动程序
- html化学式编辑器,常用的化学式编辑器有哪些?4款化学式编辑器分享
- redis java jar_Redis在java开发中使用
- 一个计算机系统配置有3台打印机,我们是局域网3站点 3台电脑共享一个打印机 3台电脑相同的系统 打印参数也设置完全一样 但是打印出来的效果却不同 请问什么原因?需要怎么设置?...
- 博客推广技巧:如何通过博客推广并实现营销
- pytorch中维度dim的理解
- IE被劫持的手动解除
- 【记录一次服务器被攻击】-[附带解决方案]