MATLAB 矩阵构造建议

Question

我需要使用 MATLAB 求解以下矩阵方程：

这个矩阵是 nxn，我可以使用以下 MATLAB 代码构造大部分矩阵：

e = ones(n,1);
A = h^(-2)*spdiags([e (h^2 - 2)*e e], [-1 0 1], n, n);

更改第一行和最后一行的好方法是什么？也许只需添加一个不错的矩阵 B ，第一行为 [ 2/h 1/h^2 0 ... 0 0 0 ]，最后一行为 [ 0 0 0 ... 0 1/h ^2 (2h + 1)/h^2] 只取 A + B。不过你会怎么做呢？

Answer 1

我认为在这种情况下最简单的方法是最好的，因为您不会修改您创建的大部分矩阵：

A(1,:)=A(1,:)+[2/h 1/h^2 zeros(1,n-2)];
A(n,:)=A(n,:)+[zeros(1,n-2) 1/h^2 2/h];

甚至替换单个元素而不是行。