01010101010.m

来自「一种新的光纤通信调制技术」· M 代码 · 共 21 行

r1=0.2;
r2=0.5455;
for k=1:100
    phi=0.1*k;
   R1(k)=( r1^2+r2^2+2*r1*r2*cos(phi) )/ ( 1+(r1*r2)^2+2*r1*r2*cos(phi));
end

r1=0.5;r2=0.5;
for k=1:100
    phi=0.1*k;
   R2(k)=( r1^2+r2^2+2*r1*r2*cos(phi) )/ ( 1+(r1*r2)^2+2*r1*r2*cos(phi));
end

r1=0.7;r2=0.7;
for k=1:100
    phi=0.1*k;
   R3(k)=( r1^2+r2^2+2*r1*r2*cos(phi) )/ ( 1+(r1*r2)^2+2*r1*r2*cos(phi));
end

phi0=0.1:0.1:10;
plot(phi0,R1,'k-',phi0,R2,'b--',phi0,R3,'r-.');