Как векторизовать этот цикл?

@qDenisq · Регистрация: 05.03.2013

Студворк — интернет-сервис помощи студентам

Всем привет.

Есть проблема с векторизацией одного цикла в программке, которая создает массив точек, находящихся между 2 прямыми на заданном отрезке. +100500 в карму тому, кто мне поможет, а то я уже голову сломал и не знаю возможно ли это?

Matlab M
1
2
3
4
5
6
7
pdata=[];
data1=[1 1.1; 2 2.2; 3 3.3; 4 4.4; 5 5.5];
data2=[1 3.3; 2 4.8; 3 6.9; 4 10.2; 5 13];
for i=1:(size(data1,1)-1)
    pdata=cat(1, pdata, [ones(floor(data2(i,2))-ceil(data1(i,2))+1,1)*data1(i,1)...
 [ceil(data1(i,2)):1:floor(data2(i,2))]']); 
end

Зосима · 05.03.2013, 11:45

А что оно должно делать?

чтот не очень ясно

@qDenisq · 05.03.2013, 12:41 **[ТС]**

Вообще в целом задача состоит в нахождении всех точек с целыми координатами (т.е. [1 3], [2 8], а не [1.3 3]) внутри треугольника с заданными вершинами. В примере я просто немного ее упростил и убрал некоторые моменты для удобства понимания. Все бы ничего, но проблема в том, что мне крайне важна скорость работы этой маленькой программы, поэтому надо по максимуму избавиться от циклов.

Зосима · 05.03.2013, 13:47

Сообщение от qDenisq

убрал некоторые моменты для удобства понимания

И стало совсем ничего не понять

объясни словами, как для домохозяйки, как из двух матриц data1 и data2 получается pdata

@qDenisq · 05.03.2013, 16:22 **[ТС]**

))) вобщем, data1 - массив точек лежащих на каком-то отрезке, как видно, взятых с шагом 1 по оси х (первый столбец мы взяли за х) С data2 то же самое - только отрезок другой (это даже может быть впринципе и не отрезок, а ломаная или кривая( суть не меняется). потом я просто для каждого х от 1 до 5 в этом случае создаю массив из точек лежащих между отрезками data1 и data2. потом объединяю и получаю один сплошной массив точек, которые лежат между этими двумя отрезками. как-то так))) на рисунке, например, я нахожу все точки в треугольнике и тут data1-массив от АВ, а data2-массив от АВС.

Зосима · 07.03.2013, 12:08

Гляди, что набросал

Нашел интересную ф-цию inpolygon и вот что вышло:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
clear, clc
 
% координаты вершин
x = [0.4 20.1 15.8];
y = [0.9 15.4 20.8];
% замыкаем многоугольник
X = [x x(1)];
Y = [y y(1)];
% массив целых значений координат
[xi yi] = meshgrid(floor(min(x)):ceil(max(x)), floor(min(y)):ceil(max(y)));
IN = inpolygon(xi,yi, X, Y);
 
plot(X,Y,'-b',xi(IN),yi(IN),'or')
grid minor
axis equal
M = [xi(IN) yi(IN)] % матрица с целочисленными кординатами внутри

Нетрудно догадаться, что можно указывать не три вершины, а больше

@qDenisq · 07.03.2013, 16:10 **[ТС]**

Ты волшебник...спасибо большое)

Новые блоги и статьи Все статьи Все блоги /
Очистка реквизитов документа при копировании Maks 09.04.2026 Алгоритм из решения ниже применим как для типовых, так и для нетиповых документов на самых различных конфигурациях. Задача: при копировании документа очищать определенные реквизиты и табличную. . .	модель ЗдравоСохранения 8. Подготовка к разному выполнению заданий anaschu 08.04.2026 https:/ / github. com/ shumilovas/ med2. git main ветка * содержимое блока дэлэй из старой модели теперь внутри зайца новой модели 8ATzM_2aurI	Блокировка документа от изменений, если он открыт у другого пользователя Maks 08.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа, разработанного в конфигурации КА2. Задача: запретить редактирование документа, если он открыт у другого пользователя. / / . . .	Система безопасности+живучести для сервера-слоя интернета (сети). Двойная привязка. Hrethgir 08.04.2026 Далее были размышления о системе безопасности. Сообщения с наклонным текстом - мои. А как нам будет можно проверить, что ссылка наша, а не подделана хулиганами, которая выбросит на другую ветку и. . .
Модель ЗдрввоСохранения 7: больше работников, больше ресурсов. anaschu 08.04.2026 работников и заданий может быть сколько угодно, но настроено всё так, что используется пока что только 20% kYBz3eJf3jQ	Дальние перспективы сервера - слоя сети с космологическим дизайном интефейса карты и логики. Hrethgir 07.04.2026 Дальнейшее ближайшее планирование вывело к размышлениям над дальними перспективами. И вот тут может быть даже будут нужны оценки специалистов, так как в дальних перспективах всё может очень сильно. . .	Горе от ума kumehtar 07.04.2026 Эта мне ментальная установка, что вот прямо сейчас, мол, мне для полного счастья не хватает (нужное вписать), и когда я этого достигну - тогда и полный кайф. Одна из самых сильных ловушек на пути. . . .	Использование значений реквизитов справочника в документе, с определенными условиями и правами Maks 07.04.2026 1. Контроль срока действия договора Алгоритм из решения ниже реализован на примере нетипового документа "ЗаявкаНаРаботу", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если. . .

Зосима 5246 / 3574 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	05.03.2013, 11:45
	А что оно должно делать? чтот не очень ясно 1

@qDenisq 0 / 0 / 0 Регистрация: 05.03.2013 Сообщений: 10
	05.03.2013, 12:41 [ТС]
	Вообще в целом задача состоит в нахождении всех точек с целыми координатами (т.е. [1 3], [2 8], а не [1.3 3]) внутри треугольника с заданными вершинами. В примере я просто немного ее упростил и убрал некоторые моменты для удобства понимания. Все бы ничего, но проблема в том, что мне крайне важна скорость работы этой маленькой программы, поэтому надо по максимуму избавиться от циклов. 0

@qDenisq 0 / 0 / 0 Регистрация: 05.03.2013 Сообщений: 10
	05.03.2013, 16:22 [ТС]
	))) вобщем, data1 - массив точек лежащих на каком-то отрезке, как видно, взятых с шагом 1 по оси х (первый столбец мы взяли за х) С data2 то же самое - только отрезок другой (это даже может быть впринципе и не отрезок, а ломаная или кривая( суть не меняется). потом я просто для каждого х от 1 до 5 в этом случае создаю массив из точек лежащих между отрезками data1 и data2. потом объединяю и получаю один сплошной массив точек, которые лежат между этими двумя отрезками. как-то так))) на рисунке, например, я нахожу все точки в треугольнике и тут data1-массив от АВ, а data2-массив от АВС. Миниатюры 0

@qDenisq 0 / 0 / 0 Регистрация: 05.03.2013 Сообщений: 10
	07.03.2013, 16:10 [ТС]
	Ты волшебник...спасибо большое) 0