工频半波积分m语言计算函数

admin

[mw_shl_code=applescript,true]function integ_rsl=halfWaveInt(A)
%实现半波积分矩形算法，输入为任意矩阵，每行为每一通道，每列为每一时刻采样值
%输出有效值
[integ_r,integ_c]=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
[/mw_shl_code]