| 查看: 4337 | 回复: 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 |
回复此楼» 猜你喜欢
河北省自然基金
已经有3人回复
-
本人42，博士刚毕业，现在找不到工作，怎么办？:(
已经有8人回复
-
没考上92，去双非就完了？其实你还有……
已经有6人回复
-
交大在职博士（哲学、社会学）
已经有5人回复
-
3,4-二羟基苯乙酮如何纯化？
已经有5人回复
-
国基评审
已经有10人回复
-
析晶
已经有5人回复
-
国自然面上和省基金B类撒花
已经有22人回复
-
2026-博士申请
已经有4人回复
-
26级硕士毕业生求博导收留
已经有4人回复
chenyuwen
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1722.1
- 散金: 110
- 红花: 6
- 帖子: 430
- 在线: 120.8小时
- 虫号: 3182070
- 注册: 2014-05-05
- 性别: GG
- 专业: 光学
|
是要展开成你写的形式，，你的有些参数不是只能取实数吗，那就将只能取实数的参数和实数方程放在一起，而虚数则与展开虚数的放在一起，。你在展开的时候要注意，一些实际可以取复数的量必须写成x+y*j的形式， 发自小木虫Android客户端 |
12楼2016-05-16 18:19:34
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