Exponenta Banner

interpolateElectricPotential

Интерполируйте электрический потенциал в электростатическом результате в произвольных пространственных местоположениях

свернуть все на странице

    Синтаксис

    Vintrp = interpolateElectricPotential(electrostaticresults,xq,yq)
    Vintrp = interpolateElectricPotential(electrostaticresults,querypoints)

    Описание

    пример

    Vintrp = interpolateElectricPotential(electrostaticresults,xq,yq) возвращает интерполированные электрические потенциальные значения в точках, заданных в xq и yq.

    пример

    Vintrp = interpolateElectricPotential(electrostaticresults,querypoints) возвращает интерполированные электрические потенциальные значения в точках, заданных в querypoints.

    Примеры

    свернуть все

    Скрипт Open Live Script

    Создайте электромагнитную модель для электростатического анализа.

    emagmodel = createpde('electromagnetic','electrostatic');

    Создайте квадратную геометрию и включайте ее в модель. Постройте геометрию с метками ребра.

    R1 = [3,4,-1,1,1,-1,1,1,-1,-1]';
g = decsg(R1, 'R1', ('R1')');
geometryFromEdges(emagmodel,g);
pdegplot(emagmodel,'EdgeLabels','on')
xlim([-1.5 1.5])
axis equal

    Figure contains an axes. The axes contains 5 objects of type line, text.

    Задайте вакуумную проницаемость в системе СИ модулей.

    emagmodel.VacuumPermittivity = 8.8541878128E-12;

    Задайте относительную проницаемость материала.

    electromagneticProperties(emagmodel,'RelativePermittivity',1);

    Примените граничные условия напряжения на ребра квадрата.

    electromagneticBC(emagmodel,'Voltage',0,'Edge',[1 3]);
electromagneticBC(emagmodel,'Voltage',1000,'Edge',[2 4]);

    Задайте плотность заряда для целой геометрии.

    electromagneticSource(emagmodel,'ChargeDensity',5E-9);

    Сгенерируйте mesh.

    generateMesh(emagmodel);

    Решите модель и постройте электрический потенциал.

    R = solve(emagmodel);
pdeplot(emagmodel,'XYData',R.ElectricPotential, ...
                  'Contour','on')
axis equal

    Figure contains an axes. The axes contains 12 objects of type patch, line.

    Интерполируйте получившийся электрический потенциал к сетке, покрывающей центральный фрагмент геометрии для x и y от -0.5 к 0.5.

    v = linspace(-0.5,0.5,51);
[X,Y] = meshgrid(v);

Vintrp = interpolateElectricPotential(R,X,Y)
    Vintrp = 2601×1

  602.2959
  616.0208
  629.0498
  641.4049
  653.0828
  664.0757
  674.4209
  684.1432
  693.2704
  701.8026
      ⋮

    Измените Vintrp и постройте получившийся электрический потенциал.

    Vintrp = reshape(Vintrp,size(X));

figure
contourf(X,Y,Vintrp)
colormap(cool)
colorbar

    Figure contains an axes. The axes contains an object of type contour.

    В качестве альтернативы можно задать сетку при помощи матрицы точек запроса.

    querypoints = [X(:),Y(:)]';
Vintrp = interpolateElectricPotential(R,querypoints);

    Входные параметры

    свернуть все

    Решение электростатической проблемы в виде ElectrostaticResults объект. Создайте electrostaticresults использование solve функция.

    Пример: electrostaticresults = solve(emagmodel)

    x- точки запроса в виде действительного массива. interpolateElectricPotential оценивает электрический потенциал в точках координаты [xq(i) yq(i)] для каждого i. Из-за этого, xq и yq должен иметь то же количество записей.

    interpolateElectricPotential преобразует точки запроса в вектор-столбцы xq(:) и yq(:). Это возвращает электрические потенциальные значения как вектор-столбец, одного размера. Чтобы гарантировать, что размерности возвращенного решения сопоставимы с размерностями точек исходного запроса, использовать reshape. Например, используйте Vintrp = reshape(Vintrp,size(xq)).

    Пример: xq = [0.5 0.5 0.75 0.75]

    Типы данных: double

    y- точки запроса в виде действительного массива. interpolateElectricPotential оценивает электрический потенциал в точках координаты [xq(i) yq(i)] для каждого i. Из-за этого, xq и yq должен иметь то же количество записей.

    interpolateElectricPotential преобразует точки запроса в вектор-столбцы xq(:) и yq(:). Это возвращает электрические потенциальные значения как вектор-столбец, одного размера. Чтобы гарантировать, что размерности возвращенного решения сопоставимы с размерностями точек исходного запроса, использовать reshape. Например, используйте Vintrp = reshape(Vintrp,size(yq)).

    Пример: yq = [1 2 0 0.5]

    Типы данных: double

    Точки запроса в виде действительной матрицы с двумя строками для 2D геометрии. interpolateElectricPotential оценивает электрический потенциал в точках координаты querypoints(:,i) для каждого i, так каждый столбец querypoints содержит точно одну 2D точку запроса.

    Пример: querypoints = [0.5 0.5 0.75 0.75; 1 2 0 0.5]

    Типы данных: double

    Выходные аргументы

    свернуть все

    Электрический потенциал в точках запроса, возвращенных как FEStruct объект со свойствами, представляющими пространственные компоненты электрического потенциала в точках запроса. Для точек запроса, которые находятся вне геометрии, Vintrp(i) = NaN. Свойства FEStruct объект только для чтения.

    Смотрите также

    | | | |

    Введенный в R2021a

    Документация Partial Differential Equation Toolbox

    Поддержка