all_zone_m.m

雷诺方程的数值解法

function [Delta]=all_zone_m()
ps_location=evalin('base','ps_location');
beta=evalin('base','beta');
n=evalin('base','n');
M=evalin('base','M');
N=evalin('base','N');
ps=[ps_location,beta];
num=size(ps);
for i=1:num(1)
    m(i)=m_in(ps(i,1),ps(i,2),ps(i,3));
end
min1=sum(m);
min2=mz(1,1,M+1);
min3=mz(N+1,1,M+1);
Dm=min1+min2-min3;
D_ratio=Dm/min1;
Delta=[min1,min2-min3,D_ratio,Dm];
end