surfview.m

模糊逻辑工具箱 模糊逻辑工具箱 模糊逻辑工具箱 模糊逻辑工具箱 模糊逻辑工具箱

        'BackgroundColor',frmColor);

    %------------------------------------
    % The STATUS text window
    labelStr='Status info go here';
    name='status';
    pos=[left bottom right-left btnHt];
    hndl=uicontrol( ...
        'Style','text', ...
        'BackgroundColor',frmColor', ...
        'HorizontalAlignment','left', ...
        'Units','pixel', ...
        'Position',pos, ...
        'Tag',name, ...
        'String',labelStr);

    surfview #update
    rotate3d on

    % Normalize all coordinates
    hndlList=findobj(figNumber,'Units','pixels');
    set(hndlList,'Units','normalized');
 
    % Uncover the figure
    set(figNumber, ...
        'Visible','on');

elseif strcmp(action,'#update');
    %====================================
    figNumber=watchon;
    oldfis=get(figNumber,'UserData');
    fis=oldfis{1};
    xaxisHndl=findobj(figNumber,'Type','uicontrol','Tag','xaxis');
    yaxisHndl=findobj(figNumber,'Type','uicontrol','Tag','yaxis');
    zaxisHndl=findobj(figNumber,'Type','uicontrol','Tag','zaxis');
    xgridHndl=findobj(figNumber,'Type','uicontrol','Tag','xgrid');
    ygridHndl=findobj(figNumber,'Type','uicontrol','Tag','ygrid');
    plotHndl=findobj(figNumber,'Type','uimenu','Tag','plottype');
    evalHndl=findobj(figNumber,'Type','uimenu','Tag','alwayseval');
    refinputHndl=findobj(figNumber,'Type','uicontrol','Tag','refinputedit');
    hndlList=[xaxisHndl yaxisHndl zaxisHndl xgridHndl ygridHndl ...
         plotHndl evalHndl refinputHndl];
    set(hndlList,'Enable','off');
    if isfield(fis, 'input')
     numInputs=length(fis.input);
    else
     numInputs=0;
    end
    if isfield(fis, 'output')
     numOutputs=length(fis.output);
    else
     numOutputs=0;
    end

    if (numInputs<1) | (numOutputs<1),
        statmsg(figNumber, ...
            'Need at least one input and one output to view output surface');
        cla     
        watchoff(figNumber)
        return
    end
    if isfield(fis, 'rule')
      numRules=length(fis.rule);
    else
      numRules=0;
    end

    if numRules<1,
        statmsg(figNumber,'Need at least one rule to view output surface');
        cla
        watchoff(figNumber)
        return
    end

    % The X-AXIS popup menu
    inLabels=getfis(fis,'inLabels');
    xAxisLabelStr=[32*ones(size(inLabels,1),1) inLabels];
    labelStr=xAxisLabelStr;
    set(xaxisHndl,'String',labelStr);

    % The Y-AXIS popup menu
    inLabels(1,:)=[];
    if ~isempty(inLabels),
        labelStr=[32*ones(size(inLabels,1)+1,1) str2mat(inLabels,'- none -')];
    else
        labelStr=' - none -';
    end
    set(yaxisHndl,'String',labelStr);

    % The Z-AXIS popup menu
    outLabels=getfis(fis,'outLabels');
    labelStr=[32*ones(size(outLabels,1),1) outLabels];
    set(zaxisHndl,'String',labelStr);

    % The REFERENCE INPUT edit field
    if numInputs==0,
        refInput=[];
    else
        Range=getfis(fis,'inRange');
        refInput=mean(Range');
    end
    set(refinputHndl,'UserData',refInput);

    set(hndlList,'Enable','on');
    surfview #refinputedit

elseif strcmp(action,'#refinputedit');
    %====================================
    figNumber=watchon;
    refInputHndl=findobj(figNumber,'Type','uicontrol','Tag','refinputedit');
    if isempty(refInputHndl)
       watchoff(figNumber);
       return;
    end
    % The edit field's UserData is where the actual reference input is stored
    oldRefInput=get(refInputHndl,'UserData');

    oldfis=get(figNumber,'UserData');
    fis=oldfis{1};
    numInputs=length(fis.input);
    inLabels=getfis(fis,'inLabels');

    % Determine the x axis
    xAxisHndl=findobj(figNumber,'Type','uicontrol','Tag','xaxis');
    xAxisValue=get(xAxisHndl,'Value');
    xAxisList=get(xAxisHndl,'String');
    xAxisStr=deblank(xAxisList(xAxisValue,:));
    xIndex=findrow(xAxisStr,inLabels);
 
    % Determine the y axis
    yAxisHndl=findobj(figNumber,'Type','uicontrol','Tag','yaxis');
    yAxisValue=get(yAxisHndl,'Value');
    yAxisList=get(yAxisHndl,'String');
    yAxisStr=deblank(yAxisList(yAxisValue,:));
    yIndex=findrow(yAxisStr,inLabels);

    newRefInput=eval(get(refInputHndl,'String'),'oldRefInput');

    if length(newRefInput)~=length(oldRefInput),
        newRefInput=oldRefInput;
    end

    % Fill in the spots that are currently occupied by NaN place-holders
    nanIndex=find(isnan(newRefInput));
    newRefInput(nanIndex)=oldRefInput(nanIndex);

    inRange=getfis(fis,'inRange');

    % Guarantee that new input falls within the accepted range for the variable
    newRefInput=max(inRange(:,1)',newRefInput);
    newRefInput=min(inRange(:,2)',newRefInput);
 
    set(refInputHndl,'UserData',newRefInput);

    % Re-insert the NaN place-holders
    if isempty(yIndex),
        newRefInput(xIndex)=NaN;
    else
        newRefInput([xIndex yIndex])=[NaN NaN];
    end

    if all(isnan(newRefInput)),
        % If it's nothing but NaNs, then don't let the user edit (or see) the field
        set(refInputHndl,'String',' ', ...
            'Enable','off')
    else
        set(refInputHndl,'String',[' ' mat2str(newRefInput,4)], ...
            'Enable','on')
    end

    % Plot right away if always-eval is checked
    evalHndl=findobj(gcf,'Type','uimenu','Tag','alwayseval');
    if strcmp(get(evalHndl,'Checked'),'on'),
        surfview #evaluate
    end
        
    watchoff(figNumber);

elseif strcmp(action,'#xaxis');
    %====================================
    % We need to prevent the x-axis and y-axis from being the same thing
    figNumber=watchon;
    xAxisHndl=get(figNumber,'CurrentObject');
    oldfis=get(figNumber,'UserData');
    fis=oldfis{1};
    numInputs=length(fis.input);
    inLabels=getfis(fis,'inLabels');
    outLabels=getfis(fis,'outLabels');
    % Determine the x axis
    xAxisValue=get(xAxisHndl,'Value');
    xAxisList=get(xAxisHndl,'String');
    xAxisStr=deblank(xAxisList(xAxisValue,:));
    xIndex=findrow(xAxisStr,inLabels);
    yAxisHndl=findobj(figNumber,'Type','uicontrol','Tag','yaxis');
    yAxisValue=get(yAxisHndl,'Value');
    inLabels(xIndex,:)=[];
    if ~isempty(inLabels),
        yAxisList=[32*ones(numInputs,1) str2mat(inLabels,'- none -')];
    else
        yAxisList=' - none -';
    end
    set(yAxisHndl,'String',yAxisList);
    set(yAxisHndl,'Value',yAxisValue);
    yAxisStr=deblank(yAxisList(yAxisValue,:));
    yIndex=findrow(yAxisStr,inLabels);

    % Now update the reference input area
    surfview #refinputedit
    watchoff(figNumber);

elseif strcmp(action,'#yaxis');
    %====================================
    % We need to prevent the x-axis and y-axis from being the same thing
    figNumber=watchon;
    % Update the reference input area
    surfview #refinputedit
    watchoff(figNumber);

elseif strcmp(action,'#grids');
    %====================================
    % Calculate immediately if necessary
    gridHndl=gco;
    % Use try-catch to avoid bogus grids
    gridNum=eval(get(gridHndl,'String'),'15');
    gridNum=abs(round(gridNum(1)));
    gridNum=max(gridNum,3);
    gridNum=min(gridNum,100);
    set(gridHndl,'String',[' ' num2str(gridNum)]);
    evalHndl=findobj(gcf,'Type','uimenu','Tag','alwayseval'); 
    if strcmp(get(evalHndl,'Checked'),'on'),
        surfview #evaluate
    end

elseif strcmp(action,'#evaluate');
    %====================================
    figNumber=get(0,'CurrentFigure');
    axHndl=get(figNumber,'CurrentAxes');
    oldfis=get(figNumber,'UserData');
    fis=oldfis{1};
    statmsg(figNumber,'Calculating new surface');

    % Determine the number of points to be evaluated
    xGridHndl=findobj(figNumber,'Type','uicontrol','Tag','xgrid');
    xGrids=eval(get(xGridHndl,'String'));
    yGridHndl=findobj(figNumber,'Type','uicontrol','Tag','ygrid');
    yGrids=eval(get(yGridHndl,'String'));

    inLabels=getfis(fis,'inLabels');
    outLabels=getfis(fis,'outLabels');

    % Determine the x axis
    xAxisHndl=findobj(figNumber,'Type','uicontrol','Tag','xaxis');
    xAxisValue=get(xAxisHndl,'Value');
    xAxisList=get(xAxisHndl,'String');
    xAxisStr=deblank(xAxisList(xAxisValue,:));
    xIndex=findrow(xAxisStr,inLabels);
 
    % Determine the y axis
    yAxisHndl=findobj(figNumber,'Type','uicontrol','Tag','yaxis');
    yAxisValue=get(yAxisHndl,'Value');
    yAxisList=get(yAxisHndl,'String');
    yAxisStr=deblank(yAxisList(yAxisValue,:));
    yIndex=findrow(yAxisStr,inLabels);

    zAxisHndl=findobj(figNumber,'Type','uicontrol','Tag','zaxis');
    zIndex=get(zAxisHndl,'Value');
    Range=getfis(fis,'inRange');

    % Retrieve the reference input from the edit field's UserData
    refInputHndl=findobj(figNumber,'Type','uicontrol','Tag','refinputedit');
    refInput=get(refInputHndl,'UserData');

    [x,y,z]=gensurf(fis,[xIndex yIndex],zIndex,[xGrids yGrids],refInput);

    % Put the results in the appropriate hiding places
    plotMenuHndl=findobj(figNumber,'Type','uimenu','Tag','plottype');
    cmapMenuHndl=findobj(figNumber,'Type','uimenu','Tag','colormap');
    if isempty(y),
%       set(axHndl,'NextPlot','replace');
%       set(figNumber,'NextPlot','replace');
        plot(x,z,'Color','blue','LineWidth',3);
        xlabel(deblank(inLabels(xIndex,:)));
        ylabel(deblank(outLabels(zIndex,:)));
        set(axHndl,'Box','on')
        set(plotMenuHndl,'Enable','off');
        set(cmapMenuHndl,'Enable','off');
        rotate3d off
        set(figNumber,'HandleVisibility','callback');
        set(axHndl,'HandleVisibility','callback');
    else
        xlabel(deblank(inLabels(xIndex,:)));
        ylabel(deblank(inLabels(yIndex,:)));
        zlabel(deblank(outLabels(zIndex,:)));
        set(axHndl,'UserData',[x y z],'Box','on')
        set(plotMenuHndl,'Enable','on');
        set(cmapMenuHndl,'Enable','on');
        rotate3d on
        surfview #plot
    end
    statmsg(figNumber,'Ready');

elseif strcmp(action,'#plot');
    figNumber=watchon;
    plotMenuHndl=findobj(figNumber,'Type','uimenu','Tag','plottype');
    plotTypeHndl=findobj(plotMenuHndl,'Checked','on');
    plotType=deblank(get(plotTypeHndl,'Label'));
    axHndl=get(figNumber,'CurrentAxes');

    % Unpack the hidden data
    xyz=get(axHndl,'UserData');
    viewData=get(axHndl,'View');
    numCols=size(xyz,2)/3;
    x=xyz(:,1:numCols);
    y=xyz(:,numCols+(1:numCols));
    z=xyz(:,2*numCols+(1:numCols));
    xStr=get(get(axHndl,'XLabel'),'String');
    yStr=get(get(axHndl,'YLabel'),'String');
    zStr=get(get(axHndl,'ZLabel'),'String');

    cla

    rotate3d on
    colorMap=get(figNumber,'Colormap');
    if all(viewData==[0 90]),
        viewData=[-37.5 30];
    end

    % Generate the plot
    if strcmp(plotType,'Surface'),
        surf(x,y,z);
    elseif strcmp(plotType,'Lit Surface'),
        surfl(x,y,z);
    elseif strcmp(plotType,'Mesh'),
        meshHndl=mesh(x,y,z);
        set(meshHndl,'LineWidth',2);
    elseif strcmp(plotType,'X Mesh'),
        meshHndl=mesh(x,y,z);
        set(meshHndl,'MeshStyle','row','FaceColor','none','LineWidth',2);
    elseif strcmp(plotType,'Y Mesh'),
        meshHndl=mesh(x,y,z);
        set(meshHndl,'MeshStyle','column','FaceColor','none','LineWidth',2);
    elseif strcmp(plotType,'Contour'),
        numContours=20;
        [c,h]=contour3(x,y,z,numContours);
        set(h,'LineWidth',3);
    elseif strcmp(plotType,'Pseudo-Color'),
        pcolor(x,y,z);
        viewData=[0 90];
        rotate3d off
    elseif strcmp(plotType,'Quiver'),
        [px,py]=gradient(z);
        quiver(x,y,px,py,2,'b')
        viewData=[0 90];
        rotate3d off
    end

    set(axHndl,'UserData',xyz,'Box','on')
    xlabel(xStr);
    ylabel(yStr);
    zlabel(zStr);
    xMin=min(min(x)); xMax=max(max(x));
    yMin=min(min(y)); yMax=max(max(y));
    zMin=min(min(z)); zMax=max(max(z));
    if zMin==zMax, zMin=-inf; zMax=inf; end;
    axis([xMin xMax yMin yMax zMin zMax])
    set(axHndl,'View',viewData)
    watchoff(figNumber);

elseif strcmp(action,'#colormap');
    figNumber=watchon;
    menuHndl=gcbo;
    newMapName=deblank(get(menuHndl,'Label'));
    if strcmp(newMapName,'Default'),
        newMap=jet(150);
        newMap=newMap(33:97,:);
    elseif strcmp(newMapName,'Blue'),
        newMap=fliplr(pink(120));
        newMap=newMap(43:107,:);
    elseif strcmp(newMapName,'Hot'),
        newMap=hot(128);
        newMap=newMap(43:107,:);
    elseif strcmp(newMapName,'HSV'),
        newMap=hsv;
    end

    colormap(newMap)
    watchoff(figNumber)

elseif strcmp(action,'#shading');
    %====================================
    figNumber=watchon;
    menuHndl=gcbo;
    shadingType=deblank(get(menuHndl,'Label'));
    surfHndl=findobj(figNumber,'Type','surface');
    if strcmp(shadingType,'Faceted'),
        set(surfHndl,'FaceColor','flat');
        set(surfHndl,'EdgeColor',mean(get(figNumber,'Colormap')));
    elseif strcmp(shadingType,'Flat'),
        set(surfHndl,'FaceColor','flat');
        set(surfHndl,'EdgeColor','none');
    elseif strcmp(shadingType,'Smooth'),
        set(surfHndl,'FaceColor','interp');
        set(surfHndl,'EdgeColor',mean(get(figNumber,'Colormap')));
    end

    watchoff(figNumber)

elseif strcmp(action,'#evaltoggle');
    %====================================
    figNumber=watchon;
    evalMenuHndl=gcbo;
    evalBtnHndl=findobj(figNumber,'Type','uicontrol','Tag','evaluate');
    if strcmp(get(evalMenuHndl,'Checked'),'on'),
        set(evalMenuHndl,'Checked','off');
        set(evalBtnHndl,'Enable','on');
    else
        set(evalMenuHndl,'Checked','on');
        set(evalBtnHndl,'Enable','off');
    end
    watchoff(figNumber)

elseif strcmp(action,'#plotselect');
    %====================================
    figNumber=watchon;
    plotSelectHndl=gcbo;
    plotMenuHndl=get(plotSelectHndl,'Parent');
    plotUnselectHndl=findobj(plotMenuHndl,'Checked','on');
    set(plotUnselectHndl,'Checked','off');
    set(plotSelectHndl,'Checked','on');
    surfview #plot
    watchoff(figNumber)

elseif strcmp(action,'#help');
    %====================================
    figNumber=watchon;
    helpwin(mfilename);
    watchoff(figNumber)
                                               
end;