matlab 矢量图不出现；解方程得出的变量值不能继续参与梯度运算

Question

matlab 为什么不显示矢量图呀？

[X,Y] = meshgrid(-6:0.1:6,0:0.1:15);    %设置绘制区域的网格
h=7.5;R=0.011;s=2;                      %对地高度选择7.5米,双回线路两导线间距4米,线径11毫米，电压取110
lamda=110000/(18*10^9*log(h/R));        %求出每根输电线的线密度
r1=sqrt((X+s).^2 +(Y-h).^2);            %r1是空间中一点到左侧输电线中心轴的距离
r2=sqrt((X-s).^2 +(Y-h).^2);            %r2是空间中一点到右侧输电线中心轴的距离
r3=sqrt((X-s).^2 +(Y-h).^2);            %r3是空间中一点到左侧镜像输电线中心轴的距离
r4=sqrt((X-s).^2 +(Y-h).^2);            %r4是空间中一点到右侧镜像输电线中心轴的距离
phi1=(18*lamda*10^9)*log(r3./r1);      %给出左侧输电线电势的表达式
phi2=(18*lamda*10^9)*log(r4./r2);      %给出右侧输电线电势的表达式
phi=phi1+phi2;                          %给出总输电线电势的表达式
[EX,EY] = gradient(-phi);               %利用数值方法求电势的负梯度，得到电场强度
contour(X,Y,phi,25);                       %绘制等势线
hold on; 
quiver(X,Y,EX,EY);
axis image;

matlab 为什么会出现这个问题，该怎么解决呀？

syms lamda;
[X,Y] = meshgrid(-3:0.4:3,-3:0.4:3);    %设置绘制区域的网格
h=7.5;R=0.011;
lamda=vpasolve([18*10^9*log(h/R)*lamda==110000],lamda);
r1=sqrt(X.^2 +(Y-h).^2);%r1是空间中一点到输电线中心轴的距离
r2=sqrt(X.^2 +(Y+h).^2);%r2是空间中一点到镜像输电线中心轴的距离
phi = (18*lamda*10^9).*log(r2./r1); %给出电势的表达式，忽略系数
[EX,EY] = gradient(-phi);         %利用数值方法求电势的负梯度，得到电场强度
contour(X,Y,phi);             %绘制等势线
hold on;                                 
quiver(X,Y,EX,EY);                   
axis image; 

Answer 1

phi是个sym 要转化成矩阵，用double()

syms lamda;
[X,Y] = meshgrid(-3:0.4:3,-3:0.4:3);    %设置绘制区域的网格
h=7.5;R=0.011;
lamda=vpasolve([18*10^9*log(h/R)*lamda==110000],lamda);
r1=sqrt(X.^2 +(Y-h).^2);%r1是空间中一点到输电线中心轴的距离
r2=sqrt(X.^2 +(Y+h).^2);%r2是空间中一点到镜像输电线中心轴的距离
phi = (18*lamda*10^9).*log(r2./r1); %给出电势的表达式，忽略系数
[EX,EY] = gradient(-double(phi));         %利用数值方法求电势的负梯度，得到电场强度
contour(X,Y,phi);             %绘制等势线
hold on;                                 
quiver(X,Y,EX,EY);                   
axis image;