关于MatLab中函数的细微差别的问题

Question

当我运行这段代码时：

dx = 4/(nx-1);
dy = 2/(ny-1);
phinew(2:ny-1,2:nx-1) = ((dy^2*(phiold(1:ny-2,2:nx-1)+phiold(3:ny,2:nx-1))+dx^2*(phiold(2:ny-1,1:nx-2)+phiold(2:ny-1,3:nx)))/(2*(dx^2+dy^2))-poissonf);

在我的代码中，这段代码运行正常，接受nx by ny矩阵的输入，并输出nx by ny矩阵。但是，当我在完全相同的位置运行以下代码时：

phinew = Smoothing(phiold,poissonf,nx,ny)

其中“平滑”函数定义为：

function [phinew,dx,dy] = Smoothing(phiold,poissonf,nx,ny)
        dx = 4/(nx-1);
        dy = 2/(ny-1);
        phinew(2:ny-1,2:nx-1) = ((dy^2*(phiold(1:ny-2,2:nx-1)+phiold(3:ny,2:nx-1))+dx^2*(phiold(2:ny-1,1:nx-2)+phiold(2:ny-1,3:nx)))/(2*(dx^2+dy^2))-poissonf);
end

此函数返回nx-1 by ny-1输入的nx by ny矩阵。

我不明白为什么会发生这样的事情。除了最后一列和最后一行完全丢失之外，输出矩阵与其应有的完全相同。我的代码是迭代的，因此需要它们具有相同的大小，所以在这个问题得到解决之前，我无法继续前进。

感谢您的时间和帮助。你们这些人是救命稻草。

Answer 1

当您在命令窗口中运行代码时，phinew已经存在。在命令窗口中，在粘贴这三行代码之前执行clear phinew，然后您会发现phinew被ny-1 nx-1了，就像您从函数中得到的一样。

如果您想强制您的函数按ny返回nx，请将phinew = zeros(nx,ny);放在函数的开头，或者将最后一列和最后一行设置为您想要的值。

编辑:回复你的评论“为什么主赋值的RHS输出一个48x48矩阵?在运行phiold命令之前，phinew和poissonf都是50x50。”

我不认为poissonf是50x50:这会导致错误Matrix dimensions must agree，因为poissonf被添加到表达式的其余部分，这是phiold的一部分，所以我将在下面忽略poissonf。

RHS总是ny-2 by nx-2，甚至在第一次迭代时也是如此。您可以通过将RHS分配给一个中间变量(例如phipiece = ... )并检查size(phipiece)来查看这一点。phinew (如果重新创建)是49x49的原因是因为它被分配给(2:ny-1,2:nx-1)，这将创建一个ny-1 by nx-1矩阵，并将第一行和第一列保留为零。

如果首先使用phinew = zeros(nx,ny);，则第一行和最后一列将保留为零。