a=str2num(get(handles.a_edit,'string')); Ts=1; N=15; eye_num=6; N_data=1000; dt=Ts/N; t=-3*Ts:dt:3*Ts; d=sign(randn(1,N_data)); dd=sigexpand(d,N); ht=sinc(t/Ts).*(cos(a*pi*t/Ts))./(1-4*a^2*t.^2/Ts^2+eps); st=conv(dd,ht); tt=-3*Ts:dt:(N_data+3)*N*dt-dt; axes(handles.axes1); cla line(tt,st); axis([0 20 -1.2 1.2]); %xlabel('t/Ts'); axes(handles.axes2); cla ss=zeros(1,eye_num*N); ttt=0:dt:eye_num*N*dt-dt; for k=3:50 ss=st(k*N+1:(k+eye_num)*N); drawnow; line(ttt,ss); hold on; end; %xlabel('t/Ts'); |
3 条回复