工频半波积分m语言计算函数
function integ_rsl=halfWaveInt(A)%实现半波积分矩形算法,输入为任意矩阵,每行为每一通道,每列为每一时刻采样值
%输出有效值
=size(A);
integ_rsl=zeros(integ_r,integ_c);
start_flag=1;
point_N=20;
half_pos=19;
for kk=1:integ_r
for ii=start_flag:integ_c
if ii<point_N
for jj=1:ii
integ_rsl(kk,ii)=integ_rsl(kk,ii)+abs(A(kk,jj));
end
else
for jj=ii-half_pos:ii
integ_rsl(kk,ii)=integ_rsl(kk,ii)+abs(A(kk,jj));
end
end
end
end
integ_rsl=integ_rsl*pi*0.025/sqrt(2);
end
页:
[1]