标题: 求大神解答。谢谢 [打印本页] 作者: yfz688 时间: 2012-8-29 21:51 标题: 求大神解答。谢谢 x0=[124.92,132.04,139.45,150.93,165.13,181.93,196.83,212.38,228.07,241.45];
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
能运算啊……你保存路径的问题
