[高手进]求助！用MATLAB计算非线性方程组，编程后输出结果显示有错，请各位高手帮帮忙

6633

收藏 2010-03-13

急！求助：以下是我解非线性方程组引用的fsolve命令的程序以及输出的错误，望各位高手帮忙！谢谢..
主程序是：
function F=myfun(x)
a=(a(1);a(2);a(3))
F=[a(1)-x(1)*normcdf((log(x(1)/a(2))+(0.025+x(2)*x(2))/2)/x(2))+a(2)*exp(-0.025)*normcdf((log(x(1)/a(2))+(0.025-x(2)*x(2))/2)/x(2));normcdf((log(x(1)/a(2))+(0.025+x(2)*x(2))/2)/x(2))*x(1)*x(2)/a(1)-a(3)];
输入命令：
a=[851169606 1059209620 0.514241];
>> x0=[5;0.4];
>> [x,fval]=fsolve(@myfun,x0,optimset('Display','iter'))
后显示的错误：
??? Error: File: d:\MATLAB7\work\myfun.m Line: 2 Column: 8
Incomplete or misformed expression or statement.

Error in ==> fsolve at 180

fuser = feval(funfcn{3},x,varargin{:});

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

aris_zzy

2010-3-13 21:36:45

你这是求解BS模型的隐含波动率什么的?
错误提示语法不完全

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2010-3-13 22:53:50

有参数最简单的方法,是用global Variables.
但global Variables要小心用.
function F=myfun(x,a)
global a
F=[a(1)-x(1)*normcdf((log(x(1)/a(2))+....

%%%in MATLAB workspace,
global a
a=[851169606 1059209620 0.514241];
x0=[5;0.4];
[x,fval]=fsolve(@myfun,x0,optimset('Display','iter'))

Optimizer appears to be converging to a minimum that is not a root:
Sum of squares of the function values is > sqrt(options.TolFun).
Try again with a new starting point.

x =
5.0000
0.4000
fval =
1.0e+008 *

8.5117
-0.0000

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

lingfei2mm

2010-3-14 10:56:57

恩，是的！我要求的就是BS模型的隐含波动率和资产价值，不知道我这程序中哪些语法出现错误，还望多多赐教，非常感谢！！！！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

lingfei2mm

2010-3-14 11:06:07

非常感谢这么详细的指导，但是我还是有点不明白这个是算出来的结果吗？好像这个结果不太对哦，是不是我赋的初值不对啊？？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2010-3-14 11:21:17

Optimizer appears to be converging to a minimum that is not a root:
Sum of squares of the function values is > sqrt(options.TolFun).
Try again with a new starting point.

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

lingfei2mm

2010-3-14 11:57:18

这个意思是不是要赋一个新的初始值啊？？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2010-3-14 12:24:44

对.意思就是要你换初始值,
初始值很重要,一定要多试几个.
有空请再检查一下你的function,
方程有没写错.

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

lingfei2mm

2010-3-14 17:02:40

按照你的解释我赋了很多初始值好像还是没有作用，现在实在是没办法了，急着求结果，就只有靠你的帮忙了，您可否帮我找出适合的初始值吗？实在是感激不尽！！！！！！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

lingfei2mm

2010-3-14 17:43:40

我赋了初值【1000000000；0.1】后显示的结果是：
Optimizer appears to be converging to a minimum that is not a root:
Sum of squares of the function values is > sqrt(options.TolFun).
Try again with a new starting point.
请问里面的初值是要变大还是应该缩小的好？？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2010-3-14 19:56:06

x(1),x(2)差异很大.
global a
a=[851169606 1059209620 0.514241];
x0=[0.01;100000000];
options=optimset('Display','iter','MaxFunEvals',15000,'MaxIter',5000);
[x,fval]=fsolve(@myfun,x0,options);
format long e
x
fval
format
%%%
%%%
5000 14951 41.344 133.887 1.55e+004 134
Maximum number of iterations reached:
increase options.MaxIter.
x =
8.516556446669145e+008
6.942190255968392e+000
fval =
1.714000699575990e-002
6.429904095747514e+000

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

lingfei2mm

2010-3-14 20:14:45

太谢谢你了！！！你太厉害了，真的非常感谢了，帮我解决了一个大难题！！这个结果就是我想要的，以后多多指教哦！！！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

lingfei2mm

2010-3-14 20:37:34

你好！！再一次的打扰你了，好像上面程序中x（2）的值不是很正确呢，算出来应该得到的是个小数才对，好像大了点，(*^__^*) 嘻嘻……，请指点一下

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

dingd

2010-3-14 22:32:32

用1stOpt，不需初值，代码：

Constant a=[851169606,1059209620,0.514241];
Function a1-x1*normalcdf((ln(x1/a2)+(0.025+x2*x2)/2)/x2)+a2*exp(-0.025)*normalcdf((ln(x1/a2)+(0.025-x2*x2)/2)/x2);
normalcdf((ln(x1/a2)+(0.025+x2*x2)/2)/x2)*x1*x2/a1-a3;

结果：
x1: 1883720899.28348
x2: 0.23333732568053

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

lingfei2mm

2010-3-15 10:42:27

谢谢哦！我想问下这个软件具体应该怎么操作呢？我没用过

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

lingfei2mm

2010-3-15 12:21:17

哈哈！我会用了，我想请教你一个问题，就是在1stopt中怎么把ecxel中的数据变化成矩阵形式？？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

aris_zzy

2010-3-15 14:06:56

KMV方程组，如果你这样编程会使得你的结果与初始值关系巨大，你可以看我那边《金融数量matlab》书上的程序，使用了一个数值技巧使得初始值基本不需要选择

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

windflakes

2011-9-11 17:21:59

请问楼上那本书的全名是什么啊？
我在网上找不到呢，我也需要KMV模型编程。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

yqx1985

2011-11-25 10:27:23

Optimizer appears to be converging to a minimum that is not a root:
Sum of squares of the function values is > sqrt(options.TolFun).
Try again with a new starting point.

如果根据那句话的意思，你仅仅换个初始值估计也没那么容易解出正确结果，它的意思说计算除了最小值，不过是局部最小值，如果你知道真实结果的范围，可以考虑在真实值附近设定初始值，否则结合全局搜索可能是一个不错的选择

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群