cyberck的博客

借用这个题目 然后用比较简单的代码来求解 %空间0到1，切了5份，时间0到1，切了十份。 clc,clear v = 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; clear clc M=[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); .