搜索
 找回密码
 立即注册
查看: 296|回复: 1

电磁学代码分享

[复制链接]
  • TA的每日心情
    开心
    2017-7-7 21:27
  • 签到天数: 8 天

    [LV.3]偶尔看看II

    3

    主题

    12

    帖子

    126

    积分

    版主

    Rank: 7Rank: 7Rank: 7

    积分
    126
    发表于 2016-12-31 23:50:08 | 显示全部楼层 |阅读模式
    本帖最后由 swiminthesky 于 2016-12-28 13:38 编辑

    王梦园
    张萌
    刘伟娜
    赵雅茹
       四个学妹的成果,非常感谢她们。打算把这个分享给大家,希望能给出更好的建议:


    1一个正电荷
    %一个正点电荷的电场线和等势线(请在“创建图形窗口”处设置断点,以观察画图过程)
    clear                                  %清除变量
    xm=2.5;                                %横坐标范围
    ym=2;                                  %纵坐标范围
    x=linspace(-xm,xm,1501);                %横坐标向量
    y=linspace(-ym,ym,1001);                %纵坐标向量
    [X,Y]=meshgrid(x,y);                   %坐标网点(矩阵)
    R1=sqrt(X.^2+Y.^2);                %一个正电荷到场点的距离
    U=1./R1;                         %计算电势
    u=-2:1:6;                            %等势线的电势向量
    figure                                 %创建图形窗口
    C=contour(X,Y,U,u,\'LineWidth\',3);      %画等势线并取等势线的坐标
    [C,h]=contour(X,Y,U,u);                %画等值线并取坐标矩阵和句柄(6.5版用)
    set(h,\'LineWidth\',3)                   %加粗曲线
    clabel(C,\'FontSize\',10)                %标记等势线的值
    hold on                                %保持图像
    plot([-xm;xm],[0;0],[0;0],[-ym;ym])    %画水平和竖直线
    plot(-1,0,\'o\',\'MarkerSize\',12) %画电荷

    [Ex,Ey]=gradient(-U,x,y);%用电势梯度求场强的两个分量
    [Ex,Ey]=gradient(-U);                  %用电势梯度求场强的两个分量
    axis equal tight                       %使坐标间隔相等
    dth=20;                               %电场线角度间隔
    th=(dth:dth:360-dth)*pi/180;           %电场线的起始角度
    r0=0.1;                                %电场线起点半径
    x0=r0*cos(th);                         %电场线的起点横坐标
    y0=r0*sin(th);                         %电场线的起点纵坐标
    streamline(X,Y,Ex,Ey,x0,y0)          %画左边电场线
    grid on                                %加网格
    title(\'一个正点电荷的电场线和等势线\',\'FontSize\',20)%显示标题
    xlabel(\'\\itx/a\',\'FontSize\',16)         %显示横坐标
    ylabel(\'\\ity/a\',\'FontSize\',16)         %显示纵坐标
    text(-xm,ym,\'电势单位:\\itkQ/a\',\'FontSize\',16)%显示电势单位
    Ex=X./R1.^3;          %用公式求场强的x分量
    Ey=Y./R1.^3;                  %用公式求场强的y分量
    h=streamline(X,Y,Ex,Ey,x0,y0);       %重画左边电场线并取句柄
    set(h,\'LineWidth\',4)                   %加粗曲线


    2一个负电荷
    %一个负点电荷的电场线和等势线(请在“创建图形窗口”处设置断点,以观察画图过程)
    clear                                  %清除变量
    xm=2.5;                                %横坐标范围
    ym=2;                                  %纵坐标范围
    x=linspace(-xm,xm,501);                %横坐标向量
    y=linspace(-ym,ym,401);                %纵坐标向量
    [X,Y]=meshgrid(x,y);                   %坐标网点(矩阵)
    R2=sqrt(X.^2+Y.^2);                %一个负电荷到场点的距离
    U=1./R2;                         %计算电势
    u=-4:0.5:4;                            %等势线的电势向量
    figure                                 %创建图形窗口
    C=contour(X,Y,U,u,\'LineWidth\',2);      %画等势线并取等势线的坐标
    [C,h]=contour(X,Y,U,u);                %画等值线并取坐标矩阵和句柄(6.5版用)
    set(h,\'LineWidth\',2)                   %加粗曲线
    clabel(C,\'FontSize\',16)                %标记等势线的值
    hold on                                %保持图像
    plot([-xm;xm],[0;0],[0;0],[-ym;ym])    %画水平和竖直线
    plot(1,0,\'o\',\'MarkerSize\',12) %画电荷

    [Ex,Ey]=gradient(-U,x,y);%用电势梯度求场强的两个分量
    [Ex,Ey]=gradient(-U);                  %用电势梯度求场强的两个分量
    axis equal tight                       %使坐标间隔相等
    dth=20;                                %电场线角度间隔
    th=(dth:dth:360-dth)*pi/180;           %电场线的起始角度
    r0=0.1;                                %电场线起点半径
    x0=r0*cos(th);                         %电场线的起点横坐标
    y0=r0*sin(th);                         %电场线的起点纵坐标
    streamline(X,Y,-Ex,-Ey,x0,y0)        %画电场线
    grid on                                %加网格
    title(\'一个负点电荷的电场线和等势线\',\'FontSize\',20)%显示标题
    xlabel(\'\\itx/a\',\'FontSize\',16)         %显示横坐标
    ylabel(\'\\ity/a\',\'FontSize\',16)         %显示纵坐标
    text(-xm,ym-0.5,\'电势单位:\\itkQ/a\',\'FontSize\',16)%显示电势单位
    Ex=-X./R2.^3;          %用公式求场强的x分量
    Ey=-Y./R2.^3;                  %用公式求场强的y分量
    h=streamline(X,Y,-Ex,-Ey,x0,y0);       %重画电场线并取句柄

    set(h,\'LineWidth\',2)                   %加粗曲线

    3X轴上两个同种正电荷
    clear                                  %清除变量
    xm=2.5;                                %横坐标范围
    ym=2;                                  %纵坐标范围
    x=linspace(-xm,xm,501);                %横坐标向量
    y=linspace(-ym,ym,401);                %纵坐标向量
    [X,Y]=meshgrid(x,y);                   %坐标网点(矩阵)
    R1=sqrt((X+1).^2+Y.^2);                %左边第一个正电荷到场点的距离
    R2=sqrt((X-1).^2+Y.^2);                %右边第二个正电荷到场点的距离
    U=1./R1+1./R2;                         %计算电势
    u=-4:0.5:4;                            %等势线的电势向量
    figure                                 %创建图形窗口
    C=contour(X,Y,U,u,\'LineWidth\',2);      %画等势线并取等势线的坐标
    [C,h]=contour(X,Y,U,u);                %画等值线并取坐标矩阵和句柄(6.5版用)
    set(h,\'LineWidth\',2)                   %加粗曲线
    clabel(C,\'FontSize\',16)                %标记等势线的值
    hold on                                %保持图像
    plot([-xm;xm],[0;0],[0;0],[-ym;ym])    %画水平和竖直线
    plot(-1,0,\'o\',1,0,\'o\',\'MarkerSize\',12) %画电荷

    [Ex,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1));%用电势梯度求场强的两个分量
    [Ex,Ey]=gradient(-U);                  %用电势梯度求场强的两个分量
    axis equal tight                       %使坐标间隔相等
    dth=20;                                %电场线角度间隔
    th=(dth:dth:360-dth)*pi/180;           %电场线的起始角度
    r0=0.1;                                %电场线起点半径
    x0=r0*cos(th);                         %电场线的起点横坐标
    y0=r0*sin(th);                         %电场线的起点纵坐标
    streamline(X,Y,Ex,Ey,x0-1,y0)          %画左边电场线(中间部分达到右边)
    streamline(X,Y,Ex,Ey,x0+1,y0)        %画右边电场线(中间部分达到左边)
    grid on                                %加网格
    title(\'X轴上等量同号正点电荷的电场线和等势线\',\'FontSize\',20)%显示标题
    xlabel(\'\\itx/a\',\'FontSize\',16)         %显示横坐标
    ylabel(\'\\ity/a\',\'FontSize\',16)         %显示纵坐标
    text(-xm,ym-0.5,\'电势单位:\\itkQ/a\',\'FontSize\',16)%显示电势单位
    Ex=(X+1)./R1.^3+(X-1)./R2.^3;          %用公式求场强的x分量
    Ey=Y./R1.^3+Y./R2.^3;                  %用公式求场强的y分量
    h=streamline(X,Y,Ex,Ey,x0-1,y0);       %重画左边电场线并取句柄
    set(h,\'LineWidth\',2)                   %加粗曲线
    h=streamline(X,Y,Ex,Ey,x0+1,y0);     %重画右边电场线并取句柄

    set(h,\'LineWidth\',2)                   %加粗曲线

    4.x轴上等量同号负点电荷
    clear                                  %清除变量
    xm=2.5;                                %横坐标范围
    ym=2;                                  %纵坐标范围
    x=linspace(-xm,xm,501);                %横坐标向量
    y=linspace(-ym,ym,401);                %纵坐标向量
    [X,Y]=meshgrid(x,y);                   %坐标网点(矩阵)
    R1=sqrt((X+1).^2+Y.^2);                %左边第一个正电荷到场点的距离
    R2=sqrt((X-1).^2+Y.^2);                %右边第二个正电荷到场点的距离
    U=1./R1+1./R2;                         %计算电势
    u=-4:0.5:4;                            %等势线的电势向量
    figure                                 %创建图形窗口
    C=contour(X,Y,U,u,\'LineWidth\',2);      %画等势线并取等势线的坐标
    [C,h]=contour(X,Y,U,u);                %画等值线并取坐标矩阵和句柄(6.5版用)
    set(h,\'LineWidth\',2)                   %加粗曲线
    clabel(C,\'FontSize\',16)                %标记等势线的值
    hold on                                %保持图像
    plot([-xm;xm],[0;0],[0;0],[-ym;ym])    %画水平和竖直线
    plot(-1,0,\'o\',1,0,\'o\',\'MarkerSize\',12) %画电荷

    [Ex,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1));%用电势梯度求场强的两个分量
    [Ex,Ey]=gradient(-U);                  %用电势梯度求场强的两个分量
    axis equal tight                       %使坐标间隔相等
    dth=20;                                %电场线角度间隔
    th=(dth:dth:360-dth)*pi/180;           %电场线的起始角度
    r0=0.1;                                %电场线起点半径
    x0=r0*cos(th);                         %电场线的起点横坐标
    y0=r0*sin(th);                         %电场线的起点纵坐标
    streamline(X,Y,Ex,Ey,x0-1,y0)          %画左边电场线(中间部分达到右边)
    streamline(X,Y,-Ex,-Ey,x0+1,y0)        %画右边电场线(中间部分达到左边)
    grid on                                %加网格
    title(\'X轴上等量同号负点电荷的电场线和等势线\',\'FontSize\',20)%显示标题
    xlabel(\'\\itx/a\',\'FontSize\',16)         %显示横坐标
    ylabel(\'\\ity/a\',\'FontSize\',16)         %显示纵坐标
    text(-xm,ym-0.5,\'电势单位:\\itkQ/a\',\'FontSize\',16)%显示电势单位
    Ex=(X+1)./R1.^3+(X-1)./R2.^3;          %用公式求场强的x分量
    Ey=Y./R1.^3+Y./R2.^3;                  %用公式求场强的y分量
    h=streamline(X,Y,Ex,Ey,x0-1,y0);       %重画左边电场线并取句柄
    set(h,\'LineWidth\',2)                   %加粗曲线
    h=streamline(X,Y,-Ex,-Ey,x0+1,y0);     %重画右边电场线并取句柄

    set(h,\'LineWidth\',2)                   %加粗曲线

    5两个异号等量电荷
    clear                                  %清除变量
    xm=2.5;                                %横坐标范围
    ym=2;                                  %纵坐标范围
    x=linspace(-xm,xm,400);                %横坐标向量
    y=linspace(-ym,ym,400);                %纵坐标向量
    [X,Y]=meshgrid(x,y);                   %坐标网点(矩阵)
    R1=sqrt((X+1).^2+Y.^2);                %左边第一个正电荷到场点的距离
    R2=sqrt((X-1).^2+Y.^2);                %右边第二个负电荷到场点的距离
    U=1./R1-1./R2;                         %计算电势
    u=-4:0.5:4;                            %等势线的电势向量
    figure                                 %创建图形窗口
    C=contour(X,Y,U,u,\'LineWidth\',2);      %画等势线并取等势线的坐标
    [C,h]=contour(X,Y,U,u);                %画等值线并取坐标矩阵和句柄(6.5版用)
    set(h,\'LineWidth\',2)                   %加粗曲线
    clabel(C,\'FontSize\',16)                %标记等势线的值
    hold on                                %保持图像
    plot([-xm;xm],[0;0],[0;0],[-ym;ym])    %画水平和竖直线
    plot(-1,0,\'o\',1,0,\'o\',\'MarkerSize\',12) %画电荷

    [Ex,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1));%用电势梯度求场强的两个分量
    [Ex,Ey]=gradient(-U);                  %用电势梯度求场强的两个分量
    axis equal tight                       %使坐标间隔相等
    dth=20;                                %电场线角度间隔
    th=(dth:dth:360-dth)*pi/180;           %电场线的起始角度
    r0=0.1;                                %电场线起点半径
    x0=r0*cos(th);                         %电场线的起点横坐标
    y0=r0*sin(th);                         %电场线的起点纵坐标
    streamline(X,Y,Ex,Ey,x0-1,y0)          %画左边电场线(中间部分达到右边)
    streamline(X,Y,-Ex,-Ey,x0+1,y0)        %画右边电场线(中间部分达到左边)
    grid on                                %加网格
    title(\'等量异号点电荷的电场线和等势线\',\'FontSize\',20)%显示标题
    xlabel(\'\\itx/a\',\'FontSize\',16)         %显示横坐标
    ylabel(\'\\ity/a\',\'FontSize\',16)         %显示纵坐标
    text(-xm,ym-0.5,\'电势单位:\\itkQ/a\',\'FontSize\',16)%显示电势单位
    Ex=(X+1)./R1.^3-(X-1)./R2.^3;          %用公式求场强的x分量
    Ey=Y./R1.^3-Y./R2.^3;                  %用公式求场强的y分量
    h=streamline(X,Y,Ex,Ey,x0-1,y0);       %重画左边电场线并取句柄
    set(h,\'LineWidth\',2)                   %加粗曲线
    h=streamline(X,Y,-Ex,-Ey,x0+1,y0);     %重画右边电场线并取句柄
    set(h,\'LineWidth\',2)                   %加粗曲线


    6.y轴上等量同号正点电荷
    clear                                  %清除变量
    xm=2.5;                                %横坐标范围
    ym=2;                                  %纵坐标范围
    x=linspace(-xm,xm,400);                %横坐标向量
    y=linspace(-ym,ym,400);                %纵坐标向量
    [X,Y]=meshgrid(x,y);                   %坐标网点(矩阵)
    R1=sqrt(X.^2+(Y+1).^2);                %下边第一个正电荷到场点的距离
    R2=sqrt(X.^2+(Y-1).^2);                %上边第一个正电荷到场点的距离
    U=1./R1+1./R2;                         %计算电势
    u=-4:0.5:4;                            %等势线的电势向量
    figure                                 %创建图形窗口
    C=contour(X,Y,U,u,\'LineWidth\',2);      %画等势线并取等势线的坐标
    %[C,h]=contour(X,Y,U,u);                %画等值线并取坐标矩阵和句柄(6.5版用)
    %set(h,\'LineWidth\',2)                   %加粗曲线
    clabel(C,\'FontSize\',16)                %标记等势线的值
    hold on                                %保持图像
    plot([-xm;xm],[0;0],[0;0],[-ym;ym])    %画水平和竖直线
    plot(0,-1,\'o\',0,1,\'o\',\'MarkerSize\',12) %画电荷

    [Ex,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1));%用电势梯度求场强的两个分量
    %[Ex,Ey]=gradient(-U);                  %用电势梯度求场强的两个分量
    axis equal tight                       %使坐标间隔相等
    dth=20;                                %电场线角度间隔
    th=(dth:dth:360-dth)*pi/180;           %电场线的起始角度
    r0=0.1;                                %电场线起点半径
    x0=r0*cos(th);                         %电场线的起点横坐标
    y0=r0*sin(th);                         %电场线的起点纵坐标
    streamline(X,Y,Ex,Ey,x0,y0+1)          %画下边电场线(中间部分达到右边)
    streamline(X,Y,Ex,Ey,x0,y0-1)        %画上边电场线(中间部分达到左边)
    grid on                                %加网格
    title(\'Y轴上等量同号正点电荷的电场线和等势线\',\'FontSize\',20)%显示标题
    xlabel(\'\\itx/a\',\'FontSize\',16)         %显示横坐标
    ylabel(\'\\ity/a\',\'FontSize\',16)         %显示纵坐标
    text(-xm,ym-0.5,\'电势单位:\\itkQ/a\',\'FontSize\',16)%显示电势单位
    Ex=X./R1.^3+X./R2.^3;          %用公式求场强的x分量
    Ey=(Y+1)./R1.^3+(Y-1)./R2.^3;                  %用公式求场强的y分量
    h=streamline(X,Y,Ex,Ey,x0,y0-1);       %重画上左边电场线并取句柄
    set(h,\'LineWidth\',2)                   %加粗曲线
    h=streamline(X,Y,Ex,Ey,x0,y0+1);     %重画下边电场线并取句柄
    set(h,\'LineWidth\',2)                   %加粗曲线

    7y轴上等量同种负电荷
    clear                                  %清除变量
    xm=2.5;                                %横坐标范围
    ym=2;                                 %纵坐标范围
    x=linspace(-xm,xm,400);                %横坐标向量
    y=linspace(-ym,ym,400);                %纵坐标向量
    [X,Y]=meshgrid(x,y);                   %坐标网点(矩阵)
    R1=sqrt(X.^2+(Y+1).^2);                %下边第一个正电荷到场点的距离
    R2=sqrt(X.^2+(Y-1).^2);                %上边第一个正电荷到场点的距离
    U=1./R1+1./R2;                         %计算电势
    u=-4:0.5:4;                            %等势线的电势向量
    figure                                 %创建图形窗口
    C=contour(X,Y,U,u,\'LineWidth\',2);      %画等势线并取等势线的坐标
    %[C,h]=contour(X,Y,U,u);                %画等值线并取坐标矩阵和句柄(6.5版用)
    %set(h,\'LineWidth\',2)                   %加粗曲线
    clabel(C,\'FontSize\',16)                %标记等势线的值
    hold on                                %保持图像
    plot([-xm;xm],[0;0],[0;0],[-ym;ym])    %画水平和竖直线
    plot(0,-1,\'o\',0,1,\'o\',\'MarkerSize\',12) %画电荷

    [Ex,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1));%用电势梯度求场强的两个分量
    %[Ex,Ey]=gradient(-U);                  %用电势梯度求场强的两个分量
    axis equal tight                       %使坐标间隔相等
    dth=20;                                %电场线角度间隔
    th=(dth:dth:360-dth)*pi/180;           %电场线的起始角度
    r0=0.1;                                %电场线起点半径
    x0=r0*cos(th);                         %电场线的起点横坐标
    y0=r0*sin(th);                         %电场线的起点纵坐标
    streamline(X,Y,Ex,Ey,x0,y0+1)          %画下边电场线(中间部分达到右边)
    streamline(X,Y,Ex,Ey,x0,y0-1)        %画上边电场线(中间部分达到左边)
    grid on                                %加网格
    title(\'Y轴上等量同号正点电荷的电场线和等势线\',\'FontSize\',20)%显示标题
    xlabel(\'\\itx/a\',\'FontSize\',16)         %显示横坐标
    ylabel(\'\\ity/a\',\'FontSize\',16)         %显示纵坐标
    text(-xm,ym-0.5,\'电势单位:\\itkQ/a\',\'FontSize\',16)%显示电势单位
    Ex=X./R1.^3+X./R2.^3;          %用公式求场强的x分量
    Ey=(Y+1)./R1.^3+(Y-1)./R2.^3;                  %用公式求场强的y分量
    h=streamline(X,Y,Ex,Ey,x0,y0-1);       %重画上左边电场线并取句柄
    set(h,\'LineWidth\',2)                   %加粗曲线
    h=streamline(X,Y,Ex,Ey,x0,y0+1);     %重画下边电场线并取句柄
    set(h,\'LineWidth\',2)                   %加粗曲线








    %三个等量正点电荷的电场线和等势线(请在“创建图形窗口”处设置断点,以观察画图过程)
    clear                                  %清除变量
    xm=2.5;                                %横坐标范围
    ym=2;                                  %纵坐标范围
    x=linspace(-xm,xm,400);                %横坐标向量
    y=linspace(-ym,ym,400);                %纵坐标向量
    [X,Y]=meshgrid(x,y);                   %坐标网点(矩阵)
    R1=sqrt((X+1).^2+Y.^2);                %左边第一个正电荷到场点的距离
    R2=sqrt((X-1).^2+Y.^2);                %右边第二个正电荷到场点的距离
    R3=sqrt(X.^2+(Y+1).^2);                %下边第三个正电荷到场点的距离
    U=1./R1+1./R2+1./R3;                   %计算电势
    u=-4:0.5:4;                            %等势线的电势向量
    figure                                 %创建图形窗口
    C=contour(X,Y,U,u,\'LineWidth\',2);      %画等势线并取等势线的坐标
    [C,h]=contour(X,Y,U,u);                %画等值线并取坐标矩阵和句柄(6.5版用)
    set(h,\'LineWidth\',2)                   %加粗曲线
    clabel(C,\'FontSize\',16)                %标记等势线的值
    hold on                                %保持图像
    plot([-xm;xm],[0;0],[0;0],[-ym;ym])    %画水平和竖直线
    plot(-1,0,\'o\',1,0,\'o\',0,-1,\'o\',\'MarkerSize\',12) %画电荷
    [Ex,Ey]=gradient(-U);                  %用电势梯度求场强的两个分量
    axis equal tight                       %使坐标间隔相等
    dth=20;                                %电场线角度间隔
    th=(dth:dth:360-dth)*pi/180;           %电场线的起始角度
    r0=0.1;                                %电场线起点半径
    x0=r0*cos(th);                         %电场线的起点横坐标
    y0=r0*sin(th);                         %电场线的起点纵坐标
    streamline(X,Y,Ex,Ey,x0-1,y0)          %画左边电场线
    streamline(X,Y,Ex,Ey,x0+1,y0)          %画右边电场线
    streamline(X,Y,Ex,Ey,x0,y0-1)          %画下边电场线
    grid on                                %加网格
    title(\'三个等量正点电荷的电场线和等势线\',\'FontSize\',20)%显示标题
    xlabel(\'\\itx/a\',\'FontSize\',16)         %显示横坐标
    ylabel(\'\\ity/a\',\'FontSize\',16)         %显示纵坐标
    text(-xm,ym-0.5,\'电势单位:\\itkQ/a\',\'FontSize\',16)%显示电势%三个等量正点电荷的电场线和等势线(请在“创建图形窗口”处设置断点,以观察画图过程)
    clear                                  %清除变量
    xm=2.5;                                %横坐标范围
    ym=2;                                  %纵坐标范围
    x=linspace(-xm,xm,400);                %横坐标向量
    y=linspace(-ym,ym,400);                %纵坐标向量
    [X,Y]=meshgrid(x,y);                   %坐标网点(矩阵)
    R1=sqrt((X+1).^2+Y.^2);                %左边第一个正电荷到场点的距离
    R2=sqrt((X-1).^2+Y.^2);                %右边第二个正电荷到场点的距离
    R3=sqrt(X.^2+(Y+1).^2);                %下边第三个正电荷到场点的距离
    U=1./R1+1./R2+1./R3;                   %计算电势
    u=-4:0.5:4;                            %等势线的电势向量
    figure                                 %创建图形窗口
    C=contour(X,Y,U,u,\'LineWidth\',2);      %画等势线并取等势线的坐标
    [C,h]=contour(X,Y,U,u);                %画等值线并取坐标矩阵和句柄(6.5版用)
    set(h,\'LineWidth\',2)                   %加粗曲线
    clabel(C,\'FontSize\',16)                %标记等势线的值
    hold on                                %保持图像
    plot([-xm;xm],[0;0],[0;0],[-ym;ym])    %画水平和竖直线
    plot(-1,0,\'o\',1,0,\'o\',0,-1,\'o\',\'MarkerSize\',12) %画电荷
    [Ex,Ey]=gradient(-U);                  %用电势梯度求场强的两个分量
    axis equal tight                       %使坐标间隔相等
    dth=20;                                %电场线角度间隔
    th=(dth:dth:360-dth)*pi/180;           %电场线的起始角度
    r0=0.1;                                %电场线起点半径
    x0=r0*cos(th);                         %电场线的起点横坐标
    y0=r0*sin(th);                         %电场线的起点纵坐标
    streamline(X,Y,Ex,Ey,x0-1,y0)          %画左边电场线
    streamline(X,Y,Ex,Ey,x0+1,y0)          %画右边电场线
    streamline(X,Y,Ex,Ey,x0,y0-1)          %画下边电场线
    grid on                                %加网格
    title(\'三个等量正点电荷的电场线和等势线\',\'FontSize\',20)%显示标题
    xlabel(\'\\itx/a\',\'FontSize\',16)         %显示横坐标
    ylabel(\'\\ity/a\',\'FontSize\',16)         %显示纵坐标
    text(-xm,ym-0.5,\'电势单位:\\itkQ/a\',\'FontSize\',16)%显示电势单位单位   
             




    来源: 电磁学代码分享
  • TA的每日心情
    开心
    2017-12-4 16:57
  • 签到天数: 10 天

    [LV.3]偶尔看看II

    0

    主题

    15

    帖子

    120

    积分

    注册会员

    Rank: 2

    积分
    120
    发表于 2017-1-19 15:27:42 | 显示全部楼层
    :handshake
    回复

    使用道具 举报

    发表回复

    您需要登录后才可以回帖 登录 | 立即注册

    本版积分规则

    1