用matlab实现这些功能。

Question

我有一系列的心电图信号1250x1双，让我们称之为"a“。我需要实现四个功能，代表特征是用来表征信号。能量、四次方、非线性能量与曲线长度

​

​

我准备实施能源与第四动力

for i=1:1250
energy = sum(a.^2,i);

power4th =  sum(a.^4,i);
end

产生两个阵列(能量和power4th)

我怎么能产生另外两个数组？让我们称他们为NonLE和CL。

Answer 1

使用矢量化而不是for循环来解决所需的所有4个公式

% generate some random numbers
a = rand(1000,1);

Energy = sum(a.^2);
Power4 = sum(a.^4);
NLEnergy = sum(-a(3:end).*a(1:end-2) + a(2:end).^2);
CurveLength = sum(a(2:end) - a(1:end-1));

.运算符允许在向量中逐元素操作。

Answer 2

实际上，我认为您可以在不使用for循环的情况下实现公式。您可以使用矩阵乘法特性。试试下面的代码：

len      = 1250;
a        = randi(10, len, 1); % // You didn' t give your vector so I generated random a.. 
Energy   = ones(1, len) * (a.^2);
power4th = ones(1, len) * (a.^4);
NonLE    = ones(1, len - 2) * ( -a(3:end) .* a(1:end-2) ) +  ones(1, len - 1) * ( a(2:end).^2 );
CL       = ones(1, len - 1) * ( a(2:end) - a(1:end-1) );

Answer 3

对于其中的三个，您并不需要一个for循环：

energy = sum(a.^2);
power_4th = sum(a.^4);
curve_length = sum(diff(a));

对于最后一个问题，您可以这样做：

nonLE = 0;
for k = 3 : length(a)
    nonLE = nonLE + a(k - 1)^2 - a(k) * a(k - 2);
end