matlab編程天然氣壓縮因子計(jì)算模型(編輯修改稿)
2025-06-19 18:26 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 ... ,,,,。... ,,,,。... ,,,,。... ,,,,。... ,,,,。... ,,,,。... ,,,,。... ,,,,。... ,,,,]。Kx=[,,,,。 ,,,,。 ,,,,。 ,,,,。 ,,,,。 ,,,,。 ,,,,。 ,,,,。 ,,,,。 ,,,,。 ,,,,。 ,,,,。 ,,,,。 ,,,,]。Z0=1。M0=0。B=0。%參數(shù)Eij、Gij以及第二維利系數(shù)的計(jì)算sum1=0。sum2=0。 for n=1:18 sum=0。 ZJCS=(T^u(n))。 for i=1:14 for j=1:14 Eij(i,j)=Ex(i,j)*(E(i)*E(j)^)。 Gij(i,j)=Gx(i,j)*(G(i)+G(j))/2。 Bij(i,j)=((Gij(i,j)+1g(n))^g(n))*((Q(i)*Q(j)+1q(n))^q(n))*(((F(i)*F(j)^)+1f(n))^f(n))*((S(i)*S(j)+1s(n))^s(n))*((W(i)*W(j)+1w(n))^w(n))。 sum=sum+x(i)*x(j)*Bij(i,j)*(K(i)*K(j)^)。 end end B=B+a(n)*ZJCS*sum。 end%K值的計(jì)算 F0=0。 G0=0。 Q0=0。 U0=0。for i=1:N F0=F0+(x(i)^2)*F(i)。 Q0=Q0+x(i)*Q(i)。 sum1=sum1+x(i)*G(i)。 sum2=sum2+x(i)*(E(i)^)。endfor i=1:(N1) for j=2:N G0=G0+x(i)*x(j)*(Gx(i,j)1)*(G(i)+G(j))。 U0=U0+((Ux(i,j)^5)1)*((E(i)*E(j))^)。 endendG0=sum1+2*G0。U0=((sum2^2)+2*U0)^。sum1=0。for i=1:N sum1=sum1+x(i)*(K(i)^)。endsum2=0。for i=1:(N1) for j=2:N sum2=sum2+x(i)*x(j)*((Kx(i,j)^5)1)*((K(i)*K(j))^)。 endendK0=(((sum1^2)+2*sum2)^)。%對(duì)比密度的計(jì)算pr=p*(K0^3)/(Z0*R*T)。 SUM1=0。SUM2=0。 %計(jì)算Cn的值，共有46個(gè)值for n=13:18 Cn=a(n)*((G0+1g(n))^g(n))*(((Q0^2)+Q0q(n))^q(n))*((F0+1f(n))^f(n))*(U0^u(n))*(T^u(n))。 SUM1=SUM1+Cn。endfor i=13:58 Cn=a(i)*((G0+1g(i))^g(i))*((Q0*Q0+Q0q(i))^q(i))*((F0+1f(i))^f(i
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1