函数输出 = poisson1D(f,Uleft,Uright,N,L) %POISSON1D求解泊松方程d2U / dX2 =F。 % u = poisson1D(f,Uleft,Uright,N,L) % % f:代表右手边的向量% Uleft:u(0) 处的狄利克雷边界条件% Uright:u(L) 处的狄利克雷边界条件% N:节点数% L:域的长度dx = L / (N - 1); uB = 零(长度(f),1); uB(2)= Uleft; uB(end-1) = 右; f = f * dx^2 - uB; b = dst(f); m = (1:length(b))'; a = b ./ (2 * ( cos(m * pi / (N - 1)) - 1 )); uSOL = idst(a); uSOL(1) = 左左; uSOL(end) = Uright; 出 = uSOL;