I need to put all the output of the function in vector
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
Hasan
am 11 Nov. 2015
Kommentiert: Hasan
am 12 Nov. 2015
function [Eg]= test_cal(L,U,P)
for t=-0.2:0.01:0;
xx=1:length(t);
x3 = 0.20*t^2 + 5.65e-5/cos(cos(8.92*t)) ...
+ 0.0003*t*cos(cos(8.922*t)) + ...
0.037*t^3*cos(cos(8.92*t)) - ...
0.031 - 0.31*t - 0.36*t^4;
x2 = 0.31*t + 0.002*sin(14.83*sin(t))^3 ...
+ 0.0012*sin(14.83*sin(t))^2 + ...
0.0008*sin(14.83*sin(t))/t - ...
0.0072 - 0.098*t^2 - ...
0.0015*cos(t)*sin(14.83*sin(t));
x1 = 0.36+ 2.46*t + 28.82*t^3 + ...
3.46*t^4 + 2.46*t^5*sin(sin(sin(99.88*t^2))) - ...
0.011*sin(sin(99.87*t^2));
x0 = 2.66 + 0.43*t*cos(30.07*t^2) + ...
0.015*t*cos(30.08*t^2)^2 - 0.48*t ...
- 0.73*t^2 - 0.097*cos(30.08*t^2) ...
- 1.19*t^3*cos(30.08*t^2);
a0 = 1.42e-10;
for tt=0.05:0.05:1
alpha0 = x0;
alpha1 = x1;
alpha3 = x3;
c = 1;
U=0;
kx = -2*pi/(3*a0):2*pi/(3*499.5*a0):2*pi/(3*a0);
for c = 1:1000
hk(c) = (exp(0)+exp(i*3*a0/2*kx(c))+exp(-i*(3*a0/2*kx(c)))); %h(k) for ky = 0.
z=sin(c);
H = [U,-alpha0*hk(c),-alpha1,-alpha3*hk(c); %H matrix
-alpha0*conj(hk(c)),U,-alpha3*conj(hk(c)),-alpha1;
-alpha1,-alpha3*hk(c),0,-alpha0*hk(c);
-alpha3*conj(hk(c)),-alpha1,-alpha0*conj(hk(c)),0];
E(:,c) = eig(H); %Solve the eigenvalue problem.
end
E3 =E(3,:);
E4 = E(4,:);
E2 = E(2,:);
Min = abs(E2-E3);
Eg = min(Min)
end
end
end
The output is like the following
Eg =
0.0069
Eg =
0.0069
Eg =
0.0069
Eg =
0.0069
Eg =
0.0069
Eg =
0.0069
Eg =
0.0069
Eg =
0.0069
So I need it as
Eg=[0.0069 0.0069 ...]
0 Kommentare
Akzeptierte Antwort
Walter Roberson
am 12 Nov. 2015
Right before the line
for tt=0.05:0.05:1
add the line
Egidx = 0;
then change the line
Eg = min(Min)
to
Egidx = Egidx + 1;
Eg(Egidx) = min(Min);
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Particle Swarm 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!