| 查看: 4194 | 回复: 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 |
» 猜你喜欢
表哥与省会女结婚,父母去帮带孩子被省会女气回家生重病了
已经有12人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有14人回复
江汉大学解明教授课题组招博士研究生/博士后
已经有3人回复
AI 太可怕了,写基金时,提出想法,直接生成的文字比自己想得深远,还有科学性
已经有11人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有11人回复
铸剑人
金虫 (正式写手)
- 应助: 80 (初中生)
- 金币: 1552.6
- 散金: 406
- 红花: 41
- 帖子: 736
- 在线: 159.5小时
- 虫号: 2517634
- 注册: 2013-06-22
- 性别: GG
- 专业: 非线性光学与量子光学
【答案】应助回帖
| Matlab tool: Nonlinear Least Squares (Curve Fitting)--Fit a Model to Complex-Valued Dat-- which can solve your problem by using two function that lsqcurvefit(Solve nonlinear curve-fitting (data-fitting) problems in least-squares sense)and lsqnonlin(Solve nonlinear least-squares (nonlinear data-fitting) problems) ,addres:http://cn.mathworks.com/help/opt ... o-complex-data.html |

16楼2016-05-18 20:49:08
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