sqrc.m

来自「基于Matlab」· M 代码 · 共 12 行

function p = sqrc(t,T,alpha)

J1=find(t==0);
t(J1)=eps;
J2=find(t==T/4/alpha);
t(J2)=t(J2)+1e-5;
J3=find(t==-T/4/alpha);
t(J3)=t(J3)-1e-5;
p=4*alpha/pi/sqrt(T)*...
  (cos((1+alpha)*pi*t/T)+T*sin((1-alpha)*pi*t/T)./(4*alpha*t))./...
  (1-(4*alpha*t/T).^2);