| 查看: 914 | 回复: 4 | ||
[求助]
（急）：用1stopt跑一段程序，解决非线性方程组的解（已有程序，需要较高版本的1stopt
|
|
程序如下：毕业论文要用，一共三段程序，求每段程序在不同参数值下的结果。即每段程序有四组解。我用的1stopt1.0解不出来，估计版本太低的缘故，希望热心肠的大神解决 Constant N=86794000000,S=5.41,H=0.6330,r=0.03,D=100000000,K=1,X=6,a=0.22,T=30/365,t1=0,sigmas=0.260893462; N=7441200000,S=20.9,H=0.6699,r=0.03,D=100000000,K=10,X=300,a=0.18,T=180/365,t1=0,sigmas=0.323578025； N=21098900000,S=9.88,H=0.5474,r=0.03,D=40000000,K=1,X=10.79,a=0.11,T=50/365,t1=0,sigmas=0.199024917； N=23565000000,S=11.08,H=0.5605,r=0.03,D=200000000,K=10,X=130,a=0.15,T=122/365,t1=0,sigmas=0.364097305； Function sigmas-(V*(N+K*D-D*K*normalcdf((log(K*V/(N*X))+r*(T-t1)+1/2*sigma*sigma*(T^(2*H)-t1^(2*H))+0.5*sigma*sigma*a^2*(T-t1))/sqrt(sigma*sigma*(T^(2*H)-t1^(2*H))+sigma*a^2*(T-t1))))/(S*N*(N+D*K)))*sigma; S*N-V+D*(1/(N+K*D))*(K*V*normalcdf(((log(K*V/(N*X))+r*(T-t1)+1/2*sigma*sigma*(T^(2*H)-t1^(2*H))+0.5*sigma*sigma*a^2*(T-t1))/sqrt(sigma*sigma*(T^(2*H)-t1^(2*H))+sigma*a^2*(T-t1))))-exp(-r*(T-t1))*N*X*normalcdf(((log(K*V/(N*X))+r*(T-t1)-1/2*sigma*sigma*(T^(2*H)-t1^(2*H))-0.5*sigma*sigma*a^2*(T-t1))/sqrt(sigma*sigma*(T^(2*H)-t1^(2*H))+sigma*a^2*(T-t1))))); W=(V-N*S)/D ; Constant N=86794000000,S=5.41,H=0.6330,r=0.03,M=100000000,L=1,X=6,T=30/365,t1=0,sigmas=0.260893462; N=7441200000,S=20.9,H=0.6699,r=0.03,M=100000000,L=10,X=30,T=180/365,t1=0,sigmas=0.323578025； N=21098900000,S=9.88,H=0.5474,r=0.03,M=40000000,L=1,X=10.79,T=50/365,t1=0,sigmas=0.199024917； N=23565000000,S=11.08,H=0.5605,r=0.03,M=200000000,L=10,X=13,T=122/365,t1=0,sigmas=0.364097305； Function sigmas-(V*(N+M*L-M*L*normalcdf((log(V/(N*X))+r*(T-t1)+1/2*sigma*sigma*(T^(2*H)-t1^(2*H)))/sqrt(sigma*sigma*(T^(2*H)-t1^(2*H)))))/(S*N*(N+M*L)))*sigma; S*N-V+M*(L/(N+L*M))*(V*normalcdf(((log(V/(N*X))+r*(T-t1)+1/2*sigma*sigma*(T^(2*H)-t1^(2*H)))/sqrt(sigma*sigma*(T^(2*H)-t1^(2*H)))))-exp(-r*(T-t1))*N*X*normalcdf(((log(V/(N*X))+r*(T-t1)-1/2*sigma*sigma*(T^(2*H)-t1^(2*H)))/sqrt(sigma*sigma*(T^(2*H)-t1^(2*H)))))); W=(V-N*S)/M ; Constant N=86794000000,S=5.41,r=0.03,M=100000000,L=1,X=6,T=30/365,t1=0,sigmas=0.260893462; N=7441200000,S=20.9,r=0.03,M=100000000,L=10,X=300,T=180/365,t1=0,sigmas=0.323578025； N=21098900000,S=9.88,r=0.03,M=40000000,L=1,X=10.79,T=50/365,t1=0,sigmas=0.199024917； N=23565000000,S=11.08,r=0.03,M=200000000,L=10,X=130,T=122/365,t1=0,sigmas=0.364097305； Function sigmas-(V*(N+M*L-M*L*normalcdf((log(L*V/(N*X))+(r+1/2*sigma*sigma)(T-t1))/sqrt(sigma*sigma*(T-t1))))/(S*N*(N+M*L)))*sigma; S*N-V+M*(1/(N+L*M))*(LV*normalcdf(((log(L*V/(N*X))+(r+1/2*sigma*sigma)*(T-t1))/sqrt(sigma*sigma*(T-t1))))-exp(-r*(T-t1))*N*X*normalcdf(((log(L*V/(N*X))+(r-1/2*sigma*sigma)*(T-t1))/sqrt(sigma*sigma*(T-t1))))); W=(V-N*S)/M ; |
回复此楼» 猜你喜欢
声子谱计算问题
已经有2人回复
-
已结束
已经有23人回复
-
物理学I论文润色/翻译怎么收费?
已经有255人回复
-
Error in routine lambda (100): wrong # or too many modes
已经有0人回复
-
祈福国自然面上
已经有57人回复
-
请问castep计算出的能带图横轴数据缩放因子怎么计算
已经有1人回复
-
Result in Physics 从见刊开始 一般多久可以被SCI收录
已经有15人回复
-
面青地会评时间？？？
已经有6人回复
-
科研求助
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助，用1stopt解方程组，两次得到的解不唯一？什么原因，大神帮忙
已经有6人回复
- sos!!1stopt解非线性方程组，一直出现常数定义错误的问题！求大神帮助！
已经有12人回复
- 求助matlab解一个非线性方程组
已经有8人回复
- matlab解非线性方程组解不开
已经有3人回复
- 求高版本1stopt，解决一个问题，微分方程组参数拟合
已经有3人回复
- 非线性方程组的解法
已经有5人回复
- 求fortran解非线性方程组
已经有4人回复
- 1stOpt程序的一个小问题
已经有5人回复
- 求助matlab 非线性方程组全部解！
已经有6人回复
- Matlab：Newton-Raphson 解线非线性性方程组的问题
已经有12人回复
- 1stopt运行程序，总是出现字符函数错误的提示，请高手帮助
已经有5人回复
- 1stopt4阶非线性常微分方程，帮小弟看一下，跪求。
已经有7人回复
- matlab求解非线性方程组，求助！
已经有6人回复
- 如何使用mathematica求解非线性微分方程组
已经有11人回复
- 1stopt问题请教(附代码)
已经有5人回复
- 哪位大哥能帮我用1stopt3.0以上版本跑一个拟合程序
已经有11人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 请教Matlab高手解多元非线性方程组中的参数
已经有6人回复
- 【求助】求助，解非线性方程组，急
已经有5人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 【求助】求一个用牛顿迭代法解非线性方程组的程序
已经有10人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15356.8
- 散金: 101
- 红花: 234
- 帖子: 3409
- 在线: 1222.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
xiaokang1213: 金币+15, ★★★★★最佳答案 2013-05-08 17:59:49
csgt0: 金币+1, 谢谢 2013-05-09 14:23:33
感谢参与，应助指数 +1
xiaokang1213: 金币+15, ★★★★★最佳答案 2013-05-08 17:59:49
csgt0: 金币+1, 谢谢 2013-05-09 14:23:33
|
1.1 v: 469550326114.98 sigma: 0.260969439328504 w: -0.0521383621322139 1.2 v: 154999232907.983 sigma: 0.336346161451454 w: -5.21844347271593 1.3 v: 208455921385.871 sigma: 0.199143912876964 w: -0.0302631841514178 1.4 v: 261533312786.74 sigma: 0.376126128056199 w: 2.1656234982059 2.1 v: 469552129840.991 sigma: 0.260958386111791 w: -0.0341014140819834 2.2 v: 155081405944.153 sigma: 0.335645389568348 w: -4.3967225825483 2.3 v: 208456252403.897 sigma: 0.199141484693512 w: -0.0219856949728495 2.4 v: 261553549485.575 sigma: 0.375994554468124 w: 2.26680235453713 3.1 v: 469553608690.895 sigma: 0.260979523693026 w: -0.0193121154194043 3.2 v: 150276581384.821 sigma: 0.347243050976798 w: -52.4449479591908 3.3 v: 208456886695.061 sigma: 0.19914820898378 w: -0.00612851218441807 3.4 v: 254021338434.329 sigma: 0.386898981212979 w: -35.3943078143853 |
2楼2013-05-08 09:23:35
xxz903
金虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 1050.6
- 散金: 190
- 红花: 4
- 帖子: 253
- 在线: 66.7小时
- 虫号: 1360336
- 注册: 2011-08-03
- 性别: GG
- 专业: 化工系统工程
感谢参与，应助指数 +1
csgt0: 应助指数-1, 非应助请勿选择应助 2013-05-09 14:23:46
csgt0: 应助指数-1, 非应助请勿选择应助 2013-05-09 14:23:46
| dingd已经给你算好了啊~ |
» 本帖已获得的红花（最新10朵）
xiaokang12133楼2013-05-08 16:11:38
4楼2013-05-08 17:59:15
5楼2013-05-08 18:00:45
20
