应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 概率论与统计 » 用Mathematica软件进行最小二乘法曲线拟合遇到问题，哪位大神帮忙看看？
333/4返回列表上一页1234下一页
查看: 3269  |  回复: 32
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

meina861010

银虫 (小有名气)
引用回帖:
19楼: Originally posted by xzczd at 2015-09-18 14:25:37
阿尔法和西塔的取值可能是复数？但是你用于拟合的原始数据又全是实的？这个模型真的没有问题吗？

如果真的没有问题，那么：
1. 这大概又是模型本身的强非线性导致的拟合失败，可能需要精心寻找合适的初值，关于 ...

你好，请问能不能把你试的拟合效果不好的代码发我一份（我写的根本就不显示拟合曲线）。我赋一下可能的初始值，刚学这个软件，写进去一堆错误，麻烦了！
赞一下 回复此楼
21楼2015-09-18 15:32:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
meina861010(feixiaolin代发): 金币+50 2015-09-21 13:08:06
引用回帖:
21楼: Originally posted by meina861010 at 2015-09-18 15:32:17
你好，请问能不能把你试的拟合效果不好的代码发我一份（我写的根本就不显示拟合曲线）。我赋一下可能的初始值，刚学这个软件，写进去一堆错误，麻烦了！...

CODE:
model = ((1/2 Exp[I θ] Tanh[r])^(x/2) Exp[-(1/2) (Abs[α]^2 + Conjugate[α]^2 Exp[I θ] Tanh[r])] (1 +
       x + x^2 + x^3 + x^4 + x^5 + x^6 + x^7 + x^8 + x^9 + x^10 + x^11 + x^12 + x^13 +
       x^14 + x^15 + x^16) (α + Conjugate[α] Exp[I θ] Tanh[r]))/(Sqrt[x! Cosh[r]] Sqrt[2 Exp[I θ] Tanh[r]]);

nlm = ComplexFit[data, model, {α, r, θ}, x(*,Method -> NMinimize*)]

Show[ListPlot[data, PlotRange -> All],
Plot[{Re@#, Im@#} &@model /. nlm // Evaluate, {x, 0, 20}, PlotRange -> All],
Frame -> True]

注意我用了上面那帖里的程序包。
赞一下 回复此楼
小木虫Mathematica版块已毁（当然原本也不咋的），建议大家前往百度贴吧或Stackexchange。
22楼2015-09-18 15:55:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)
引用回帖:
21楼: Originally posted by meina861010 at 2015-09-18 15:32:17
你好，请问能不能把你试的拟合效果不好的代码发我一份（我写的根本就不显示拟合曲线）。我赋一下可能的初始值，刚学这个软件，写进去一堆错误，麻烦了！...

还有，用这个方法拟合出来的曲线其实又是在复数域上的，但是虚部非常小（注意观察上面那段代码执行出来的图象上的红线），可以认为只是误差。
赞一下 回复此楼
小木虫Mathematica版块已毁（当然原本也不咋的），建议大家前往百度贴吧或Stackexchange。
23楼2015-09-18 15:56:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meina861010

银虫 (小有名气)
引用回帖:
23楼: Originally posted by xzczd at 2015-09-18 15:56:52
还有，用这个方法拟合出来的曲线其实又是在复数域上的，但是虚部非常小（注意观察上面那段代码执行出来的图象上的红线），可以认为只是误差。...

谢谢啊，运行之后，还是遇到一样的问题：不显示拟合曲线，如附件红色区域所示。难道是我软件的问题吗，但是画其它的可以啊？
用Mathematica软件进行最小二乘法曲线拟合遇到问题，哪位大神帮忙看看？
22.jpg
赞一下 回复此楼
24楼2015-09-18 16:45:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meina861010

银虫 (小有名气)
引用回帖:
23楼: Originally posted by xzczd at 2015-09-18 15:56:52
还有，用这个方法拟合出来的曲线其实又是在复数域上的，但是虚部非常小（注意观察上面那段代码执行出来的图象上的红线），可以认为只是误差。...

没有显示你运行出的红线，
赞一下 回复此楼
25楼2015-09-18 16:47:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)
引用回帖:
25楼: Originally posted by meina861010 at 2015-09-18 16:47:13
没有显示你运行出的红线，...

……你就不能认真点读我链的那个帖子吗？在Stackexchange的代码已经不是最新的了，最新的程序包在github上面（链接也在那帖里），语法和这边的已经有所不同了。还有，你不觉得你nlm那行的输出有点奇怪吗？如果不知道错在哪里，那么，请仔细读一读Re，Im，ReplaceAll(也就是那个“/.”)的自带帮助。总之，你先从github上把最新版本的程序包下载下来。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

meina861010
小木虫Mathematica版块已毁（当然原本也不咋的），建议大家前往百度贴吧或Stackexchange。
26楼2015-09-18 19:28:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meina861010

银虫 (小有名气)
送红花一朵
引用回帖:
26楼: Originally posted by xzczd at 2015-09-18 19:28:54
……你就不能认真点读我链的那个帖子吗？在Stackexchange的代码已经不是最新的了，最新的程序包在github上面（链接也在那帖里），语法和这边的已经有所不同了。还有，你不觉得你nlm那行的输出有点奇怪吗？如果不知 ...

请息怒，弱弱的问一句最新版本的程序包是TransformedFit.nb[ 1407 lines (1381 sloc) ]这个吧。所有代码都是copy的你的，除了程序包。请理解一个还没入门新手的蠢笨吧
如果可以的话，能不能把你编好的能拟合出曲线的nb文件给俺，省得老麻烦你
赞一下 回复此楼
27楼2015-09-18 22:42:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meina861010

银虫 (小有名气)
引用回帖:
27楼: Originally posted by meina861010 at 2015-09-18 22:42:58
请息怒，弱弱的问一句最新版本的程序包是TransformedFit.nb这个吧。所有代码都是copy的你的，除了程序包。请理解一个还没入门新手的蠢笨吧
如果可以的话，能不能把你编好的能拟合出曲线的nb文件给俺，省 ...

可以看着你的正确的代码学习
赞一下 回复此楼
28楼2015-09-18 22:44:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meina861010

银虫 (小有名气)
引用回帖:
26楼: Originally posted by xzczd at 2015-09-18 19:28:54
……你就不能认真点读我链的那个帖子吗？在Stackexchange的代码已经不是最新的了，最新的程序包在github上面（链接也在那帖里），语法和这边的已经有所不同了。还有，你不觉得你nlm那行的输出有点奇怪吗？如果不知 ...

最新的程序包及运行代码nb文件，让您老确认一下是不是这个样子，运行文件见附件。
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : TransformedFit.m
    • 2015-09-19 10:35:13, 6.4 K
29楼2015-09-19 10:36:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meina861010

银虫 (小有名气)
引用回帖:
26楼: Originally posted by xzczd at 2015-09-18 19:28:54
……你就不能认真点读我链的那个帖子吗？在Stackexchange的代码已经不是最新的了，最新的程序包在github上面（链接也在那帖里），语法和这边的已经有所不同了。还有，你不觉得你nlm那行的输出有点奇怪吗？如果不知 ...

上传错了，这个是正确版本。
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : TransformedFit.nb
    • 2015-09-19 10:40:59, 104.65 K
30楼2015-09-19 10:41:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 meina861010 的主题更新
333/4返回列表上一页1234下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人考085602 化学工程 专硕 +7 不知道叫什么！ 2026-03-15 8/400 2026-03-15 20:11 by 棒棒球手
[考研] 283求调剂 +6 小楼。 2026-03-12 9/450 2026-03-14 22:15 by cfxczx2015
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 331求调剂（0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-09 6/300 2026-03-14 02:14 by JourneyLucky
[考研] 求调剂，一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 材料与化工（0856）304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[硕博家园] 深圳大学硕士招生（2026秋，传感器方向，仅录取第一志愿） +4 xujiaoszu 2026-03-11 7/350 2026-03-13 17:28 by xujiaoszu
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 282分材料专业求调剂院校 +18 枫桥ZL 2026-03-09 25/1250 2026-03-13 10:47 by 白夜悠长
[考研] 0857 资源与环境 285分 +6 未名考生 2026-03-09 6/300 2026-03-11 21:08 by 30660438
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
信息提示
关闭
请填处理意见
关闭