应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab拟合模型参数
51/1返回列表
查看: 2137  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lxyy

银虫 (小有名气)

[求助] matlab拟合模型参数

CA，CB,CC,CW,CR每个时间点t下都有具体数值，要怎么用matlab编程啊
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 新建 Microsoft Word 文档.doc
    • 2012-03-23 11:05:51, 20.5 K

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-03-23 11:07:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lxyy

银虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2012-03-23 15:34:22:
这种微分方程拟合不是Matlab的强项，用1stOpt试试！

可是我不怎么会
回复此楼
4楼2012-03-23 19:15:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
lxyy: 金币+10 2012-03-23 19:15:43
这种微分方程拟合不是Matlab的强项，用1stOpt试试！
赞一下 回复此楼
2楼2012-03-23 15:34:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
dbb627: 金币+2, 谢谢应助！ 2012-03-23 17:16:55
lxyy: 金币+60, 如果可以，剩下的金币会再给你，呵呵 2012-03-23 19:17:04
function IsothermTR
clear all
clc
global keq1 keq2 t0 c0

keq1 = ?;         % 平衡常数k1
keq2 = ?;   % 平衡常数k2

t0 = [0.5, 1.0, 1.5, 2.0, 2.5, 3.0];
c0 = [?, ? ?]        % A C W 组成的初始浓度
ci=[0.184         0.2443        0.1348;
    0.4457        0.3214        0.1733;
    0.7023        0.4183        0.2713;
    0.7758        0.4789        0.1056;
    0.3365        0.7035        0.1951;
    0.7946        0.8864        0.2235];
% cA cC cR 对应时间数值

k0=[2 2 3];   % k1+, k2, k3+ 初始向值
lb = [0 0 0];
ub = [1000  1000 500];  % 上下限

% 使用函数lsqnonlin()进行参数估计
[k,resnorm,residual,exitflag] = lsqnonlin(@ObjFunc,k0,lb,ub,[],ci);      

k1plus=k(1);
k1minus=keq1/k1plus;
k2=k(2);
k3plus=k(3);
k3minus=keq1/k3plus;

% ------------------------------------------------------------------
function f = ObjFunc(k,ci)    % 目标函数
global t0  c0

[t,c_cal] = ode45(@Euqations,t0,c0,[],k);
f= c_cal-ci;

% ------------------------------------------------------------------
function dcdt = Euqations(t,c,k)
global keq1 keq2

k1plus=k(1);
k1minus=keq1/k1plus;
k2=k(2);
k3plus=k(3);
k3minus=keq1/k3plus;

cA=c(1);
cC=c(2);
cR=c(3);
cB=?;                 % cB物料衡算表达
cW=?;                % cW物料衡算表达

dcAdt = k1plus*cA*cB-k1minus*cC+k2*cA*cW;
dccdt = -k1plus*cA*cB+k1minus*cC+k3plus*cB*cR-k3minus*cC*cW;
dcRdt = k2*cA*cW-k3plus*cB*cR+k3minus*cC*cW;
dcdt = [dcAdt; dcCdt;dcRdt];
赞一下 回复此楼
3楼2012-03-23 17:13:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lxyy

银虫 (小有名气)
引用回帖:
3楼: Originally posted by change0618 at 2012-03-23 17:13:43:
function IsothermTR
clear all
clc
global keq1 keq2 t0 c0

keq1 = ?;         % 平衡常数k1
keq2 = ?;   % 平衡常数k2

t0 = ;
c0 =         % A C W 组成的初始浓度
ci=;
% cA cC cR 对应时间数值

k0=;   % k ...

我试试看，谢谢
回复此楼
5楼2012-03-23 19:16:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 317求调剂 +5 申子申申 2026-03-19 10/500 2026-03-20 15:58 by 蔡诚
[考研] 广西大学家禽遗传育种课题组2026年硕士招生（接收计算机专业调剂） +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 281求调剂（0805） +14 烟汐忆海 2026-03-16 25/1250 2026-03-20 15:47 by yuncha
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +21 rare12345 2026-03-18 21/1050 2026-03-20 14:31 by 无懈可击111
[考研] 295材料求调剂，一志愿武汉理工085601专硕 +4 Charlieyq 2026-03-19 4/200 2026-03-20 14:26 by 无懈可击111
[考研] 材料学硕297已过四六级求调剂推荐 +6 adaie 2026-03-19 6/300 2026-03-20 14:12 by 星空星月
[考研] 本人考085602 化学工程 专硕 +18 不知道叫什么！ 2026-03-15 20/1000 2026-03-20 13:52 by danney002
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 085700资源与环境308求调剂 +10 墨墨漠 2026-03-18 11/550 2026-03-20 11:57 by 学员8dgXkO
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 一志愿华中科技大学，080502，354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 330求调剂 +3 小材化本科 2026-03-18 3/150 2026-03-18 21:55 by 无懈可击111
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
信息提示
关闭
请填处理意见
关闭