matlab解发酵matlab动力学仿真实例方程

来源：蜘蛛抓取(WebSpider) 时间：2016-11-24 00:46 标签： 发酵动力学
请问谁知道用matlab解二阶动力学方程嘛,求指导
如果你想求时域内的结果推荐用状态空间的方法求解,当然也可以采用龙哥库塔法求解.如果是频域,则需要通过傅里叶变换,则需要有这方面的基础知识.　　用状态空间方法求解时域响应,请看实例文章“基于Matlab_Simulink的多自由度机械振动系统仿真”,百度一下就能下载找到,介绍很详细.　　例如这个例子,如下图：　　启动Matlab／Simulink之后,新建一模型文件打开一空白的编辑窗口．采用Continuous库中的State—Space模块进行建模,选择信号源模块组中的正弦信号模块为输入信号.搭建模型如下：　　状态方程系数ABCD参数如下：　　m1=100;m2=100;m3=100;m4=100;　　k1=50;k2=50;k3=50;k4=50;　　C1=100;C2=150;C3=100;C4=150;　　M=[m1 0 0 0;0 m2 0 0;0 0 m3 0;0 0 0 m4];　　C=[C1+C2 -C2 0 0;-C2 C2+C3 -C3 0;0 -C3 C3+C4 -C4;0 0 -C4 C4];　　K=[k1+k2 -k2 0 0;-k2 k2+k3 -k3 0;0 -k3 k3+k4 -k4;0 0 -k4 k4];　　A=cat(1,cat(2,zeros(4,4),eye(4)),cat(2,-inv(M)*K,-inv(M)*C));　　G=eye(4);B=cat(1,zeros(4,4),-inv(M)*G);　　C0=cat(2,eye(4),zeros(4,4));D=zeros(4,4);　　双击正弦输入模块,在弹出的对话框中输入：振幅设置为向量[10 15 20 25],输入频率为10,其余采用默认值．双击State—space模块,设置参数：A为A、B为B、C为Co、D为D．初始条件设为0．&&&&仿真就可以得到x1至x4&&&&希望对你有帮助
谢谢哈，我就是用ode45求解的，但是结果是nan，不知道哪里错了，要不下次麻烦你帮我看一下程序吧
可以，你再向我提一个问题，把你有问题的代码贴上来。给你查一下
哦哦，我现在用手机登陆的，电脑不在身边，下次发给你看一下可以么
好的。没问题
为您推荐：
其他类似问题
扫描下载二维码小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
求matlab编程求动力学方程的方法
求动力学方程，动力学模型为-dC/dt=kC^n*H2^m
C-反应物浓度，H2-氢气浓度，n、m为反应级数
下面是小弟编的程序，用积分法求解，不知道那错了，求解。
万分感谢！
改了一下，还是有问题
好棒，能否把matlab程序发我一份，小弟想学一学，万分感谢！
能否指正下我编的程序，后面有很多数据要处理，谢谢
dingd用的是1stopt软件，不是MATLAB。
如果你想找人改你的MATLAB代码，请把代码以文本的形式发上来，而不是图片格式。
这是matlab文件
哦哦，谢谢！我用matlab不知道出错在哪，文件已上传
请问去哪可以下到1stopt软件？这个软件我在网上搜了好久没有安装包
你是做微分方程参数拟合，得用这个版本：http://emuch.net/bbs/viewthread.php?tid=8050245
1stOpt软件常见问题与解答：http://emuch.net/bbs/viewthread.php?tid=9176191
好的，谢谢！请问我这个问题用matlab好解决吗
不过这种问题，MATLAB效果不一定有1stOpt好。
能否帮忙看一下我的matlab程序，我做了俩，一个是积分法求方程（文件已上传），一个是微分法，微分法比较简单，但是拟合效果太差，积分法我不会，刚用matlab，有些东西不是很理解，万分感谢
你应该参照的是黄华江那本书上的代码。
可成功运行的MATLAB代码，可参考：http://emuch.net/bbs/viewthread.php?tid=6425538
是啊，我用您推荐的1stopt，怎么算不出来，您看问题出在哪？麻烦老师了
文本形式的代码
哦哦，是我复制的问题吗，我直接粘过去的，不是输的
Title &Type your title here&;
Parameters b1,b2,b3;
Variable t,CA,H2;
ODEFunction CA'=-b1*CA^b2*H2^b3;
0.2& && &&&0.291& && &&&0.023
0.3& && &&&0.186& && &&&0.015
0.4& && &&&0.134& && &&&0.015
0.5& && &&&0.112& && &&&0.008
0.6& && &&&0.101& && &&&0.007
请问您用的1stopt是哪个版本，在网上我下载了一个，按照下面的输入算不出来，是什么原因呢，谢谢！
Title &Type your title here&;
Parameters b1,b2,b3;
Variable t,CA,H2;
ODEFunction CA'=-b1*CA^b2*H2^b3;
0.2& && &&&0.291& && &&&0.023
0.3& && &&&0.186& && &&&0.015
0.4& && &&&0.134& && &&&0.015
0.5& && &&&0.112& && &&&0.008
0.6& && &&&0.101& && &&&0.007
是的吧，用的Auto2Fit,依然是解决不了,发愁
变量CA和H2肯定是有函数关系的，通过关系式变化，消掉一个变量，不然原方程就不是常微分，而是偏微分了。
如果CA和H2都要保留，那就得构成常微分方程组，就是说此时还少了H2' 的表达式。
这个得购买吧，网上可以下到6.0的吗，不然只能求助matlab了
跟用什么软件没关系。
CA H2 和 时间 t的函数关系，必须得先搞清楚。
哦哦，谢谢！就是求解微分方程必须得把各变量关系搞清楚吧，不过H2只是个反应条件，没啥具体的关系。我看的许多人先把方程左边的微分数值求出来，再求右边的参数，这样就不是求解微分方程了，不用知道Ca和H2的关系了，不过这样精度不高啊。
网上买的书，黄华江的《实用化工计算机模拟》，带电子案例，淘宝有售
您好，我下载的1stopt为啥没法用啊，您的是花钱买的吗
Title &HDS1&;
Parameter k,n;
Variable t,CA;
ODEFunction CA'=-k*CA^n;
0.4& && & 1.04936
0.5& && & 0.88105
0.66667& &0.83847
1& && && &0.80326
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
matlab求解动力学 方程中的参数
实验数据位（t,c）=(0,6.94),(2,6.4),(4,6.3),(8,6.2),(48,6.1),动力学方程为：dc/dt=-k*(c0-c)^(2/3)*(c-c~)^p，其中，c0指初始浓度（6.9），c~指平衡浓度（6.1）。
如何求方程中的k与p.文献中给出的方法是先积分求出一个含有k、p的方程，然后拟合求出。如何写程序呢？？谢谢大家帮忙了
你的思路是对的，关键是微分方程比较复杂，应先积分后去对数吧？
不需要积分。
你的目的是求 K和 P。
两边取log后， log(a)=log(K)+ (2/3)long(c0-c)+P log(c-c~)
令 y=log(a)-(2/3) log(C0-C); x=log(C-C~);
就是一般线性方程 y = m* x + b
其中 m=P; b=log(K)
如果数据点是近似在一条直线上不就求出m、b了？
你好，那个a=/(t-t0) 或者= c(t+1)-c（t）,这个是因为啥？原方程两边取对数后方程左侧是dc/dt,这个怎么取对数？
我明白你的意思了。但是dc/dt可以用两个相邻点的纵坐标差减去横坐标差来代替？这样岂不是误差很大啊？有没有积分可以算出来的，
这样得出的结果和别人积分所得的差别很大，拟合出来的相关系数相差很大。呵呵，不过也是一种方法，呵呵，谢谢了。
对你&&这个不懂哎，1stOpt是什么？是matlab里的一种算法吗？？
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研您的位置： &
巧解机器人动力学方程——MATLAB语言应用一例
优质期刊推荐

我要回帖

更多关于 发酵动力学 的文章

 

随机推荐