Question

急急急~matlab 解方程组！！！求code!!

L*sigma_L= A*sigma_A*N(d1)
where d1= (Log(A/B)+(r+1/2*(sigma_A)^2)*T)/(sigma_A*sqrt(T))；

求A 和 sigma_A, N()是standard normal cdf, 其他所有值都是输入值，随便输入什么值都行，只要能解出这组方程，求matlab的code, 急急急！！！跪谢！！！

Answer 1

1
First, write a file that computes F, the values of the equations at x.
Save this function file as myfun.m somewhere on your MATLAB path.
M文件：
function F = myfun(x)
% 输入参数
L = 10;
sigma_L = 1;
d1 = 1.6;
B = 5;
T = 3;
r = 0.05;
A = x(1);
sigma_A = x(2);
F = [L*sigma_L - A*sigma_A*normcdf(d1,0,1);(log(A/B)+(r+1/2*(sigma_A)^2)*T)/(sigma_A*sqrt(T)) - d1];
2 Next, set up the initial point and call fsolve:
命令窗口：
x0 = [8; 1]; % Make a starting guess at the solution
[x,fval] = fsolve(@myfun,x0); % Call solver
A = x(1)
sigma_A = x(2)
3 fsolve finds an answer:
计算结果：
A = 6.2081
sigma_A = 1.7042

追问

朋友 我错了，我少写了一个方程式，正确写法为：
L= A*N(d1) - B*exp(-r*T)*N(d1-sigma_A*sqrt(T))；
L*sigma_L= A*sigma_A*N(d1)；
where d1= (Log(A/B)+(r+1/2*(sigma_A)^2)*T)/(sigma_A*sqrt(T))；
还是求A 和sigma_A~大神 我这有一组数 特别大。不知道你能不能解出来：L=15692*10^9,
B=5156*10^9, T=1/12, r=0.2, sigma_L=908344, 这数据结果我真心急！！如果可以解出来 不知道能留个q，以后好交流，跪谢了先！！

追答

这组数据中的sigma_L过大，是你算错了，还是写错了？