求助——生化反应动力学建模

huangzishan

一道作业:请高手解决下,或者指点迷津也行.
      详细内容请下载附件>

[ 本帖最后由 huright 于 2007-6-6 13:41 编辑 ]

huangzishan

我的 QQ:85388057   有什么不清楚的地方可以找我聊下.
      大家帮帮忙啊!!!

sbstoneface

应该可以用神经网络吧

huangzishan

不要下了就跑了啊,解决下啊,各位大哥:'(
   吾将感激涕淋

longe329

呵呵，有没有人做出来啊，做出来了希望代码放出来哦，我也是做优化的，比楼主这个问题还复杂，也很期待这个问题的答案哦，
谢谢了，希望大人能够做出来
顺便帮楼主顶了

xjzuo

这种题目自己不先动手做做，恐怕即使别人知道，也没有那么多时间去帮你完成。

rocwoods

干脆出点血得了。舍不得孩子套不到狼。大家的时间也是时间啊

huangzishan

老大啊,我都研究好久了,就是连基本思路都没有啊!只知道大概要怎么进行!
     我们老师什么都没说啊!他自己思路都有点混乱!说的乱七八糟的,不会从我们的观点深入!!

huangzishan

问题老没得到解决,有没有诚心一点的,不要老是利益为先,动不动谈价钱!
      浮一下!

eight

原帖由 huangzishan 于 2007-6-7 22:44 发表

                               
登录/注册后可看大图

问题老没得到解决,有没有诚心一点的,不要老是利益为先,动不动谈价钱!
      浮一下!

切勿动怒，建议先看看本版各置顶贴，特别是“动手者”、“聚宝盆”等

huangzishan

function df=m1(t,f)
%OBJFUN Summary of this function goes here
%   Detailed explanation goes here
% 全局变量
global Um;
global Ks;
global Yxs;
global m;
% 转换因变量
% f(1):Cx,f(2):Cs
% dCx/dt :df(1) ,dCs/dt :df(2) ,
% ***************
% 建立 2*1 阶矩阵
df=zeros(2,1);
df(1)=Um*f(1)*f(2)/(Ks+f(2));
df(2)=(-1)*df(1)/Yxs+m*f(1);
function data=objfun(k)
%OBJFUN Summary of this function goes here
%   Detailed explanation goes here
% 全局参数
global Um;
global Ks;
global Yxs;
global m;
Um=k(1);
Ks=k(2);
Yxs=k(3);
m=k(4);
% 时间
tspan=[3 11 19 27 35 43 51 59 67 75 83];
% 瞬时细胞浓度实验值
Cx=[0.020,0.029,0.041,0.059,0.084,0.12,0.18,0.27,0.41,0.66,1.23];
% 瞬时底物浓度实验值
Cs=[30.00,29.81,29.55,29.16,28.61,27.80,26.61,24.82,22.04,17.42,8.35];
% 瞬时浓度值计算
% Cx,Cs初值
f0=[0.020 30.00];   
%options
[t,f]=ode45(@m1,tspan,f0);
% 计算值与实验值的残差平方和 data
test=[Cx;Cs]';
temp1=f-test;
temp2=temp1./test;
data=sum(sum(temp2.^2));
function y=ansfun(k)
% 全局参数
% global Um;
% global Ks;
% global m;
% global Yxs;
% 全局参数的数值范围
Um=k(1);
Ks=k(2);
Yxs=k(3);
m=k(4);
%
format long e;
% 求解过程参数设置
% 初值以及精度
Lb=[1e0 1e2 1e-4 1e-4];
% 范围设置
Ub=[Um Ks Yxs m];
% 参数设置
options=optimset('largescale','off','display','final','tolx',1e-8,'tolfun',1e-8,'MaxFunEvals',6000);
% 求解过程
[p,fval,exitflag,output]=fmincon(@objfun,Lb,[],[],[],[],Lb,Ub,[],options);
y=p;
disp('fval=');
disp(fval);
disp('exitflag=');
disp(exitflag);
disp('output=');
disp(output);
disp('revelent error=');
error=gridfun(p);
disp(error);
         main()
k=[1e8 1e8 0.3 4];
ansfun(k);

huangzishan

我过几天还会用其他的方法的!!!
     大家期待吧!