PS
# RectifiersDPV.m4
cct_init(svg_font(Times,11bp__))

 tmax = 3.5
 period = tmax/1.75
 eamp = 0.5
 omega = twopi_/period
 ed = eamp/6
 r = 0.95

#                           For DC voltages
 ted = asin(ed/eamp)/twopi_*period
 npts = 40; sum = 0
 for i=0 to npts do {
   t = ted+(period/2-2*ted)*i/npts
   sum += (eamp*cos(omega*t-pi_/2) - ed)*r }

define(`Sine',`Cosine(eamp,omega,`$1',-pi_/2)')
define(`axis',`line thick 0.4 right tmax chop -tmax*0.04 chop 0')

Threephase: [ linethick_(1)
 Ai: axis ; "N`'svg_sub(o)" ljust
 move right 0.3 from Ai.end
 EOrig: Ai.start+(tmax*0.04,0)
 E1: sinusoid(eamp,omega,-pi_/2,-ted,tmax) with .Origin at EOrig
 E2: sinusoid(eamp,omega,-pi_/2-twopi_/3,-ted,tmax) with .Origin at EOrig
 E3: sinusoid(eamp,omega,-pi_/2+twopi_/3,-ted,tmax) with .Origin at EOrig
 "e`'svg_sub(1)" at EOrig+(period/3,eamp*0.8) ljust above
 "e`'svg_sub(2)" at EOrig+(2*period/3,eamp*0.8) ljust above
 "e`'svg_sub(3)" at EOrig+(period,eamp*0.8) ljust above

 Ad: axis with .end at Ai.end+(0,-eamp*2.1); "N`'svg_sub(o)" ljust
 DOrig: Ad.start+(tmax*0.04,0)
 Id: axis with .end at Ad.end+(0,-eamp*1.2); "0" ljust
 IOrig: Id.start+(tmax*0.04,0)

 for i=1 to 11 by 2 do {
   tl = max(0,period/2-i*period/6)
   tr = min(period/2,tmax-(i*period/6-period/2))
   sinusoid(eamp,omega,-pi_/2,tl,tr,thick 0.4 dashed) \
     with .Origin at DOrig+(i*period/6-period/2,0)
   tl = max(period/12,-(i*period/6-period/2+period/12) )
   tr = period/2-period/12
   S1: sinusoid(eamp*r,omega,-pi_/2,tl,tr) \
     with .Origin at DOrig+(i*period/6-period/2,0)+(0,-ed*r); round(at S1.End)
   S2: sinusoid(eamp*r,omega,-pi_/2,tl,tr) \
     with .Origin at IOrig+(i*period/6-period/2,0)+(0,-ed*r); round(at S2.End)
   }
 "e`'svg_sub(1)" at DOrig+(period/3,eamp*0.8) ljust above
 "e`'svg_sub(2)" at DOrig+(2*period/3,eamp*0.8) ljust above
 "e`'svg_sub(3)" at DOrig+(period,eamp*0.8) ljust above

 "e`'svg_sub(L)" at DOrig+(period/4,eamp*0.8) below
 "i`'svg_sub(L)" at IOrig+(period/6,eamp*0.6) rjust above

 vEd = r*eamp*cos(pi_/12)/(5*pi_/12) - r*ed
 VEd: axis with .end at Ad.end+(0,vEd); "E`'svg_sub(D)" ljust
 IEd: axis with .end at Id.end+(0,vEd); "I`'svg_sub(D)" ljust
 ]

Cct3:[
 No: dot; {"N`'svg_sub(o)" at last [].s below rjust }
 { inductor(to Here+(Rect_(dimen_,60)),W); E2: Here; "e`'svg_sub(2)" ljust }
 { inductor(to Here+(Rect_(dimen_,-60)),W); E3: Here; "e`'svg_sub(3)" ljust }
 { inductor(from No-(dimen_,0) to No,W); E1: last line.start
    "e`'svg_sub(1)" at E1 above }
 D3: diode(down_ from E3); b_current(i`'svg_sub(3),rjust)
   llabel(,,D`'svg_sub(3))
 line from E1 to (E1.x,D3.start.y+dimen_/3) then right_ (E3.x+E1.x)/2-E1.x \
   then down_ dimen_/3
 D1: diode(down_); b_current(i`'svg_sub(1),rjust); llabel(,,D`'svg_sub(1))
 line from E2 left E2.x-E1.x+dimen_/3 then down E2.y-E3.y \
   then to 2 between D3.start and D1.start
 D2: diode(down_); b_current(i`'svg_sub(2),rjust); llabel(,,D`'svg_sub(2))
 { line from D2.end to D3.end }
 dot(at D1.end)
 arrowline(down elen_/2); rlabel(,i`'svg_sub(L))
 resistor(right_ elen_); rlabel(+,e`'svg_sub(L),-); llabel(,R)
 line to (Here,No) then to No
 ] with .e at Threephase.w-(0.5,0)
 move left 0.2 from Cct3.w

Fullwave: [
 linethick_(1)
 Ai: axis ; "N`'svg_sub(o)" ljust
 EOrig: Ai.start+(tmax*0.05,0)
 E1: sinusoid(eamp,omega,-pi_/2,-ted,tmax) with .Origin at EOrig
 "e`'svg_sub(1)" at EOrig+(period/3,eamp*0.8) ljust above
 E2: sinusoid(eamp,omega, pi_/2,-ted,tmax) with .Origin at EOrig
 "e`'svg_sub(2)" at EOrig+(period*5/6,eamp*0.8) ljust above

 Ad: axis with .end at Ai.end+(0,-eamp*2.1); "N`'svg_sub(o)" ljust
 DOrig: Ad.start+(tmax*0.05,0)
 Id: axis with .end at Ad.end+(0,-eamp*1.2); "0" ljust
 IOrig: Id.start+(tmax*0.05,0)

 for i=0 to 3 do {
   if i<3 then { tf = period/2 } else { tf = Ad.end.x-(DOrig.x+3*period/2) }
   sinusoid(eamp,omega,-pi_/2,0,tf,thick 0.4 dashed) \
     with .Origin at DOrig+(i*period/2,0)
   line from DOrig+(i*period/2-ted,0) right 2*ted
   line from IOrig+(i*period/2-ted,0) right 2*ted
   if i<3 then {tf = period/2-ted } else {tf = Ad.end.x-(DOrig.x+3*period/2) }
   Ed: sinusoid(eamp*r,omega,-pi_/2,ted,tf) \
     with .Origin at DOrig+(i*period/2,-ed*r)
   Ed: sinusoid(eamp*r,omega,-pi_/2,ted,tf) \
     with .Origin at IOrig+(i*period/2,-ed*r)
   }
 "e`'svg_sub(1)" at DOrig+(period/6,eamp*0.9) rjust above
 "e`'svg_sub(2)" at DOrig+(period/6+period/2,eamp*0.9) rjust above
 "e`'svg_sub(L)" at DOrig+(period/4,eamp*0.8) below
 "i`'svg_sub(L)" at IOrig+(period/6,eamp*0.6) rjust above

 vEd = sum/(period/2/(period/2-2*ted)*npts)
 VEd: axis with .end at Ad.end+(0,vEd); "E`'svg_sub(D)" ljust
 IEd: axis with .end at Id.end+(0,vEd); "I`'svg_sub(D)" ljust
 ] with .sw at Threephase.nw+(0,0.5)

Cct2:[
 T: transformer(right_ elen_,R,8,AW,8)
 line up elen_/2 from T.P1
 gap(right_ elen_)
 { sinusoid(0.08,omega*8,-pi_/2,0,twopi_/(omega*8)) with .c at last [] }
 line to T.P2
 line from T.S1 to (T.P1,T.S1) then down elen_/2
 D1: diode(down_ elen_*3/4); rlabel(,D`'svg_sub(1))
 line from T.S2 to (T.P2,T.S2) then to (T.P2,D1.start)
 D2: diode(down_ elen_*3/4); rlabel(,D`'svg_sub(2))
 corner; line to D1.end; corner
 dot(at 0.5 between D1.end and D2.end)
 arrowline(down elen_/2); rlabel(,i`'svg_sub(L))
 resistor(right_ elen_); rlabel(+,e`'svg_sub(L),-); llabel(,R)
 line to (Here,T.TS)-(0,elen_/2) then left Here.x-T.TS.x then to T.TS
 {"N`'svg_sub(o)" at T.TS-(0,elen_/2) below }
 thinlines_
 move to T.TS-(0,0.1)
 { arrow to (T.P1,Here) chop 0.05 "e`'svg_sub(1)" below }
 { arrow to (T.P2,Here) chop 0.05 "e`'svg_sub(2)" below }
 linethick_(1)
 ] at (Cct3,Fullwave)

Halfwave: [
 linethick = 1
 Ai: axis ; "N`'svg_sub(o)" ljust
 EOrig: Ai.start+(tmax*0.05,0)
 Ei: sinusoid(eamp,omega,-pi_/2,0,tmax) with .Origin at EOrig
 "e`'svg_sub(i)" at EOrig+(period/3,eamp*0.8) ljust above

 Ad: axis with .end at Ai.end+(0,-eamp*1.75); "N`'svg_sub(o)" ljust
 DOrig: Ad.start+(tmax*0.05,0)
 sinusoid(eamp,omega,-pi_/2,0,period/2,thick 0.4 dashed) with .Origin at DOrig
 sinusoid(eamp,omega,-pi_/2,0,period/2,thick 0.4 dashed) \
   with .Origin at DOrig+(period,0)
 "e`'svg_sub(i)" at DOrig+(period/5,eamp*0.9) rjust above

 line from Ad.start to DOrig+(ted,0)
 Ed: sinusoid(eamp*r,omega,-pi_/2,ted,period/2-ted) \
   with .Origin at DOrig+(0,-ed*r)
 "e`'svg_sub(L)" at DOrig+(period/3,eamp*0.6) below
 line from Ed.End right period/2+2*ted
 Ed2: sinusoid(eamp*r,omega,-pi_/2,ted,period/2-ted) \
   with .Origin at DOrig+(period,-ed*r)
 line from Ed2.End to Ad.end

 vEd = sum/(period/(period/2-2*ted)*npts)
 VEd: axis with .end at Ad.end+(0,vEd); "E`'svg_sub(D)" ljust above

 Id: axis with .end at Ad.end+(0,-eamp*1.0); "0" ljust
 IOrig: Id.start+(tmax*0.05,0)
 line from Id.start to IOrig+(ted,0)
 Sd: sinusoid(eamp*r,omega,-pi_/2,ted,period/2-ted) \
   with .Origin at IOrig+(0,-ed*r)
 "i`'svg_sub(L)" at IOrig+(period/3,eamp*0.7) above
 line from Sd.End right period/2+2*ted
 Sd2: sinusoid(eamp*r,omega,-pi_/2,ted,period/2-ted) \
   with .Origin at IOrig+(period,-ed*r)
 line from Sd2.End to Id.end
 IEd: axis with .end at Id.end+(0,vEd); "I`'svg_sub(D)" ljust above
 ] with .sw at Fullwave.nw+(0,0.5)

Cct1:[ diode(down_ elen_); b_current(i`'svg_sub(L),rjust,O,E); rlabel(D)
 resistor(right_ elen_); rlabel(+,e`'svg_sub(L),-); llabel(,R)
 line up elen_
 { "N`'svg_sub(o)" ljust below }
 gap(left_ elen_); rlabel(-,e`'svg_sub(i),+)
 ] at (Cct3,Halfwave)

command "</g>" # end font
PE