QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5357|回复: 3
打印 上一主题 下一主题

subs赋值太慢,大家有遇到没,进来看下啊。

[复制链接]
字体大小: 正常 放大
cpf933        

5

主题

3

听众

13

积分

升级  8.42%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2008-6-6 10:36 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
subs赋值太慢,大家有遇到没,进来看下啊。

function G=myfun(x)
syms x1 x2 x3 x4
t=[5 10 15 20 30 45 60 90 120];
p=[1;2;3;5;10;20;50;100];
i=[1.306 1.036 0.880 0.723 0.564 0.435 0.356 0.262 0.209
1.714 1.356 1.127 0.960 0.749 0.570 0.470 0.357 0.291
1.952 1.544 1.272 1.097 0.856 0.650 0.536 0.413 0.338
2.253 1.780 1.454 1.270 0.992 0.749 0.620 0.483 0.398
2.660 2.101 1.701 1.504 1.177 0.885 0.733 0.579 0.480
3.068 2.422 1.949 1.739 1.361 1.021 0.847 0.674 0.561
3.607 2.846 2.276 2.049 1.605 1.200 0.997 0.800 0.669
4.014 3.166 2.523 2.284 1.789 1.336 1.110 0.896 0.750];
for m=1:8
for n=1:9
F(m,n)=(x1.*(1+x2.*log(p(m,1)))/(t(1,n)+x3).^x4-i(m,n)).^2;
end
end
Q=sum(sum(F));
G=subs(Q,[x1,x2,x3,x4],[x(1),x(2),x(3),x(4)])
上面是我定义的m函数
[x,fval]=fminsearch ('myfun', [15,1,10,1])
然后用fminsearch求极小值,发现subs函数赋值太慢,有没有高手帮忙解决下啊
上面的程序可以运行就是很慢。
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
liwenhui        

70

主题

65

听众

5193

积分

独孤求败

  • TA的每日心情
    擦汗
    2018-4-26 23:29
  • 签到天数: 1502 天

    [LV.Master]伴坛终老

    自我介绍
    紫薇软剑,三十岁前所用,误伤义士不祥,乃弃之深谷。 重剑无锋,大巧不工。四十岁前恃之横行天下。 四十岁后,不滞于物,草木竹石均可为剑。自此精修,渐进至无剑胜有剑之境。

    社区QQ达人 邮箱绑定达人 发帖功臣 元老勋章 新人进步奖 风雨历程奖 最具活力勋章

    群组计量经济学之性

    群组LINGO

    回复

    使用道具 举报

    zcc21        

    0

    主题

    4

    听众

    89

    积分

    升级  88.42%

    该用户从未签到

    回复

    使用道具 举报

    723816893        

    1

    主题

    2

    听众

    145

    积分

    升级  22.5%

    该用户从未签到

    新人进步奖

    我觉得不用定义x1,x2,x3,x4直接用用x(1),x(2),x(3),x(4)就可,速度还行吧,一眨眼
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2025-11-25 04:16 , Processed in 0.652828 second(s), 70 queries .

    回顶部