Matlab绘制电场线向量,为什么全是点?我写的程序function [m] = f(z1,z2)%定义函数,z1,z2分别为q1,q3的z坐标v=-2:0.1:2;[x,y]=meshgrid(v,v);%做出定义域q=1.6*10^(-19);%电子电量,略小ep=8.85*10^(-12);q1=-q;q2=-q;q3=2*q;%代
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/06 12:02:21
![Matlab绘制电场线向量,为什么全是点?我写的程序function [m] = f(z1,z2)%定义函数,z1,z2分别为q1,q3的z坐标v=-2:0.1:2;[x,y]=meshgrid(v,v);%做出定义域q=1.6*10^(-19);%电子电量,略小ep=8.85*10^(-12);q1=-q;q2=-q;q3=2*q;%代](/uploads/image/z/9908166-30-6.jpg?t=Matlab%E7%BB%98%E5%88%B6%E7%94%B5%E5%9C%BA%E7%BA%BF%E5%90%91%E9%87%8F%2C%E4%B8%BA%E4%BB%80%E4%B9%88%E5%85%A8%E6%98%AF%E7%82%B9%3F%E6%88%91%E5%86%99%E7%9A%84%E7%A8%8B%E5%BA%8Ffunction+%5Bm%5D+%3D+f%28z1%2Cz2%29%25%E5%AE%9A%E4%B9%89%E5%87%BD%E6%95%B0%2Cz1%2Cz2%E5%88%86%E5%88%AB%E4%B8%BAq1%2Cq3%E7%9A%84z%E5%9D%90%E6%A0%87v%3D-2%3A0.1%3A2%3B%5Bx%2Cy%5D%3Dmeshgrid%28v%2Cv%29%3B%25%E5%81%9A%E5%87%BA%E5%AE%9A%E4%B9%89%E5%9F%9Fq%3D1.6%2A10%5E%28-19%29%3B%25%E7%94%B5%E5%AD%90%E7%94%B5%E9%87%8F%2C%E7%95%A5%E5%B0%8Fep%3D8.85%2A10%5E%28-12%29%3Bq1%3D-q%3Bq2%3D-q%3Bq3%3D2%2Aq%3B%25%E4%BB%A3)
Matlab绘制电场线向量,为什么全是点?我写的程序function [m] = f(z1,z2)%定义函数,z1,z2分别为q1,q3的z坐标v=-2:0.1:2;[x,y]=meshgrid(v,v);%做出定义域q=1.6*10^(-19);%电子电量,略小ep=8.85*10^(-12);q1=-q;q2=-q;q3=2*q;%代
Matlab绘制电场线向量,为什么全是点?
我写的程序
function [m] = f(z1,z2)%定义函数,z1,z2分别为q1,q3的z坐标
v=-2:0.1:2;
[x,y]=meshgrid(v,v);%做出定义域
q=1.6*10^(-19);%电子电量,略小
ep=8.85*10^(-12);
q1=-q;q2=-q;q3=2*q;%代入题目要求
plot(0,z1,'*',0,0,'*',0,z2,'*');%画出电子位置,用*表示
hold on
phi=4.*pi.*ep.*(q1./(sqrt(x.^2+(y-z1).^2))+q2./(sqrt(x.^2+y.^2))+q3./(sqrt(x.^2+(y-z2).^2)));%代入电势phi的计算公式
[dx,dy]=gradient(phi);%用gradient函数计算phi的梯度
contour(x,y,phi);%画出等势线
quiver(x,y,dx,dy,'MaxHeadSize',1.2);%画出电场向量
colormap hsv
hold off
end%结束函数
为什么没有向量的箭头?
Matlab绘制电场线向量,为什么全是点?我写的程序function [m] = f(z1,z2)%定义函数,z1,z2分别为q1,q3的z坐标v=-2:0.1:2;[x,y]=meshgrid(v,v);%做出定义域q=1.6*10^(-19);%电子电量,略小ep=8.85*10^(-12);q1=-q;q2=-q;q3=2*q;%代
程序没有什么错,只是画线时数据要归一化才能显示出来
function [m] = f(z1,z2)%定义函数,z1,z2分别为q1,q3的z坐标
v=-2:0.1:2;
[x,y]=meshgrid(v,v);%做出定义域
q=1.6*10^(-19);%电子电量,略小
ep=8.85*10^(-12);
q1=-q;q2=-q;q3=2*q;%代入题目要求
plot(0,z1,'*',0,0,'*',0,z2,'*');%画出电子位置,用*表示
hold on
phi=4.*pi.*ep.*(q1./(sqrt(x.^2+(y-z1).^2))+q2./(sqrt(x.^2+y.^2))+q3./(sqrt(x.^2+(y-z2).^2)));%代入电势phi的计算公式
phiM=4*pi*ep*(q1/(sqrt(2^2+(2-z1)^2))+q2/(sqrt(2^2+2^2))+q3/(sqrt(2^2+(2-z2)^2)));%代入电势phi的计算公式
dv = linspace(-phiM,phiM,10);
contour(x,y,phi,dv,'k-');%画出等势线
[dx,dy]=gradient(phi);
d=sqrt(dx.^2+dy.^2);
dx=dx./d;
dy=dy./d;
hold on;
quiver(x,y,dx,dy);
colormap hsv
hold off
end%结束函数
>>f(1,-1)