求助帖

X0=[3810 3826 3864 3961 3995 3958 4034 4122 4136 4084 4195 4287 4217 4294 4398 4415 4394 4527 4476 4477 4442 4480 4299 4229 4112 4206 4334 4401 4376 4378 4309 4283 4418 4446 4529 4658 4814 4911 4942 4620 4612 4829 4911 4910 4947 5023 5132 5114 5106 5122 5166 5063 4887 4968 4785 4478 4576 4690 4528 4193 4053 4277 4054 3913 3687 3776 3727 3507 3709 3878 3970 3924 3806 3823 3957 3992 4018 4026 4124 4071 3726 3663 3789 3706 3664 3623 3757 3695 3662 3744 3928 3928 3886 3955 3965];
max1=max(X0)
X00=X0+3*max1;
n=length(X00);
for i=1:n
leijia(i)=sum(X00(1:i));
end
for i=1:n-1
junzhi(i,1)=-(leijia(i)+leijia(i+1))/2;%
end
junzhi=[junzhi ones(n-1,1)];%均值序列
y=X00(2:n);
yn=y';
a=(junzhi'*junzhi)^(-1)*junzhi'*yn%待估参数
max1 =
5166
警告: 矩阵为奇异工作精度。
a =
NaN
NaN
NaN
请问为什么会警告啊？？