求助求助。。半周积分算法仿真

求助大佬，①为什么我这个程序在设置谐波和基波不同初相角的时候，偶次谐波的误差会比奇次谐波误差大(c和d是直流分量，默认为0)。。。

②然后用矩形近似、梯形法则计算出来的电压电流误差会不一样。。。
搞了3天了不知道为啥好难受