本人新手小白,因论文需求
使用下面的代码求解
clc
clear
Ex=[92.5,80,70,60,27.5;92.5,80,70,60,27.5;95,85,75,65,30;92.5,80,70,60,27.5;92.5,80,70,60,27.5;92.5,80,70,60,27.5;92.5,80,70,60,27.5;95,85,75,65,30]
En=[2.5,1.67,1.67,1.67,9;2.5,1.67,1.67,1.67,9;1.67,1.67,1.67,1.67,10;2.5,1.67,1.67,1.67,9;
2.5,1.67,1.67,1.67,9;2.5,1.67,1.67,1.67,9;2.5,1.67,1.67,1.67,9;1.67,1.67,1.67,1.67,10]
He=[0.4,0.3,0.4,0.4,0.1;0.3,0.3,0.4,0.2,0.1;0.5,0.2,0.1,0.3,0.1;0.3,0.3,0.4,0.4,0.2;0.3,0.3,0.4,0.4,0.2;0.4,0.4,0.5,0.3,0.4;0.4,0.4,0.4,0.3,0.1;0.4,0.3,0.1,0.5,0.1;0.5,0.2,0.1,0.3,0.4]
Xi=[76,76,76,76,76;83,83,83,83,83;68,68,68,68,68;82,82,82,82,82;67,67,67,67,67;76,76,76,76,76;69,69,69,69,69;79,79,79,79,79]
W=[0.0625,0.0296,0.0330,0.0231,0.0200,0.0628,0.0645,0.0330]
Z=zeros(1,6)
for m=1:8
ui0=zeros(1,1000)'
for n=1:5
En_=En(m,n)
Ex_=Ex(m,n)
He_=He(m,n)
xi_=xi(m,n)
Eni=normrnd(En_,He_,1,1000)'
ui1=exp(-((xi_.*ones(1,1000)'-Ex_.*ones(1,1000)').^2)./(2.*Eni.^2))
ui0=[ui0,ui1]
end
Z=[Z;ui0]
end
Z(:,1:1)=[]
Z(1:1,:)=[]
Z_1=zeros(1,5)
for j=i:1000:8000
Z_1=Z(j:j,:)
Z0=[Z0;Z1]
end
Z0(1:1,:)=[]
D0=W*Z0
D=[D;D0]
end
D(1:1,:)=[]
R=ones(1,1)
For i=1:1000
sum=0
sum2=0
For j=1:5
sum1=D(i,j)*F(j)
sum=sum+sum1
sum2=sum2+D(i,j)
end
r=sum/sum2
R={R,r]
end
R(1)=[]
Erx=mean(R)
Ern=std(R,0)
g=Ern/Erx
disp('关联度矩阵Z');disp(Z_1);
disp('综合评判向量D');disp(D(1:1,:));
disp('多次运算得到的综合评判分数的期望值Erx');disp(Erx);
disp('多次运算得到的综合评判分数的熵Ern');disp(Ern);
disp('置信度因子g');disp(g);
出现未定义函数或变量“x”
求各位大佬指点一下,以及看一下是否还存在其他的问题,感激不尽!