先谢了哈
global steps L R Q A B C D E W
lamda=input('lamda=');
L=input('L=');
R=input('R=');
Q=input('Q=');
A=1-2*L/Q;
>> B=2*L*(1-L/Q);
>> C=-[2/R+2/Q*(1-2*L/R)];
>> D=-[2*L/R-(1-2*L/R)*(1-2*L/Q)];
>>E=B.*(D-A)./2*(1-A.*D);
>> W=[lamda*B*[4-(A+D)^2]^1/2/2*pi*(1-A*D)]^1/2;
>> ezplot('w^2/W^2-e^2/E^2=1',[-2,2,-1.5,1.5])
global steps L R Q A B C D E W
lamda=input('lamda=');
L=input('L=');
R=input('R=');
Q=input('Q=');
A=1-2*L/Q;
>> B=2*L*(1-L/Q);
>> C=-[2/R+2/Q*(1-2*L/R)];
>> D=-[2*L/R-(1-2*L/R)*(1-2*L/Q)];
>>E=B.*(D-A)./2*(1-A.*D);
>> W=[lamda*B*[4-(A+D)^2]^1/2/2*pi*(1-A*D)]^1/2;
>> ezplot('w^2/W^2-e^2/E^2=1',[-2,2,-1.5,1.5])