谁能帮我写2个matlab画图程序,物理的,像点电荷的电场线什么的,随便2个就行,急
发布网友
发布时间:2022-04-23 12:30
我来回答
共1个回答
热心网友
时间:2023-10-13 10:51
1. 等量异号点电荷电场中电势分布
[x,y]=meshgrid(-2:0.1:2,-2:0.1:2);
%以0.1为步长建立平面数据网格
z=1./sqrt((x-1).^2+y.^2+0.01)... %写出电势表达式
-1./sqrt((x+1).^2+y.^2+0.01);
[px,py]=gradient(z);
%求电势在x,y方向的梯度即电场强度
contour(x,y,z,[-12,-8,-5,-3,-1,... %画出等势线
-0.5,-0.1,0.1,0.5,1,3,5,8,12])
hold on %作图控制
quiver(x,y,px,py,'k') %画出各点上电场的大小和方向
2. 等量同号点电荷的电场线的绘制
下面是写微分方程的函数文件:
function ydot=dcx1fun(t,y,flag,p1,p2)
%p1,p2是参量,表示电量
ydot=[p1*(y(1)+2)/(sqrt((y(1)+2).^2+y(2).^2).^3)+...
p2*(y(1)-2)/(sqrt((y(1)-2).^2+y(2).^2).^3);
%dx/dt=Ex
p1*y(2)/(sqrt((y(1)+2).^2+y(2).^2).^3)+...
p2*y(2)/(sqrt((y(1)-2).^2+y(2).^2).^3)];
%dy/dt=Ey
编写好函数文件后,命名为dcx1fun.m存在当前路径下,然后开始编写解微分方程的主程序dcx1.m:
p1=10; p2=10; %点电荷所带电量
axis([-5,5,-5,5]); %设定坐标轴范围 -5≤x≤5,-5≤y≤5
hold on %图形控制,不可擦除模式
plot(2,0,'*r'); plot(-2,0,'*r') %绘制两源电荷
a=(pi/24):pi/12:(2*pi-pi/24);
%圆周上电场线起点所对应的角度
b=0.1*cos(a);c=0.1*sin(a);
%电场线起点所对应的相对坐标
b1=-2+b;b2=2+b; %把起点圆周的圆心放置在源电荷处
b0=[b1 b2]; c0=[c c]; %初始条件,所有电场线的起点
%的横、纵坐标构成了矢量b0和c0
for i=1:48 %循环求解48次微分方程
[t,y]=ode45('dcx1fun',[0:0.05:40],[b0(i),c0(i)],[ ],p1,p2);
%调用ode45求解,对应一个初条件(起点),求解出一条电场线
plot(y(:,1),y(:,2),'b') %绘制出此条电场线
end %结束循环,共绘制出48条电场线
参考于《Matlab 在基础物理学中的应用》
http://bnucourse.bnu.e.cn/course/physics/05/jcwlxshyjy.pdf
matlab描绘电场线论文咋写,还有两个点和周围的电场线描绘要立体的好看的...
电偶极子的电场,把我很早以前的一个小程序 clear all;clc;q=2e-6;k=9e9;a=1.5;b=-1.5;x=-5:0.3:5;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.^...
随机(正弦)振动
正弦振动多用于找出产品设计或包装设计的脆弱点。看在哪一个具体频率点响应最大(共振点);正弦振动在任一瞬间只包含一种频率的振动,而随机振动在任一瞬间包含频谱范围内的各种频率的振动。由于随机振动包含频谱内所有的频率,所以样品上的共...
如果两个点电荷电量不等,如何用MATLAB作出他们之间的电场线(请分同号...
正负点电荷 [X,Y] = meshgrid(-2:.2:2);Z = X.*exp(-X.^2 - Y.^2);[DX,DY] = gradient(Z,.2,.2);contour(X,Y,Z)hold on quiver(X,Y,DX,DY)colormap hsv hold on;plot(0.7,0,'w-','MarkerSize',8,'Color','b');plot(0.7,0,'wo','MarkerSize',8,'Color','...
matlab绘制两个点电荷的三维电场线和等势线
x1=-10.5; y1=0;%点电荷1坐标 x2=10.5; y2=0;%点电荷2坐标 r1=sqrt((x-x1).^2+(y-y1).^2); p1=Q1./r1;%电势的表达式 r2=sqrt((x-x2).^2+(y-y2).^2); p2=Q2./r2;z=p1+p2;figure;contour(x,y,z,v);%二维曲面等势线绘图 axis equal;figure;mesh(x,y,z);%...
如何根据公式T=16E(V-E)/V^2*exp[-2a/h √ 2m(V-E)] ,用matlab作图?
a=(pi/24):pi/12:(2*pi-pi/24); %圆周上电场线起点所对应的角度 b=0.1*cos(a);c=0.1*sin(a); %电场线起点所对应的相对坐标 b1=-2+b;b2=2+b; %把起点圆周的圆心放置在源电荷处 b0=[b1 b2]; c0=[c c]; %初始条件,所有电场线的起点 的横、纵坐标构成了矢量b0 和...
如何上好高三物理一轮复习课
2.物理过程的分析要注意细节,要善于找出两个相关过程的连接点(临界点)对于一个复杂的物理问题,首先要根据题目所描述的情景建立正确的物理模型,然后对物理过程进行分析,对于多过程的物理问题,学生一定要注意分析物理过程的细节,弄清各个过程的运动特点及相关联系,找出相关过程之间的物理量之间的关系,做到了这一点,也就...