plotpsd文档
这个情节时间序列的功率谱密度函数使用周期图函数。需要Matlab信号处理工具箱。
内容
语法
plotpsd (y, Fs) plotpsd (y, x) plotpsd (…, LineProperty LineValue) plotpsd(…,计算lnx) plotpsd (…, db) plotpsd(…,“λ”)h = plotpsd (…)
描述
plotpsd (y, Fs)图1 d阵列的功率谱y在采样频率Fs使用周期图函数。采样频率Fs必须是一个标量。
plotpsd (y, x)阴谋的功率谱y一个独立的变量引用x。这个语法需要x和y的长度和相等x必须是等距的,monotically增加。对于时间序列,x可能的单位是时间;进行空间分析x可能的单位长度。
plotpsd (…, LineProperty LineValue)指定绘制的线条样式和任何组合LineSpec属性(例如,“颜色”、“r”,“线宽”,2等等)。
plotpsd(…,“计算lnx”)specfies一semilogx情节。
plotpsd(…,“分贝”)情节功率谱在分贝。
“λ”plotpsd (…)标签水平轴作为波长而不是默认频率。注意,这个语法假设λ= 1 / f。
h = plotpsd (…)返回一个句柄h绘制的图形对象。
示例1:火车吹口哨
使用内置的火车吹口哨,情节开始绘图上下文的时间序列:
负载火车t =(0:长度(y) 1) / Fs;情节(t, y)盒子从包含“时间(s)”
如果你有耳机,你可以玩,火车口哨是这样的:
soundsc (y, Fs)
火车上信号的功率谱是这样的:
plotpsd (y, Fs)包含的频率(赫兹)
这使得火车的三个角很清楚!不喜欢默认的蓝色线条吗?画一个胖红线相反:
plotpsd (y, Fs,“颜色”,“红色”,“线宽”4)包含的频率(赫兹)
希望看到作为一个红色线绘制在分贝?
plotpsd (y, Fs,“米”,“数据库”)包含的频率(赫兹)
这里有一个黑线与分贝垂直规模和对数尺度在水平方向上。
plotpsd (y, Fs,“k”,“数据库”,“计算lnx”)轴紧包含的频率(赫兹)
假设你有一些测量y与一些时间向量t你不想去通过努力计算采样率。如果是这种情况,只需输入t而不是Fs:
plotpsd (y, t,”乙:“,“计算lnx”)包含的频率(赫兹)
示例2:海冰范围
加载这个海冰范围的数据带有CDT,并且只使用数据从1988年之后,因为在那之前的东西不是在日常决议:
负载seaice_extent.mat%一些示例数据,CDT%指数自1989年以来的日期:印第安纳州= t > datetime(1989年,1,1);图绘制(t(印第安纳州),extent_N(印第安纳州))ylabel北半球的海冰范围(10 ^ 6公里^ 2)”盒子从%消除了丑,contraining框架轴紧
显然有一定的周期性。我们可以使用机灵做散点图的数据相对于每年的儒略日:
jday机灵= (t);散射(jday(印第安纳州),extent_N(印第安纳州),10日datenum (t(印第安纳州))“填充”)cb = cbdate (“yyyy”);% colorbar格式为日期集(cb,“ydir”,“反向”)%翻转colorbar轴轴紧ylabel北半球的海冰范围(10 ^ 6公里^ 2)”包含年的一天
从上面的两个图,我们可以预期有能量的海冰范围时间序列1 /年频率除了一些长期趋势应表示为低频宽带能量。
画出功率谱密度作为时间的函数(在这里,λ),假设每年365.25条平行样品的采样频率:
图plotpsd (extent_N(印第安纳州),365.25,“λ”甘氨胆酸)组(,“xscale”,“日志”,“yscale”,“日志”)轴紧包含的时间(年)
正如所料,1年的高峰时期脱颖而出。但也有一些其他的小山峰,特别是在6个月和4个月的时间。使用vline告诉他们:
vline ((1:12) / 12,“——”,“颜色”rgb (“灰色”)yl = ylim;% y情节的极限文本(12/12,yl (2),“12”,“绿色”,“高级”)文本(6/12,yl (2),“6”,“绿色”,“高级”)文本(4/12,yl (2),' 4(个月),“绿色”,“高级”)
作者信息
这个函数和支持文档写的金宝app乍得a。格林德克萨斯大学奥斯汀分校的地球物理研究所(UTIG), 2015年10月。适用于2019年的气候数据的工具箱。