设有一长直矩形a2a接地金属槽上盖板对地绝缘ф1.docxVIP

设有一长直矩形(a*2a)接地金属槽，上盖板对地绝缘 4)1=10，侧壁，底4)2=0,求槽内的电位 分布。 一.有限差分法 取步长h为1, a为10,贝U长边与短边上的节点数分别为 21和11。具体采用简单迭代 法，等位线取30条，代码和计算结果如下： function [n, v1] = FDM(a, h)%h 为步长 px = 2 * a / h + 1; % 长边上的点数 py = a / h + 1; % 短边上的点数 v1 = ones(py, px); % 设置一 10*20 的矩阵 v1(1, :) = ones(1, px) * 10; % 初始化上边界 v1(py, :) = zeros(1, px); % 初始化下边界 for i = 1 : py % 初始化左右边界 v1(i, 1) = 0; v1(i, px) = 0; end v2 = v1; maxt = 1; t = 0; n = 0; while(maxt 1e-6) maxt = 0; n = n + 1; for i = 2 : py - 1 for j = 2 : px - 1 v2(i, j) = 0.25 * (v1(i, j + 1) + v1(i + 1, j) + v1(i - 1, j) + v1(i, j - 1)); % 计算边界内的值 t = abs(v2(i, j) - v1(i, j)); if(t maxt) % 误差判断 maxt = t; end end end v1 = v2; end axis([0, 20, 0, 10, 0, 10]); contour(v1, 30); end 计算的结果如下: 0 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 0 0 4.9451 6.867 7.7425 8.198 8.4589 8.6174 8.716 8.776 8.8083 8.8186 8.8083 8.776 8.716 8.6174 8.4589 8.198 7.7425 6.867 4.9451 0 0 2.9134 4.7805 5.9051 6.5907 7.4707 7.0201 7.2947 7.2947 7.0201 7.4707 6.5907 7.5795 7.6388 7.6577 7.6388 7.5795 5.9051 4.7805 2.9134 0 0 1.9281 3.4365 4.5066 5.2395 5.7362 6.0708 6.2924 6.4325 6.5098 6.5345 6.5098 6.4325 1.9281 0 6.2924 6.0708 5.7362 5.2395 4.5066 3.4365 0 1.3623 2.531 3.4451 4.1245 4.6144 4.9597 5.1958 5.3482 5.4334 5.4608 5.4334 5.3482 1.3623 0 5.1958 4.9597 4.6144 4.1245 3.4451 2.531 0 0.99018 1.88 2.6186 3.1988 3.6373 3.9578 4.1829 4.331 4.4148 4.4419 4.4148 4.331 0.99018 0 4.1829 3.9578 3.6373 3.1988 2.6186 1.88 0 0.71836 1.3804 1.9502 2.415 2.7783 3.0513 3.2472 3.378 3.4528 3.4772 3.4528 3.378 0.71836 0 3.2472 3.0513 2.7783 2.415 1.9502 1.3804 0 0.50289 0.973 1.3868 1.7328 2.0095 2.2218 2.3765 2.4812 2.5414 2.5611 2.5414 2.4812 0.50289 0 2.3765 2.2218 2.0095 1.7328 1.3868 0.973 0 0.32019 0.62194 0.89112 1 .1197 1.3054 1.4497 1.5561 1.6286 1.6706 1.6844 1.6706 1.6286 0.32019 0 1.5561 1.4497 1.3054 1.1197 0.89112 0.62194 0 0.15592 0.30348 0.43605 0.54959 0.6426 0.71545 0.76952 0.80655 0.82807 0.83513 0.82807 0.80655 0.76952 0.71545 0.6426 0.155