pvt.m

来自「《实用化工计算机模拟：MATLAB在化学工程中的应用 》这本书光盘里的程序~」· M 代码 · 共 24 行

function PVT
% 已知P和T，计算V
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2003/02/26 $

clear all
clc

P = 9.33;   % atm
T = 300.2;  % K
n = 2;      % mol
a = 4.17;
b = 0.0371;
R = 0.08206;
V0 = n*R*T/P
[V,fval] = fzero(@PVTeq,V0,[],P,T,n,a,b,R)

% ------------------------------------------------------------------
function f = PVTeq(V,P,T,n,a,b,R)
f = (P + a*n^2/V^2) * (V-n*b) - n*R*T;