output1.m

关于高层建筑强风条件下风压系数计算

function varargout = Output1(varargin)
global  flnSaveMRI  intmeth  anal  WD  WS  rf  mem_list  mem_out 
global  MRI
global  mri  anal  recurrence  AnL_A  AnL_B  Qs
% Much of this content is automatically generated by the 'guide' template.
% This file is slightly different in that it operates independently from the others.
% The four plots and peak wind effect values are shown initially using the create functions
% and the variable 'mem_out'. This information is updated within the Output1.fig when changes 
% are made to the member number (Out_11) or the mean recurrence interval (Out_12). Both of 
% these use the callback function for Out_11. The variables affected by a particular callback 
% function are listed just below as global variables.
%
% OUTPUT1 M-file for Output1.fig
%      OUTPUT1, by itself, creates a new OUTPUT1 or raises the existing
%      singleton*.
%
%      H = OUTPUT1 returns the handle to a new OUTPUT1 or the handle to
%      the existing singleton*.
%
%      OUTPUT1('Property','Value',...) creates a new OUTPUT1 using the
%      given property value pairs. Unrecognized properties are passed via
%      varargin to Output1_OpeningFcn.  This calling syntax produces a
%      warning when there is an existing singleton*.
%
%      OUTPUT1('CALLBACK') and OUTPUT1('CALLBACK',hObject,...) call the
%      local function named CALLBACK in OUTPUT1.M with the given input
%      arguments.
%
%      *See GUI Options on GUIDE's Tools menu.  Choose "GUI allows only one
%      instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help Output1

% Last Modified by GUIDE v2.5 03-Nov-2005 08:44:19

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @Output1_OpeningFcn, ...
                   'gui_OutputFcn',  @Output1_OutputFcn, ...
                   'gui_LayoutFcn',  [], ...
                   'gui_Callback',   []);
if nargin & isstr(varargin{1})
    gui_State.gui_Callback = str2func(varargin{1});
end

if nargout
    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
    gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% --- Executes just before Output1 is made visible.
function Output1_OpeningFcn(hObject, eventdata, handles, varargin)
global  flnSaveMRI  intmeth  anal  WD  WS  rf  mem_list  mem_out  
global  MRI
global  mri  anal  recurrence  AnL_A  AnL_B  Qs
handles.output = hObject;
guidata(hObject, handles);
mem_out        = mem_list(1);
load ([flnSaveMRI  '_' num2str(mem_out)])
mem_pop        = num2str(mem_list(1));        % Member no. for popup window (Out_11) 
for i=2:length(mem_list)
    mem_pop = ([ mem_pop '|' num2str(mem_list(i)) ]);
end

if isempty(mri)
    mri = 100;  % Set a default mri (100 years) if mri is not specified
end

set(handles.Out_11,'string', mem_pop);
set(handles.Out_12,'string', mri);
set(handles.Out_13,'string', recurrence);

mmm=find(MRI==mri);
 
set(handles.Out_14, 'string', [num2str(round(100*Mem_Bij_Obs_max_mat(mmm,1) ) / 100) ',' ...
    num2str(round(100*Mem_Bij_Obs_max_mat(mmm,2) ) / 100) ',' num2str(round(100*Mem_Bij_Obs_max_mat(mmm,3) ) / 100)]);
set(handles.Out_15, 'string', [num2str(round(100*Mem_Bij_Obs_min_mat(mmm,1) ) / 100) ',' ...
  num2str(round(100*Mem_Bij_Obs_min_mat(mmm,2) ) / 100) ',' num2str(round(100*Mem_Bij_Obs_min_mat(mmm,3) ) / 100)]);
if Qs==1
    set(handles.Out_16, 'string', [num2str(round(100*Mem_Bij_Est_max_mat(mmm,1) ) / 100) ',' ...
        num2str(round(100*Mem_Bij_Est_max_mat(mmm,2) ) / 100) ',' num2str(round(100*Mem_Bij_Est_max_mat(mmm,3) ) / 100)]); 
    set(handles.Out_17, 'string', [num2str(round(100*Mem_Bij_Est_min_mat(mmm,1) ) / 100) ',' ...
        num2str(round(100*Mem_Bij_Est_min_mat(mmm,2) ) / 100) ',' num2str(round(100*Mem_Bij_Est_min_mat(mmm,3) ) / 100)]); 
else 
    set(handles.Out_16, 'string', ' ')
    set(handles.Out_17, 'string', ' ')
end 

if strcmp(anal, 'A')==1
    set(handles.Out_18, 'string', 'Applied (wind) and Inertial Loads' ) 
else
    set(handles.Out_18, 'string', 'Applied (wind) Loads Only' ) 
end
Out_11_Callback(hObject, eventdata, handles)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% --- Outputs from this function are returned to the command line.
function varargout = Output1_OutputFcn(hObject, eventdata, handles)
varargout{1} = handles.output;

%%%%%%%%%%%%%%%%%  FUNCTIONS  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function Out_11_Callback(hObject, eventdata, handles)
global  flnSaveMRI  intmeth  anal  WD  WS  rf  mem_list  mem_out  
global MRI
global mri  anal  recurrence  AnL_A  AnL_B  Qs
mem_no   = get(handles.Out_11,'value');
mem_out  = mem_list(mem_no);
mri_str  = get(handles.Out_12,'string');
mri      = str2num(mri_str);
load ([flnSaveMRI  '_' num2str(mem_out)])
rankr     = 999/recurrence/mri;
mmm      = find(MRI==mri);

symb_list_obs = {'b--','b','b:'};
symb_list_est = {'r--','r','r:'};
legend_list_obs = {'Start Joint (I) (Obs)', 'Center Node (Obs)', 'End Joint (J) (Obs)'};
legend_list_obs_est = {'Start Joint (I) (Obs)', 'Start Joint (I) (Est)', 'Center Node (Obs)', 'Center Node (Est)',...
    'End Joint (J) (Obs)', 'End Joint (J) (Est)'};
for axpos_index = 1:3
    if Qs == 1
        subplot(handles.Out_plot_11)
        plot(MRI, Mem_Bij_Obs_max_mat(:,axpos_index), symb_list_obs{axpos_index}, MRI,...
            Mem_Bij_Est_max_mat(:,axpos_index), symb_list_est{axpos_index}, 'LineWidth', 1.5);
        grid('on'); axis([ 0  max(MRI)  0  1.4]); set(handles.Out_plot_11, 'YTickLabel',...
            [0:0.2:1.4]); ylabel('Peak Wind Effect');
        hold on;
        if axpos_index == 3
            legend(legend_list_obs_est);
        end

        subplot(handles.Out_plot_12)
        plot(MRI, Mem_Bij_Obs_min_mat(:,axpos_index), symb_list_obs{axpos_index}, MRI,...
            Mem_Bij_Est_min_mat(:,axpos_index), symb_list_est{axpos_index}, 'LineWidth', 1.5);
        grid('on'); axis([ 0  max(MRI)  0  1.4]); set(handles.Out_plot_12, 'YTickLabel',...
            [0:0.2:1.4]); xlabel('MRI (years)'); ylabel('Peak Wind Effect');
        hold on;
        if axpos_index == 3
            legend(legend_list_obs_est);
        end   
    else           
        subplot(handles.Out_plot_11)
        plot(MRI, Mem_Bij_Obs_max_mat(:,axpos_index), symb_list_obs{axpos_index}, 'LineWidth', 1.5);
        grid('on'); axis([ 0  max(MRI)  0  1.4]); set(handles.Out_plot_11, 'YTickLabel',...
            [0:0.2:1.4]); ylabel('Peak Wind Effect');
        hold on;
        if axpos_index == 3
            legend(legend_list_obs);
        end
       
        subplot(handles.Out_plot_12)
        plot(MRI, Mem_Bij_Obs_min_mat(:,axpos_index),symb_list_obs{axpos_index}, 'LineWidth', 1.5);
        grid('on'); axis([ 0  max(MRI)  0  1.4]); set(handles.Out_plot_12, 'YTickLabel',...
            [0:0.2:1.4]); xlabel('MRI (years)'); ylabel('Peak Wind Effect');
        hold on;
        if axpos_inde x == 3
            legend(legend_list_obs);
        end
    end
end

for axpos_index = 1:3
    if Qs == 1
        subplot(handles.Out_plot_11)
        plot(mri, Mem_Bij_Obs_max_mat(mmm,axpos_index), 'bo', mri, Mem_Bij_Est_max_mat(mmm,axpos_index), 'ro');
        grid('on'); axis([ 0  max(MRI)  0  1.4]); set(handles.Out_plot_11, 'YTickLabel',...
            [0:0.2:1.4]); ylabel('Peak Wind Effect');
        hold on;
        if axpos_index < 3
            hold on;
        else
            hold off;
        end

        subplot(handles.Out_plot_12)
        plot(mri, Mem_Bij_Obs_min_mat(mmm,axpos_index), 'bo', mri, Mem_Bij_Est_min_mat(mmm,axpos_index), 'ro');
        grid('on'); axis([ 0  max(MRI)  0  1.4]); set(handles.Out_plot_12, 'YTickLabel',...
            [0:0.2:1.4]); ylabel('Peak Wind Effect');
        if axpos_index < 3
            hold on;
        else
            hold off;
        end
    else
        subplot(handles.Out_plot_11)
        plot(mri,Mem_Bij_Obs_max_mat(mmm,axpos_index), 'bo');
        grid('on'); axis([ 0  max(MRI)  0  1.4]); set(handles.Out_plot_11, 'YTickLabel',...
            [0:0.2:1.4]); ylabel('Peak Wind Effect');
        if axpos_index < 3
            hold on;
        else
            hold off;
        end

        subplot(handles.Out_plot_12)
        plot(mri, Mem_Bij_Obs_min_mat(mmm,axpos_index), 'bo' );
        grid('on'); axis([ 0  max(MRI)  0  1.4]); set(handles.Out_plot_12, 'YTickLabel',...
            [0:0.2:1.4]); xlabel('MRI (years)'); ylabel('Peak Wind Effect');
        if axpos_index < 3
            hold on;
        else
            hold off;
        end
    end
end

set(handles.Out_12, 'string', mri);
set(handles.Out_14, 'string', [num2str(round(100*Mem_Bij_Obs_max_mat(mmm,1) ) / 100) ',' ...
    num2str(round(100*Mem_Bij_Obs_max_mat(mmm,2) ) / 100) ',' num2str(round(100*Mem_Bij_Obs_max_mat(mmm,3) ) / 100)]);
    
set(handles.Out_15, 'string', [num2str(round(100*Mem_Bij_Obs_min_mat(mmm,1) ) / 100) ',' ...
  num2str(round(100*Mem_Bij_Obs_min_mat(mmm,2) ) / 100) ',' num2str(round(100*Mem_Bij_Obs_min_mat(mmm,3) ) / 100)]); 
if Qs==1
    set(handles.Out_16, 'string', [num2str(round(100*Mem_Bij_Est_max_mat(mmm,1) ) / 100) ',' ...
        num2str(round(100*Mem_Bij_Est_max_mat(mmm,2) ) / 100) ',' num2str(round(100*Mem_Bij_Est_max_mat(mmm,3) ) / 100)]);
    set(handles.Out_17, 'string', [num2str(round(100*Mem_Bij_Est_min_mat(mmm,1) ) / 100) ',' ...
        num2str(round(100*Mem_Bij_Est_min_mat(mmm,2) ) / 100) ',' num2str(round(100*Mem_Bij_Est_min_mat(mmm,3) ) / 100)]);
else
    set(handles.Out_16, 'string', ' ')
    set(handles.Out_17, 'string', ' ')
end

function Out_Close_Callback(hObject, eventdata, handles)
close(Output1)
Page_Seven