24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3170  |  回复: 32

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朵)

小木虫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 的主题更新
信息提示
请填处理意见