Question

matlab 复杂的复数取实部

这是我的m文件
clc
clear all
syms a p real

C=0.0001*10^(-6)
w=2*pi*5*10^(6)
u=4*pi*10^(-7)
LP=0*10^(-3)

L=10*10^(-3)
A2=1.0*10^(6)
u2=1
Rdc=1.7
A1=0*10^(6)

B1=w*u*A1
B2=w*u*A2*u2

a1=sqrt(a^2+B1*i)
a2=sqrt(a^2+B2*i)

%%%Coils of finite cross section

R1=6.545*10^(-3)
L1=L+0.2*10^(-3)
R2=7.0*10^(-3)
L2=L1+0.53*10^(-3)
N=10

%%%%%%%%%
J1=besselj(1,p)
I=int(p*J1,p,a*R1,a*R2)

k=a1-a2/u2
r1=a+a1
s=a2/u2+a1
M=(r1*k+(a-a1)*s*exp(2*C*a1))/((a-a1)*k+r1*s*exp(2*C*a1))*(exp(-2*a*(L1+LP))+exp(-2*a*(L1+LP))-2*exp(-a*(L2+L1+2*LP)))
G=2*exp(-a*(L2-L1))-2
simple([real(M),imag(M)])

ReZ1=int((1/a^5)*(I^2)*(2*(L2-L1)+(G+real(M))/a),a,0,100)

其中I的定积分解不出来怎么办，还有M分不开实部和虚部，最后导致结果ReZ1算不出.

请高手帮忙啊！！！！

Answer 1

不要声明syms real, change'real ' into other name of your syms , pl.: real1 . Cuz real used for MATLAB function real() , you know.