Question

matlab 利用有限差分法解偏微分方程 矩阵out of memory

利用有限差分法解偏微分方程遇到的问题：由于采用以下格式，所以矩阵很大而很容易超出矩阵的最大容量。现nr的取值只能去到20，想将nr取值越大越好，最少达到50，请问能有什么办法。当nr取值为25时的报警为
??? Error using ==>zeros
Out of memory. Type HELP MEMORY for your options.

Error in ==> FDM_forming0623_arbitrary_point_load at 18
a=sparse(zeros(n,n));

%程序
nr=20;nct=nr;nt=nr;n=nr*nct*nt;
a=sparse(zeros(n,n));
b=sparse(zeros(n,1));
r=sparse(zeros(n,1));ct=r;t=r;
for i=1:nr
for j=1:nct
for m=1:nt
……
end
end
end
通过一系列for循环将数据写入矩阵a和b；
%方程求解
u=a\b;

Answer 1

这个真好像没有办法，matlab矩阵太大了就是不行！

Answer 2

算法不能修改嘛?重新设计一下算法吧。