Question

matlab程序

syms c1 c2 a c3 r z;
w=a^2*( c1*(0.0065^2-r^2)^2 +c2*( 0.0065^2-r^2)^3)
dw1=diff(w,'r')
dw2=diff(w,'a')
dw11=diff(dw1,'r')
dw22=diff(dw2,'a')
dw12=diff(dw1,'a')
X1=- dw11*z
X2=- (dw1/r)*z+ (dw22/(r*r))*z
X6=(dw12/r+dw12-dw2/r)*z
X=[X1;X2;0;0;0;X6]
syms Y1 Y2 Y3 Y4 Y5 Y6;
D=[86.74,-8.6,27.497,1.058,0,0
-8.6,130.74,-4.8,-1.850,0,0
27.497,-4.8,96.62,13.45,0,0
1.058,-1.85,13.45,45.39,0,0
0,0,0,0,67.47,7.61
0,0,0,0,7.61,30.35]*10^9
Y=D*X
v=1/2*(X1*Y(1)+X2*Y(2)+X6*Y(6))
v1=int(v,z,-0.0002,0.0002)
v2=int(v1,a,0,2*pi)
V=int(v2,r,0,0.2)+c1*10^5*int(int(a^2*(0.0065^2-r^2)^2,r,0,0.2),a,0,2*pi)+c2*10^5*int(int(a^2*(0.0065^2-r^2)^3,r,0,0.2),a,0,2*pi)
E1=diff(V,c1)
E2=diff(V,c2)
[c1,c2]=solve('E1=0','E2=0')

运行完会得到

E1 =

piecewise([36893488147419103232*Im(c1) + 1558749874228457*Im(c2) = 0 and 36893488147419103232*Re(c1) + 1558749874228457*Re(c2) = 0, (47288355635982428642814465285912802442923095608451552988320968500784650061211171682417274775803*pi*c1)/13339248680138825712795377473000846984696702233481792977345513677711592534520627200000000000000000。。。。。。。。。很长一串。。。这个结果是什么意思？

这样解出的c1c2都是零，结果不对，

Answer 1