关于matlab仿真的一个程序这个错在哪里挖,function pre = GM11(x)%x,建模数据序列x = [x(1),x(2),...,x(n)]% pre,由GM11模型得到的预测值,即由x估计出的x(n+1)n = length(x);x1 = cumsum(x);z1 = (x1(1:n-1)+x1(2:n))/2;Y = x(2:l

Question

关于matlab仿真的一个程序
这个错在哪里挖,
function pre = GM11(x)
%x,建模数据序列x = [x(1),x(2),...,x(n)]
% pre,由GM11模型得到的预测值,即由x估计出的x(n+1)
n = length(x);
x1 = cumsum(x);
z1 = (x1(1:n-1)+x1(2:n))/2;
Y = x(2:length(x))';
B = [-z1',ones(n-1,1)];
A = (B'*B)\B'*Y;%GM(1,1)模型的a,b参数组成的向量.
beta = A(2)/(1+.5*A(1));
alpha = A(1)/(1+.5*A(1));
%初始序列x0(k)模拟后的随时间t变化的序列
f = @(t)(beta-(alpha)*x(1))*exp(-(A(1)*(t-2)));
pre = f(n+1);
主程序为GNP=[326.82;410.22;455.86;500.72;598.79;709.00;863.23;1084.33;1395.12;1616.03;1810.49;2011.77;2174.97;2479.32];
>> lGNP=length(GNP);
>> pre=arrayfun(@(k) GM11(GNP(lGNP-k:end)),3,10)
为什么运行不出来啊,郁闷死.希望大虾指点啊,应该咋改啊,
最后一行打错了,应该是pre=arrayfun(@(k) GM11(GNP(lGNP-k:end)),3:10)

Answer 1

请问pre=arrayfun(@(k)是什么意思?
怎么找不到arrayfun()函数?