错误:尝试引用非结构数组的字段.(使用全局变量)

Question

我正在尝试修改其他人使用全局变量的代码，并不断收到我无法弄清楚的错误。这是我第一次使用全局变量，所以我怀疑它与此有关。 我收到此错误：

Attempt to reference field of non-structure array.
Error in likelihood (line 35)
X=ModelInfo.X;

这是我运行的代码：

clearvars -global
clearvars
global ModelInfo
%Pull in data from MySQL database
ModelInfo.X=.....   %76x2 double type
ModelInfo.y=.....   %76x1 double type
addpath('C:/.../Sampling Plans');
addpath('C:/.../Constructing a Surrogate');
%Number of variables
k=size(ModelInfo.X,2);
%Number of sample points
n=size(ModelInfo.X,1);
%Set upper and lower bounds for search of log theta
UpperTheta=ones(1,k).*2;
LowerTheta=ones(1,k).*-3;
%Run GA search of likelihood
[ModelInfo.Theta,MinNegLnLikelihood]=ga(@likelihood,k,[],[],[],[],LowerTheta,UpperTheta);

。它引用了一个外部函数可能性（由其他人编写）：

function [NegLnLike,Psi,U]=likelihood(x)
global ModelInfo
X=ModelInfo.X;
y=ModelInfo.y;
theta=10.^x;
p=2; 
n=size(X,1);
one=ones(n,1);
% Pre-allocate memory
Psi=zeros(n,n);
% Build upper half of correlation matrix
for i=1:n
    for j=i+1:n
        Psi(i,j)=exp(-sum(theta.*abs(X(i,:)-X(j,:)).^p)); 
    end
end
% Add upper and lower halves and diagonal of ones plus 
% small number to reduce ill-conditioning
Psi=Psi+Psi'+eye(n)+eye(n).*eps; 
% Cholesky factorisation
[U,p]=chol(Psi);
% Use penalty if ill-conditioned
if p>0
    NegLnLike=1e4;
else
    % Sum lns of diagonal to find ln(abs(det(Psi)))
    LnDetPsi=2*sum(log(abs(diag(U))));
    % Use back-substitution of Cholesky instead of inverse
    mu=(one'*(U(U'y)))/(one'*(U(U'one)));
    SigmaSqr=((y-one*mu)'*(U(U'(y-one*mu))))/n;
    NegLnLike=-1*(-(n/2)*log(SigmaSqr) - 0.5*LnDetPsi);
end

更新：

想通了问题所在。 我在开始时将 ModelInfo 声明为全局变量，然后在%Pull in data from MySQL database部分中，我创建了一个数据库连接，然后使用 clearvars -except 命令清除了一些我们不再需要的变量，但我没有在该列表中包含 ModelInfo。 之后，我将数据放入模型信息中。 由于它是全局的，我想它删除了本地副本或全局副本（不太确定），但它仍然出现在工作区的变量列表中，并且isstruct(ModelInfo)仍然返回 1......所以很明显，除了错误之外，没有任何迹象表明存在问题。

Answer 1

如果这确实是您运行的代码：

%Pull in data from MySQL database
ModelInfo.X=.....   %76x2 double type
ModelInfo.y=.....   %76x1 double type

我怀疑这与全局变量有任何关系。你应该有一些某种数据，而不是.....的 5 个点，以使任何代码正常工作（因为上面的注释指示"从...中提取数据"）。如果它不属于同一类，我也会尝试以注释建议的形式制作这些数据，即使其double和大小

。

例如，尝试查看将这些行替换为此代码是否可以使您的代码正常工作

ModelInfo.X=rand(76,2)   %76x2 double type
ModelInfo.y=rand(76,1)   %76x1 double type