Question

我哪里出错了？问题在于： R = sqrt（bsxfun（@ minus，XX，（XX）'）。^ 2 + bsxfun（@ minus，YY，（YY）'）。^ 2）; < / p>

EO = 8.8541e-12; %eps0
A2 = 1.0e-2; %2a
N = 100; %num of subareas in a plate
M = sqrt(N); %num of subareas in one axis
DX = A2/M; % 2b
DY = DX; %2b
DL = DX; %2b

%  SECOND, CALCULATE THE ELEMENS OF THE COEFFICIENT MATRIX A
% Write all the subareas centers coordinates to X and Y 
%  SECOND, CALCULATE THE ELEMENS OF THE COEFFICIENT MATRIX A
% Write all the subareas centers coordinates to X and Y 
x = linspace(DL*0.5,DL*(M-0.5),M);
y = x;
[XX,YY]=meshgrid(x,y);
%%L = NaN(N,N);
**R = sqrt(bsxfun(@minus,XX,(XX)').^2+bsxfun(@minus,YY,(YY)').^2);**
idx_diagL = find(eye(N)~=0);
idx_not_diagL = find(eye(N)==0);
L(idx_not_diagL) = DL^2./(4.*pi*EO*R(mod(idx_not_diagL,10),floor(idx_not_diagL/N)));
L(idx_diagL)  = DL*0.8814/(pi*EO); %ln(1+sqrt(2)= 0.8814

Answer 1

问题是不您说的是，但在以下行中：

L(idx_not_diagL) = DL^2./(4.*pi*EO*R(mod(idx_not_diagL,10),floor(idx_not_diagL/N)));

简而言之，您的问题是您为矩阵R指定零索引，但在MATLAB中零索引是非法的（它们从1开始，而不是0！）。

现在，你在哪里获得零指数？您可以使用以下行和列下标索引到R(..., ...)：

mod(idx_not_diagL, 10)

和

floor(idx_not_diagL / N))

偶尔会得到零值。

作为修复，我建议对行索引使用以下内容：

mod(idx_not_diagL, 10) + 1

对于列索引，还有另一个问题。 floor(...)表达式的范围是0到99.一旦你修复它（我不确定你想要实现什么），它就可以工作。

Matlab-Subscript索引必须是实数正整数或逻辑

1 个答案: