lxzq.m

信号实验常用的simulink模型和mfile,可直接在matlan下运行。

    time_line=lxzq_DAT(7);
    axes(get(time_line,'parent'));	% set the right axes

    % Obtain coordinates of mouse click location in axes units
    pt=get(gca,'currentpoint');
    x=pt(1,1);
    y=pt(1,2);


% find closest vertex of line to mouse click location (call it fixed_x, fixed_y)

    line_x=get(time_line,'XData');
    line_y=get(time_line,'YData');
    units_str = get(gca,'units');  % save normalized state
    set(gca,'units','pixels');     % distance must be in pixels
    p=get(gca,'pos');
    xa=get(gca,'xlim');
    ya=get(gca,'ylim');
    
    dist=((line_x-x)*p(3)/(xa(2)-xa(1))).^2 + ...
         ((line_y-y)*p(4)/(ya(2)-ya(1))).^2;
    [temp,i]=min(dist);
    fixed_x=line_x(i);
    fixed_y=line_y(i);
    set(time_line,'LineStyle','-');

    lxzq_DAT(5)=fixed_x;
    lxzq_DAT(6)=fixed_y;

    set(gca,'units',units_str );
    set(gcf,'WindowButtonMotionFcn', 'lxzq(''move'')');
    set(gcf,'WindowButtonUpFcn', 'lxzq(''up'')');
    % set(gcf,'userdata',u);
end

elseif strcmp(action,'move'),
    % u = get(gcf,'userdata');
    freq=lxzq_DAT(1);
    amp=lxzq_DAT(2);
    N=lxzq_DAT(3);
    M=lxzq_DAT(4);
    fixed_x=lxzq_DAT(5);
    fixed_y=lxzq_DAT(6);
    time_line=lxzq_DAT(7);
    freq_line=lxzq_DAT(8);
    freq_field=lxzq_DAT(9);
    popup=lxzq_DAT(10);
    t=lxzq_DAT(13:13+N-1);
    window=lxzq_DAT(13+N:13+N+N-1);
    ratio = lxzq_DAT(13+N+N);

    pt=get(gca,'currentpoint');
    x=pt(1,1);
    y=pt(1,2);

    amp1=y/fixed_y*amp;
    if (abs(amp1)>1.0),
       amp1=1.0*sign(amp1);
    end;
    if (abs(amp1)<0.05),
       amp1=0.05*sign(amp1);
    end;
    if (amp1 == 0),
        amp1=0.05;
    end;
    freq1=fixed_x/x*freq;
    if freq1<0, freq1=0.01; elseif freq1>100, freq1=100; end;
    val = get(popup,'Value');
    if (val == 1),
        f=amp1*sin(freq1*t*2*pi);
     elseif (val == 2),   % square wave1
        tt=freq*t*2*pi;
        tmp=rem(tt,2*pi);
        f=amp*(2*rem((tt<0)+(tmp>2*ratio*pi | tmp<-2*ratio*pi)+1,2)-1);
      elseif (val == 3),   % square wave2
        tt=freq1*t*2*pi;
        tmp=rem(tt,2*pi);
        f=amp1*(2*rem((tt<0)+(tmp>2*ratio*pi | tmp<-2*ratio*pi)+1,2)-1);
        f=(f>0).*f;
    elseif (val == 4),   % sawtooth
        tt=freq1*t*2*pi;
        f=amp1*((tt < 0) + rem(tt,2*pi)/2/pi - .5)*2;
      elseif (val == 5),   % self-defined waveform
	f=eval(fun_str);
    end;
    set(time_line,'YData',f);
    F=fft(window.*f,2*M);
    F=F(1:M);
    w=(0:M-1)*pi/M;
    FF=abs(F);
    FF=FF*2;
    FF(1)=FF(1)/2;
    FF(M)=FF(M)/2;
    FF=FF/length(f);
    set(freq_line,'YData',FF);
    set(freq_field,'String',num2str(freq1));

elseif strcmp(action,'up'),
    pt=get(gca,'currentpoint');
    x=pt(1,1);
    y=pt(1,2);

    set(gcf,'WindowButtonMotionFcn','');
    set(gcf,'WindowButtonUpFcn','');

    % u=get(gcf,'userdata');
    freq=lxzq_DAT(1);
    amp=lxzq_DAT(2);
    fixed_x=lxzq_DAT(5);
    fixed_y=lxzq_DAT(6);

    amp1=y/fixed_y*amp;
    if (abs(amp1)>1.0),
       amp1=1.0*sign(amp1);
    end;
    if (abs(amp1)<0.05),
       amp1=0.05*sign(amp1);
    end;
    freq1=fixed_x/x*freq;
    if freq1<0, freq1=0.01; elseif freq1>100, freq1=100; end;

    set(lxzq_DAT(8),'linestyle','-','markersize',8);
    lxzq_DAT(1)=freq1;  % set amplitude and frequency
    lxzq_DAT(2)=amp1;
    % set(gcf,'userdata',u);
    lxzq('redraw');

elseif strcmp(action,'done'),
    % close the figure window that is showing the window fnction:
    % u = get(gcf,'userdata');
    show_window=findobj('Tag','window_window');
    if ( ~isempty(show_window) & lxzq_DAT(11)~=-1 )
    	close(lxzq_DAT(11));
    end;
    close(gcf);
    clear global lxzq_DAT
    clear global ADDIT_DAT
    clear global fun_str
    clear show_window

elseif strcmp(action,'redraw'),
    % recomputes time and frequency waveforms and updates display
    % u = get(gcf,'userdata');
    freq=lxzq_DAT(1);
    amp=lxzq_DAT(2);
    N=lxzq_DAT(3);
    M=lxzq_DAT(4);
    time_line=lxzq_DAT(7);
    freq_line=lxzq_DAT(8);
    freq_field=lxzq_DAT(9);
    popup=lxzq_DAT(10);
    t=lxzq_DAT(13:13+N-1);
    ratio = lxzq_DAT(13+N+N);

    window=lxzq_DAT(13+N:13+N+N-1);
    val = get(popup,'Value');
    if (val == 1),
        f=amp*sin(freq*t*2*pi);
  elseif (val == 2),   % square wave1
        tt=freq*t*2*pi;
        tmp=rem(tt,2*pi);
        f=amp*(2*rem((tt<0)+(tmp>2*ratio*pi | tmp<-2*ratio*pi)+1,2)-1);
    elseif (val == 3),   % square wave2
        tt=freq*t*2*pi;
        tmp=rem(tt,2*pi);
        f=amp*(2*rem((tt<0)+(tmp>2*ratio*pi | tmp<-2*ratio*pi)+1,2)-1);
         f=(f>0).*f;
       
    elseif (val == 4),   % sawtooth
        tt=freq*t*2*pi;
        f=amp*((tt < 0) + rem(tt,2*pi)/2/pi - .5)*2;
      elseif (val == 5),   % self-defined waveform
	f=eval(fun_str);
    end;
    set(time_line,'YData',f);
    F=fft(window.*f,2*M);
    F=F(1:M);
    w=(0:M-1)*pi/M;
    FF=abs(F);
    FF=FF*2;
    FF(1)=FF(1)/2;
    FF(M)=FF(M)/2;
    FF=FF/length(f);
    set(freq_line,'YData',FF);
    set(freq_field,'String',num2str(freq));

    drawnow;

elseif strcmp(action,'setfreq'),
    x = str2num(get(lxzq_DAT(9),'string'));
    if isempty(x),   % handle the non-numeric case
        set(lxzq_DAT(9),'string',num2str(lxzq_DAT(1)));
    else,
	if ( x>=0 & x<=100 ),
	        lxzq_DAT(1)=x;
	else,
		set(lxzq_DAT(9),'string',num2str(lxzq_DAT(1)));
	end;
    end;

elseif strcmp(action,'setzkb2'),
    N=lxzq_DAT(3);
    r = get(zkb2,'value');
    lxzq_DAT(13+N+N) = r;
    set(zkb3,'string',num2str(r));
    lxzq('redraw');

elseif strcmp(action,'setzkb3'),
    N=lxzq_DAT(3);
    r = str2num(get(zkb3,'string'));
    if isempty(r),
       set(zkb3,'string',num2str(lxzq_DAT(13+N+N)));
    else,
       if ( r>=0.05 & r<=0.95 ),
          lxzq_DAT(13+N+N) = r;
          set(zkb2,'value',r);
          lxzq('redraw');
       else,
          set(zkb3,'string',num2str(lxzq_DAT(13+N+N)));
       end;
    end;



elseif strcmp(action,'info'),
    ttlStr = '连续周期信号的频谱分析'; 

    hlpStr1=Str2mat( ...                                              
        '                                             ',...
         ' 我们看到，在演示窗口中，上方的图绘出了一个  ',...
         ' 连续时间周期信号的时域波形，用鼠标单击并拖动',...
         ' 信号波形上的任意一点就可以改变信号的周期与幅',...
         ' 值；左边的频率编辑框给出了信号的频率值，也可',...
         ' 以通过用鼠标选中此编辑框，用键盘输入频率值。',...       
         ' 下方的图给出了该连续周期信号经傅立叶变换后的',...
         ' 频谱密度的图形即信号的频域波形。            ',...                                           
         ' 频域信号是通过连续信号的采样点经快速傅立叶  ',...
         ' 变换算法(FFT)而得到的。采样间隔由傅立叶变换 ',...
         ' 的点数决定。                                ',...       
         '                                             ' ,...
         ' 注意：下图中，频率范围是从0 到100 Hz。处于  ',...
         ' 负频域的 DFT 是正频域 DFT 的镜象，故在此只  ',...
         ' 画出了正频域的 DFT。抽样频率是 200Hz，这意  ',...
         ' 味着“奈奎斯特”频率是100Hz。在高于“奈奎斯 ',...
         ' 特”频率的频率点处的 DFT 与在低频（负频域） ',...
         ' 处的DFT相同。                               ' ,...                                     
         '                                             ');

    hlpStr2=Str2mat( ...                                               
        '                                              ',...
         '用鼠标选中输入信号的图标，单击该图标就给出一个',...
         '下拉式菜单，从中选择信号即可改变输入的信号。  '  ,...
         '自定义信号要求输入一个自选的周期信号。        ',...
         '                                              ' ,... 
         ' 菜单“窗口”让你选择一个窗口函数，在对信号   ',...
         ' 进行DFT 以前，这个窗口将与信号波形相乘——   ',...
         ' 即“开窗截取”。                             '   ,... 
         ' 检查框“显示窗口函数”若被选中，则在显示信   ',...
         ' 号波形和其DFT的同时，也显示出窗口函数的形状  '  ,...
         '                                              '  ,...
         ' 用鼠标单击信息窗口，就可看到给出的信息       '  );


       helpwin(ttlStr, hlpStr1, hlpStr2);

end