Pregunta

He estado trabajando en un proyecto que involucra un problema de origen inverso conocido dentro del campo de onda electromagnética. El problema que tengo es ese; Tengo que definir 3 puntos en un espacio 2D. Estos puntos deben tener Ax, y coordinar, por supuesto y un valor que definirá su corriente. Como esto:

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

También tengo que definir un círculo alrededor de esto y dividirlo en 200 puntos. Como el primer punto sería; decir R=2 ; B1(2,0) ;B50(0,2);B100(-2,0) y así.

Ahora realmente estoy teniendo dificultades para definir un espacio en Matlab y rodearlo. Entonces, lo que estoy pidiendo es que me ayude a definir un espacio 2D y lo haga como lo describí. ¡Gracias por cualquier ayuda chicos!

¿Fue útil?

Solución

Este tipo de código puede ser usado. Mirar grid En el editor de variables.

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

Tienes que definir una cuadrícula lo suficientemente grande para tu necesidad.

Otros consejos

Aquí hay un ejemplo completo que podría ser de ayuda:

%# 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

Licenciado bajo: CC-BY-SA con atribución
No afiliado a StackOverflow
scroll top