MATLAB绘制电磁场举例:
clc;close all;clear all;warning off;%清除变量
rand(‘seed’, 100);
randn(‘seed’, 100);
format long g;
m=12
for k=1:m
for j=1:m
if k==1
V(j,k)=1;
elseif((j==1)|(j==m)|(k==m))
V(j,k)=0;
else
V(j,k)=0.5;
end
end
end
cha=0.01;
delta=0;
n=0;
while(1)
n=n+1;
for k=2:m-1
for j=2:m-1
Vnew(j,k)=1/4*(V(j+1,k)+V(j-1,k)+V(j,k+1)+V(j,k-1));
d=abs((Vnew(j,k)-V(j,k))/V(j,k));
if d>delta
delta=d;
end
V(j,k)=Vnew(j,k);
end
end
if delta<cha
break;
end
if(n>100)
break;
end
delta=0.;
end
k=1:m;
j=1:m;
% 绘图
figure;
[DX,DY] = gradient(V,.4,.4);
hold on;
quiver(k,j,DX,DY,2)
hold off;
k=1:m;
j=1:m;
[DX,DY]=gradient(V,.4,.4)
A=(DX.^2+DY.^2).^0.5;
[DA,DB]=gradient(A,.4,.4);
hold on;
quiver(k,j,DA,DB,2)
hold off;
程序结果:
原文地址:https://blog.csdn.net/corn1949/article/details/135930594
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。
如若转载,请注明出处:http://www.7code.cn/show_64679.html
如若内容造成侵权/违法违规/事实不符,请联系代码007邮箱:suwngjj01@126.com进行投诉反馈,一经查实,立即删除!