原程序太长,我先截取了一段重新定义了一下发觉还是求不出来程序如下。
先定义了一些数值:
syms cm11 cm21
f=50
w=6.28*f
N=6
Ng=500
C=220*(10^-9)
Ltp=18*(10^-6)
Udc=25
k1=0.6
k2=0.6
kf=5*N^2
kt=k2*kf^0.5
tn=0.00002
M1p=kt*Ltp
M1s=N*kt*Ltp
M2p=N*k1*Ltp
M2s=N^2*k1*Ltp
Lf=kf*Ltp
Lts=N^2*Ltp
p5=(kf+N*kt+N*kt*k1-kt^2)/(N^2+N*kt+N*kt*k1-N^2*k1^2)
p1=N*kt*C*Ltp
p2=N*kt*k1*C*Ltp
p3=N^2*C*Ltp-N^2*k1^2*C*Ltp
p4=kt+p5*N*k1
w1=220/(p3*p5-p1-p2)
w2=(1-p5)/(p3*p5-p1-p2)
w3=(p4*Udc)/(p3*p5-p1-p2)
T1=-w3/w2
T2=(w*w1)/(w2^0.5(w^2-w2))
T3=w1/(w2-w^2)
T4=w3/w2
t1=(2*N*kt