Построение графика по трем точкам

@puhlak · Регистрация: 31.01.2013

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

Коллеги, добрый вечер.

Просьба подсказать.
Есть файл с данными вида:

1 1 5
1 2 7
1 3 5
1 4 1
2 1 4
2 2 5
2 3 1
2 4 7

Или 3 файла раздельных

Кликните здесь для просмотра всего текста

Файл 1:
1
1
..
2

Файл 2:
1
2
3
4
..
4

Файл 3:
5
7
5
..
7

Можно любой сделать.
Нужно построить 3х-мерный график (X,Y,Z)

plot3(x,y,z) рисует проволочный график, неинформативно.
Хотелось бы нарисовать красивый, как поверхность. Где высокие значения - это красный цвет, низкие - зеленый. Или наоборот.
Как это сделать?

Привел x,y,z к матрицам и теперь при вызове surf(x,y,z) вижу:

Matrix dimensions must agree.

При этом:

Code
1
2
3
4
5
6
7
8
9
10
11
>> size(x)
ans =
        3721           1
 
>> size(y)
ans =
        3721           1
 
>> size(z)
ans =
        3721           1

@puhlak · 31.01.2013, 20:18 **[ТС]**

Вот пример данных.
RTS1 - это в 1 файле 3 столбца.
RTS3 - это 3 файла по одному столбцу.

@puhlak · 31.01.2013, 20:36 **[ТС]**

Пробовал:
Как в этой теме Построение поверхности по трем точкам
Получается неверно.

Пробовал:

Matlab M
1
2
3
4
5
6
fid = fopen('MyFile.txt','r'); 
a = fscanf(fid, '%f %f %а', [3 inf]);
a = a'; 
fclose(fid); 
 
surf(a);

Опять же неправильно рисует.

@puhlak · 31.01.2013, 20:45 **[ТС]**

Вот если на 2д. Должно быть аналогично в 3д:

Зосима · 01.02.2013, 10:49

Надо было глянуть еще здесь: Построить трехмерный график по точкам
Там в конце как раз я и предложил алгоритм преобразования трехстолбцовой матрицы (первый файл) в нужный формат. В данном случае получим:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
clear;
clc;
% считываем файл 
M = load('export2_matlab.csv');
% раскидываем столбцы 
X = M(:,1);
Y = M(:,2);
Z = M(:,3);
% формируем векторы Y и X
fst_y = find(Y == Y(1));
len_y = fst_y(2)-1; % длинна вектора Y
Y = Y(1:len_y);
 
len_x = length(fst_y); % длинна вектора X
X = X(fst_y+1);
 
% формируем матрицу Z
Z = reshape(Z,len_y,len_x); 
 
surf(X,Y,Z,'FaceAlpha',0.5)
shading interp
xlabel('X')
ylabel('Y')

Результат *трах-тибидох-тибидох*

@puhlak · 01.02.2013, 11:12 **[ТС]**

Большое спасибо!

Добавлено через 14 минут
Коллега, пара вопросов по коду:

fst_y = find(Y == Y(1)); - "ind = find(X) locates all nonzero elements of array X, and returns the linear indices of those elements in vector ind." - подскажите, а зачем искать эти линейные индексы, не совсем понимаю. Что-то связано с необычным хранением данных в матлабе?

len_y = fst_y(2)-1 - Из числа линейных индексов вычитаем 1? Зачем? Не понял.

Y = Y(1:len_y); - это сформировали массив y - просто с 1 до конца массива?

X = X(fst_y+1); - не понял, почему +1

Заранее спасибо!

Зосима · 01.02.2013, 12:01

Дружочек, для начала вспомни структуру трехстолбцовой матрицы:
$https://www.cyberforum.ru/cgi-bin/latex.cgi?\begin{pmatrix}X & Y & Z\\1 & 1 & 522\\ 1 & 2 & 541\\ 1 & 3 & 154\\ 1 & 4 & 542\\ 2 & 1 & 151\\ 2 & 2 & 121\\ 2 & 3 & 115\\ 2 & 4 & 113\\ 3 & 1 & 721\\ 3 & 2 & 152\\ 3 & 3 & 421\\ 3 & 4 & 564\end{pmatrix}$
Теперь смотри второй столбец (Y).
Строка fst_y = find(Y == Y(1)); находит номера позиций, где значения Y равны первому значению, т.е. начинают повторяться.
Первым элементом массива fst_y будет 1, второй fst_y(2) - номер элемента где Y опять будет равен Y(1), а вот длинна неповторяющихся позиций Y равна как раз таки:
len_y = fst_y(2) - 1;
Y = Y(1:len_y); - берем из столбика Y первые неповторяющиеся значения и сохраняем все это опять таки в Y

X = X(fst_y+1) - поидее можно было без прибавления 1

@n6 · 20.02.2013, 13:35

А не подскажете, почему у меня после выполнения кода, представленного выше, матлаб выдает ошибку: "Attempted to access fst_y(2); index out of bounds because numel(fst_y)=1." ?

@puhlak · 20.02.2013, 13:36 **[ТС]**

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

А не подскажете, почему у меня после выполнения кода, представленного выше, матлаб выдает ошибку: "Attempted to access fst_y(2); index out of bounds because numel(fst_y)=1." ?

Скорее всего неверные данные.

@n6 · 20.02.2013, 17:23

А как представить данные в нужной форме? У меня они в экселе, сведены в таблицу из трех столбцов.

@puhlak · 20.02.2013, 17:47 **[ТС]**

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

А как представить данные в нужной форме? У меня они в экселе, сведены в таблицу из трех столбцов.

Я использовал файлы .csv
Смотрите в аттаче в сообщениях выше.

@n6 · 26.02.2013, 11:36

Да я тоже такие файлы использую, но выдает вот ту ошибку. Там какие-то конкретные требования есть к файлу?

@puhlak · 26.02.2013, 11:41 **[ТС]**

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

Да я тоже такие файлы использую, но выдает вот ту ошибку. Там какие-то конкретные требования есть к файлу?

разделитель десятичный - точка.
Разбить по колонкам.
Вроде все.

Зосима · 26.02.2013, 12:14

n6, ты бы свои выложил, разобрались бы быстрее

или ты шпиён иносраный?!

@n6 · 26.02.2013, 14:02

так давно выложил..
Зосима, не притворяйся как будто не видел мою другую тему

построение поверхности

прямая ссылка - файл

Зосима · 26.02.2013, 14:44

n6, забыл...

* ошибочка вышла *

Тьфу ты! Файлы с одинаковыми именами

нарисовал старый
А в твоем случае ситуация хитрее: X - имеет всего одно значение!

о какой поверхности может идти рэчь?

Вот, держи кривульку в пространстве

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
clear all
clc
 
A = load('export2_matlab2.csv');
X = A(:,1);
Y = A(:,2);
Z = A(:,3);
plot3(X,Y,Z,'r')
grid on
xlabel('X')
ylabel('Y')
zlabel('Z')
% облет камеры
for i = -50:50
    view(i,30)
    pause(0.001)
end

@n6 · 26.02.2013, 15:00

Спасибо) Но это не полный файл. В итоге Х у меня будет равен 1,2,4 и т.д., и для каждого значения будет похожее множество значений в других столбцах. Значит, проблема в том, что здесь одинаковые значения Х? ОК, попробую дополнить файл и построить поверхность.

Зосима · 26.02.2013, 15:03

Кроме того, эту кривульку можно скопировать вдоль оси Х и тогда можно получить, нечто похожее на поверхность

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
clear all
clc
 
% считываем файл 
M = load('export2_matlab2.csv');
% раскидываем столбцы 
a = unique(M(:,1));
b = unique(M(:,2));
z = unique(M(:,3));
 
[X Y] = meshgrid(a-2:0.5:a+2, b);
Z = z*ones(1,size(X,2));
 
surf(X,Y,Z,'FaceAlpha',0.5)
xlabel('X')
ylabel('Y')

@n6 · 26.02.2013, 15:25

Спасибо! Пригодится. Но все же хотелось бы разобраться с поверхностью.
Добавляю другие значения Х, матлаб выдает ошибку:

Error using reshape
To RESHAPE the number of elements must not change.

А это с чем может быть связано?

Зосима · 26.02.2013, 15:49

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

А это с чем может быть связано?

Связано с тем, что ты пытаешься вектор впихнуть матрицу несоответствующего размера

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

Добавляю другие значения Х

Куда ты их добавляешь?

Новые блоги и статьи Все статьи Все блоги /
AkelPad-скрипты, структуры, и немного лирики.. testuser2 05.04.2026 Такая программа, как AkelPad существует уже давно, и также давно существуют скрипты под нее. Тем не менее, прога живет, периодически что-то не спеша дополняется, улучшается. Что меня в первую очередь. . .	Отображение реквизитов в документе по условию и контроль их заполнения Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеСпецтехники", разработанного в конфигурации КА2. Данный документ берёт данные из другого нетипового документа. . .	Фото всей Земли с борта корабля Orion миссии Artemis II kumehtar 04.04.2026 Это первое подобное фото сделанное человеком за 50 лет. Снимок называют новым вариантом легендарной фотографии «The Blue Marble» 1972 года, сделанной с борта корабля «Аполлон-17». Новое фото. . .	Вывод диалогового окна перед закрытием, если документ не проведён Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать программный контроль на предмет проведения документа. . .
Программный контроль заполнения реквизитов табличной части документа Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: 1. Реализовать контроль заполнения реквизита. . .	wmic не является внутренней или внешней командой Maks 02.04.2026 Решение: DISM / Online / Add-Capability / CapabilityName:WMIC~~~~ Отсюда: https:/ / winitpro. ru/ index. php/ 2025/ 02/ 14/ komanda-wmic-ne-naydena/	Программная установка даты и запрет ее изменения Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: при создании документов установить период списания автоматически. . .	Вывод данных в справочнике через динамический список Maks 01.04.2026 Реализация из решения ниже выполнена на примере нетипового справочника "Спецтехника" разработанного в конфигурации КА2. Задача: вывести данные из ТЧ нетипового документа. . .

@puhlak 1 / 1 / 0 Регистрация: 31.01.2013 Сообщений: 20
	31.01.2013, 20:45 [ТС]
	Вот если на 2д. Должно быть аналогично в 3д: Изображения 0

@puhlak 1 / 1 / 0 Регистрация: 31.01.2013 Сообщений: 20
	01.02.2013, 11:12 [ТС]
	Большое спасибо! Добавлено через 14 минут Коллега, пара вопросов по коду: fst_y = find(Y == Y(1)); - "ind = find(X) locates all nonzero elements of array X, and returns the linear indices of those elements in vector ind." - подскажите, а зачем искать эти линейные индексы, не совсем понимаю. Что-то связано с необычным хранением данных в матлабе? len_y = fst_y(2)-1 - Из числа линейных индексов вычитаем 1? Зачем? Не понял. Y = Y(1:len_y); - это сформировали массив y - просто с 1 до конца массива? X = X(fst_y+1); - не понял, почему +1 Заранее спасибо! 0

Зосима 5246 / 3574 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	01.02.2013, 12:01
	Дружочек, для начала вспомни структуру трехстолбцовой матрицы: $https://www.cyberforum.ru/cgi-bin/latex.cgi?\begin{pmatrix}X & Y & Z\\1 & 1 & 522\\ 1 & 2 & 541\\ 1 & 3 & 154\\ 1 & 4 & 542\\ 2 & 1 & 151\\ 2 & 2 & 121\\ 2 & 3 & 115\\ 2 & 4 & 113\\ 3 & 1 & 721\\ 3 & 2 & 152\\ 3 & 3 & 421\\ 3 & 4 & 564\end{pmatrix}$ Теперь смотри второй столбец (Y). Строка fst_y = find(Y == Y(1)); находит номера позиций, где значения Y равны первому значению, т.е. начинают повторяться. Первым элементом массива fst_y будет 1, второй fst_y(2) - номер элемента где Y опять будет равен Y(1), а вот длинна неповторяющихся позиций Y равна как раз таки: len_y = fst_y(2) - 1; Y = Y(1:len_y); - берем из столбика Y первые неповторяющиеся значения и сохраняем все это опять таки в Y X = X(fst_y+1) - поидее можно было без прибавления 1 1

@n6 0 / 0 / 0 Регистрация: 06.02.2013 Сообщений: 15
	20.02.2013, 13:35
	А не подскажете, почему у меня после выполнения кода, представленного выше, матлаб выдает ошибку: "Attempted to access fst_y(2); index out of bounds because numel(fst_y)=1." ? 0

@n6 0 / 0 / 0 Регистрация: 06.02.2013 Сообщений: 15
	20.02.2013, 17:23
	А как представить данные в нужной форме? У меня они в экселе, сведены в таблицу из трех столбцов. 0

@n6 0 / 0 / 0 Регистрация: 06.02.2013 Сообщений: 15
	26.02.2013, 11:36
	Да я тоже такие файлы использую, но выдает вот ту ошибку. Там какие-то конкретные требования есть к файлу? 0

Зосима 5246 / 3574 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	26.02.2013, 12:14
	n6, ты бы свои выложил, разобрались бы быстрее или ты шпиён иносраный?! 0

@n6 0 / 0 / 0 Регистрация: 06.02.2013 Сообщений: 15
	26.02.2013, 14:02
	так давно выложил.. Зосима, не притворяйся как будто не видел мою другую тему построение поверхности прямая ссылка - файл 0

@n6 0 / 0 / 0 Регистрация: 06.02.2013 Сообщений: 15
	26.02.2013, 15:00
	Спасибо) Но это не полный файл. В итоге Х у меня будет равен 1,2,4 и т.д., и для каждого значения будет похожее множество значений в других столбцах. Значит, проблема в том, что здесь одинаковые значения Х? ОК, попробую дополнить файл и построить поверхность. 0