xcty2021 2022-05-27 17:42 采纳率: 66.7%
浏览 23
已结题

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

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

img

[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 为什么会出现这个问题,该怎么解决呀?

img

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; 

  • 写回答

1条回答 默认 最新

  • 胸毛男 2022-05-27 19:22
    关注

    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; 
     
    

    img

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 系统已结题 6月4日
  • 已采纳回答 5月27日
  • 创建了问题 5月27日

悬赏问题

  • ¥15 企业资源规划ERP沙盘模拟
  • ¥15 树莓派控制机械臂传输命令报错,显示摄像头不存在
  • ¥15 前端echarts坐标轴问题
  • ¥15 CMFCPropertyPage
  • ¥15 ad5933的I2C
  • ¥15 请问RTX4060的笔记本电脑可以训练yolov5模型吗?
  • ¥15 数学建模求思路及代码
  • ¥50 silvaco GaN HEMT有栅极场板的击穿电压仿真问题
  • ¥15 谁会P4语言啊,我想请教一下
  • ¥15 这个怎么改成直流激励源给加热电阻提供5a电流呀