风流怎闲却 发表于 2015-12-11 22:27
楼主,最近我也在算权益成本,能分享一下你的最终的代码吗?万分感谢,邮箱chenyipu_work@163.com,跪谢跪谢 ...
function x=Problem_1()
clear all
% x是1408个r的解向量
% 把函数Problem_1()文件和ProblemData_1数据保存放在桌面或者指定的matlab搜索文件夹(Current Folder设置)下
% 在Command Window里面输入x=Problem_1()即可得到x
clc
load ProblemData_1
x=zeros(1,length(P));
% 初始值
r0=0.1;
for i=1:length(P)
x(i)=fsolve(@(r)myfun(r,P(i),B(i,:),FROE(i,:)),r0);
end
x=x';
end
% 求解的r的函数
function y = myfun(r,P,B,FROE)
F=B(1);
for i=1:13
if i<13
F=F+(FROE(i)-r)*B(i)/((1+r)^i);
else
F=F+(FROE(i)-r)*B(i)/(r*(1+r)^12);
end
end