声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2988|回复: 2

[1stopt] 请问 1stopt怎么保存每次迭代结果

[复制链接]
发表于 2012-2-3 22:15 | 显示全部楼层 |阅读模式

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

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

x
问题:下面是迭代的最终结果,我想保存每次迭代结果,比如迭代第一次、第二次、第三次……一直到1000次时,参数V和sigma的
值(我想看看他们的变化情况)。而且像看看目标函数的变化情况,请问应该怎么弄!

源代码:

Constant N=100, M=90, k=1, X=100, T=3, ti=0, S=80, sigmas=0.75, r=0.0448, H=0.88;
Function
sigmas-(V*(((N+k*M-M*k*normcdf(((log(k*V/(N*X))+r*(T-ti)+1/2*sigma*sigma*(T-ti+T^(2*H)-ti^(2*H)))/(sigma*sqrt((T-ti+T^(2*H)-ti^(2*H)))))))/(N*(N+k*M))))/S)*sigma;
S*N-V+M*((1/(N+k*M))*(k*V*normcdf(((log(k*V/(N*X))+r*(T-ti)+1/2*sigma*sigma*(T-ti+T^(2*H)-ti^(2*H)))/(sigma*sqrt((T-ti+T^(2*H)-ti^(2*H))))))-exp(-r*(T-ti))*N*X*normcdf(((log(k*V/(N*X))+r*(T-ti)+1/2*sigma*sigma*(T-ti+T^(2*H)-ti^(2*H)))/(sigma*sqrt((T-ti+T^(2*H)-ti^(2*H)))))  -(sigma*sqrt((T-ti+T^(2*H)-ti^(2*H))))  )));
W= (V-N*S)/M;



====== 结果 ======
迭代数: 1000
计算用时(时:分:秒:毫秒): 00:00:38:250
计算中止原因: 达到最大迭代数
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
函数表达式 1: 0.75-(v*(((100+1*90-90*1*normcdf(((log(1*v/(100*100))+0.0448*(3-0)+1/2*sigma*sigma*(3-0+3^(2*0.88)-0^(2
            *0.88)))/(sigma*sqrt((3-0+3^(2*0.88)-0^(2*0.88)))))))/(100*(100+1*90))))/80)*sigma
         2: 80*100-v+90*((1/(100+1*90))*(1*v*normcdf(((log(1*v/(100*100))+0.0448*(3-0)+1/2*sigma*sigma*(3-0+3^(2
            *0.88)-0^(2*0.88)))/(sigma*sqrt((3-0+3^(2*0.88)-0^(2*0.88))))))-exp(-0.0448*(3-0))*100*100*normcdf(((log(1
            *v/(100*100))+0.0448*(3-0)+1/2*sigma*sigma*(3-0+3^(2*0.88)-0^(2*0.88)))/(sigma*sqrt((3-0+3^(2*0.88)
            -0^(2*0.88)))))  -(sigma*sqrt((3-0+3^(2*0.88)-0^(2*0.88))))  )))
         3: w-( (v-100*80)/90)
目标函数值: 0.360296522833016
v: 14455.0009295106
sigma: 1.13615526367969
w: 71.7222326442722
====== 计算结束 ======
回复
分享到:

使用道具 举报

 楼主| 发表于 2012-2-13 08:50 | 显示全部楼层
再次顶起 请大家帮忙看看啊 谢谢
 楼主| 发表于 2012-3-10 09:39 | 显示全部楼层
顶起来啊 大家帮忙啊
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-5-1 11:26 , Processed in 0.053351 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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