通过网盘分享的文件:电流场模拟静电场实验实验报告模板
链接: https://pan.baidu.com/s/1b4NQgKrpdKKGpjqQaT-B5w?pwd=dcfk
提取码: dcfk
r=[5.25.215.195.205.225.185.205.235.205.19;...4.214.204.224.194.204.184.204.214.194.20;...3.403.413.393.403.423.383.403.413.403.39;...2.812.802.822.792.802.812.802.792.802.80;...2.42.412.422.392.402.382.402.412.402.4;...1.911.901.891.901.911.881.901.911.901.89;...1.601.611.591.601.611.591.601.611.601.59];theta=[03060105135180210270300340;...03060105135180210270300340;...03060105135180210270300340;...03060105135180210270300340;...03060105135180210270300340;...03060105135180210270300340;...03060105135180210270300340];tongzhoudianlan(r,theta)% 假设测量点存储在一个 N×2 的矩阵中,每一行是 (x, y)% points = [x1, y1; x2, y2; ...]; % 替换为你的测量点 points = [0, 3.4; 1, 3.8; 2,4.2;3,4.8;4,5.3;5,5.7;9,5.8;10,5.3;11,4.8;12,4.2;...0,6.6;1,6.6;2,6.8;3,6.9;4,7;5,7.1;9,7;10,6.9;11,6.8;12,6.6;...0,9;1,8.9;2,8.8;3,8.7;4,8.6;5,8.6;9,7;10,8.4;11,8.5;12,8.5;...0,11.6;1,11.4;2,11.0;3,10.7;4,10.3;5,10;9,9.8;10,10;11,10.2;12,10.6;...2,15;3,13.4;4,12.3;5,11.4;6,10.8;7.4,10.6;9,10.9;10,11.4;11,12.3;12,13.4];% 6,16;6,12;7,11.4;8,11.5;9,12.2;9,16;9.6,13.8;5.6,12.5;7.2,11.4;8.2,11.6;];%% 提取 x 和 y 坐标x = points(:, 1);y = points(:, 2);% 假设每个点对应的电势值为 V(需要根据实际情况提供)v=[3:8];V=[ones(10,1)*3;ones(10,1)*4;ones(10,1)*5;ones(10,1)*6;ones(10,1)*7];%;ones(10,1)*8% V = [V1; V2; ...]; % 替换为你的电势值% 创建网格点[xq, yq] = meshgrid(linspace(0, 16, 10), linspace(0, 16, 10));% 插值生成平滑的电势场Vq = griddata(x, y, V, xq, yq, 'cubic');% 计算电势场的梯度(电场)[Ex, Ey] = gradient(-Vq); % 电场是电势的负梯度% 绘制等势线figure;% contour(xq, yq, Vq, 20); % 20 表示绘制 20 条等势线hold on;x1=reshape(x,10,5);y1=reshape(y,10,5);scatter(x, y, 'filled'); % 绘制测量点% % 绘制电场线quiver(xq, yq, Ex, Ey); % 绘制电场线plot(x1(:,1),y1(:,1))plot(x1(:,2),y1(:,2))plot(x1(:,3),y1(:,3))plot(x1(:,4),y1(:,4))plot(x1(:,5),y1(:,5))% contour(xq, yq, Vq, 7); % 叠加等势线xlabel('x');ylabel('y');title('电场线与等势线');colorbar;axis equal tightbox on grid on
夜雨聆风