abcd2h.m

来自「《Home Networking Basis》（《家庭网络基础》）的配套MATL」· M 代码 · 共 12 行

function [ht,h]=abcd2h(a,b,c,d);
r=100;
sz=length(a);
h=2*r*ones(sz,1)./(r*(c*r+d)+a*r+b);
phs=unwrap(angle(h));
phsjt=round(phs(sz)/pi);
phslt=phs(sz-250);
phs=[phs(1:sz-250)' spline([sz-250 sz],[phslt phsjt*pi],[sz-250+1:sz])]';
hma=abs(h).*exp(j*phs);
h(1:sz)=real(hma(1:sz))+sqrt(-1)*imag(hma(1:sz));
h(sz+1:2*sz-1)=real(hma(sz:-1:2))-sqrt(-1)*imag(hma(sz:-1:2));
ht=real(ifft(h));