📄 yj.m
字号:
function yj(action,in1,in2);
%yj 系统特性演示 - 1 : 一阶系统
% Illustrating some basic signal processing concepts, such as
% sampling, aliasing, and windowing.
% possible actions:
% 'start'
% 'down'
% 'move'
% 'up'
% 'redraw'
% 'done'
% 'setR'
% 'setwindow'
% 'showwin'
% 'setC'
set(0,'DefaultFigureColor',[.7 .7 .7]);
if nargin<1,
action='start';
end;
global yj_DAT
%global ADDIT_DAT % this is for WINDOW pop-up button
global fun_str
if strcmp(action,'start'),
%====================================
% Graphics initialization
oldFigNumber = watchon;
figNumber = figure;
set(gcf, ...
'NumberTitle','off', ...
'Name','一阶系统特性', ...
'backingstore','off',...
'Units','normalized');
%====================================
% Information for all buttons
labelColor=192/255*[1 1 1];
top=0.95;
bottom=0.05;
left=0.75;
yInitLabelPos=0.90;
left = 0.78;
labelWid=0.18;
labelHt=0.05;
btnWid = 0.18;
btnHt=0.07;
% Spacing between the label and the button for the same command
btnOffset=0.003;
% Spacing between the button and the next command's label
spacing=0.05;
%====================================
% The CONSOLE frame
frmBorder=0.02;
yPos=0.05-frmBorder;
frmPos=[left-frmBorder yPos btnWid+2*frmBorder 0.9+2*frmBorder];
h=uicontrol( ...
'Style','frame', ...
'Units','normalized', ...
'Position',frmPos, ...
'BackgroundColor',[0.5 0.5 0.5]);
%====================================
% The SIGNAL command popup button
btnNumber=1;
yLabelPos=top-(btnNumber-1)*(btnHt+labelHt+spacing);
labelStr='一阶系统 ';
% Generic label information
labelPos=[left yLabelPos-labelHt labelWid labelHt];
uicontrol( ...
'Style','text', ...
'Units','normalized', ...
'Position',labelPos, ...
'BackgroundColor',labelColor, ...
'HorizontalAlignment','left', ...
'String',labelStr);
btnPos=[left yLabelPos-labelHt-btnHt-btnOffset btnWid btnHt];
popup=uicontrol('Style','Popup','String','一阶低通|一阶高通',...
'Position', btnPos, ...
'BackgroundColor','w', ...
'ForegroundColor','k', ...
'Units','normalized',...
'CallBack',[...
'global yj_DAT;',...
'popup=yj_DAT(11);',...
'val = get(popup,''Value'');',...
'if val<=2',...
' yj(''redraw'');',...
'else',...
' getstr;',...
'end']);
%====================================
% The WINDOW command popup button
winHndl=0;
% Generic label information
% Generic popup button information
%====================================
% The FUNDAMENTAL editable text box
btnNumber=2;
yLabelPos=top-(btnNumber-1)*(btnHt+labelHt+spacing);
labelPos=[left yLabelPos-labelHt-0.10 labelWid+0.02 labelHt];
R_text = uicontrol( ...
'Style','text', ...
'Position', labelPos, ...
'Units','normalized', ...
'BackgroundColor',[0.5 0.5 0.5], ...
'ForegroundColor','k', ...
'HorizontalAlignment','left', ...
'String','回路电阻(R)：');
btnPos=[left+0.02 yLabelPos-labelHt-btnHt-btnOffset-0.10 ...
0.4*btnWid+frmBorder btnHt];
R_field = uicontrol( ...
'Style','edit', ...
'BackgroundColor','w', ...
'ForegroundColor','k', ...
'Position', btnPos, ...
'Units','normalized', ...
'HorizontalAlignment','left', ...
'String','10',...
'CallBack','yj(''setR''); yj(''redraw'')');
%====================================
% The FUNDAMENTAL editable text box
btnNumber=3;
yLabelPos=top-(btnNumber-1)*(btnHt+labelHt+spacing);
labelPos=[left yLabelPos-labelHt-0.15 labelWid labelHt];
C_text = uicontrol( ...
'Style','text', ...
'Position', labelPos, ...
'Units','normalized', ...
'BackgroundColor',[0.5 0.5 0.5], ...
'ForegroundColor','k', ...
'HorizontalAlignment','left', ...
'String','回路电容（C)：');
btnPos=[left+0.02 yLabelPos-labelHt-btnHt-btnOffset-0.15 ...
0.4*btnWid+frmBorder btnHt];
C_field = uicontrol( ...
'Style','edit', ...
'BackgroundColor','w', ...
'ForegroundColor','k', ...
'Position', btnPos, ...
'Units','normalized', ...
'HorizontalAlignment','left', ...
'String','1',...
'CallBack','yj(''setC''); yj(''redraw'')');
% The MoXing button
% val = get(popup,'Value');
%if val==1
% plot([0 2],[2,2]);
% end
%========================================
% The CLOSE button
done_button=uicontrol('Style','Pushbutton', ...
'Position',[left bottom btnWid 2*labelHt], ...
'Units','normalized','Callback',...
'yj(''done'')','String','关闭');
%====================================
% Create intial signal
N=50; % number of samples
% N=21;
C=1;
R=10; % hertz
t=R*C;
w=logspace(-4,3);
M=50; % length of fft
window=ones(N,1); % use window to lower side-lobes in the freq. domain
% (makes peaks wider)
ax_tu=axes('Position',[.12 .05 .6 .2],'XLim',...
[0,10],'YLim',[-5 5]);
ylabel('幅频特性');
%ax_win=axes('Position',[left yLabelPos-labelHt-0.1 labelWid 2*labelHt]);
% time domain
val = get(popup,'Value');
if (val == 1),
num=[0 1];den=[t 1];
[m,p]=bode(num,den,w);
p=p*pi./180;
plot([0 2 2 2],[2 2 1.5 2.5],'k',[2 3],[2.5 2.5],'k',[2 3],[1.5 1.5],'k',[3 3],[2.5 1.5],'k',[3 7],[2,2],'k',[6 6 5.5 6.5],[2 0 0 0],'k',[5.5 6.5],[-1 -1],'k',[6 6 7 0],[-1 -3 -3 -3],'k');
axis([-2 8 -5 5]);
hold off
ylabel('一阶系统模型');
text('position',[2.4 0.4],'string','R');
text('position',[5 -.5],'string','C');
elseif (val == 2), % high pass
num=[t 0];den=[t 1];
[m,p]=bode(num,den,w);
p=p*pi./180;
tu=plot([0 2 2 2],[2 2 1.5 2.5],'k',[2.5 2.5],[1.5 2.5],'k',[2.5 5],[2 2],'k',[5 5],[2 1],'k',[4.5 5.5],[1,1],'k',[4.5 4.5 5.5 5.5],[1 0 1 0],'k',[4.5 5.5],[0 0],'k',[5 5 5 0],[0 -1 -1 -1 -1],'k');
axis([-2 8 -5 5]);
ylabel('一阶系统模型');
hold off
end;
ax_pha=axes('Position',[.12 .36 .6 .22],'XLim',[1e-4 1e4],'YLim',[-2.5 2.5]);
pha_line= semilogx(w,p, '-','markersize',8,'EraseMode','xor');
axis([1e-4 1e4 -2.5 2]);
grid;
ylabel('相频特性');
xlabel('频率 (Hz)');
ax_amp=axes('Position',[.12 .67 .6 .26],'XLim',[1e-4 1e4],'YLim',[-60 10]);
amp_line= semilogx(w,20*log10(m),'-','markersize',8,'EraseMode','xor');
axis([1e-4 1e4 -60 10]);
%ax_ax=axes('Position',[.52 .18 .6 .3]);
% (set to xor mode to prevent re-rendering, that is, for speed)
grid;
ylabel('幅频特性');
% xlabel('频率 (Hz)');
text(0.10, 22.25,'一阶系统的波特图');
%set(amp_line,'ButtonDownFcn','yj(''down'')');
%set(pha_line,'ButtonDownFcn','yj(''down'')');
yj_DAT = [R; C; N; M; C_field; 0; 0; ...
amp_line; pha_line; R_field; popup; -1; gcf; ...
w(:); window(:); 0];
%ADDIT_DAT = winHndl;
watchoff(oldFigNumber);
elseif strcmp(action,'down'),
popup=yj_DAT(11);
val = get(popup,'value');
if val<=2
% assumes that a line was clicked
amp_line=yj_DAT(8);
axes(get(amp_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(amp_line,'XData');
line_y=get(amp_line,'YData');
units_str = get(gca,'units'); % save normalized state
set(gca,'units','pixels'); % distance must be in pixels
pp=get(gca,'pos');
xa=get(gca,'xlim');
ya=get(gca,'ylim');
dist=((line_x-x)*pp(3)/(xa(2)-xa(1))).^2 + ...
((line_y-y)*pp(4)/(ya(2)-ya(1))).^2;
[temp,i]=min(dist);
fixed_x=line_x(i);
fixed_y=line_y(i);
set(amp_line,'LineStyle',':');
yj_DAT(6)=fixed_x;
yj_DAT(7)=fixed_y;
set(gca,'units',units_str );
set(gcf,'WindowButtonMotionFcn', 'yj(''move'')');
set(gcf,'WindowButtonUpFcn', 'yj(''up'')');
% set(gcf,'userdata',u);
end
elseif strcmp(action,'move'),
% u = get(gcf,'userdata');
R=yj_DAT(1);
C=yj_DAT(2);
N=yj_DAT(3);
M=yj_DAT(4);
C_field=yj_DAT(5);
fixed_x=yj_DAT(6);
fixed_y=yj_DAT(7);
amp_line=yj_DAT(8);
pha_line=yj_DAT(9);
R_field=yj_DAT(10);
popup=yj_DAT(11);
w=yj_DAT(14:14+N-1);
window=yj_DAT(14+N:14+N+N-1);
pt=get(gca,'currentpoint');
x=pt(1,1);
y=pt(1,2);
val = get(popup,'Value');
t=R*C;
if (val == 1),
num=[0 1];den=[t 1];
[m,p]=bode(num,den,w);
p=p*pi./180;
elseif (val == 2), % high pass
num=[t 0];den=[t 1];
[m,p]=bode(num,den,w);
p=p*pi./180;
end;
set(amp_line,'YData',m);
set(pha_line,'YData',p);
set(R_field,'String',num2str(R));
set(C_field,'String',num2str(C));
%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');
R=yj_DAT(1);
C=yj_DAT(2);
fixed_x=yj_DAT(6);
fixed_y=yj_DAT(7);
set(yj_DAT(8),'linestyle','-','markersize',8);
yj_DAT(1)=R; % set amplitude and frequency
yj_DAT(2)=C;
% set(gcf,'userdata',u);
yj('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) & yj_DAT(12)~=-1 )
close(yj_DAT(12));
end;
close(gcf);
clear global yj_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');
R=yj_DAT(1);
C=yj_DAT(2);
N=yj_DAT(3);
M=yj_DAT(4);
C_field=yj_DAT(5);
amp_line=yj_DAT(8);
pha_line=yj_DAT(9);
R_field=yj_DAT(10);
popup=yj_DAT(11);
w=yj_DAT(14:14+N-1);
window=yj_DAT(14+N:14+N+N-1);
val = get(popup,'Value');
ax_tu=axes('Position',[.12 .05 .6 .2],'XLim',...
[0,10],'YLim',[-5 5]);
t=R*C;
if (val == 1),
num=[0 1];den=[t 1];
[m,p]=bode(num,den,w);
p=p*pi./180;
plot([0 2 2 2],[2 2 1.5 2.5],'k',[2 3],[2.5 2.5],'k',[2 3],[1.5 1.5],'k',[3 3],[2.5 1.5],'k',[3 7],[2,2],'k',[6 6 5.5 6.5],[2 0 0 0],'k',[5.5 6.5],[-1 -1],'k',[6 6 7 0],[-1 -3 -3 -3],'k');
axis([-2 8 -5 5]);
hold off
ylabel('一阶系统模型');
text('position',[2.4 0.4],'string','R');
text('position',[5 -.5],'string','C');
elseif (val == 2), % high pass
num=[t 0];den=[t 1];
[m,p]=bode(num,den,w);
p=p*pi./180;
plot([0 2 2 2],[2 2 0.5 3.5],'k',[2.3 2.3],[0.5 3.5],'k',[2.3 6],[2 2],'k',[5 5],[2 1.2],'k',[4.9 5.1],[1.2,1.2],'k',[4.9 4.9],[1.2 -1.2], 'k',[ 5.1 5.1],[1.2 -1.2],'k',[4.9 5.1],[-1.2 -1.2],'k',[5 5 6 0],[-1.2 -3 -3 -3],'k');
axis([-2 8 -5 5]);
ylabel('一阶系统模型');
hold off
text('position',[2.4 0.4],'string','C');
text('position',[4.4 -.5],'string','R');
end;
set(amp_line,'YData', 20*log10(m));
% amp_line= semilogx(w,20*log10(m));
%grid
set(pha_line,'YData',p);
set(R_field,'String',num2str(R));
set(C_field,'String',num2str(C));
drawnow;
elseif strcmp(action,'setR'),
x = str2num(get(yj_DAT(10),'String'));
if isempty(x), % handle the non-numeric case
set(yj_DAT(10),'String',num2str(yj_DAT(1)));
else
yj_DAT(1)=x;
end;
elseif strcmp(action,'setC'),
x = str2num(get(yj_DAT(5),'String'));
if isempty(x), % handle the non-numeric case
set(yj_DAT(10),'String',num2str(yj_DAT(2)));
else
yj_DAT(2)=x;
end;
%set(0,'DefaultFigureColor',[.5 .5 .5])
end
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -