pmdswjtu 发表于 2012-2-7 22:52

灰色代码

clc,clear
x0=;
n=length(x0);
x1=cumsum(x0)
a_x0=diff(x0);
a_x0=
for i=2:n
z(i)=0.5*(x1(i)+x1(i-1));
end
B=[-x0(2:end)',-z(2:end)',ones(n-1,1)];
Y=a_x0(2:end)';
u=B\Y
x=dsolve('D2x+a1*Dx+a2*x=b','x(0)=c1,x(5)=c2');
x=subs(x,{'a1','a2','b','c1','c2'},{u(1),u(2),u(3),x1(1),x1(5)});
yuce=subs(x,'t',0:n-1);
digits(6),x=vpa(x)
x0_hat=
epsilon=x0-x0_hat
delta=abs(epsilon./x0)

alair006 发表于 2012-2-8 08:19

nuanyang 发表于 2012-2-8 08:34

这是干什么用的？

alair006 发表于 2012-2-8 10:39

467857726de 发表于 2012-8-27 19:56

不错哦  谢谢楼主

空木葬花 发表于 2014-3-9 15:17

非常感谢楼主的福利！
页: [1]
查看完整版本: 灰色代码