西电电磁场大作业

最新推荐文章于 2024-07-09 22:51:38 发布

weijunshaomanhh

最新推荐文章于 2024-07-09 22:51:38 发布

阅读量457

点赞数 11

文章标签： matlab 信息与通信

本文链接：https://blog.csdn.net/m0_56162460/article/details/134895290

版权

具体作业题目就是标题

梯度的仿真

x=linspace(-2, 2, 10);          %在-2到2取25个点
y=linspace(-2, 2, 10);          %在-2到2取25个点
z=linspace(-2, 2, 10);
[xx, yy, zz]=meshgrid(x, y, z);     %生成网格采样点
r = sqrt( xx.^2 + yy.^2 +zz.^2);        %生成矩阵Z
% h=contour(xx,yy,zz, 12);        %以12个等高线层级绘制矩阵Z的等高线图
% clabel(h);                        %写等高线的值
[dx, dy, dz]=gradient(r,0.2,0.2,0.2);       %求梯度
grid on
hold on                         %做下一幅图时保持原来图像
quiver3(xx,yy,zz,dx,dy,dz);           %画矢量图箭头
axis equal;                     %等比例显示

散度（二维）的仿真

syms x y z real                         %定义符号变量
F = [ cos(x+2*y), sin(x-2*y) ];         %定义函数F
g = divergence(F,[x y]);                 %求函数F的散度，符号形式
divF=matlabFunction(g);                 %将散度转换为函数形式
x = linspace(-2.5,2.5,20);
y = linspace(-2.5,2.5,20);
[X,Y]=meshgrid(x,y);                    %定义网格
Fx=cos(X+2*Y);                          %F的x分量
Fy=sin(X-2*Y);                          %F的y分量
div_num=divF(X,Y);                      %散度的数值形式
pcolor(X,Y,div_num);                    %绘制散度
shading interp;                         %差值
colorbar;                               %绘制色条
hold on;                                %保持绘图模式打开
quiver(X,Y,Fx,Fy,'k','linewidth',1);    %绘制箭头图

旋度（二维）的仿真

syms x y z real                         %定义符号变量
F = [ cos(x+2*y), sin(x-2*y) ];         %定义函数F
G = curl([F,0],[x y z]);                 %计算F的旋度，并赋予G       
curlF=matlabFunction(G(3));             %将G的z分量，赋予curlF
x=linspace(-2.5,2.5,20);    
[X,Y]=meshgrid(x,x);                    %定义网格
Fx=cos(X+2*Y);                          %计算F的x分量
Fy=sin(X-2*Y);                          %计算F的y分量
rot=curlF(X,Y);                         %计算旋度的值
pcolor(X,Y,rot);                        %绘制旋度
shading interp;                         %颜色做插值，平滑过渡
colorbar;                               %绘制色条，可以在图片窗口对图片的颜色进行设置
hold on;                                %保持模式打开，保持图片
quiver(X,Y,Fx,Fy,'k','linewidth',1);    %绘制箭头图，并设置颜色为黑色，线宽为1