24小时热门版块排行榜    

查看: 4233  |  回复: 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的回帖

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

chenyuwen

金虫 (正式写手)

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

发自小木虫Android客户端

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

2楼2016-05-14 23:33:48
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 初始318分求调剂(有工作经验) +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 278求调剂 +6 烟火先于春 2026-03-17 6/300 2026-03-21 01:57 by JourneyLucky
[考研] 307求调剂 +10 冷笙123 2026-03-17 10/500 2026-03-21 01:54 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 一志愿华中农业071010,总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 260求调剂 +3 朱芷琳 2026-03-20 3/150 2026-03-20 20:35 by 学员8dgXkO
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
信息提示
请填处理意见