ftprpplt.m

很多matlab的源代码

          m1=(1)./(1+j*2*pi*(f_axs+k)); m2=(1)./(1+j*2*pi*(f_axs-k));
          pha = angle(m1.*m1+m2.*m2);
            set(lin1,'xdata',t,'ydata',sig,...
            'linestyle','-','erasemode','back','color','y');
           set(lin3,'xdata',f_axs,'ydata',pha,...
            'linestyle','-','erasemode','back','color','c');
           drawnow
        end
     end
  end



elseif sig == 5

   if prop == 1 %SCALING
     if disp == 1 %mag
      t=-8:0.02:8;
      sig=sinc(t);sig = sig.*sig;
      f_axs=-1:0.02:1;
      spec = urect(f_axs);
      set(axs1,'xlim',[-8,8],'ylim',[-0.5,1.5]);
      set(axs2,'xlim',[-1,1],'ylim',[0,3]);

      for k=1:0.05:3
          sig = sinc(t/k);mag = k*urect(k*f_axs);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
          set(lin3,'xdata',f_axs,'ydata',mag,...
          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    pause(1)
      for k=3:-0.05:1
          sig = sinc(t/k);mag = k*urect(k*f_axs);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
          set(lin3,'xdata',f_axs,'ydata',mag,...
          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
     elseif disp == 2 %phase
      t=-2:0.01:2;
      sig=sinc(t);
      f_axs=-1:0.02:1;
      spec = urect(f_axs);
      set(axs1,'xlim',[-2,2],'ylim',[-0.5,1.5]);
      set(axs2,'xlim',[-1,1],'ylim',[-3.5,3.5]);

      for k=1:0.05:3
          sig = sinc(t/k);mag = k*urect(k*f_axs);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
%          set(lin3,'xdata',f_axs,'ydata',mag,...
%          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    pause(1)
      for k=3:-0.05:1
          sig = sinc(t/k);mag = k*urect(k*f_axs);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
%          set(lin3,'xdata',f_axs,'ydata',mag,...
%          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    end

    elseif prop == 2 %DELAY
     if disp == 1
      t=-8:0.01:8;
      sig=sinc(t);
      f_axs=-1:0.02:1;
      spec = urect(f_axs);
      set(axs1,'xlim',[-8,8],'ylim',[-0.5,1.5]);
      set(axs2,'xlim',[-1,1],'ylim',[0,1.5]);

      for k=0:-0.05:-1
          sig = sinc(t-k);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
%          set(lin3,'xdata',f_axs,'ydata',k*sinc(k*f_axs),...
%          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    pause(1)
      for k=-1:0.05:1
          sig = sinc(t-k);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
%          set(lin3,'xdata',f_axs,'ydata',k*sinc(k*f_axs),...
%          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    pause(1)
      for k=1:-0.05:0
          sig = sinc(t-k);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
%          set(lin3,'xdata',f_axs,'ydata',k*sinc(k*f_axs),...
%          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
     elseif disp == 2
      t=-5:0.01:5;
      sig=sinc(t);
      f_axs=-1:0.02:1;
      spec = urect(f_axs);
      set(axs1,'xlim',[-5,5],'ylim',[-0.5,1.5]);
      set(axs2,'xlim',[-1,1],'ylim',[-3.5,3.5]);

      for k=0:-0.05:-1
          sig = sinc(t-k);
          set(lin1,'xdata',t,'ydata',sig,...
         'linestyle','-','erasemode','back','color','y');
         set(lin3,'xdata',f_axs,'ydata',-2*pi*f_axs*k,...
         'linestyle','-','erasemode','back','color','c');
         drawnow
      end
    pause(1)
      for k=-1:0.05:1
          sig = sinc(t-k);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
          set(lin3,'xdata',f_axs,'ydata',-2*pi*f_axs*k,...
         'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    pause(1)
      for k=1:-0.05:0
          sig = sinc(t-k);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
          set(lin3,'xdata',f_axs,'ydata',-2*pi*f_axs*k,...
          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    end

    elseif prop == 3 %MODULATION
     if disp == 1 %Mag
        t=-3:0.01:3;
        sig=sinc(t);
        f_axs=-6:0.01:6;
        spec = urect(f_axs);
        set(axs1,'xlim',[-3,3],'ylim',[-1.2,1.2]);
        set(axs2,'xlim',[-6,6],'ylim',[0,1.2]);

        for k=0:0.1:5
          sig = sinc(t).*cos(2*pi*k*t);
          mag=0.5*(urect(f_axs-k)+urect(f_axs+k));
            set(lin1,'xdata',t,'ydata',sig,...
            'linestyle','-','erasemode','back','color','y');
           set(lin3,'xdata',f_axs,'ydata',mag,...
           'linestyle','-','erasemode','back','color','c');
           drawnow
        end
      pause(1)
        for k=5:-0.1:0
          sig = sinc(t).*cos(2*pi*k*t);
          mag=0.5*(urect(f_axs-k)+urect(f_axs+k));
            set(lin1,'xdata',t,'ydata',sig,...
            'linestyle','-','erasemode','back','color','y');
           set(lin3,'xdata',f_axs,'ydata',mag,...
            'linestyle','-','erasemode','back','color','c');
           drawnow
        end
     elseif disp == 2 %phase
        t=-2:0.01:2;
        sig=sinc(t);
        f_axs=-6:0.01:6;
        spec = urect(f_axs);
        set(axs1,'xlim',[-2,2],'ylim',[-1.2,1.2]);
        set(axs2,'xlim',[-3,3],'ylim',[-3.5,3.5]);

        for k=0:0.1:5
            sig = sinc(t).*cos(2*pi*k*t);
            set(lin1,'xdata',t,'ydata',sig,...
            'linestyle','-','erasemode','back','color','y');
%           set(lin3,'xdata',f_axs,'ydata',0.5*(sinc(f_axs-k)+sinc(f_axs+k)),...
%           'linestyle','-','erasemode','back','color','c');
           drawnow
        end
      pause(1)
        for k=5:-0.1:0
            sig = sinc(t).*cos(2*pi*k*t);
            set(lin1,'xdata',t,'ydata',sig,...
            'linestyle','-','erasemode','back','color','y');
%           set(lin3,'xdata',f_axs,'ydata',0.5*(sinc(f_axs-k)+sinc(f_axs+k)),...
%            'linestyle','-','erasemode','back','color','c');
           drawnow
        end
     end
  end




elseif sig == 6

   if prop == 1 %SCALING
     if disp == 1 %mag
      t=-8:0.02:8;
      sig=sinc(t);sig = sig.*sig;
      f_axs=-2:0.02:2;
      spec = tri(f_axs);
      set(axs1,'xlim',[-8,8],'ylim',[0,1.2]);
      set(axs2,'xlim',[-2,2],'ylim',[0,3]);

      for k=1:0.05:3
          sig = sinc(t/k);sig=sig.*sig;mag = k*tri(k*f_axs);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
          set(lin3,'xdata',f_axs,'ydata',mag,...
          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    pause(1)
      for k=3:-0.05:1
          sig = sinc(t/k);sig=sig.*sig;mag = k*tri(k*f_axs);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
          set(lin3,'xdata',f_axs,'ydata',mag,...
          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
     elseif disp == 2 %phase
      t=-8:0.02:8;
      sig=sinc(t);sig=sig.*sig;
      f_axs=-2:0.02:2;
      spec = tri(f_axs);
      set(axs1,'xlim',[-8,8],'ylim',[0,1.2]);
      set(axs2,'xlim',[-2,2],'ylim',[-3.5,3.5]);

      for k=1:0.05:3
          sig = sinc(t/k);sig=sig.*sig;mag = k*tri(k*f_axs);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
%          set(lin3,'xdata',f_axs,'ydata',mag,...
%          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    pause(1)
      for k=3:-0.05:1
          sig = sinc(t/k);sig=sig.*sig;mag = k*tri(k*f_axs);
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
%          set(lin3,'xdata',f_axs,'ydata',mag,...
%          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    end

    elseif prop == 2 %DELAY
     if disp == 1
      t=-3:0.01:3;
      sig=sinc(t);sig = sig.*sig;
      f_axs=-2:0.02:2;
      spec = tri(f_axs);
      set(axs1,'xlim',[-3,3],'ylim',[0,1.2]);
      set(axs2,'xlim',[-2,2],'ylim',[0,1.5]);

      for k=0:-0.05:-1
          sig = sinc(t-k);sig = sig.*sig;
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
%          set(lin3,'xdata',f_axs,'ydata',k*sinc(k*f_axs),...
%          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    pause(1)
      for k=-1:0.05:1
          sig = sinc(t-k);sig = sig.*sig;
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
%          set(lin3,'xdata',f_axs,'ydata',k*sinc(k*f_axs),...
%          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    pause(1)
      for k=1:-0.05:0
          sig = sinc(t-k);sig = sig.*sig;
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
%          set(lin3,'xdata',f_axs,'ydata',k*sinc(k*f_axs),...
%          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
     elseif disp == 2
      t=-5:0.01:5;
      sig=sinc(t);sig=sig.*sig;
      f_axs=-2:0.02:2;
      spec = tri(f_axs);
      set(axs1,'xlim',[-5,5],'ylim',[0,1.2]);
      set(axs2,'xlim',[-2,2],'ylim',[-3.5,3.5]);

      for k=0:-0.05:-1
          sig = sinc(t-k);sig = sig.*sig;
          set(lin1,'xdata',t,'ydata',sig,...
         'linestyle','-','erasemode','back','color','y');
         set(lin3,'xdata',f_axs,'ydata',-2*pi*f_axs*k,...
         'linestyle','-','erasemode','back','color','c');
         drawnow
      end
    pause(1)
      for k=-1:0.05:1
          sig = sinc(t-k);sig = sig.*sig;
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
          set(lin3,'xdata',f_axs,'ydata',-2*pi*f_axs*k,...
         'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    pause(1)
      for k=1:-0.05:0
          sig = sinc(t-k);sig = sig.*sig;
          set(lin1,'xdata',t,'ydata',sig,...
          'linestyle','-','erasemode','back','color','y');
          set(lin3,'xdata',f_axs,'ydata',-2*pi*f_axs*k,...
          'linestyle','-','erasemode','back','color','c');
          drawnow
      end
    end

    elseif prop == 3 %MODULATION
     if disp == 1 %Mag
        t=-3:0.01:3;
        sig=sinc(t);sig = sig.*sig;
        f_axs=-6:0.01:6;
        spec = tri(f_axs);
        set(axs1,'xlim',[-3,3],'ylim',[-1.2,1.2]);
        set(axs2,'xlim',[-6,6],'ylim',[0,1.2]);

        for k=0:0.1:6
          sig = sinc(t);sig=sig.*sig.*cos(2*pi*k*t);
          mag=0.5*(tri(f_axs-k)+tri(f_axs+k));
            set(lin1,'xdata',t,'ydata',sig,...
            'linestyle','-','erasemode','back','color','y');
           set(lin3,'xdata',f_axs,'ydata',mag,...
           'linestyle','-','erasemode','back','color','c');
           drawnow
        end
      pause(1)
        for k=6:-0.1:0
          sig = sinc(t);sig=sig.*sig.*cos(2*pi*k*t);
          mag=0.5*(tri(f_axs-k)+tri(f_axs+k));
            set(lin1,'xdata',t,'ydata',sig,...
            'linestyle','-','erasemode','back','color','y');
           set(lin3,'xdata',f_axs,'ydata',mag,...
            'linestyle','-','erasemode','back','color','c');
           drawnow
        end
     elseif disp == 2 %phase
        t=-2:0.01:2;
        sig=sinc(t);sig=sig.*sig;
        f_axs=-6:0.01:6;
        spec = tri(f_axs);
        set(axs1,'xlim',[-2,2],'ylim',[-1.2,1.2]);
        set(axs2,'xlim',[-6,6],'ylim',[-3.5,3.5]);

        for k=0:0.1:6
            sig = sinc(t);sig=sig.*sig.*cos(2*pi*k*t);
            set(lin1,'xdata',t,'ydata',sig,...
            'linestyle','-','erasemode','back','color','y');
%           set(lin3,'xdata',f_axs,'ydata',0.5*(sinc(f_axs-k)+sinc(f_axs+k)),...
%           'linestyle','-','erasemode','back','color','c');
           drawnow
        end
      pause(1)
        for k=6:-0.1:0
            sig = sinc(t);sig=sig.*sig.*cos(2*pi*k*t);
            set(lin1,'xdata',t,'ydata',sig,...
            'linestyle','-','erasemode','back','color','y');
%           set(lin3,'xdata',f_axs,'ydata',0.5*(sinc(f_axs-k)+sinc(f_axs+k)),...
%            'linestyle','-','erasemode','back','color','c');
           drawnow
        end
     end
  end



end