请问如何得到这样的曲线族？

20wangz

请问该如何画出这样的曲线族？

有两个自变量，pai和pait，由它们可导出中间变量tau3，并可进一步导得中间变量a0，最后得到结果w.
（由tau3导出a0的方法为a0=-Rg*tau3*log(pait)/cp，Rg和cp为已知常量。）

现在由其它的条件得知，有EH1=a0*Cw/(CH1*(tauH3-tau3))，其中Cw,CH1,tauH3为已知常量。如果我想得到EH1分别取０,０.３,０.６,０.９和１.０时以pai为橫坐标，以w为纵坐标的曲线族，我该怎么办呢？


极为感谢！！！

（附上我的不完整的程序，主要是推导出的结果。请注意，程序是不完整的。）

clear all;
clc;
%
%
k=1.4;m=(k-1)/k;Cw=1.0;
Rg=287.05;cp=1004.68;
CH=1.2;CL=1.2;CH1=1.2;
CHmin=min(Cw,CH);CLmin=min(CL,Cw);
tauH1=4;tauL1=1; tauH3=4.67;
etaC=0.9;etaT=0.9;
EH=0.9;EL=EH;ER=EH;
  PAI=1:0.1:20;
  EEH1=[0,0.３,0.6,0.9,1];
  gong=zeros(1,length(EEH1));
  w=zeros(length(EEH1),length(PAI));
%   
%
for jj=1:length(EEH1)
     EH1=EEH1(jj);
     ii=0;
     for pai=PAI
         ii=ii+1;
         x=pai^m;y=pait^m;
         a1=(-1+x*y)*etaT;
         a2=-1+x+etaC;
         a3=1-x+(-1+x)*ER-etaC;
         a4=x*y-a1;
         a5=(-1+x)*x*y+a1*a3;
         tau3=(x*y*(-a2*CLmin*Cw*EL*tauL1+CHmin*EH*(a2*CLmin*EL*(-ER*tauH1+tauL1)+Cw*tauH1*((-1+...
             x)*ER-etaC))))/(a5*Cw^2+a3*a4*CHmin*Cw*EH-a2*CLmin*EL*(a4*(Cw+CHmin*EH*(-1+ER))+a1*Cw*ER));
         a0=-Rg*tau3*log(pait);
         gong(ii)=(a0*(a5*Cw^2+a3*a4*CHmin*Cw*EH-a2*CLmin*EL*(a4*(Cw+CHmin*EH*(-1+ER))+a1*Cw*ER))+a5*Cw*...
             (CHmin*EH*tauH1+CLmin*EL*tauL1)+CLmin*EL*(a4*CHmin*EH*(etaC*tauH1+a3*tauL1)...
             -a2*(-x*y*(CHmin*EH+CLmin*EL)*tauL1+a4*(CHmin*EH*tauH1+CLmin*EL*tauL1)+a1*(CHmin*EH*ER*tauH1+...
             CLmin*EL*tauL1))))/(a5*Cw^2+a3*a4*CHmin*Cw*EH-a2*CLmin*EL*(a4*(Cw+CHmin*EH*(-1+ER))+a1*Cw*ER));        
     end
     w(jj,:)=gong;     
end
plot(PAI,w)

lin2009

"想得到EH1分别取０,０.３,０.６,０.９和１.０"应先推导出对应pait的值，再画，如下。
clear all;
clc;
%
k = 1.4;
m = (k - 1)/k;
Cw = 1.0;
Rg = 287.05;
cp = 1004.68;
CH = 1.2;
CL = 1.2;
CH1 = 1.2;
CHmin = min(Cw, CH);
CLmin = min(CL, Cw);
tauH1 = 4;
tauL1 = 1;
tauH3 = 4.67;
etaC = 0.9;
etaT = 0.9;
EH = 0.9;
EL = EH;
ER = EH;
PAI = 1:0.1:20;
EEH1 = [ 0, 0.3, 0.6, 0.9, 1 ];
gong = zeros(1, length(EEH1));
w = zeros(length(EEH1), length(PAI));
%
%
for pait = 3:0.01:3.05 %仅说明问题，不代表实际值。
                  %"由EH1＝０,０.３,０.６,０.９和１.０"推导对应pait的值
   for jj = 1 : length(EEH1)
        EH1 = EEH1(jj);
        ii = 0;
        for pai = PAI
            ii = ii + 1;
            x = pai^m;
            y = pait^m;
            a1 = ( -1 + x*y)*etaT;
            a2 = -1 + x + etaC;
            a3 = 1 - x + ( -1 + x)*ER - etaC;
            a4 = x*y - a1;
            a5 = ( -1 + x)*x*y + a1*a3;
            tau3 = (x*y*( -a2*CLmin*Cw*EL*tauL1 + CHmin*EH*(a2*CLmin*EL*( -ER*tauH1 + tauL1) + Cw*tauH1*(( -1 + ...
                x)*ER - etaC))))/(a5*Cw^2 + a3*a4*CHmin*Cw*EH - a2*CLmin*EL*(a4*(Cw + CHmin*EH*( -1 + ER)) + a1*Cw*ER));
            a0 = -Rg*tau3*log(pait);
            gong(ii) = (a0*(a5*Cw^2 + a3*a4*CHmin*Cw*EH - a2*CLmin*EL*(a4*(Cw + CHmin*EH*( -1 + ER)) + a1*Cw*ER)) + a5*Cw*...
                (CHmin*EH*tauH1 + CLmin*EL*tauL1) + CLmin*EL*(a4*CHmin*EH*(etaC*tauH1 + a3*tauL1)...
                -a2*( -x*y*(CHmin*EH + CLmin*EL)*tauL1 + a4*(CHmin*EH*tauH1 + CLmin*EL*tauL1) + a1*(CHmin*EH*ER*tauH1 + ...
                CLmin*EL*tauL1))))/(a5*Cw^2 + a3*a4*CHmin*Cw*EH - a2*CLmin*EL*(a4*(Cw + CHmin*EH*( -1 + ER)) + a1*Cw*ER));
        end
        w(jj, :) = gong;
    end
    plot(PAI, w);
    hold on;
end

20wangz

2# lin2009

的确，要得到要求的图需要得到pait，可问题的关键是这一步是不可能直接完成的，因为我的提问给出的计算式中，EH1=a0*Cw/(CH1*(tauH3-tau3))，这里面又出现了tau3，而这个tau3又是pai和pait的函数。