residual.m

MDPSAS工具箱是马里兰大学开发的

function [A,g1,g1fine,Twfine,Tqfine] = residual(A)

% residual.m Nonlinear equation model object residual update.
%
% INPUT/OUTPUT PARAMETERS
%        A : an object derived from the neqmodel class

unpack(A)

g1 = e_wq*F*(Tq^4-Tw^4) + e_ws*(2-F)*(1-Tw^4) ...
    + k_w/alf*(Tq-Tw) + q_w*u;
g1 = Dt*Tw - g1;
g1 = setbval(g1,Tw-Tw0,'t','min');

g2 = e_qw*F*(Tw^4-Tq^4) + e_qs*(1-F)*(1-Tq^4) ...
    + nu_q*(1-Tq) + q_q*u;
g2 = Dt*Tq - g2;

periflag = 0;
if periflag
    Tdiff = getbval(Tq,'t','min')-getbval(Tq,'t','max');
    g2 = setbval(g2,Tdiff,'t','min');
else
    g2 = setbval(g2,Tq-Tq0,'t','min');
end

A = set(A,'resid',{g1,g2});

if nargout > 1
    tfine = quadgrid('slab',100,'t',[25 35]);
    tfine = relquad(tfine,0);
    Dtfine = linearoperator(tfine,'d','t');

    Twfine = interp(Tw,tfine);
    Tqfine = interp(Tq,tfine);
    ufine  = interp(u,tfine);

    g1fine = e_wq*F*(Tqfine^4-Twfine^4) + e_ws*(2-F)*(1-Twfine^4) ...
        + k_w/alf*(Tqfine-Twfine) + q_w*ufine;
    g1fine = Dtfine*Twfine - g1fine;
end