yfz688 发表于 2012-8-29 21:51

求大神解答。谢谢

x0=;
x1=ones(1,10);x11=ones(1,10);x01=ones(1,10);
s=0;t1=0:1:9;t=1:1:20;
for i=1:10
    s=s+x0(i);
  x1(i)=s;
end
B=ones(9,2);
for i=1:9
    B(i)=-0.5*(x1(i+1)+x1(i));
end
y=x0(2:end)';
u=inv(B'*B)*B'*y;
for i=0:9
    x11(i+1)=(x1(1)-u(2)/u(1))*exp(-u(1)*i)+u(1)/u(2);
end
x01(1)=x0(1);
for i=2:10
    x01(i)=x11(i)-x11(i-1);
end
% jingdujianyanyuyuce
x=0;
for i=1:10
    E(i)=x0(i)-x01(i);
    e(i)=E(i)/x0(i);
    x=x+x0(i);
end
E
e
x=0.1*x
s=0;s1=0;
for i=1:10
    s1=s1+(x0(i)-x)^2;
    s=s+E(i);
end
s1=sqrt(0.1*s1)
E1=1/9*(s-E(1))
s=0;
for i=2:10
    s=s+(E(i)-E1)^2;
end
s2=sqrt(1/9*s)
c=s2/s1
% yuceweilaishinian
for i=10:19
    x11(i+1)=(x1(1)-u(2)/u(1))*exp(-u(1)*i)+u(1)/u(2);
end
for i=11:20
    x01(i)=x11(i)-x11(i-1);
end
x01
plot(t,x01,'--')
hold on
plot(t1,x0,'o',t1,x01(1:10),'-')
legend('yuceshuju','yuanshuju','niheshuju',2)
输出是找不到指定模块。请问是什么原因

hzs2012 发表于 2012-8-29 23:33

在我的R2011a上可以运行。。。。。。

秋leaves 发表于 2012-8-30 09:04

能运算啊……你保存路径的问题

darker50 发表于 2012-8-30 09:23

   可能是版本的问题,换个版本试试!

yfz688 发表于 2012-8-30 09:46

秋leaves 发表于 2012-8-30 09:04 static/image/common/back.gif
能运算啊……你保存路径的问题

好吧。谢了

yfz688 发表于 2012-8-30 09:47

darker50 发表于 2012-8-30 09:23 static/image/common/back.gif
可能是版本的问题,换个版本试试!

好吧。我换个版本
页: [1]
查看完整版本: 求大神解答。谢谢