exn642.m

来自「基天子matlab的双缝干涉从单缝(线光源)发出的单色光射到与之平行的双缝上」· M 代码 · 共 24 行

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第六章例6－4－2程序
% 电位分布及电场计算
% 西安电子科技大学出版社出版  陈怀琛编著 2007年5月
%
fprintf('输入电位分布方程 V(x,y) \n');
fprintf('例如:  log(x.^2 + y.^2) \n');
V = input(': ','s');   % 读入字符串 V(x,y)
NGrid = 20;            % 绘图的网格线数
xMax = 5;              % 绘图区从 x= -xMax 到 x= xMax
xPlot = linspace( -xMax, xMax, NGrid); % 绘图取的x值
%yPlot = xPlot; %linspace( -yMax, yMax, NGrid); % 
[x,y]=meshgrid(xPlot);
VPlot=eval(V);
[ExPlot, EyPlot] = gradient(-VPlot);		%电场是电位的负梯度
clf; subplot(1,2,1),meshc(VPlot);  % 画含等高线的三维曲面
xlabel('x');  ylabel('y'); zlabel('电位');
% 规定等高线图的范围及比例
subplot(1,2,2), axis([-xMax xMax -xMax xMax]);    % 建立第二子图
cs = contour(x,y,VPlot); % 画等高线
clabel(cs); hold on;          % 在等高线图上加上编号
%@ 在等高线图上加上电场方向
quiver(x,y,ExPlot,EyPlot);  % 画电场 E 的箭头图
xlabel('x');  ylabel('y');hold off;