ok，我帮你看看

帮你修改好了，这回应该没问题了。你那个程序有问题，写得乱七八糟的。或者直接粘来的可能作者没有把输入和使用的注意事项弄清楚。看明白花了好久，还不如自己写一个呢。以后还是尽量自己写吧，粘现成的真的不放心。
我改成脚本文件了，你可以自己改成函数文件。


hope helpful~

具体的错误，其实他solve at each node 的那个part和put的边界条件都有点问题，只是这两个问题正好抵消了（通俗的说），put弄出来就对了，call就错了。

sptprice = zeros(100,1);
optprice = zeros(100,1);
for(i = 1:100)
     sptprice(i,1) = i;
     optprice(i,1) = finDiffImplicit(50,sptprice(i,1),0.05,0.3,0:1:100,0:0.001:1,'CALL');
end
plot(sptprice,optprice,'-')
grid on

lz没搞明白FDM的原理，stock价格1到100，网格上界还是100，这是不可以的。
正确的做法是设网格上界S_max,使得（0,S_max）以非常大的概率覆盖住stock，
把上界定到300试试

别忘了加好评哦，亲

太感谢了！！！
按你说的 我把上届改成300 图像果然变正常了
但是另外5张本来就正常的图 我当时写的上届就是100 是不是我也需要把这些的上届都改成300才好呢？

我试着把6张图的上届都定为300，发现这回显式看涨和显式看跌的图像都出不来，但隐式看涨、隐式看跌、CN看涨、CN看跌的图像都出来了，而且这4张图中，看跌的2张图和上届定为100的几乎没区别，看涨的2张图的option price的最高点由原先的不到50变成了超过50.
为什么这次显式的会出问题呢？
是不是我不需要把所有的上届都改为300，只需把隐式改为300，显式和CN仍然定为100呢？

显式的程序   如下
function oPrice = finDiffExplicit(X,S0,r,sig,Svec,tvec,oType)
% Function to calculate the price of a vanilla European
% Put or Call option using the explicit finite difference method
%
% oPrice = finDiffExplicit(X,r,sig,Svec,tvec,oType)
%
% Inputs: X - strike
%       : S0 - stock price
%       : r - risk free interest rate
%       : sig - volatility
%       : Svec - Vector of stock prices (i.e. grid points)
%       : tvec - Vector of times (i.e. grid points)
%       : oType - must be 'PUT' or 'CALL'.
%
% Output: oPrice - the option price
%
% Notes: This code focuses on details of the implementation of the
%        explicit finite difference scheme.
%        It does not contain any programatic essentials such as error
%        checking.
%        It does not allow for optional/default input arguments.
%        It is not optimized for memory efficiency, speed or
%        use of sparse matrces.

% Author: Phil Goddard (phil@goddardconsulting.ca)
% Date  : Q4, 2007

% Get the number of grid points
M = length(Svec)-1;
N = length(tvec)-1;
% Get the grid sizes (assuming equi-spaced points)
dt = tvec(2)-tvec(1);

% Calculate the coefficients
% To do this we need a vector of j points
j = 1:M-1;
sig2 = sig*sig;
j2 = j.*j;
aj = 0.5*dt*(sig2*j2-r*j);
bj = 1-dt*(sig2*j2+r);
cj = 0.5*dt*(sig2*j2+r*j);

% Pre-allocate the output
price(1:M+1,1:N+1) = nan;

% Specify the boundary conditions
switch oType
    case 'CALL'
        % Specify the expiry time boundary condition
        price(:,end) = max(Svec-X,0);
        % Put in the minimum and maximum price boundary conditions
        % assuming that the largest value in the Svec is
        % chosen so that the following is true for all time
        price(1,:) = 0;
        price(end,:) = (Svec(end)-X)*exp(-r*tvec(end:-1:1));
    case 'PUT'
        % Specify the expiry time boundary condition
        price(:,end) = max(X-Svec,0);
        % Put in the minimum and maximum price boundary conditions
        % assuming that the largest value in the Svec is
        % chosen so that the following is true for all time
        price(1,:) = (X-Svec(end))*exp(-r*tvec(end:-1:1));
        price(end,:) = 0;
end

% Form the tridiagonal matrix
A = diag(bj);  % Diagonal terms
A(2:M:end) = aj(2:end); % terms below the diagonal
A(M:M:end) = cj(1:end-1); % terms above the diagonal

% Calculate the price at all interior nodes
offsetConstants = [aj(1); cj(end)];
for i = N:-1:1
    price(2:end-1,i) = A*price(2:end-1,i+1);
    % Offset the first and last terms
    price([2 end-1],i) = price([2 end-1],i) + ...
        offsetConstants.*price([1 end],i+1);
end

% Calculate the option price
oPrice = interp1(Svec,price(:,1),S0);

画图程序  如下
>> sptprice = zeros(100,1);
optprice = zeros(100,1);
for(i = 1:100)
     sptprice(i,1) = i;
     optprice(i,1) = finDiffExplicit(50,sptprice(i,1),0.05,0.3,0:1:300,0:0.001:1,'CALL');
end
plot(sptprice,optprice,'-')
grid on

>> sptprice = zeros(100,1);
optprice = zeros(100,1);
for(i = 1:100)
     sptprice(i,1) = i;
     optprice(i,1) = finDiffExplicit(50,sptprice(i,1),0.05,0.3,0:1:300,0:0.001:1,'PUT');
end
plot(sptprice,optprice,'-')
grid on

真的非常感谢！！！
主要是我实在没本事自己写程序啊 所以只好粘现成的了

恩 必须的亲

理论上是这样的

数值计算的内部运作本来就很微妙，我一时也想不起来原因。
出结果也许有巧合的成分。
“显式差分格式”不稳定，不一定收敛到真实解，不一定算出结果。

做期权计算，建议使用稳定的方法，比如“隐式差分”。如果有文献
说明了一种方法不稳定，最好不要用。

谢谢
那我还是把隐式和CN的上届都定位300，然后显示就还是100吧，至少这样图像都是正常的

再弱弱的问一句， 此时 S_max=300，那M和N分别等于多少呢？

和上界大小无关，只取决于你划分的格子的密度。这个你自己设定。N是时间维度的密度，跟股票就更加没有关系了。再说一遍，这个程序其实最要命的地方是他最后一部分写得有问题，画图的问题倒是其次，我用我自己的方法画出来也跟你差不多。所以你不用拼命去改S的上界，我帮你改的股价是10，上界是15都没有问题。楼主一直无视我的话，我很怀疑你有没有运行过我给你的东西。隐式方法是不可能画出不稳定的诡异图像的，所以不是画图的问题，是程序本来就写错了。

亲啊 我不是无视你啊 我真的有认真看啊而且非常感谢你 只是我实在是matlab白痴 所以你之前给我的脚本 我不知道该怎么才能改成函数文件啊 不会用啊 能不能麻烦你一步一步告诉我该怎么做呢
如果你觉得这个程序实在不好 要是你有你的程序 能成功弄出价格 能不能告诉我呢
或者你直接做 成功以后 直接把matlab程序和生成的图像发给我

能告诉我你的qq吗 方便联系 要是我的matlab部分真的搞定了 我愿意支付酬劳表示感谢
真的很谢谢你啊亲

严格地讲，普通欧式期权的PDE在“价格”这一维度上没有指定边界条件，本质是一个“开放问题”。如果直接解“开放问题”，很麻烦，PDE数值解中有专门的研究。
设定一个上界是技术上的考虑，方便实施，换个角度看，相当于转化成了一个上界非常高的“向上敲出期权”。这种“向上敲出期权”的极端状况就是普通欧式期权

恩，其实没那么复杂,从hedge的角度来说deep in the money的option就可以看成一个forward，这样就好算了。pde的边界条件其实是隐含股价是0和无穷的只是一般都不写。有限差分就用很大的数代替无穷咯~

请问 你知道怎样用matlab算出B-S PDE的准确值并画出曲线吗？我想用FDM算出的近似值和准确值做比较 然后算出误差并作图   谢谢！！！

很感谢你那么耐心的解答我的问题 如果我的论文搞定了 我愿意支付酬劳表示感谢！

我也在做相关论文，至今FDM的原理还不是很清楚

感谢解答，也解决了我的问题

第一次接触，学习学习。

层主你好，请问你的程序中的Boundary是什么意思？

边界条件的意思