function varargout = untitled(varargin)
% UNTITLED MATLAB code for untitled.fig
%      UNTITLED, by itself, creates a new UNTITLED or raises the existing
%      singleton*.
%
%      H = UNTITLED returns the handle to a new UNTITLED or the handle to
%      the existing singleton*.
%
%      UNTITLED(\'CALLBACK\',hObject,eventData,handles,...) calls the local
%      function named CALLBACK in UNTITLED.M with the given input arguments.
%
%      UNTITLED(\'Property\',\'Value\',...) creates a new UNTITLED or raises the
%      existing singleton*.  Starting from the left, property value pairs are
%      applied to the GUI before untitled_OpeningFcn gets called.  An
%      unrecognized property name or invalid value makes property application
%      stop.  All inputs are passed to untitled_OpeningFcn via varargin.
%
%      *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 untitled

% Last Modified by GUIDE v2.5 28-Aug-2015 09:45:23

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct(\'gui_Name\',       mfilename, ...
                   \'gui_Singleton\',  gui_Singleton, ...
                   \'gui_OpeningFcn\', @untitled_OpeningFcn, ...
                   \'gui_OutputFcn\',  @untitled_OutputFcn, ...
                   \'gui_LayoutFcn\',  [] , ...
                   \'gui_Callback\',   []);
if nargin && ischar(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 untitled is made visible.
function untitled_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% varargin   command line arguments to untitled (see VARARGIN)

% Choose default command line output for untitled
handles.output = hObject;

% Update handles structure
guidata(hObject, handles);

% This sets up the initial plot - only do when we are invisible
% so window can get raised using untitled.
if strcmp(get(hObject,\'Visible\'),\'off\')
  %  plot(rand(5));
end

% UIWAIT makes untitled wait for user response (see UIRESUME)
% uiwait(handles.figure1);


% --- Outputs from this function are returned to the command line.
function varargout = untitled_OutputFcn(hObject, eventdata, handles)
% varargout  cell array for returning output args (see VARARGOUT);
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure
varargout{1} = handles.output;

% --- Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
axes(handles.axes1);
cla;

popup_sel_index = get(handles.popupmenu1, \'Value\');
switch popup_sel_index
    case 1
        plot(rand(5));
    case 2
        plot(sin(1:0.01:25.99));
    case 3
        bar(1:.5:10);
    case 4
        plot(membrane);
    case 5
        surf(peaks);
end


% --------------------------------------------------------------------
function FileMenu_Callback(hObject, eventdata, handles)
% hObject    handle to FileMenu (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)


% --------------------------------------------------------------------
function OpenMenuItem_Callback(hObject, eventdata, handles)
% hObject    handle to OpenMenuItem (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
file = uigetfile(\'*.fig\');
if ~isequal(file, 0)
    open(file);
end

% --------------------------------------------------------------------
function PrintMenuItem_Callback(hObject, eventdata, handles)
% hObject    handle to PrintMenuItem (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
printdlg(handles.figure1)

% --------------------------------------------------------------------
function CloseMenuItem_Callback(hObject, eventdata, handles)
% hObject    handle to CloseMenuItem (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
selection = questdlg([\'Close \' get(handles.figure1,\'Name\') \'?\'],...
                     [\'Close \' get(handles.figure1,\'Name\') \'...\'],...
                     \'Yes\',\'No\',\'Yes\');
if strcmp(selection,\'No\')
    return;
end

delete(handles.figure1)


% --- Executes on selection change in popupmenu1.
function popupmenu1_Callback(hObject, eventdata, handles)
% hObject    handle to popupmenu1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: contents = get(hObject,\'String\') returns popupmenu1 contents as cell array
%        contents{get(hObject,\'Value\')} returns selected item from popupmenu1


%  Executes during object creation, after setting all properties.
function popupmenu1_CreateFcn(hObject, eventdata, handles)
% hObject    handle to popupmenu1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: popupmenu controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,\'BackgroundColor\'), get(0,\'defaultUicontrolBackgroundColor\'))
     set(hObject,\'BackgroundColor\',\'white\');
end

set(hObject, \'String\', {\'plot(rand(5))\', \'plot(sin(1:0.01:25))\', \'bar(1:.5:10)\', \'plot(membrane)\', \'surf(peaks)\'});


% --- Executes on button press in pushbutton4.
%新增一行按钮触发事件
function pushbutton4_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton4 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
data = get(handles.uitable4,\'data\');
 [a b] = size(data);
 data(a+1,:) = {\'\'};
 set(handles.uitable4, \'data\', data);



function edit1_Callback(hObject, eventdata, handles)
% hObject    handle to edit1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,\'String\') returns contents of edit1 as text
%        str2double(get(hObject,\'String\')) returns contents of edit1 as a double


% --- Executes during object creation, after setting all properties.
function edit1_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,\'BackgroundColor\'), get(0,\'defaultUicontrolBackgroundColor\'))
    set(hObject,\'BackgroundColor\',\'white\');
end



function edit2_Callback(hObject, eventdata, handles)
% hObject    handle to edit2 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,\'String\') returns contents of edit2 as text
%        str2double(get(hObject,\'String\')) returns contents of edit2 as a double


% --- Executes during object creation, after setting all properties.
function edit2_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit2 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,\'BackgroundColor\'), get(0,\'defaultUicontrolBackgroundColor\'))
    set(hObject,\'BackgroundColor\',\'white\');
end


% --- Executes on button press in pushbutton5.
%计算按钮点击触发事件
function pushbutton5_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton5 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
data = get(handles.uitable4,\'data\');
v1=str2double(data(1,1));
v2=str2double(data(2,1));
p1=str2double(data(1,2));
p2=str2double(data(2,2));
w=0;
w1=0;
q=0;
qt=0;
q1=0;
Cp=1.050;
Cv=0.762
m=Cp/Cv;
%每个节点T和T1的数量关系relCur
relCur=0.0;
%第一个节点计算及画图start
if(v1==v2)
    w1=0;
    q1=Cv*((p2/p1)-1);
    relCur=p2/p1;
    axes(handles.axes1);
    cla;
    x=0:pi/200:v1;
    plot([v1,v2],[p1,p2],\'LineWidth\',2);
    hold on;
else
    axes(handles.axes1);
    cla;
    n=(log10(p1)-log10(p2))/(log10(v2)-log10(v1));
    k1=p1*(v1^n);
    k2=p2*(v2^n);
    if(round(k1)==round(k2) | (k1==k2))
        if((n-1)~=0)
            q1=(n-m)/(n-1)*Cv*((v1^(n-1))/(v2^(n-1)) -1);
        else
            q1=0.287*log(v2 /v1);
        end
        relCur=(v1^(n-1))/(v2^(n-1));
        if(v1<v2)
           fun=@(v) (((v1.^n)*p1)./(v.^n));
            w1=quad(fun,v1,v2); 
           
            x=v1:pi/200:v2;
            y=k1./(x.^n);
            plot(x,y,\'LineWidth\',2);
           
            hold on;
           
        else
            fun=@(v) (((v2.^n)*p2)./(v.^n));
            w1=-quad(fun,v2,v1);
            x=v2:pi/200:v1;
            y=k1./(x.^n);
            plot(x,y,\'LineWidth\',2);
            
            hold on;
        end
        
    else
        disp(\'k1<>k2\');
    end
    
end
q=q+q1;
if(q1>=0)
    qt=qt+q1;
end
%第一个节点计算及画图end
w=w+w1;
%遍历2-n-1个节点，进行计算及绘制曲线 end
 [a b] = size(data);
 
 i=2;
 axes(handles.axes1);

 vv2=0.0;
 
 pp2=0.0;
 while i<=a-1
     qq1=0;
     ww1=0;
     vv1=str2double(data(i,1));
     pp1=str2double(data(i,2));
     vv2=str2double(data(i+1,1));
     pp2=str2double(data(i+1,2));
     if(vv1==vv2)
        qq1=Cv*((pp2/pp1)-1)*relCur;
        relCur=relCur*(pp2/pp1);
        plot([vv1,vv2],[pp1,pp2],\'LineWidth\',2);
        hold on;
     else
        n=(log10(pp1)-log10(pp2))/(log10(vv2)-log10(vv1));
        kk1=pp1*(vv1^n);
        kk2=pp2*(vv2^n);
        if(round(kk1)==round(kk2) | (kk1==kk2))
            if((n-1)~=0)
                qq1=(n-m)/(n-1)*Cv*((vv1^(n-1))/(vv2^(n-1)) -1)*relCur;
            else
                qq1=0.287*log(vv2/vv1)*relCur;
            end
             relCur=relCur*((vv1^(n-1))/(vv2^(n-1)));
            if(vv1<vv2)
             fun=@(v) (((vv1.^n)*pp1)./(v.^n));
             ww1=quad(fun,vv1,vv2);
               % syms v;
             % ww1=int(((vv1^n)*pp1)/(v^n),v,vv1,vv2);
             x=vv1:pi/200:vv2;
             y=kk1./(x.^n);
             plot(x,y,\'LineWidth\',2);
             hold on;
            else
               fun=@(v) (((vv2.^n)*pp2)./(v.^n));
               ww1=-quad(fun,vv2,vv1);
               x=vv2:pi/200:vv1;
               y=kk1./(x.^n);
                plot(x,y,\'LineWidth\',2);
                hold on; 
            end
            w=w+ww1;
        else
            disp(\'kk1<>kk2\');
        end
       
     end 
     i=i+1;
     q=q+qq1;
     if(qq1>=0)
         qt=qt+qq1;
     end
 end
 %遍历2-n-1个节点，进行计算及绘制曲线 end
 %最后一个节点的计算及绘图处理 start
 qqq1=0;
 if(vv2==v1)
     qqq1=Cv*((p1/pp2)-1)*relCur;
     plot([v1,v1],[pp2,p1],\'LineWidth\',2);
 else
     www1=0;
     n=(log10(pp2)-log10(p1))/(log10(v1)-log10(vv2));
     kkk1=pp2*(vv2^n);
     kkk2=p1*(v1^n);
     if(round(kkk1)==round(kkk2) | (kkk1==kkk2))
         if((n-1)~=0)
            qqq1=(n-m)/(n-1)*Cv*((vv2^(n-1))/(v1^(n-1)) -1)*relCur;
         else
            qqq1=0.287*log(v1/vv2)*relCur; 
         end
         if(vv2<v1)
             %syms v;
              %www1=int(((vv2^n)*pp2)/(v^n),v,vv2,v1);
              fun=@(v) (((vv2.^n)*pp2)./(v.^n));
             www1=quad(fun,vv2,v1);
         x=vv2:pi/200:v1;
         y=kkk1./(x.^n);
         plot(x,y,\'LineWidth\',2);
         else
              fun=@(v) (((v1.^n)*p1)./(v.^n));
             www1=-quad(fun,v1,vv2);
            x=v1:pi/200:vv2;
            y=kkk1./(x.^n);
            plot(x,y,\'LineWidth\',2); 
         end
         w=w+www1;
     else
         disp(\'kkk1<>kkk2\');
     end
     
 end
 q=q+qqq1;
 if(qqq1>=0)
 qt=qt+qqq1;
 end
  %最后一个节点的计算及绘图处理 end
  %GUI界面赋值
  j=1;
 
  [a b] = size(data);
  maxV=0.0;
  maxP=0.0;
  while j<=a
      if(str2double(data(j,1))>maxV)
          maxV=str2double(data(j,1));
      else
      end
      if(str2double(data(j,2))>maxP)
          maxP=str2double(data(j,2));
      else
      end
      j=j+1;
  end
   xlim([0,round(maxV*1.1)]);
   ylim([0,round(maxP*1.1)]);
 set(handles.text13,\'String\',num2str(q/qt));
 set(handles.text11,\'String\',[num2str(q),\'T1\']);
 set(handles.text9, \'String\', num2str(w));
disp(\'刘晓勇\');


% --- Executes on button press in pushbutton6.
%删除一行按钮触发事件
function pushbutton6_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton6 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
data = get(handles.uitable4,\'data\');
 [a b] = size(data);
 if a<=3
     
 else
     data(end,:)=[];
     set(handles.uitable4, \'data\', data);
 end