高分求用MATLAB编译的高数或大学物理题(40~50行左右)
来源:学生作业帮 编辑:搜搜考试网作业帮 分类:综合作业 时间:2024/06/26 06:05:35
高分求用MATLAB编译的高数或大学物理题(40~50行左右)
![高分求用MATLAB编译的高数或大学物理题(40~50行左右)](/uploads/image/z/16334718-6-8.jpg?t=%E9%AB%98%E5%88%86%E6%B1%82%E7%94%A8MATLAB%E7%BC%96%E8%AF%91%E7%9A%84%E9%AB%98%E6%95%B0%E6%88%96%E5%A4%A7%E5%AD%A6%E7%89%A9%E7%90%86%E9%A2%98%EF%BC%8840%7E50%E8%A1%8C%E5%B7%A6%E5%8F%B3%EF%BC%89)
】研究偶极子(Dipole)的电势(Electric potential)和电场强度(Electric field density)
clear;
clf;
q=2e-6;
k=9e9;
a=1.5;
b=-1.5;
x=-6:0.6:6;
y=x;
[X,Y]=meshgrid(x,y);
rp=sqrt((X-a).^2+(Y-b).^2);
rm=sqrt((X+a).^2+(Y+b).^2);
V=q*k*(1./rp-1./rm);
[Ex,Ey]=gradient(-V);
AE=sqrt(Ex.^2+Ey.^2);
Ex=Ex./AE;Ey=Ey./AE;
cv=linspace(min(min(V)),max(max(V)),49);
contourf(X,Y,V,cv,'k-')
%axis('square')
title('\fontname{隶书}\fontsize{22}偶极子的场'),hold on
quiver(X,Y,Ex,Ey,0.7)
plot(a,b,'wo',a,b,'w+')
plot(-a,-b,'wo',-a,-b,'w-')
xlabel('x');ylabel('y'),hold off
clear;
clf;
q=2e-6;
k=9e9;
a=1.5;
b=-1.5;
x=-6:0.6:6;
y=x;
[X,Y]=meshgrid(x,y);
rp=sqrt((X-a).^2+(Y-b).^2);
rm=sqrt((X+a).^2+(Y+b).^2);
V=q*k*(1./rp-1./rm);
[Ex,Ey]=gradient(-V);
AE=sqrt(Ex.^2+Ey.^2);
Ex=Ex./AE;Ey=Ey./AE;
cv=linspace(min(min(V)),max(max(V)),49);
contourf(X,Y,V,cv,'k-')
%axis('square')
title('\fontname{隶书}\fontsize{22}偶极子的场'),hold on
quiver(X,Y,Ex,Ey,0.7)
plot(a,b,'wo',a,b,'w+')
plot(-a,-b,'wo',-a,-b,'w-')
xlabel('x');ylabel('y'),hold off