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:31
if nargin == 0 % LAUNCH GUI
fig = openfig(mfilename); % no 'reuse'
% set(fig,'Color',get(0,'defaultUicontrolBackgroundColor'));
handles = guihandles(fig);
figure(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'),handles.rfna_plot); % 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))');
guidata(fig, handles);
if nargout > 0
varargout{1} = fig;
end
elseif 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);
end
end
function 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 real and imaginary format
set(handles.aux_ro,'string',sprintf(' Mag= %8.6f Phase=%6.1f ',20*log10(sqrt(data.real(index)^2+data.imag(index)^2)),(180/pi)* atan2(data.imag(index),data.real(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'}
set(handles.line,'xdata', data.real, 'ydata',data.imag,'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', 10*log10(data.real.^2+data.imag.^2),'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,In5)
% function out1 = smith_chart(Action,In1,In2,In3,In4)
% Action = 'draw'
% In1 = [rvalues];
% In2 = [xvalues];
% In3 = Zo of line (optional)
% In4 = position (optional)
% In5 = parent figure handle
NSEGSc = 63;
if strcmp(Action,'draw')
if nargin <=3
Zo = 50;
else
Zo = In3;
end;
figure(In5)
if nargin <=4
out1=axes('box','off','xgrid','off','ygrid','off','zgrid','off',...
'color','none','nextplot','add','visible','off','parent',In5);
else
out1=axes('Units','pixels','position',In4,'box','off','xgrid','off','ygrid','off','zgrid','off',...
'color','none','nextplot','add','visible','off', 'parent',In5);
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;
