求助——生化反应动力学建模
一道作业:请高手解决下,或者指点迷津也行.详细内容请下载附件>
[ 本帖最后由 huright 于 2007-6-6 13:41 编辑 ]
回复 #2 huangzishan 的帖子
我的 QQ:85388057 有什么不清楚的地方可以找我聊下.大家帮帮忙啊!!!
RE:一道极难的题目
应该可以用神经网络吧我再说下....
不要下了就跑了啊,解决下啊,各位大哥:'(吾将感激涕淋 呵呵，有没有人做出来啊，做出来了希望代码放出来哦，我也是做优化的，比楼主这个问题还复杂，也很期待这个问题的答案哦，
谢谢了，希望大人能够做出来
顺便帮楼主顶了 这种题目自己不先动手做做，恐怕即使别人知道，也没有那么多时间去帮你完成。 干脆出点血得了。舍不得孩子套不到狼。大家的时间也是时间啊
回复 #6 xjzuo 的帖子
老大啊,我都研究好久了,就是连基本思路都没有啊!只知道大概要怎么进行!我们老师什么都没说啊!他自己思路都有点混乱!说的乱七八糟的,不会从我们的观点深入!!
回复 #1 huangzishan 的帖子
问题老没得到解决,有没有诚心一点的,不要老是利益为先,动不动谈价钱!浮一下! 原帖由 huangzishan 于 2007-6-7 22:44 发表 http://www.chinavib.com/forum/images/common/back.gif
问题老没得到解决,有没有诚心一点的,不要老是利益为先,动不动谈价钱!
浮一下!
切勿动怒，建议先看看本版各置顶贴，特别是“动手者”、“聚宝盆”等
我自己来!求人不如求己!!!
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=;
% 瞬时细胞浓度实验值
Cx=;
% 瞬时底物浓度实验值
Cs=;
% 瞬时浓度值计算
% Cx,Cs初值
f0=;
%options
=ode45(@m1,tspan,f0);
% 计算值与实验值的残差平方和 data
test=';
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=;
% 范围设置
Ub=;
% 参数设置
options=optimset('largescale','off','display','final','tolx',1e-8,'tolfun',1e-8,'MaxFunEvals',6000);
% 求解过程
=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=;
ansfun(k);
回复 #12 huangzishan 的帖子
我过几天还会用其他的方法的!!!大家期待吧!
页:
[1]