| 查看: 4528 | 回复: 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 |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
两块石头
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有8人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
chenyuwen
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1722.1
- 散金: 110
- 红花: 6
- 帖子: 430
- 在线: 120.8小时
- 虫号: 3182070
- 注册: 2014-05-05
- 性别: GG
- 专业: 光学
15楼2016-05-16 18:30:32
chenyuwen
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1722.1
- 散金: 110
- 红花: 6
- 帖子: 430
- 在线: 120.8小时
- 虫号: 3182070
- 注册: 2014-05-05
- 性别: GG
- 专业: 光学
meifangh2楼2016-05-14 23:33:48
3楼2016-05-15 01:59:52
chenyuwen
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1722.1
- 散金: 110
- 红花: 6
- 帖子: 430
- 在线: 120.8小时
- 虫号: 3182070
- 注册: 2014-05-05
- 性别: GG
- 专业: 光学
4楼2016-05-15 09:24:42
