Effective-Medium-Approximation--EMA-/PD.m at main · xiongyuup/Effective-Medium-Approximation--EMA- · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
function PS=PD(kesei1,kesei2,f1)


% 输出解
% a=-2-2*(1-f1)/f1;
% b=2*(1-f1)*kesei1/f1+kesei2-kesei1;
% c=(1-f1)/f1*kesei1*kesei2;

% x=-2-2*f1/(1-f1);
% y=f1/(1-f1)*(2*kesei1-kesei2)+2*kesei2-kesei1;
% z=(f1/(1-f1+1))*kesei1*kesei2;
% solution = roots([x, y, z]);


l=2*(1-f1)+2*f1;
m=(1-f1)*kesei1-2*(1-f1)*kesei2-2*f1*kesei1+f1*kesei2;
n=-f1*kesei1*kesei2-(1-f1)*kesei1*kesei2;

% y=f1/(1-f1)*(2*kesei1-kesei2)+2*kesei2-kesei1;
% z=(f1/(1-f1+1))*kesei1*kesei2;
% solution = roots([x, y, z]);

% 输出解
% disp('方程的解为：');
% disp(solution);
% solution = roots([a ,b, c]);
% disp(solution);
solution = roots([l ,m, n]);
if real(solution(1))>0
    PS=solution(1);
else
    PS=solution(2);
end
% disp(solution);