矩形波を合成するのはサイン波より何倍も難しいのですが
ガンバってるのですね。



後半の話、
ごめんなさい。
(1 - (-0.37))/2
です。

サンプリング周波数/4 の時に k=0.5 となり、 IIRである>>792の式は
out+=k*(in+old -2*out) -> out = out+(in+old -2*out)/2 = (in+out)/2
で >>791のFIRフィルタとなります。