untitled.m

这是关于数字处理信号的MATLAB源代码

t=200;
N=4000;
deltat=t/N;
n=1:N;
ta=n*deltat;
zp=N/2*deltat;
x=1./((ta-zp).*(ta-zp)+1);
 y=imag(hilbert(x));
plot(ta-zp,x),hold on;
plot(ta-zp,y);
axis([-3,3,-2,2]);
grid;
z=0.0;
for i=1:N 
  z=z+x(i)*y(i)*deltat;
end;