Question

我要使用以下代码

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)，但无论如何我都无法获得数字。你能帮我吗？

Answer 1

您的代码在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}

无法获得变量gamma1的解，这是怎么回事？

1 个答案: