creating for loop in fzero
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
fartash2020
le 12 Nov 2017
Commenté : Star Strider
le 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 commentaires
Réponse acceptée
Star Strider
le 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 commentaires
Star Strider
le 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.
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Optimization dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!