fdiff4.m

来自「薛定宇高等应用数学问题的MATLAB求解例子源码」· M 代码 · 共 10 行

   function [t,yd]=fdiff3(S,gam)
   tt=S.breaks; M=[ones(3,1); 0]; t=tt(2:end-1);
   aa=[3:-1:0]'; aa1=gamma(aa'+1)./gamma(aa'+1-gam);
   for i=1:length(t)
       kk=find(tt<t(i)); if length(kk)>0, k=kk(end); end
%       yd(i)=S.coefs(k,:)*[[(t(i)-tt(k))*ones(3,1);t(i)].^aa];
       %       yd(i)=(S.coefs(k,:).*aa1)*[(t(i)-tt(k))*M.^(aa-gam)];
       yd(i)=(S.coefs(k,:).*aa1)*[[(t(i)-tt(k))*ones(3,1); t(i)].^[aa-gam]];
   end