Vectorization of a function
显示 更早的评论
I have dotted all the variable that are defined as a vector, but matlab gives me an error.
Do you see where is my problem?
function [C] = BlackScholesCall(S,K,t,r,sigma)
% Calculates the price of a call option
% INPUT S 1x1 ... Current stock price (underlying)
% K 1x1 ... Strike price
% t 1x1 ... Time to maturity
% r 1x1 ... Risk-free interest rate
% sigma 1x1 ... standard deviation (volatility of the underlying)
% OUTPUT C 1x1 ... The price of a call option
% USAGE BlackScholesCall(S,K,t,r,sigma)
d1=(log(S/K.))+(r+(1/2)*sigma^2)*t)/(sigma*sqrt(t));
d2=d1-sigma*sqrt(t.);
C=(S*normcdf(d1))-(K.*(exp(-r*.t))*normcdf(d2));
end
Code to call the function:
S = 22
K = 20:25
t = 0.1:0.1:0.6
r = 0.02
sigma = 0.25
C = BlackScholesCall(S, K, t, r, sigma)
Thank you very much!
2 个评论
Geoff Hayes
2019-9-30
编辑：Geoff Hayes
2019-9-30
Francesco - please copy and paste the full error message to this question. In the function header, you have
t 1x1 ... Time to maturity
which implies that t is a scalar...but in your input to this function, you define t as an array
t = 0.1:0.1:0.6
Which should it be - a scalar or an array? Perhaps this is the problem...you are passing in an array but the code is expecting a scalar? Are you the author of BlackScholesCall?
Francesco Rossi
2019-9-30
采纳的回答
the cyclist
2019-9-30
编辑：the cyclist
2019-9-30
The problem is that you don't seem to have a grasp on how the "." syntax actually works. (Sorry if that sounds harsh!)
For example, you have
log(S/K.)
I'm not sure what you intend there -- just stating that K is a vector? -- but that's just not how it works. You don't just "dot the vectors". Oversimplifying a bit -- you dot the operations not the variables.
Your code is just a bit tricky to fix, because you are also trying to use two vectors of different lengths, and I'm guessing you actually want all combinations of K and t to generate results. So, you actually have two-dimensional input. Is that right?
2 个评论
the cyclist
2019-9-30
编辑：the cyclist
2019-9-30
Oh, maybe not as tricky as I thought. Are the values of K and t paired, such that each pair of inputs gives one value of C? In that case, this should work:
d1=(log(S./K))+(r+(1/2)*sigma^2)*t./(sigma*sqrt(t));
d2=d1-sigma*sqrt(t);
C=(S*normcdf(d1))-(K.*(exp(-r*t)).*normcdf(d2));
in place of what you had. Notice how I moved some of your dots to be associated with vector operations, not the variables themselves.
the cyclist
2019-9-30
Still easier than I expected if you do not want K and t paired, but there just happened to be the same number of inputs for each.
Due to "implicit expansion", you can send in a row vector for K, and a column vector of t, and get all the combinations. Just change your input to this:
t = (0.1:0.1:0.6)'
(in addition to making the changes I suggested to your code).
更多回答（1 个）
meghannmarie
2019-9-30
1 个投票
I think you have some of your dot operators wrong:
d1 = (log(S./K) +(r+(1/2)*sigma^2)*t)/(sigma*sqrt(t));
d2 = d1-sigma*sqrt(t);
C = (S*normcdf(d1)) - (K.*(exp(-1.*t)).*normcdf(d2));
类别
在 帮助中心 和 File Exchange 中查找有关 Matrix Indexing 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
