function varargout = rfna_plot(varargin)% RFNA_PLOT Application M-file for rfna_plot.fig%    FIG = RFNA_PLOT launch rfna_plot GUI.%    RFNA_PLOT('callback_name', ...) invoke the named callback.% Last Modified by GUIDE v2.0 07-May-2002 18:01:31if nargin == 0  % LAUNCH GUI	fig = openfig(mfilename);  % no 'reuse'    	% set(fig,'Color',get(0,'defaultUicontrolBackgroundColor'));    handles = guihandles(fig);   	% Generate a structure of handles to pass to callbacks, and store it.     handles.Zo=50;    handles.smith = smith_chart('draw',...                                [5 10 20 30 50 75 100 150 200 300 500]/handles.Zo,...                                [5 10 20 30 50 75 100 150 200 300 500]/handles.Zo,...                                handles.Zo,...                                get(handles.axes1,'position'));  % place on top of existing axis                                    handles.polar = line('color',[0 1 0],'parent',handles.smith,'visible','off','xdata',[],'ydata',[]);                                 set(handles.axes1,'color','none'); % this eventually ends up on top and must be transparent to show                                        % smith chart coordinates below                                  vis='off';    set(handles.smith,'visible',vis);    set(get(handles.smith,'children'),'visible',vis);    handles.line    = line('color',[0 1 0],'parent',handles.axes1,'linewidth',2);    handles.overlay = line('color',[0 1 1],'parent',handles.axes1,'visible','off');            CLWc   = 20;   % Cursor Label width    CROWc  = 70;   % Cursor readout width    MKWc   = 38;   % Mark button width     CURDXc = 8;    HTXTc  = 20;           y  = 40;       % position in pixels    x  = 20;                         positions = [x+0.2*CURDXc+MKWc                             y          CLWc   HTXTc; % x label                       x+0.2*CURDXc+MKWc                             y-HTXTc    CLWc   HTXTc; % y label                       x+CLWc+0.2*CURDXc+MKWc                        y          CROWc  HTXTc; % x readout                       x+CLWc+CROWc+CURDXc+0.2*CURDXc+MKWc           y          CROWc  HTXTc; % x expand                       x+CLWc+0.2*CURDXc+MKWc                        y-HTXTc    CROWc  HTXTc; % y readout                       x+CLWc+CROWc+CURDXc+0.2*CURDXc+MKWc           y-HTXTc    CROWc  HTXTc; % y expand                       x                                             y          HTXTc  HTXTc; % peak button                       x+HTXTc                                       y          HTXTc  HTXTc; % valley button                       x+0.1*CURDXc                                  y-HTXTc    MKWc   HTXTc];% mark button                       curcol=[0 0  1];                                 colors    =  [1 1 1;   % label color                    1 1 0;   % readout color (zeros means track line color)                    0 0 1;   % expansion box color                    1 1 0;   % mark line color                    1 0 0;   % cursor color line 1                    0 0 1];  %                         format_s  = ['%8.5f';'%8.5f'];  % formats for x and y related readouts      font_size = 12;      handles.cursor_id = add_cur(handles.axes1,'init',positions,colors,['X:';'Y:'],...                                             ['x';'+'],font_size,format_s,'on');                     add_cur(handles.cursor_id,'set','axis_cb','rfna_plot(''axis_synch'',gcbo,[],guidata(gcbo))');                    handles	  guidata(fig, handles);	if nargout > 0		varargout{1} = fig;	endelseif ischar(varargin{1}) % INVOKE NAMED SUBFUNCTION OR CALLBACK	try		if (nargout)			[varargout{1:nargout}] = feval(varargin{:}); % FEVAL switchyard		else			feval(varargin{:}); % FEVAL switchyard		end	catch		disp(lasterr);	endendfunction varargout = plot_data(h, data, handles, varargin)         handles.data = data;         switch data.name            case {'s11','s22'}                % create a list for the display popup                s={'Smith','Zin','vswr','return loss'}            case {'s21','s12'}                s={'dB mag','phase','polar'}                         end;            set(handles.display,'string',s,'value',1);         guidata(handles.rfna_plot, handles);         display_Callback(h,[],handles,varargin);  function varargout = axis_synch(h,data, handles, varargin)            set(handles.smith,'xlim',get(handles.axes1,'xlim'),'ylim',get(handles.axes1,'ylim'));                         function  varargout = smith_cursor(h,data, handles, varargin)         % cursor for smith chart aux readout         [xy,index] = add_cur(handles.cursor_id,'get','position');         [number,h_line]=add_cur(handles.cursor_id,'get','active_line');         data = get(h_line,'userdata');  % need for freq info          x=xy(1); y=xy(2);         z  = handles.Zo*(1+(x+j*y))/(1-(x+j*y));         rs = real(z);         xs = imag(z);         index=min(length(data.freq),index);         f  = data.freq(index);         set(handles.aux_ro,'string',sprintf(' F=%8.2f MHz R=%4.1f ohms  X=%4.1f ohms',f,rs,xs));         mag=abs(x+j*y);         set(handles.aux_r2,'string',sprintf(' Radius= %5.3f Angle=%6.1f  VSWR=%4.2f', mag, (180/pi)*atan2(y,x), (1+mag)/(1-mag)));         function  varargout = xfer_cursor(h,data, handles, varargin)         % cursor for xfer function (s12,s21) aux readout         [xy,index] = add_cur(handles.cursor_id,'get','position');         [number,h_line]=add_cur(handles.cursor_id,'get','active_line');         data = get(h_line,'userdata');  % this is in dB and phase format         set(handles.aux_ro,'string',sprintf(' Mag= %8.6f   Phase=%6.1f ',10^(data.mag(index)/20), data.ph(index)));         set(handles.aux_r2,'string','');% --------------------------------------------------------------------function varargout = Overlay_Callback(h, eventdata, handles, varargin)        if strcmp(get(handles.overlay,'visible'),'on')            set(handles.overlay,'visible','off');        else             set(handles.overlay,'visible','on','xdata',...                 get(handles.line,'xdata'),'ydata',get(handles.line,'ydata'),...                 'userdata',get(handles.line,'userdata'));         end;% --------------------------------------------------------------------function varargout = display_Callback(h, eventdata, handles, varargin)                data=handles.data;                switch data.name                case {'s11','s22'}                % convert mag phase to real & imag                    [real_part,imag_part]=dbphase2xy(data.mag,data.ph);                                        set(handles.line,'xdata', real_part, 'ydata',imag_part,'userdata',data);                    set(handles.axes1,'xlim',get(handles.smith,'xlim'),'ylim',get(handles.smith,'ylim'),...                                     'xgrid','off','ygrid','off','xcolor',[0 0 0],'ycolor',[0 0 0]);                    set(get(handles.axes1,'xlabel'),'string','Real');set(get(handles.axes1,'ylabel'),'string','Imag');                    set(handles.smith,'visible','on');                    set(get(handles.smith,'children'),'visible','on');                    add_cur(handles.cursor_id,'set','move_cb','rfna_plot(''smith_cursor'',gcbo,[],guidata(gcbo))');                case  {'s12','s21'}                    set(handles.smith,'visible','off');                    set(get(handles.smith,'children'),'visible','off');                    set(handles.aux_r2,'visible','off');                    set(handles.line,'xdata', data.freq, 'ydata', data.mag,'userdata',data);                    set(handles.axes1,'xlim',[data.freq(1),data.freq(end)],'ylimmode','auto',...                                      'xgrid','on','ygrid','on','xcolor',[1 1 1],'ycolor',[1 1 1]);                     set(get(handles.axes1,'xlabel'),'string','MHz');set(get(handles.axes1,'ylabel'),'string','dB');                    add_cur(handles.cursor_id,'set','move_cb','rfna_plot(''xfer_cursor'',gcbo,[],guidata(gcbo))');                end;                                function out1 = smith_chart(Action,In1,In2,In3,In4) % function out1 = smith_chart(Action,In1,In2,In3,In4)% Action = 'draw'%    In1 = [rvalues];%    In2 = [xvalues];%    In3 = Zo of line (optional)%    In4 = position   (optional)    NSEGSc = 63;    if strcmp(Action,'draw')        if nargin <=3            Zo = 50;        else            Zo = In3;        end;                if nargin <=4           out1=axes('box','off','xgrid','off','ygrid','off','zgrid','off',...                     'color','none','nextplot','add','visible','off');        else           out1=axes('Units','pixels','position',In4,'box','off','xgrid','off','ygrid','off','zgrid','off',...                     'color','none','nextplot','add','visible','off');        end;        set(out1,'xlim',[-1.1 1.1],'ylim',[-1.1,1.1]);        xlabel('');        ylabel('');        theta = linspace(-pi, pi, NSEGSc);                for r=(In1)           z     =((r/(r+1) + 1/(r+1)*exp(j*theta)));           line('xdata',real(z),'ydata',imag(z),'color',[0.5 0.5 0]);           xpos = (r-1)/(r+1);            h=text(xpos, 0, ftoa('%4w',r*Zo));            % h=text(xpos, 0, sprintf('%4.1f',r*Zo));           set(h, 'VerticalAlignment', 'top', 'HorizontalAlignment', 'right','color',[1 1 1]);        end;                % draw unit circle , and imag==0 line         z = exp(j*theta);        line('xdata',real(z),'ydata',imag(z),'color',[1 1 0]);        line('xdata',[-1 1],'ydata',[0 0],'color',[1 1 0]);                v = linspace(0,10,NSEGSc);        r = [ 0 v.^2];                for x = In2            z = r+j*x*ones(1,NSEGSc+1);            g = (z-1)./(z+1);            line('xdata',real(g),'ydata', imag(g), 'color',[0 0.5 0.5]);            line('xdata',real(g),'ydata',-imag(g), 'color',[0 0.5 0.5]);                        g= ((j*x-1)/(j*x+1));            xpos = real(g);            ypos = imag(g);                        s = ftoa('%4w',Zo*x);            h=text([xpos xpos], [ypos -ypos], [' j' s ; '-j' s]);            set(h(1),'VerticalAlignment', 'bottom','color',[1 1 1]);            set(h(2),'VerticalAlignment', 'top','color',[1 1 1]);            if xpos == 0                set(h, 'Horizontalalignment', 'center');            elseif xpos < 0               set(h, 'Horizontalalignment', 'right');            end        end;        rmin = min(In2);        rmax = max(In2);        line('xdata',[(rmin-1)/(rmin+1),(rmax-1)/(rmax+1)],'ydata',[0 0],'color',[1 1 0]);    else        end;        % end;