| 查看: 4290 | 回复: 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 |
回复此楼» 猜你喜欢
085410 273求调剂
已经有9人回复
-
297求调剂
已经有16人回复
-
291求调剂
已经有11人回复
-
计算机22408 281分，求调剂
已经有7人回复
-
316求调剂
已经有5人回复
-
一志愿郑州大学 22408 305分求调剂
已经有5人回复
-
299求调剂
已经有8人回复
-
312求调剂
已经有6人回复
-
280求调剂
已经有10人回复
-
280求调剂
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教matlab非线性拟合结果不唯一的问题
已经有4人回复
- 请问水射流的波动，可以用什么函数形式拟合吗？
已经有0人回复
- 如何用MATLAB拟合这样一个四次曲线？
已经有1人回复
- 有关MATLAB的问题...
已经有110人回复
- matlab 拟合问题求助
已经有6人回复
- 求助Matlab拟合曲线方程
已经有2人回复
- matlab如何在一个图里拟合多条直线
已经有5人回复
- 求助用matlab 多元非线性拟合2
已经有21人回复
- 一组非线性数据拟合
已经有6人回复
- 如何用matlab实现分峰拟合？
已经有18人回复
- 求拟合一组数据
已经有11人回复
- 如何用matlab进行指数拟合【已搜无重复】
已经有2人回复
- 想拟合一条简单曲线，请教一下大家~
已经有6人回复
- matlab中fminbnd的求助
已经有2人回复
- 如何用matlab来拟合一组球坐标数据，方程应该是纺锤体的方程。。
已经有2人回复
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
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2016-05-15 13:44:32
6楼2016-05-15 22:17:27
|
就靠你了啊！我折腾了近一个月了，一直没有思路。 我手头已知数据有OT,DO的测试数据，复数的。 S0，T0，a,都是已知的。 拟合参数为tv，tc，phi_a,phi_v,phi_c必须拟合出来是实数。 大神你忙的话就给个思路吧..真的不是很多拟合这东西 http://muchongimg.xmcimg.com/data/b ... _1463322789_412.png  QQ截图20160515223016.png |
7楼2016-05-15 22:35:43
chenyuwen
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1722.1
- 散金: 110
- 红花: 6
- 帖子: 430
- 在线: 120.8小时
- 虫号: 3182070
- 注册: 2014-05-05
- 性别: GG
- 专业: 光学
8楼2016-05-15 23:45:50
最怜宵舞
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1711.4
- 散金: 114
- 红花: 16
- 沙发: 6
- 帖子: 420
- 在线: 24小时
- 虫号: 4171587
- 注册: 2015-10-25
- 性别: MM
- 专业: 计算数学与科学工程计算
9楼2016-05-15 23:54:34
10楼2016-05-16 14:32:40
5