lphase.m

这是matlab中的一些有用并且常见的原代码

clear,close all
M=input('M= ');n=input('n= ');
L1=input('L1= ');
if isempty(L1) L1=3;end
L2=input('L2= ');
if isempty(L2) L2=3;end
dw=2*pi/16/M;
k=1:16*M; w=k*dw;
H1=L1*exp(-n*j*k*dw);
H2=L2*exp(-(M-n-1)*j*k*dw);
H=H1+H2;
x1(k)=real(H1);y1(k)=imag(H1);
x2(k)=real(H);y2(k)=imag(H);
A=diff(unwrap(angle([H1',H2',H'])))*180/2/pi;
figure(1)
%subplot(2,1,1),plot(w,abs(H));
%subplot(2,1,2),plot(w,angle(H));
%figure(2)
for k=1:1000
    if k==1
        h1=line([0,x1(1),x2(1)],[0,y1(1),y2(1)],'Linewidth',3);
        h2=line([0,x2(1)],[0,y2(1)],'Linewidth',3,'color','r');
        set(h1,'erasemode','xor');
        set(h2,'erasemode','xor');
        axis([-6,6,-6,6]),axis square
        pause
    else
        set(h1,'Xdata',[0,x1(k),x2(k)],'Ydata',[0,y1(k),y2(k)])
        set(h2,'Xdata',[0,x2(k)],'Ydata',[0,y2(k)])
        fft(rand(2^18,1));
        A(k-1,:),pause
        drawnow
    end,
end