利用偏微分方程学习偏微分方程

clear;
N = 7; M = 7; T = 10;
c = 2; sigma = 2;
dat = DataGen(N,M,T,c,sigma);
[a,f] = learner(dat);

%%%% we check the solution a going with sigma %%%%
for sigma = 0.5:0.1:3
	dat = DataGen(N,M,T,c,sigma);
	[a,f] = learner(dat);
	disp(sum(sum(a)));
end

clear;
N = 7; M = 7; T = 10;
c = 2; sigma = 2;
dat = DataGen(N,M,T,c,sigma);
[a,f] = learner(dat);

%%%% we check the solution f going with c %%%%
for c = 1:0.05:2
	dat = DataGen(N,M,T,c,sigma);
	[a,f] = learner(dat);
	disp(sum(f));
end

function ret = DataGen(N,M,T,c,sigma)
% This function is to generate three-dimension data matrix of the PDE
% \frac{
   \partial Q}{
   \partial t} = - c(x,y)Q(x,y,t) -\int
% p(i,j;x,y)Q(i,j,t)didj
% we set c(x,y)\equiv c = 2 and p N(0,\sigma)
lambda = 3