p_wav.m

来自「心电信号产生程序 仿真人体正常状态下的心电信号」· M 代码 · 共 17 行

function [pwav]=p_wav(x,a_pwav,d_pwav,t_pwav,li)
l=li;
a=a_pwav;
x=x+t_pwav;
b=(2*l)/d_pwav;
n=100;
p1=1/l;
p2=0;
for i = 1:n
    harm1=(((sin((pi/(2*b))*(b-(2*i))))/(b-(2*i))+(sin((pi/(2*b))*(b+(2*i))))/(b+(2*i)))*(2/pi))*cos((i*pi*x)/l);             
    p2=p2+harm1;
end
pwav1=p1+p2;
pwav=a*pwav1;