cyberck的博客

借用这个题目然后用比较简单的代码来求解%空间0到1，切了5份，时间0到1，切了十份。clc,clearv = 1;dx = 0.2;x = 0:dx:1;dt = 0.1;nx = 6;nt = 11;t = 0:dt:1;b = v/(2*dx*dx);c = b;a = 1/dt+b+c;Uold = exp(x);Unew = ones(1,nx); Unew(1) = exp(t(2));Unew(nx) = exp(1+t(2));%其实最关键的就是得

tic;clearclcM=[10,20,40,80,160,320,640];%x的步数K=M; %时间t的步数for p=1:length(M)hx=1/M(p);ht=1/K(p);r=ht/hx^2; %网格比x=0:hx:1;t=0:ht:1;numerical=zeros(M(p)+1,K(p)+1);numerical(:,1)=exp(x); %初始值numerical(1,:)=exp(t); %边值numerical(M(p)+1,:)=exp(t+1); .