[问答] GLS模型非线性求解 [推广有奖]

0关注
0粉丝

高中生

75%

还不是VIP/贵宾

-

0%

威望: 0 级
论坛币: 0 个
通用积分: 1.0000
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 198 点
帖子: 3
精华: 0
在线时间: 73 小时
注册时间: 2009-4-7
最后登录: 2021-12-12

楼主

wah79 发表于 2016-4-8 14:51:52 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求助：

本人matlab菜鸟，在用fsolve程序解非线性方程时遇到问题，具体程序、数据以及所用方程如附件。
程序以及提示错误如下，百思不得其解，求解答！！

function x=Problem_1()
clear all
%  x是100个r的解向量
%  把函数Problem_1()文件和ProblemData_1数据保存放在桌面或者指定的matlab搜索文件夹(Current Folder设置)下
%  在Command Window里面输入x=Problem_1()即可得到x
clc
load matlab1111
x=zeros(1,length(P));
%  初始值
r0=0;
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:12
if i<12
F=F+(FROE(i)-r)*B(i)/((1+r)^i);
else
   F=F+(FROE(i)-r)*B(i)/(r*(1+r)^i);
end
end
y=P-F;
%  求解y等于0
end

提示错误：
错误使用 trustnleqn (line 28)
Objective function is returning undefined values at initial point. FSOLVE cannot continue.
出错 fsolve (line 366)
[x,FVAL,JACOB,EXITFLAG,OUTPUT,msgData]=...
出错 Problem_1 (line 12)
x(i)=fsolve(@(r)myfun(r,P(i),B(i,:),FROE(i,:)),r0);