how can I make a GUI for this script?

Question

0 个投票

clc,clear
close all
disp('Formation Water Properties');
S=input('Salinity(Weight percent solids)= ');
T=input('The temperature in degrees F = ');
P=input('Pressure in psia = ');
Pb=input('Bubble point Pressure = ');
z=input('Solution Gas Water Ratio at Pb ');
  fprintf('%5s  %5s  %5s  %5s       %s     %5s\n','P','Rsw','Bw','D','cw','Uw');
for i=100:100:P
    if i<Pb
    P(i)=i;
  A0=8.15839; A1=-6.1226510^-2; A2=1.91663*10^-4; A3=-2.1654*10^-7;
  A= A0+A1*T+A2*T^2+A3*T^3;
  B0=1.01021*10^-2; B1=-7.44241*10^-5; B2=3.0555*10^-7; B3=-2.94883*10^-10;
  B=B0+B1*T+B2*T^2+B3*T^3;
  C0=-9.02505; C1=0.130237; C2=-8.53425*10^-4; C3=2.34122*10^-6; 
  C4=-2.37049*10^-9;
  C = (C0 + C1*T + C2*T^2 + C3*T^3 + C4*T^4)*(10^-7);
  Rsw(i)=A + B* P(i) + C*P(i);
    Rswb(i)=Rsw(i)*10^(-0.0840655*S*T^-.285854);
    X=-1.0001*(10^-2)+1.3339*(10^-4)*T+5.50654*(10^-7)*T^2;
    Y(i)=-1.95301*(10^-9)*P(i)*T -1.72834*(10^-13)*P(i)^2*T-3.58922*(10^-7)*P(i)-2.25341*(10^-10)*P(i)^2;
    Bw(i)=(1+Y(i))*(1+X);
    mi=62.368+0.438603*S + 1.60074*(10^-3)*S^2;
    Aw(i)=mi/Bw(i);
    A1 = 7.033; A2 = 0.5415; A3 = -537.0; A4 = 403300;
    cw(i)= 1 / (A1*P(i) + A2*S + A3*T + A4);
    A0=109.574; A1=-8.40564; A2=0.313314; A3=8.72213*10^-3;
    A = A0+ A1*S + A2*S^2 + A3*S^3;
    B0=-1.12166; B1=2.63951*10^-2; B2=-6.79461*10^-4; B3=-5.47119*10^-5; 
    B4=1.55586*10^-6;
    B = B0 + B1*S + B2*S^2+B3*S^3+B4*S^4;
    Uo=A*T^B;
    Uw(i)=Uo*(0.9994+4.0295e-05*P(i)+3.1062e-09*P(i)^2);
    fprintf('%5d  %5.2f  %5.2f  %5.2f  %5d  %5.3f\n',P(i),Rswb(i),Bw(i),Aw(i),cw(i),Uw(i));
      end
      if i>=Pb
             P(i)=i;
             A0=8.15839; A1=-6.1226510^-2; A2=1.91663*10^-4; A3=-2.1654*10^-7;
    A= A0+A1*T+A2*T^2+A3*T^3;
    B0=1.01021*10^-2; B1=-7.44241*10^-5; B2=3.0555*10^-7; B3=-2.94883*10^-10;
    B=B0+B1*T+B2*T^2+B3*T^3;
    C0=-9.02505; C1=0.130237; C2=-8.53425*10^-4; C3=2.34122*10^-6; 
    C4=-2.37049*10^-9;
    C = (C0 + C1*T + C2*T^2 + C3*T^3 + C4*T^4)*(10^-7);
    Rsw(i)=A + B* P(i) + C*P(i);
    Rswb(i)=Rsw(i)*10^(-0.0840655*S*T^-.285854);
     X=-1.0001*(10^-2)+1.3339*(10^-4)*T+5.50654*(10^-7)*T^2;
    Y(i)=-1.95301*(10^-9)*P(i)*T -1.72834*(10^-13)*P(i)^2*T-3.58922*(10^-7)*P(i)-2.25341*(10^-10)*P(i)^2;
    Bw(i)=(1+Y(i))*(1+X);
    mi=62.368+0.438603*S + 1.60074*(10^-3)*S^2;
    Aw(i)=mi/Bw(i);
    A1 = 7.033; A2 = 0.5415; A3 = -537.0; A4 = 403300;
    cw(i)= 1 / (A1*P(i) + A2*S + A3*T + A4);
    A0=109.574; A1=-8.40564; A2=0.313314; A3=8.72213*10^-3;
    A = A0+ A1*S + A2*S^2 + A3*S^3;
    B0=-1.12166; B1=2.63951*10^-2; B2=-6.79461*10^-4; B3=-5.47119*10^-5; 
    B4=1.55586*10^-6;
    B = B0 + B1*S + B2*S^2+B3*S^3+B4*S^4;
    Uo=A*T^B;
    Uw(i)=Uo*(0.9994+4.0295e-05*P(i)+3.1062e-09*P(i)^2);
    fprintf('%5d  %5.2f  %5.2f  %5.2f  %5d  %5.3f\n',P(i),Rswb(i),Bw(i),Aw(i),cw(i),Uw(i));
        end 
    end
    i=100:100:P;
    figure('name','p& Rsw ');
    plot(P(i), Rswb(i));
    xlabel('Pressure');
    ylabel('Rsw, scf/STB');
    title('p& Rsw ');
      figure('name','p& Bw ');
    plot(P(i), Bw(i));
    xlabel('Pressure');
    ylabel('Bw, res bbl/STB');
    title('p& Bw ');
    figure('name','p& D ');
    plot(P(i), Aw(i));
    xlabel('Pressure');
    ylabel('Dw, lb/cu ft');
    title('p& D ');
    figure('name','p&Cw');
    plot(P(i), cw(i));
    xlabel('Pressure');
    ylabel('Cw, 1/psi');
    title('p& Cw ');
    figure('name','p&Uw');
    plot(P(i), Uw(i));
    xlabel('Pressure');
    ylabel('UW,cp');
    title('p& Uw ');