24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 4314  |  回复: 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

金虫 (正式写手)

引用回帖:
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的回帖
查看全部 16 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302分求调剂 +10 凡语祈愿 2026-04-08 11/550 2026-04-14 16:50 by jiangguiquan11
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[考研] 279学硕食品专业求调剂院校 20+6 孤独的狼爱吃羊 2026-04-12 28/1400 2026-04-14 15:44 by zs92450
[考研] 071000生物学调剂求助 +18 zzzzwww 2026-04-09 21/1050 2026-04-14 15:39 by zs92450
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 求调剂 +16 雪逢冬 2026-04-10 16/800 2026-04-14 14:27 by 逆水乘风
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
[考研] 一志愿厦大生物学332求调剂 +11 池池池池池池 2026-04-08 11/550 2026-04-13 14:10 by 科研论
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[硕博家园] 新一代电子信息294求调剂 不挑学校 +7 Ytyt11 2026-04-09 8/400 2026-04-12 16:57 by ajpv风雷
[考研] 电气专硕320求调剂 +6 小麻子111 2026-04-10 6/300 2026-04-12 10:54 by lemon6009
[考研] 求调剂,262机械专硕 +8 嗯yyl 2026-04-08 8/400 2026-04-12 02:31 by 秋豆菜芽
[考研] 22408调剂315分 +3 zhuangyan123 2026-04-09 3/150 2026-04-12 00:25 by 蓝云思雨
[考研] 求调剂 +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文 +13 Emmy~ 2026-04-09 13/650 2026-04-09 14:47 by only周
[考研] 招收有机化学、化工,药学,食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
信息提示
请填处理意见