求助matlab时域分析程序-峰值、均方根、裕度、峭度

x=0:0.1:2*pi

y=sin(x) %信号

ma = max(y) %最大值

mi = min(y) %最小值

me = mean(y) %平均值

pk = ma-mi %峰-峰值

av = mean(abs(y)) %绝对值的平均值(整流平均值)

va = var(y) %方差

st = std(y) %标准差

ku = kurtosis(y) %峭度

rm = rms(y) %均方根

S = rm/av %波形因子

C = pk/rm %峰值因子

Kr = sum(y.^4)/sqrt(sum(y.^2)) %峭度因子

I = pk/av %脉冲因子

xr = mean(sqrt(abs(y)))^2

L = pk/xr %裕度因子

刚好自己要用，整理了下。

如果数据都放在a.txt记事本中。

代码如下：（如果a是向量这按下面的，如果是矩阵，需算两次）

load 'a.txt'

amin=min(a)%%最小值 若为矩阵，则amin=min(min(a))

amax=max(a)%%最大值 同上

aavg=mean(a)%%平均值 同上

avar=var(a)

plot(a)

use function 【skewness】

e.g.

X = randn([5 4])

y = skewness(X)

y =

-0.29330.04820.27350.4641

For vectors, skewness(x) is the skewness of the elements of x. For matrices, skewness(X) is a row vector containing the sample skewness of each column.

For more information of this function you can write 【help skewness】at command line.....

---