MATLAB program to compute the potential in a square region by solving Laplace’s equation using the Finite Difference (FD) approximation method.using the “surf” and “contour” MATLAB commands

7 次查看（过去 30 天）

Sai Hitesh Gorantla 2020-6-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/543290-matlab-program-to-compute-the-potential-in-a-square-region-by-solving-laplace-s-equation-using-the-f

评论： KSSV 2020-6-7

My Code:

clear all;
[X,Y] = meshgrid(linspace(0,1,20));
V = zeros(size(X));
[M,N] = size(V);
V(2:M-1,1) = zeros(M-2,1);
V(2:M-1,N) = 100;
V(1,2:N-1) = 40;
V(M,2:N-1) = 10;
for r = 1:10000
    V(2:M-1,2:N-1) = 0.25*(V(3:M,2:N-1) + V(1:M-2,2:N-1) + V(2:M-1,2:N-1) + V(4:M-1,2:N-1));
end
surf(X,Y,V);

Getting Error:

>> contour(X,Y,V)

Unrecognized function or variable 'X'.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Sai Hitesh Gorantla 2020-6-6

please respond

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

KSSV 2020-6-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/543290-matlab-program-to-compute-the-potential-in-a-square-region-by-solving-laplace-s-equation-using-the-f#answer_447500

在 MATLAB Online 中打开

Just rectifies the error you said:

clc; clear all ;
clear all;
[X,Y] = meshgrid(linspace(0,1,20));
V = zeros(size(X));
[M,N] = size(V);
V(2:M-1,1) = zeros(M-2,1);
V(2:M-1,N) = 100;
V(1,2:N-1) = 40;
V(M,2:N-1) = 10;
for r = 1:10000
    V(2:M-1,2:N-1) = 0.25*(V(3:M,2:N-1) + V(1:M-2,2:N-1) ....
                     + V(2:M-1,2:N-1) + V(2:M-1,2:N-1));
end
surf(X,Y,V);