chap4_7mended3.m

利用MATLAB解决<偏微分方程数值算法>中的一些问题,如插值求解

function[E]=chap4_7mended3(h,tao)   %习题4.7
 a=1;s=a*tao/h;
for k=1:1/tao+1
  for i=1:1/h+1
      U(k,i)=0;
  end
end
for i=1:1/h+1
    X(1,i)=(i-1)*h;U(1,i)=0;
end
for i=1:1/h+1
    X(1,i)=(i-1)*h;U(2,i)=tao*X(1,i)+U(1,i);
end
for k=1:1/tao+1
    T(1,k)=(k-1)*tao;U(k,1)=0;U(k,1/h+1)=sin(T(1,k));
end
for k=3:1/tao+1
    for i=2:1/h
    U(k,i)=2*(1-s^2)*U(k-1,i)+s^2*U(k-1,i-1)+s^2*U(k-1,i+1)-U(k-2,i)+tao^2*(T(1,k-1)^2-X(1,i)^2)*sin(X(1,i)*T(1,k-1));
end
end
k=1/tao+1
  for i=1:1/h+1
      W(k,i)=sin(0.5*T(1,k));E(k,i)=abs(W(k,i)-U(k,i));
  end
  plot(X,E(k,i));
  hold on;