evalvar.m

来自「MDPSAS工具箱是马里兰大学开发的」· M 代码 · 共 25 行

function A = evalvar(A)

% INPUT/OUTPUT PARAMETERS
%        A : an object derived from the naemodel class
%
% Copyright (c) by Raymond A. Adomaitis, 1998-2003

unpack(A)

r = get(S,'qp');
r = r{1};

theta = [-pi:0.01:pi];
for i = 1:N
   qlamp(i,1) = (h_lamp*Q_lamp/8/pi^2) ...
   *sum( (h_lamp^2+r(i)^2+r_lamp^2-2*r_lamp*r(i)*cos(theta)).^(-3/2) )*0.01;
   end

qlamp = 0.70*qlamp; % 0.3 is absorbed by quartz window
qlamp = u*qlamp;    % Operate at a fraction of u

% qlamp = Q - scalarfield(S,qlamp);
Q = scalarfield(S,qlamp);
A = set(A,'var',Q,'Q');