| 查看: 4122 | 回复: 15 | ||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | ||
[求助]
matlab 拟合一个复数函数的系数。求思路,感觉没法解决 已有2人参与
|
||
|
function [OT DT] = hemomodel(omega, tc, tv, phi_a, phi_c, phi_v) S0 = 0.98; T0 = 2.3; alpha = 0.8; O = S0*T0/(phi_a+phi_c+phi_v) ... *(phi_a + phi_c*(exp(-alpha*tc+1j*omega*tc)-1)./(-alpha*tc+1j*omega*tc) ... + phi_v*(exp(1j*omega*tv)-1)./(1j*omega*tv).*exp(-alpha*tc+1j*omega*tc)); T = T0/(phi_a+phi_c+phi_v) ... *(phi_a + phi_c*(exp(1j*omega*tc)-1)./(1j*omega*tc) ... + phi_v*(exp(1j*omega*tv)-1)./(1j*omega*tv).*exp(1j*omega*tc)); D = T - O; OT = O./T; DT = D./T; (初始值的话。可以尝试这个。tc=0.5,tv=1,phi_a=0.005,phi_v=0.005,phi_c=0.015) 上面的是一个复数函数。我有实验数据OT,DT。现在需要的是通过已知的OT,DT数据进行拟合函数中系数的tc,tv,phi_a,phi_v,phi_c. 拟合出来的系数必须是实数。 我是从一片文献中看到的论文,他们是使用matlab中的函数lsqcurvefit进行拟合的。我本人用此函数拟合的时候遇到了2个问题。 1.lsqcurvefit函数对复数函数拟合,出现的复数。此函数的系数是有物理意义的,不可能是复数,我不知道怎么调试... 2.lsqcurvefit函数一次只能拟合一条函数。而OT,DT显然是相互关联的,如果分开你和OT,DT的话,拟合的参数结果不一样。 小弟拜托大家了,我自己拟合了很久了.还是不行啊!!@laosam280 |
» 猜你喜欢
求个博导看看
已经有16人回复
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗?
已经有8人回复
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
上海工程技术大学【激光智能制造】课题组招收硕士
已经有5人回复
求助院士们,这个如何合成呀
已经有4人回复
临港实验室与上科大联培博士招生1名
已经有9人回复
需要合成515-64-0,50g,能接单的留言
已经有4人回复
自荐读博
已经有4人回复
写了一篇“相变储能技术在冷库中应用”的论文,论文内容以实验为主,投什么期刊合适?
已经有6人回复
带资进组求博导收留
已经有10人回复
chenyuwen
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1736.1
- 散金: 110
- 红花: 6
- 帖子: 429
- 在线: 120.5小时
- 虫号: 3182070
- 注册: 2014-05-05
- 性别: GG
- 专业: 光学
15楼2016-05-16 18:30:32
chenyuwen
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1736.1
- 散金: 110
- 红花: 6
- 帖子: 429
- 在线: 120.5小时
- 虫号: 3182070
- 注册: 2014-05-05
- 性别: GG
- 专业: 光学
2楼2016-05-14 23:33:48
3楼2016-05-15 01:59:52
chenyuwen
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1736.1
- 散金: 110
- 红花: 6
- 帖子: 429
- 在线: 120.5小时
- 虫号: 3182070
- 注册: 2014-05-05
- 性别: GG
- 专业: 光学
4楼2016-05-15 09:24:42







回复此楼
meifangh