creating for loop in fzero
2 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
fartash2020
il 12 Nov 2017
Commentato: Star Strider
il 28 Nov 2017
Dear all, I want to plot "z" vs. "kappa" where kappa is changing from 1 to 0. this is my @ff function:
function y = f(x)
y = (1/pi)*((0.1*x)-(((sqrt((1+x)./2)).*exp(1).^-(6*(sqrt(2*0.1*(1+x))))))*(kappa*(sqrt((1-x)/(1+x)))-(sqrt((1+x)/(1-x)))));
after that using fzero I want to find "z" for every single kappa in the region in order to get the plotted data:
fun = @ff; % function
x0 = 0; % initial interval
z = fzero(fun,x0);
now how can I create a loop to change kappa from 1 to 0 and get "z" for each of them on my plot. many Thanks in Advance,
0 Commenti
Risposta accettata
Star Strider
il 12 Nov 2017
Try this:
f = @(x,kappa) (1/pi)*((0.1*x)-(((sqrt((1+x)./2)).*exp(-(6*(sqrt(2*0.1*(1+x))))))).*(kappa*(sqrt((1-x)./(1+x)))-(sqrt((1+x)./(1-x)))));
N = 10; % Length Of ‘kappa’ Vector
kappa = linspace(0, 1, N);
for k1 = 1:N
z(k1) = fzero(@(x) f(x, kappa(k1)), 1E-3);
end
figure(1)
plot(kappa, z)
xlabel('\kappa')
ylabel('z(\kappa)')
grid
4 Commenti
Star Strider
il 28 Nov 2017
Try this:
N = 100;
kd = 6;
kappa = 0.1;
lambda = linspace(0.05,0.25, N);
Delta_E = linspace(0.1,0.35, N);
[Lm,DE] = meshgrid(lambda, Delta_E);
f = @(x,lambda,Delta_E) ((Delta_E)+(lambda*x)-(((sqrt((1+x)./2)).*exp(-(kd*(sqrt(2*lambda*(1+x))))))).*(kappa*(sqrt((1-x)./(1+x)))-(sqrt((1+x)./(1-x)))));
for k1 = 1:length(lambda)
for k2 = 1:length(Delta_E)
z(k1,k2) = fzero(@(x) f(x,lambda(k1),Delta_E(k2)), rand);
end
end
figure(1)
mesh(lambda, Delta_E, z)
You might find it worthwhile to plot the absolute value of ‘z’, because the results are complex when they are not NaN, although the plot itself reveals little.
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Optimization in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!