Использование системы координат в Matlab - PullRequest
0 голосов
/ 18 октября 2011

Я работал над проектом, связанным с проблемой обратного источника, известной в поле электромагнитных волн. У меня проблема в том, что; Я должен определить 3 точки в 2D-пространстве. Эти точки должны иметь координаты x, y и значение, которое будет определять их текущее значение. Как это:

A1(2,3)=1
A2(2,-2)=2
and so on.

Также я должен определить круг вокруг этого и разделить его на 200 пунктов. Как первый пункт будет; скажем R=2 ; B1(2,0) ;B50(0,2);B100(-2,0) и т. д.

Теперь мне действительно трудно определить пространство в MATLAB и обвести его. Так что я прошу помочь мне определить 2D-пространство и сделать это так, как я описал. Спасибо за любую помощь, ребята!

Ответы [ 2 ]

0 голосов
/ 27 октября 2011

Вот полный пример, который может помочь:

%# points
num = 3;
P = [2 3; 2 -2; -1 1];          %# 2D points coords
R = [2.5 3 3];                  %# radii of circles around points

%# compute circle points
theta = linspace(0,2*pi,20)';   %'
unitCircle = [cos(theta) sin(theta)];
C = zeros(numel(theta),2,num);
for i=1:num
    C(:,:,i) = bsxfun(@plus, R(i).*unitCircle, P(i,:));
end

%# prepare plot
xlims = [-6 6]; ylims = [-6 6];
line([xlims nan 0 0],[0 0 nan ylims], ...
    'LineWidth',2, 'Color',[.2 .2 .2])
axis square, grid on
set(gca, 'XLim',xlims, 'YLim',ylims, ...
    'XTick',xlims(1):xlims(2), 'YTick',xlims(1):xlims(2))
title('Cartesian Coordinate System')
xlabel('x-coords'), ylabel('y-coords')

hold on

%# plot centers
plot(P(:,1), P(:,2), ...
    'LineStyle','none', 'Marker','o', 'Color','m')
str = num2str((1:num)','A%d');    %'
text(P(:,1), P(:,2), , str, ...
    'HorizontalAlignment','left', 'VerticalAlignment','bottom')

%# plot circles
clr = lines(num);
h = zeros(num,1);
for i=1:num
    h(i) = plot(C(:,1,i), C(:,2,i), ...
        'LineStyle','-', 'Marker','.', 'Color',clr(i,:));
end
str = num2str((1:num)','Circle %d');    %'
legend(h, str, 'Location','SW')

hold off

screenshot

0 голосов
/ 18 октября 2011

Этот вид кода может быть использован.Посмотрите на grid в редакторе переменных.

grid = zeros(50, 50);
R = 10;
angles = (1:200)/2/pi;
x = cos(angles)*R;
y = sin(angles)*R;

center = [25 20];

for n=1:length(angles)   
    grid(center(1)+1+round(x(n)), center(2)+1+round(y(n))) = 1;
end

Вы должны определить сетку, достаточно большую для ваших нужд.

...