clear?all
current=[10087500-50-87-100-87-500508710087500-50-198-243-148...
????-57-25571982431485725-57-198-243-148-57-2557198];
n=length（current）;
%濾波器，去除三次諧波和直流分量
for?i=3:n
????y（i-2）=（1/（3^0.5））*（current（i）-current（i-2））;
end
%兩點(diǎn)乘積算法
for?i=1:n
????I1（i）=（0.5*（y（i）^2+y（i+3）^2））^0.5;
????if?