函数 u = poisson1Dneumann(F,x0,xEnd) %POISSON1DNUEMANN用Neumann求解一维泊松方程d2U / dX2 = F % 边界条件 dUdX = 0 在 X = 0 和 X = L。 % u = poisson1Dneumann(F,x0,xEnd) % % u:解向量% F:右侧向量% x0:域的起始坐标。 % xEnd:域的结束坐标。 % 检查兼容性xInt = linspace(x0,xEnd,length(F)); fInt = trapz(xInt,F); 如果 (fInt > 0.0001) || (fInt < -0.0001) disp('不满足兼容条件'); 结尾% 解决方案N = 长度(F); dx = (xEnd - x0) / (N - 1);
b = dct(F); m = (0:length(b)-1)'; a
2022-05-22 17:17:09
2KB
matlab
1