24小时热门版块排行榜    

查看: 2890  |  回复: 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: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的回帖
查看全部 4 个回答

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

银虫 (正式写手)

引用回帖:
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 2026国自然函评费到账 +7 羊腰板 2026-08-21 7/350 2026-08-21 21:28 by iwuli
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +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
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 人气不行了 +7 fansofjerry 2026-08-21 7/350 2026-08-21 14:13 by Vivilian
[基金申请] 感觉是下周放榜了 +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
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
信息提示
请填处理意见