关于MATLAB中fsolve函数的问题》

Question

%在MATLAB中有一个自带的函数blsprice.若已知s,k,r,t,v的值，则可以直接求出c,p的值。即blsprice(s,k,r,t,v)=[c,p].
%现在我欲解决的问题是，在已知k,r,t,c,p的情况下，求s,v的值。我编写的代码如下，运行中出错了，请大家帮忙看看哪里出错了。非常感谢！

%m文件：
function a=myfun(s,v)
k=340
r=0.08
t=150/365
c=29.8299
p=8.8335
a=@(s,v)(blsprice(s,k,r,t,v)-[c,p]);

%主函数：
x=fsolve('myfun',[340,0.5]);

错误使用fsolve
FSOLVE requires all values returned by user functions to be of data type double.

Answer 1

m文件修改如下：

function a=myfun(x)
k=340;
r=0.08;
t=150/365;
c=29.8299;
p=8.8335;
a=blsprice(x(1),k,r,t,x(2))-[c,p];

这样从语法上来说，不存在错误了；但是运行结果似乎找不到解，这就取决于：

①这个模型是否有解，② 你给的初值是否在解的附近。

追问

您好！感谢您的帮助。图片上是blsprice函数所求解的方程组。它以s,k,r,t,v为输入，以c,p为输出。现在我需要以c,p,k,r,t为输入，求出s,v，d1,d2.关键是求s,v.这里有4个方程，4个未知数，理论上应该是可以求出来的吧？带入的初值[340,0.5]也与正确值[350,02]比较接近。请问看了这个方程组之后您是否能帮我想办法写出正确的代码？

追答

你把[350,0.2]代入计算可知，这并不是方程的解，你检查你的数据

x=[350,0.2];
k=340;
r=0.08;
t=150/365;
c=29.8299;
p=8.8335;
a=blsprice(x(1),k,r,t,x(2))-[c,p]

a =

    -4.72878122259601e-005          20.9963527121878

追问

可是这么算就是正确的啊。而且我不利用这个自带的blsprice函数已经解决了问题，用这个函数倒没办法解决。这也是我想不明白的问题。