求教关于PDE求解偏微分方程

follwind

小弟做的是热传导偏微分方程求解，用pdepe求解，按照matlab自带的例子编了一段程序，有两个偏微分方程，[c,f,s] = pdexpde(x,t,T,DTDx,Ti)，也就是s有两种情况，一种是不带内热源，一种是带内热源，所处时间段不同则偏微分方程不同，现在令我苦恼的是我把初始温度设成924时没有任何问题，能求解完成，但是我把初始温度设置成950或是更高就不行了，本来一共是求解14段，结果求到中途就停止了，不同的初始温度中止的时间也不一样，温度越高中止的时间越早，下面是求解中止时matlab命令窗口显示的错误信息，我整了好几天也没能整明白问题出在哪里，想哭啊，还请各位高人不吝赐教！！
     非常谢谢！！！
Warning: Failure at t=4.794292e+001.  Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.703275e-013) at time t.
(Type "warning off MATLAB:ode15s:IntegrationTolNotMet" to suppress this warning.)
In D:\MATLAB\toolbox\matlab\funfun\ode15s.m at line 658
  In D:\MATLAB\toolbox\matlab\funfun\pdepe.m at line 287
  In E:\follwind\TEMPSIM.m at line 85??? Subscripted assignment dimension mismatch.
Error in ==> D:\MATLAB\toolbox\matlab\funfun\pdepe.mOn line 291  ==>   sol(:,:,i) = y(:,i:npde:end);
Error in ==> E:\follwind\TEMPSIM.m
On line 85  ==>         sol=pdepe(m,@pdexpde,@pdexic,@pdexbc,x,t,[],Ti);

[ 本帖最后由 xinyuxf 于 2007-7-12 15:05 编辑 ]