johnl
2010-05-30, 13:18
matlab基础太差，求各位高手帮忙解决。
theta=[0.0485 0.0459 0.0720 0.0686 0.1504];
V31=[1.13 1.01 1.07 1.04 1.17];
V32=[1.01 0.94 0.89 0.92 1.05]
z=[0 1];
初值 k0 = [1 1 1 1 1];x0=[0.0565 0.3017 0.5357 0.0564 0.1355]；
需要先求解动力学方程，然后采用最小二乘lsqnonlin估计。
dxdz = ...
(theta*V31*k(1)*x(3))
(theta*V32*k(2)*x(4))
(theta*(k(1)+k(2)+k(3)+k(5))*x(3))
(theta*V31*k(3)*x(3)-theta*k(4)*x(4)
(theta*V32*k(5)*x(3)+theta*V31*k(4)*x(4))]
结果为在反应器出口处测得的结果，即在Z=1处测得的结果
结果(五组)如下：结果与theta(i),V31(i),V32(I)向对应
exp=
[ 1 0.134 0.377 0.162 0.09 0.217
2 0.132 0.390 0.281 0.076 0.168
3 0.142 0.406 0.201 0.086 0.163
4 0.152 0.422 0.243 0.08 0.152
5 0.102 0.329 0.113 0.121 0.286]
我的邮箱[email protected]
theta=[0.0485 0.0459 0.0720 0.0686 0.1504];
V31=[1.13 1.01 1.07 1.04 1.17];
V32=[1.01 0.94 0.89 0.92 1.05]
z=[0 1];
初值 k0 = [1 1 1 1 1];x0=[0.0565 0.3017 0.5357 0.0564 0.1355]；
需要先求解动力学方程，然后采用最小二乘lsqnonlin估计。
dxdz = ...
(theta*V31*k(1)*x(3))
(theta*V32*k(2)*x(4))
(theta*(k(1)+k(2)+k(3)+k(5))*x(3))
(theta*V31*k(3)*x(3)-theta*k(4)*x(4)
(theta*V32*k(5)*x(3)+theta*V31*k(4)*x(4))]
结果为在反应器出口处测得的结果，即在Z=1处测得的结果
结果(五组)如下：结果与theta(i),V31(i),V32(I)向对应
exp=
[ 1 0.134 0.377 0.162 0.09 0.217
2 0.132 0.390 0.281 0.076 0.168
3 0.142 0.406 0.201 0.086 0.163
4 0.152 0.422 0.243 0.08 0.152
5 0.102 0.329 0.113 0.121 0.286]
我的邮箱[email protected]