请教为什么输出数组不受我子程序设定的最大最小值限制，值传不到主程序里
function ex%--------------------------------------------------------------------------
L = 50;
s1 = 0.5;
s2 = 0;
T = 5;
qr = 0.1699;
f = 0.4432;
a = 0.0526;
n = 1.3924;
ks = 31.6;
aa=0.0037;
bb=14.029;
x = linspace(0,L,50);
t = linspace(0,T,25);
options=odeset('RelTol',1e-4,'AbsTol',1e-4,'NormControl','off','InitialStep',1e-7)
sol = pdepe(0,@unsatpde,@unsatic,@unsatbc,x,t,options,s1,s2,qr,f,a,n,ks);
u1=real(sol(:,:,1))
figure
surf(x,t,u1)
% -------------------------------------------------------------------------
function = unsatpde(x,t,u,DuDx,s1,s2,qr,f,a,n,ks)
= sedprop(u,qr,f,a,n,ks);
c=;
f=;
s = ;
% -------------------------------------------------------------------------
function u0 = unsatic(x,s1,s2,qr,f,a,n,ks)
u0 = ;
% -------------------------------------------------------------------------
function = unsatbc(xl,ul,xr,ur,t,s1,s2,qr,f,a,n,ks)
pl = [-0.5];
ql = [-1];
pr = ;
qr = ;
%-------------------------------------------------------------------------
function = sedprop(u,qr,f,a,n,ks)
m = 1-1/n;
if u<qr
c=1e-20;
k=0;
u=qr;
Dth=aa*exp(bb*qr);
elseif u>f
c=1e-20;
k=ks;
u=f;
Dth=aa*exp(bb*f);
else
k=ks*((u-qr)/(f-qr))^0.5*(1-(1-((u-qr)/(f-qr))^(1/m))^m)^2;
Dth=aa*exp(bb*u);
end
我想通过子程序控制u的范围在之间，即，但是u的输出结果很多不在这个范围之间，是怎么回事？
如：
输出结果：
Columns 49 through 50
0.4432 0
0.8439 0
0.8176 0
0.7902 0
0.7748 0
0.7645 0
0.7491 0
0.7385 0
0.7302 0
0.7236 0
0.7168 0
0.7078 0
0.7037 0
0.7016 0
0.7017 0
0.7053 0
0.7054 0
0.7036 0
0.7064 0
0.7033 0
0.7060 0
0.7040 0
0.7058 0
0.7036 0
0.7058 0 子程序传递有问题？这个问题怎么解决呢？请教大家 也没有人理我:@(，提个建议也好啊
页:
[1]