24小时热门版块排行榜    

查看: 2886  |  回复: 3

QiChyi

银虫 (初入文坛)

[求助] 求助Matlab lsqcurve函数曲线拟合

本人编写一个利用lsqcurvefit函数拟合曲线的函数,代码如下:
总是提示错误:Optimization terminated: norm of the current step is less than OPTIONS.TolX.
不晓得是怎么回事,请求高人指点!(注:由于种种原因,仅能用Matlab编程,其他软件不采用)
w=2*pi*[0.00099 0.001 0.002 0.004 0.010 0.021544 0.0464 0.100 0.2154 0.4641 1 2.154 4.64 10 20 40 70 110 222 446 1000];
ydata=[5.5981 5.5924 5.3308 5.2180 5.1315 5.0785 5.0289 5 5.1115 5.0825 5.0579 5.0401 5.0175 5.0271 5.0151 5.005 5.0002 4.9950 4.9906 4.9862 4.9875];
ep0=[3.3 0.1 -4.8 4 7e-4 5e-11 0.68 0.78 5];

ept=inline('ep(1)*w.^(-ep(2))+real(ep(3)./(1+(j*w*ep(5)).^(1-ep(7)))+ep(4)./(1+(j*w*ep(6)).^(1-ep(8))))+ep(9)','ep','w');
ep=lsqcurvefit(ept,ep0,w,ydata);
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mycastle

银虫 (正式写手)


jjdg(金币+1): 感谢参与 2011-09-08 09:59:52
可能是OPTIONS.TolX的值的原因吧.
TolX - Termination tolerance on X [ positive scalar ]
是不是ep0中的猜值5e-11选得不合适?
希言自然~
2楼2011-09-08 00:35:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mycastle

银虫 (正式写手)

【答案】应助回帖


jjdg(金币+1): 感谢参与 2011-09-08 09:59:59
QiChyi(金币+5): 谢谢 mycastle,再尝试一哈 2011-09-08 10:16:24
可以输入optimset命令
设置X tolerance项。默认的是1.0e-6

或者吧那个猜值改的大一点试试。
希言自然~
3楼2011-09-08 00:40:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mycastle

银虫 (正式写手)

引用回帖:
3楼: Originally posted by mycastle at 2011-09-08 00:40:23:
可以输入optimset命令
设置X tolerance项。默认的是1.0e-6

或者吧那个猜值改的大一点试试。

不好意思,打错了

可以输入optimtool命令
设置X tolerance项。默认的是1.0e-6

或者吧那个猜值改的大一点试试
希言自然~
4楼2011-09-08 00:42:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 QiChyi 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金发布提前给出明确的时间点 +7 kulium 2026-08-21 8/400 2026-08-21 18:38 by 陈晨晨陈啊
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 14/700 2026-08-21 16:51 by applepetty
[基金申请] 人气不行了 +7 fansofjerry 2026-08-21 7/350 2026-08-21 14:13 by Vivilian
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +15 医学老男孩 2026-08-20 18/900 2026-08-21 13:46 by 医学老男孩
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 只有每年这种时候来逛逛小木虫 +23 yaoyewhu2008 2026-08-20 25/1250 2026-08-21 12:19 by Poppy1104
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 13/650 2026-08-21 08:30 by 北京莱茵编辑
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
信息提示
请填处理意见