设为首页收藏本站联系我们
切换到宽版

声振论坛

 找回密码
 我要加入

QQ登录

只需一步，快速开始

声振论坛»论坛 › 算法编程 › 科学分析软件 › 多元函数拟合的软件和程序
返回列表 发新帖
查看: 3306|回复: 6

[其它软件] 多元函数拟合的软件和程序

[复制链接]
wlxydgq
发表于 2010-3-2 11:39 | 显示全部楼层 |阅读模式

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要 登录 才可以下载或查看，没有账号？我要加入

x
:@Q ，我现在需要拟合一个五元函数。我用的是1stopt 1.5版的，尚不说拟合好坏，关键是拟合的速度太慢，需要好多天。
大家都是用什么软件拟合多元函数的参数呢，哪一个速度会快点呢？希望大家推荐一二，如果方便的话，给一个现成的程序，我将不胜感激。
注：我用过maple和Fortran，这两个软件可以做到吗？怎么写程序呢？
回复
分享到：

使用道具 举报

adawin
发表于 2010-3-3 17:33 | 显示全部楼层
您好，1stopt这么慢啊，我还真想买呢，能和您交流下吗，我也用matlab，也许能帮到你
回复 支持 反对

使用道具 举报

dingd
发表于 2010-3-3 21:48 | 显示全部楼层
原帖由 wlxydgq 于 2010-3-2 11:39 发表

                               
登录/注册后可看大图

:@Q ，我现在需要拟合一个五元函数。我用的是1stopt 1.5版的，尚不说拟合好坏，关键是拟合的速度太慢，需要好多天。
大家都是用什么软件拟合多元函数的参数呢，哪一个速度会快点呢？希望大家推荐一二，如果方便的话 ...


最好把问题贴上来才好发现哪儿有问题！
回复 支持 反对

使用道具 举报

adawin
发表于 2010-3-4 19:30 | 显示全部楼层
主任，我的问题算写清楚了么，麻烦解答下啦，谢谢
回复 支持 反对

使用道具 举报

wlxydgq
 楼主| 发表于 2010-3-5 11:35 | 显示全部楼层
计算程序没有问题，只是运算速度慢。
就是想找一个运算快点的方法。
回复 支持 反对

使用道具 举报

wlxydgq
 楼主| 发表于 2010-3-5 11:41 | 显示全部楼层

回复 板凳 dingd 的帖子

程序如下：
Title "Type your title here";
Parameters p51,p52,p53,p54,p55,p56,p57,p58,p59,p510,p511,p61,p62,p63,p64,p65,p66,p67,p68,p69,p610,p611,p71,p72,p73,p74,p75,p76,p77,p78,p79,p710,p711;
Constant R=83.1451,Tc=647.096,p50=-10097698.7020424,p60=242005341.091914,p70=-1926309306.61686;
Variable V,pLG,VG,VL,T,y;
Function   y=R*T*(1/V-(1/10)*((10*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1)))*R*T*VG^6+
60*ln(pLG*VL/(R*T))*R*T*VL^6*VG^6-60*ln(pLG*VG/(R*T))*R*T*VL^6*VG^6-(10*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+
p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1)))*R*T*VL^6-80*pLG*VL^7*VG^6+80*pLG*VG^7*VL^6+10*VL^5*R*T*VG^7-10*VL^7*R*T*VG^5-10*VL^5*pLG*VG^8+10*VL^8*VG^5*pLG-
10*VL^8*R*T*VG^4+10*VL^4*R*T*VG^8+10*VL^9*VG^4*pLG-10*VL^4*pLG*VG^9+10*VL^5*R*T*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+
p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))*VG^2+10*VL^5*VG*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+
p77*(Tc^7/T^7-1))+10*VL^5*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))*VG^3-
10*VL^2*VG^4*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))-10*VL^3*VG^4*R*T*(p60+
p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))+10*VL^4*VG^3*R*T*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+
p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))+10*VL^4*VG^2*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+
p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))-10*VL*VG^5*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+
p77*(Tc^7/T^7-1))-10*VL^3*VG^5*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))-
10*VL^2*VG^5*R*T*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))+(5*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-
1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1)))*R*T*VL^2*VG^6-(8*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-
1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))*R*T*VL^6*VG-(5*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-
1)+p57*(Tc^7/T^7-1)))*R*T*VL^6*VG^2+(8*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+
p67*(Tc^7/T^7-1)))*R*T*VL*VG^6)/(VG^4*VL^4*R*T*(-VL^3+3*VL^2*VG+VG^3-3*VL*VG^2)*V^2)+(1/10)*(-(8*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+
p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))*R*T*VL^6*VG^2+40*VL^5*R*T*VG^8-40*VL^8*R*T*VG^5+80*pLG*VG^8*VL^6-40*VL^5*pLG*VG^9-80*VL^8*VG^6*pLG+40*VL^9*VG^5*pLG+
60*ln(pLG*VL/(R*T))*R*T*VL^6*VG^7-60*ln(pLG*VG/(R*T))*R*T*VL^6*VG^7+60*ln(pLG*VL/(R*T))*R*T*VL^7*VG^6-60*ln(pLG*VG/(R*T))*R*T*VL^7*VG^6+(10*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-
1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1)))*R*T*VG^7-(10*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+
p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1)))*R*T*VL^7+(5*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+
p57*(Tc^7/T^7-1)))*R*T*VL^2*VG^7-(5*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1)))*R*T*VL^7*VG^2-
10*VG*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))*R*T*VL^6+40*VL^5*R*T*(p60+p61*(Tc/T-1)+
p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))*VG^3+40*VL^5*VG^2*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-
1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))+40*VL^5*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+
p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))*VG^4+10*VL*VG^6*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-
1))+5*VL^3*VG^6*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))+8*VL^2*VG^6*R*T*(p60+p61*(Tc/T-
1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))-40*VL^2*VG^5*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+
p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))-40*VL^4*VG^5*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+
p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))-40*VL^3*VG^5*R*T*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+
p67*(Tc^7/T^7-1))-(5*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1)))*R*T*VL^6*VG^3+(8*(p60+
p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))*R*T*VL*VG^7-(8*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+
p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))*R*T*VL^7*VG)/(VL^4*VG^4*R*T*(-VL^3+3*VL^2*VG+VG^3-3*VL*VG^2)*V^3)-(1/10)*((15*(p50+
p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1)))*R*T*VL^2*VG^6+(18*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+
p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))*R*T*VL*VG^6+(20*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+
p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1)))*R*T*VG^6-(15*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+
p57*(Tc^7/T^7-1)))*R*T*VL^6*VG^2-(18*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))*R*T*VL^6*VG-
(20*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1)))*R*T*VL^6+60*ln(pLG*VL/(R*T))*R*T*VL^6*VG^6-
60*ln(pLG*VG/(R*T))*R*T*VL^6*VG^6-90*pLG*VL^7*VG^6+90*pLG*VG^7*VL^6+30*VL^5*R*T*VG^7-30*VL^7*R*T*VG^5-30*VL^5*pLG*VG^8+30*VL^8*VG^5*pLG+30*VL^5*R*T*(p60+p61*(Tc/T-1)+
p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))*VG^2+30*VL^5*VG*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-
1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))+30*VL^5*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+
p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))*VG^3-30*VL*VG^5*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-
1))-30*VL^3*VG^5*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))-30*VL^2*VG^5*R*T*(p60+
p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))/(R*T*VL^3*VG^3*(-VL^3+3*VL^2*VG+VG^3-3*VL*VG^2)*V^4)+
(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))/V^5+(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+
p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))/V^6+(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+
p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))/V^7);
Data;
V,pLG,VG,VL,T,y=
45421047.33        0.006116548        3710987.122        18.01900658        273.16        0.0005
22709507.93        0.006116548        3710987.122        18.01900658        273.16        0.001
15138991.54        0.006116548        3710987.122        18.01900658        273.16        0.0015
11353729.76        0.006116548        3710987.122        18.01900658        273.16        0.002
9082568.923        0.006116548        3710987.122        18.01900658        273.16        0.0025
7568457.791        0.006116548        3710987.122        18.01900658        273.16        0.003
6486945.84        0.006116548        3710987.122        18.01900658        273.16        0.0035
5675807.806        0.006116548        3710987.122        18.01900658        273.16        0.004
5044918.542        0.006116548        3710987.122        18.01900658        273.16        0.0045
4540202.959        0.006116548        3710987.122        18.01900658        273.16        0.005
4127249.636        0.006116548        3710987.122        18.01900658        273.16        0.0055
3783117.628        0.006116548        3710987.122        18.01900658        273.16        0.006
。。。。。。。。。。。。。。。。。。。
回复 支持 反对

使用道具 举报

wlxydgq
 楼主| 发表于 2010-3-5 11:42 | 显示全部楼层
参数有6列，1000多行，就不意义列举了
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-5-9 20:44 , Processed in 0.097804 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表