无法获得变量gamma1的解,这是怎么回事?

时间:2019-10-17 13:19:03

标签: maple

我要使用以下代码

f := sqrt(7200*8200)*10^6: 
Ig := 0, 9330; Ug := 13.134: 
f_gr := 10^10:
r_e := .8: 
L_e := .5*10^(-9): 
C_e := 2.5*10^(-12): 
Cka := .4*10^(-12): 
Ckp := .8*10^(-12): 
L_b := .3*10^(-9): 
I_0 := I*f*Ig*(cos(.4*f/f_gr)+I*sin(.4*f/f_gr))/(f_gr*gamma1):
I_e := I_0+Ig:
U_e := I_e*(r_e+I*(2*Pi*f*L_e)):
U_p := -I*(1-gamma1)*I_0/(2*Pi*f*C_e):
Uska := Ug+U_p:
Iska := I*(2*Pi*f*Cka*Uska):
Irb := I_0+Iska:
Urb := 2.8*Irb:
Uskp := Urb+Uska:
Iskp := I*(2*Pi*f*Ckp*Uskp):
Irk := Uskp/773.5:
Ib := Irb+Iskp+Irk:
Ulb := I*(2*Pi*f*L_b*Ib):
Uv := U_e+Ulb+Urb+U_p:
Ik := Ig-Iska-Iskp-Irk:
Uk := Ug-U_e:

求解gamma1的Re(Uk)*Re(Ik)+Im(Uk)*Im(Ik)=0.186(我代表虚数单位)。但这是行不通的。我尝试使用蛮力,将一些数字分配给gamma1,然后找到Re(Uk)*Re(Ik)+Im(Uk)*Im(Ik),但无论如何我都无法获得数字。你能帮我吗?

1 个答案:

答案 0 :(得分:2)

您的代码在Ig的定义中有明显的错误,您已为其分配了表达式序列0, 9330

也许您打算使用Ig:=0.9330

restart;
f := sqrt(7200*8200)*10^6: 
Ig := 0.9330:
Ug := 13.134: 
f_gr := 10^10:
r_e := .8: 
L_e := .5*10^(-9): 
C_e := 2.5*10^(-12): 
Cka := .4*10^(-12): 
Ckp := .8*10^(-12): 
L_b := .3*10^(-9): 
I_0 := I*f*Ig*(cos(.4*f/f_gr)+I*sin(.4*f/f_gr))/(f_gr*gamma1):
I_e := I_0+Ig:
U_e := I_e*(r_e+I*(2*Pi*f*L_e)):
U_p := -I*(1-gamma1)*I_0/(2*Pi*f*C_e):
Uska := Ug+U_p:
Iska := I*(2*Pi*f*Cka*Uska):
Irb := I_0+Iska:
Urb := 2.8*Irb:
Uskp := Urb+Uska:
Iskp := I*(2*Pi*f*Ckp*Uskp):
Irk := Uskp/773.5:
Ib := Irb+Iskp+Irk:
Ulb := I*(2*Pi*f*L_b*Ib):
Uv := U_e+Ulb+Urb+U_p:
Ik := Ig-Iska-Iskp-Irk:
Uk := Ug-U_e:

现在,

solve(Re(Uk)*Re(Ik)+Im(Uk)*Im(Ik)=0.186, {gamma1});

       {gamma1 = -0.9296926790 - 0.6316639400 I}