[matlab]余弦分频电路

function p=cos_divide(h,divide)pt=h;pt=(pt-min(pt))/(max(pt)-min(pt));pt=pt*2-1;pt=acos(pt);p=pt;phase_compensate=0;if length(pt)>1    for n=2:length(p)        if pt(n)<pt(n-1)            p(n)=2*pi-pt(n);        end        p(n)=p(n)+phase_compensate;        if p(n-1)-p(n)>pi            p(n)=p(n)+2*pi;            phase_compensate=phase_compensate+2*pi;        end        if p(n)>2*pi*divide            p(n)=p(n)-2*pi*divide;            phase_compensate=phase_compensate-2*pi*divide;        end    endendp=cos(p/divide);