24小时热门版块排行榜    

Znn3bq.jpeg
查看: 4287  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 271求调剂 +7 2261744733 2026-04-11 8/400 2026-04-11 17:04 by zchychina
[考研] 366求调剂 +7 不知名的小卅 2026-04-11 7/350 2026-04-11 16:45 by PeggyPeng17
[考研] 调剂求助 +6 果然有我 2026-04-11 7/350 2026-04-11 16:22 by 明月此时有
[考研] 材料考研调剂 +22 云木达达 2026-04-11 24/1200 2026-04-11 15:44 by wp06
[考研] 296求调剂 +13 汪!?! 2026-04-08 14/700 2026-04-11 10:51 by 逆水乘风
[考研] 283求调剂 +22 那个噜子 2026-04-09 22/1100 2026-04-11 10:41 by 逆水乘风
[考研] 268分085602化学工程调剂 +27 月照花林。 2026-04-09 27/1350 2026-04-11 10:40 by maddjdld
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 广东省 085601 329分求调剂 +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[考研] 277 数一104,学硕,求调剂 +19 瓶子PZ 2026-04-09 21/1050 2026-04-11 09:39 by 逆水乘风
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 301求调剂 +5 149. 2026-04-10 5/250 2026-04-10 15:45 by 柴小白
[考研] 277求调剂 +19 倪建设 2026-04-06 19/950 2026-04-10 09:24 by guosr9609
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 材料考研求调剂总分280 +30 mkjlz1 2026-04-06 35/1750 2026-04-08 21:25 by cyh—315
[考研] 一志愿吉大化学327求调剂 +12 王王白石 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[考研] 259求调剂 +5 就爱吃土豆呀呀 2026-04-07 5/250 2026-04-07 22:40 by JourneyLucky
[论文投稿] Decision: Revise for Editor还会送审吗 100+3 CccccccccFD 2026-04-04 5/250 2026-04-07 10:58 by 北京莱茵润色
[考研] 材料调剂 +5 小刘同学吖吖 2026-04-06 5/250 2026-04-06 18:34 by sherry_1901
[考研] 298求调剂 +7 manman511 2026-04-05 7/350 2026-04-05 10:29 by 唐沐儿
信息提示
请填处理意见