Find values that put equation to zero
13 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
MARCO FORTI
am 9 Apr. 2021
Kommentiert: MARCO FORTI
am 10 Apr. 2021
In this model:
where all the variables are of dymension TX1, I need to find all the (T) values of gamma that put the equation to zero.
I tryed with the following code, but the program returns me all NAN, how can I fix the code to get it?
j=1; %increment level. The loop will work for j>1 as well.
for i=1:j:T;
fnc = @(gamma,C,S) 1./(1+gamma).*(C(i)-(S(i)*gamma)) - M(i); % Model function
rn = @(gamma) norm(M - fnc(gamma,C,S)); % Residual norm
gamma0 = [15]; % Inizializzazione
[gamma, ResNorm] = fzero(rn, gamma0);
end
0 Kommentare
Akzeptierte Antwort
Matt J
am 10 Apr. 2021
Bearbeitet: Matt J
am 10 Apr. 2021
gamma=nan(T,1);
j=1; %increment level. The loop will work for j>1 as well.
for i=1:j:T;
ci=C(i); si=S(i); mi=M(i);
fnc = @(gamma) 1./(1+gamma).*(ci-(si*gamma)) - mi;
gamma0 = 15;
[gamma(i), ResNorm] = fzero(fcn, gamma0);
end
3 Kommentare
Matt J
am 10 Apr. 2021
You may be initializing at a bad point. This works:
T=5;
[C,S,M]=deal(rand(T,1), rand(T,1), rand(T,1));
gamma=nan(T,1);
j=1; %increment level. The loop will work for j>1 as well.
for i=1:j:T;
ci=C(i); si=S(i); mi=M(i);
fnc = @(gamma) 1./(1+gamma).*(ci-(si*gamma)) - mi;
gamma0 = 15;
[gamma(i), ResNorm] = fzero(fnc, 0);
end
gamma
Weitere Antworten (1)
Alan Stevens
am 9 Apr. 2021
Bearbeitet: Alan Stevens
am 9 Apr. 2021
If your expression is equal to zero, then you can rearrange it to get gamma directly. A little algebra gives
gamma = (C-M)./(S+M);
No need for fzero at all.
Siehe auch
Kategorien
Mehr zu Gamma Functions finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!