24小时热门版块排行榜    

查看: 4238  |  回复: 15

meifangh

新虫 (小有名气)

[求助] 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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenyuwen

金虫 (正式写手)

为什么不把实部虚部分开,然后去拟合?

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

2楼2016-05-14 23:33:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meifangh

新虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by chenyuwen at 2016-05-14 23:33:48
为什么不把实部虚部分开,然后去拟合?

大神能不能给个具体的思路

发自小木虫Android客户端
3楼2016-05-15 01:59:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenyuwen

金虫 (正式写手)

那你最好能把函数的具体形式写出来,用类似公式编辑器

发自小木虫Android客户端
4楼2016-05-15 09:24:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
试试1stOpt吧,可以直接进行复数拟和,夜不需要猜初值。

[ 发自手机版 http://muchong.com/3g ]
5楼2016-05-15 13:44:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meifangh

新虫 (小有名气)

引用回帖:
5楼: Originally posted by dingd at 2016-05-15 13:44:32
试试1stOpt吧,可以直接进行复数拟和,夜不需要猜初值。

1stopt高端版本没有啊,好像旧版本没有!!!复数拟合
6楼2016-05-15 22:17:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meifangh

新虫 (小有名气)

引用回帖:
4楼: Originally posted by chenyuwen at 2016-05-15 09:24:42
那你最好能把函数的具体形式写出来,用类似公式编辑器

就靠你了啊!我折腾了近一个月了,一直没有思路。
我手头已知数据有OT,DO的测试数据,复数的。
S0,T0,a,都是已知的。
拟合参数为tv,tc,phi_a,phi_v,phi_c必须拟合出来是实数。
大神你忙的话就给个思路吧..真的不是很多拟合这东西

http://muchongimg.xmcimg.com/data/b ... _1463322789_412.png
matlab 拟合一个复数函数的系数。求思路,感觉没法解决
QQ截图20160515223016.png

7楼2016-05-15 22:35:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenyuwen

金虫 (正式写手)

引用回帖:
7楼: Originally posted by meifangh at 2016-05-15 22:35:43
就靠你了啊!我折腾了近一个月了,一直没有思路。
我手头已知数据有OT,DO的测试数据,复数的。
S0,T0,a,都是已知的。
拟合参数为tv,tc,phi_a,phi_v,phi_c必须拟合出来是实数。
大神你忙的话就给个思路吧.. ...

将你所有公式中的复数展开,分成虚部和
实部,这样一个复数公式就变成两个实数公式,然后对每个实数公式根据你的数据分别拟合,,

发自小木虫Android客户端
8楼2016-05-15 23:45:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

最怜宵舞

木虫 (正式写手)

9楼2016-05-15 23:54:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meifangh

新虫 (小有名气)

引用回帖:
8楼: Originally posted by chenyuwen at 2016-05-15 23:45:50
将你所有公式中的复数展开,分成虚部和
实部,这样一个复数公式就变成两个实数公式,然后对每个实数公式根据你的数据分别拟合,,
...

1.就是分开拟合实部和虚部,怎么看待呢?实部拟合出来的,跟虚部拟合出来的数据,之间什么关系?
2.拟合出来的数据怎么组合?
3.分别实部,虚部拟合,拟合出来的系数值不一样...怎么办?
10楼2016-05-16 14:32:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 meifangh 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286求调剂 +3 Faune 2026-03-21 3/150 2026-03-21 11:14 by ms629
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 307求调剂 +10 冷笙123 2026-03-17 10/500 2026-03-21 01:54 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 一志愿北京化工大学0703化学318分,有科研经历,求调剂 +4 一瓶苯甲酸 2026-03-14 4/200 2026-03-20 20:36 by fen_rao
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
信息提示
请填处理意见