matlab求解高次方程

楼主，最近我也在算权益成本，能分享一下你的最终的代码吗？万分感谢，邮箱chenyipu_work@163.com，跪谢跪谢

tongleye 发表于 2015-11-29 21:43
需要！！谢谢楼主！

您好，您有权益成本计算的matlab代码吗？可以转发给我一份吗？邮箱chenyipu_work@163.com

你好，可以请教一下你后面几期B是怎么算的吗？然后可以分享一下代码吗？邮箱lizhaohui903@163.com

楼主好，我之前计算过，出现了好多负数。现在很是需要计算出权益资本成本，不知能否分享下你的计算代码呢？如果可以，劳烦发下代码，邮箱：sophia.travel@hotmail.com,谢谢！！！

楼主好，能够把您最后的计算代码发给我看一下吗？谢谢啦，837536374@qq.com

fantuanxiaot 发表于 2014-11-18 14:38
这是解答

这个是matlab的程序还是sas程序啊？

楼主，求分享，lx_5988@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

风流怎闲却 发表于 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

KJ001 发表于 2015-12-2 10:34
楼主，能否分享下计算的代码？

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