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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求教origin拟合简单的渗流方程（又称逾渗方程），做复合材料的应该会遇到吧
212/3返回列表上一页123下一页
查看: 7874  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 计算强帖 ，点击这里进行查看

glazio

铁虫 (小有名气)
引用回帖:
7楼: Originally posted by dbb627 at 2011-11-25 23:22:39:
是可以改善的
[code]

A=[0.1001   3.5e06
0.1002      3.3e06
0.11        2.9e05
0.12        9.0e04
0.15        1.5e04
0.2        3.3e03
0.3        7.1e02
0.4        2.8e02
0.5        1.5e0 ...

这几天网络故障，好不容易盗网上线。。。

现在还剩一个问题，即
1. 我的目的是解决拟合参数共享的问题，才投入Origin的怀抱的。请问Matlab目前有解决这个问题的方案么？如果有的话请提示一二。
赞一下 回复此楼
11楼2011-11-28 22:24:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

拟合参数共享？不太明白，简而言之还有其他方程要用到上面参数吗，这个可以自己构建一个目标函数，然后再对其优化，用最小二乘或者遗传算法求参数就可以了，这个还是比较容易实现的。见帖http://muchong.com/bbs/viewthread.php?tid=3751370&page=1#pid3104538
赞一下 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
12楼2011-11-28 22:28:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

glazio

铁虫 (小有名气)
你给的帖子里的情况好像并不是我说的共享参数拟合。下面举个例子说明一下。

假如说在数据图中有两组实验数据集Y1和Y2，可以用某个含有三个参数a，b，c的方程y=f(x)描述。普通的拟合一般是用F(x)对数据集Y1和Y2分开拟合，这样会得到两组不同的参数(a1，b1,c1)和(a2，b2,c2)。
若假设数据集Y1和Y2的参数a是相同的，那么就应该使用共享参数拟合了，这样拟合的结果就应该是(a1，b1,c1)和(a1，b2,c2)。
赞一下 回复此楼
13楼2011-11-29 22:38:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

glazio(金币+10): dbb627版主不愧是matlab神人，所有问题都得到回答，非常感谢！ 2011-11-30 08:57:21
cenwanglai(计算强帖+1): 你的回复都挺好。不过我不熟悉，没有仔细看明白。就这里给一个EPI吧~ 2011-12-22 20:09:16
引用回帖:
13楼: Originally posted by glazio at 2011-11-29 22:38:41:
你给的帖子里的情况好像并不是我说的共享参数拟合。下面举个例子说明一下。

假如说在数据图中有两组实验数据集Y1和Y2，可以用某个含有三个参数a，b，c的方程y=f(x)描述。普通的拟合一般是用F(x)对数据集Y1和Y2 ...

这个没有问题
matlab可以做的
见下面的例子
CODE:
x=rand(1,7);
y1=2*x+3*sin(x)+6*x.^2;
y2=2*x+3*cos(x)+2*x.^2;%参数 2 3 6 2前面2 3共用
F=@(p,x)[p(1)*x+p(2)*sin(x)+p(3)*x.^2;p(1)*x+p(2)*cos(x)+p(4)*x.^2];
p = lsqcurvefit(F, [1 1 1 1], x,[y1;y2])

Local minimum found.

Optimization completed because the size of the gradient is less than
the default value of the function tolerance.




p =

    2.0000    3.0000    6.0000    2.0000
赞一下(1人) 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
14楼2011-11-30 00:26:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
共享参数拟合问题，1stOpt应该是目前最高效也最易使用的工具了，建议试试！
赞一下 回复此楼
15楼2011-12-04 21:50:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

715211229

木虫 (正式写手)
引用回帖:
3楼: Originally posted by dbb627 at 2011-11-25 20:14:26:
origin不常用，matlab可以拟合出较好的结果
代码如下
[code]


A=[0.1001       3.5e06
0.1002       3.3e06
0.11         2.9e05
0.12        9.0e04
0.15        1.5e04
0.2        3.3e03
0.3     ...

我用matlab cftool工具箱，怎么拟合无法进行？
赞一下 回复此楼
我是蜗牛
16楼2011-12-20 06:01:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)
呵呵，中文的matlab，很少见，这个有个参数cftool工具箱可能调不了，就是数据的归一化。
赞一下 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
17楼2011-12-20 11:28:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)
引用回帖:
16楼: Originally posted by 715211229 at 2011-12-20 06:01:11:
我用matlab cftool工具箱，怎么拟合无法进行？

呵呵，中文的matlab，很少见，你将归一化勾上，在调整初值试试。还有如果初值不合适，也肯能没有合适的结果。

[ Last edited by dbb627 on 2011-12-20 at 11:32 ]
赞一下 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
18楼2011-12-20 11:29:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

715211229

木虫 (正式写手)
引用回帖:
18楼: Originally posted by dbb627 at 2011-12-20 11:29:15:
呵呵，中文的matlab，很少见，你将归一化勾上，在调整初值试试。还有如果初值不合适，也肯能没有合适的结果。
[ Last edited by dbb627 on 2011-12-20 at 11:32 ]

不知道哪里出错了

matlab技术论坛提供的一个界面汉化包，最高到r200a
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : MATLAB 2009a(7.8.1)汉化包.rar
    • 2011-12-20 19:56:05, 1.36 M
我是蜗牛
19楼2011-12-20 19:56:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)
你用命令行运行能出结果吗？
赞一下 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
20楼2011-12-20 20:28:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 glazio 的主题更新
212/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否退出参与的面上项目解除限项 +11 koalala 2026-08-24 13/650 2026-08-24 10:49 by anjeeshine
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 建议基金发布提前给出明确的时间点 +11 kulium 2026-08-21 14/700 2026-08-24 02:12 by 流流伤
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭