Question about fsolve

Question

Tommy 2012-6-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/41667-question-about-fsolve

After reading the answer of my last question, I change my code about fsolve, but there is still notification of: No solution found.

fsolve stopped because the problem appears to be locally singular.

Would anyone could check my code and give me some advises plz? Thanks a lot.:

%%function for omega of every time period for OI scheme
function f=function_1_5_3(x,Z)
r0=.02;%interest rate
sigma=.15;%vatality rate of risky asset
mu0=.06;%drift rate of risky asset
gamma=5;%risk aversion rate
M=10000;%number of trajectories
N=55;%time period
T=55;%total time period
R=40;%time of retirement
dt=T/N;%each time period
t=1:dt:T;
omega=x;
Rf=exp(r0);%riskless reture
load OI_a_matrix;%set a grid of at
Rs=exp(mu0+sigma*Z);%risky market return
a=20*mat(:,N-2);
a_1=20*mat(:,N-1);
W=((a.*(Rf+omega*(Rs-Rf))-a_1).^(-gamma)).*(Rs-Rf)/M;
f=W;

and

clear
M=10000;
x0=0.4;
Z=randn(M,1);
x=fsolve(@(x)function_1_5_3(x,Z),x0)

in which the loaded matrix is a random matrix which saved at first:

mat=rand(M,N)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Sargondjani 2012-6-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/41667-question-about-fsolve#answer_51666

first you should try yourself to solve the equation when there is just 1 line (1 a and a_1). see what's happening... and then increase to a couple of lines.

and it seem there is no (serial) correlation between a and a_1, so why not just use their expected values??? (im not a statics expert, but if they are not correlated it think that is a fair thing to do)

i will try to help with limited time, but then i first need to understand better what you want to do in words: what is the equation? do you have documentation on it? because i still find the equation obscure. i especially want to know what 'a' and a_1 and omega are (in words)

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Tommy 2012-6-25

Do you have an email address instead of the contact method on Mathworks. I do this on the basis of a paper, and I send the original paper to you via email attachment if you would provide me your email address. Thank you.

请先登录，再进行评论。

Answer 2

Sargondjani 2012-6-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/41667-question-about-fsolve#answer_51362

first of all you should delete the /M in the objective, because this only scales your problem in a bad way. another problem could arise when Rs-Rf is very close to zero, because then the objective will get very close to zero, no matter what omega is (and will thus appear singular)

does your problem solve when you have just one random element in Z, a and a_1?

actually im also a bit confused in what you are trying to do:

- before you had an omega for every element and then there is a solution for every element, but now with only one omega for all elements, there is no solution W(:)=0 (minimizing least squares seems more appropriate in that case)
-also i see time variables, but you dont do anything with (R,dt,N,t,T)
-what are 'a' and 'a_1', because I am surprised that thoser random numbers enther the objective as well (the random market return i can understand)

3 个评论
显示 1更早的评论隐藏 1更早的评论

Tommy 2012-6-22

Actually I only have the random element in Z. for others, I want to solve an individual optimal question by envelop method, so I need a and a_1 as a group of random numbers. I change something of my code of function 1_5_3 as following:

%% function for omega of every time period for OI scheme

function f=function_1_5_3(x,Z)

r0=.02;%interest rate

sigma=.15;%vatality rate of risky asset

mu0=.06;%drift rate of risky asset

gamma=5;%risk aversion rate

M=10000;%number of trajectories

N=55;%time period

T=55;%total time period

R=40;%time of retirement

dt=T/N;%each time period

t=1:dt:T;

omega=x;

Rf=exp(r0);%riskless reture

load OI_a_matrix;%set a grid of at

Rs=exp(mu0+sigma*Z);%risky market return

a=20*mat(:,N-2);

a_1=20*mat(:,N-1);

W=((a.*(Rf+omega*(Rs-Rf))-a_1).^(-gamma)).*(Rs-Rf);%regard as function 4

f=mean(W);

But there is still some problem which shows:

Equation solved, fsolve stalled.

fsolve stopped because the relative size of the current step is less than the

default value of the step size tolerance squared and the vector of function values

is near zero as measured by the default value of the function tolerance.

Thank you very much!

Sargondjani 2012-6-23

This is not a problem. Like it says the equation is solved! (any further improvements are smaller than the stopping criteria)

but i still wonder if you got what you want.... if you have properly defined the problem

what do you mean with "I want to solve an individual optimal question" because now there is still only one value for omega. and a and a_1 ARE random numbers...

... and still nothing happens with your time variables

Tommy 2012-6-25

Actually, the dt has no meaning in this program. It exists only because that I just paste the parameters from another program which also solve relevant question. a and a_1 is required for random numbers, since I want use these groups of random members in controlled interval to evaluate omega. The problem I paste about fsolve stalled is really a problem, since I could not get what I need. Ever I give a start point x0, I get an answer really closed to x0. I try to change TolX and TolFun, but they did not work. Would you please tell me what should I do? Hereby I provide the mathematical equation about omega I want to solve: 0=exp(-sigma)*Et((a(t)*(Rf+omega*(Rs-Rf))-a(t+1))^(-gamma)*(Rs-Rf)), where Et stands for the expectation at time t.Thank you very much.

请先登录，再进行评论。

Question about fsolve

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

更多回答（1 个）

3 个评论
显示 1更早的评论隐藏 1更早的评论

另请参阅

类别

标签

产品

Community Treasure Hunt

Question about fsolve

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

更多回答（1 个）

3 个评论 显示 1更早的评论隐藏 1更早的评论

另请参阅

类别

标签

产品

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

3 个评论
显示 1更早的评论隐藏 1更早的评论