weilinhy 发表于 2012-2-3 22:15

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

问题:下面是迭代的最终结果,我想保存每次迭代结果,比如迭代第一次、第二次、第三次……一直到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
====== 计算结束 ======

weilinhy 发表于 2012-2-13 08:50

再次顶起 请大家帮忙看看啊 谢谢

weilinhy 发表于 2012-3-10 09:39

顶起来啊 大家帮忙啊
页: [1]
查看完整版本: 请问 1stopt怎么保存每次迭代结果