Подсчитать среднее значение по массивам из файла

@SpaceQuester · Регистрация: 19.01.2016

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

Друзья!

Помогите освежить в чистом C как считать массивы данных из файла и подсчитать среднее значение. Подзабыл я к сожалению.

Итак в файл результатов сваливается такая информация: "t, theta(1), theta(2), ..., theta(N-1), theta(N)", это столбцы в файле разделенные табуляцией. Число нодов N мы знаем. Сколько строчек по времени t мы изначально не знаем, тут уж как получится.
Но вытащить данные мы должны как массивы: t[0..size], theta(1)[0..size], theta(2)[0..size], ... , theta(N-1)[0..size], theta(N)[0..size].

Затем нужно подсчитать $https://www.cyberforum.ru/cgi-bin/latex.cgi?\rho$ параметр. Он считается как: $https://www.cyberforum.ru/cgi-bin/latex.cgi?\rho_{total} = \frac{1}{N} \sqrt{\sum_{i=0}^N{sin(\theta_i)^2} + \sum_{i=0}^N{cos(\theta_i)^2}}$ причем значения $https://www.cyberforum.ru/cgi-bin/latex.cgi?\rho$ нужно собрать во все моменты времени.
Так у нас получится массив функций $https://www.cyberforum.ru/cgi-bin/latex.cgi?\rho$ в каждый момент времени (это значит от 0 до size, где size - длинна всех массивов, она очевидно одинаковая). Далее мы усредняем по второй половине временного массива, т.е. искомая величина наша будет: $https://www.cyberforum.ru/cgi-bin/latex.cgi?\langle \rho_{total} \rangle = \frac{1}{size/2} \sum_{j = size/2}^{size}{\rho_{total}[j]}$ , где j - это индекс когда мы по времени в массиве идём.

Аналогично будут считаться другое $https://www.cyberforum.ru/cgi-bin/latex.cgi?\rho_{glial}$ :
$https://www.cyberforum.ru/cgi-bin/latex.cgi?\rho_{glial} = \frac{1}{N/2} \sqrt{\sum_{i=0}^{N/2}{sin(\theta_i)^2} + \sum_{i=0}^{N/2}{cos(\theta_i)^2}}$
$https://www.cyberforum.ru/cgi-bin/latex.cgi?\langle \rho_{glial} \rangle = \frac{1}{size/2} \sum_{j = size/2}^{size}{\rho_{glial}[j]}$

И наконец $https://www.cyberforum.ru/cgi-bin/latex.cgi?\rho_{neuron}$ :
$https://www.cyberforum.ru/cgi-bin/latex.cgi?\rho_{neuron} = \frac{1}{N/2} \sqrt{\sum_{i=N/2}^N{sin(\theta_i)^2} + \sum_{i=N/2}^N{cos(\theta_i)^2}}$
$https://www.cyberforum.ru/cgi-bin/latex.cgi?\langle \rho_{neuron} \rangle = \frac{1}{size/2} \sum_{j = size/2}^{size}{\rho_{neuron}[j]}$

На выходе нужно отдать три значения: $https://www.cyberforum.ru/cgi-bin/latex.cgi?\rho$ , $https://www.cyberforum.ru/cgi-bin/latex.cgi?\rho_{glial}$ и $https://www.cyberforum.ru/cgi-bin/latex.cgi?\rho_{neuron}$ :

В MATLAB это считается вот так:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
filename4 = 'results.txt';
path4 = '/Users/sergey/Solver';
Data = importdata( [ path4 '/' filename4 ] );
time = Data(:, 1);
theta = Data(:, 2:end);
 
rx_total = (1/(size(theta, 2)))*sum(cos(theta(:, :)')); % Sum of mean x-part of theta
ry_total = (1/(size(theta, 2)))*sum(sin(theta(:, :)')); % Sum of mean y-part of theta
r_total = sqrt(rx_total.*rx_total + ry_total.*ry_total);
r_total_mean = mean(r_total(:, round(size(r_total, 2)*0.5):end));
 
rx_wired_glial = (1/(size(theta(:, 1:end/2), 2)))*sum(cos(theta(:, 1:end/2)')); % Sum of mean x-part of theta
ry_wired_glial = (1/(size(theta(:, 1:end/2), 2)))*sum(sin(theta(:, 1:end/2)')); % Sum of mean y-part of theta
r_wired_glial = sqrt(rx_wired_glial.*rx_wired_glial + ry_wired_glial.*ry_wired_glial);
r_wired_glial_mean = mean(r_wired_glial(:, round(size(r_wired_glial, 2)*0.5):end))
 
rx_random_neuron = (1/(size(theta(:, end/2+1:end), 2)))*sum(cos(theta(:, end/2+1:end)')); % Sum of mean x-part of theta
ry_random_neuron = (1/(size(theta(:, end/2+1:end), 2)))*sum(sin(theta(:, end/2+1:end)')); % Sum of mean y-part of theta
r_random_neuron = sqrt(rx_random_neuron.*rx_random_neuron + ry_random_neuron.*ry_random_neuron);
r_random_neuron_mean = mean(r_random_neuron(:, round(size(r_random_neuron, 2)*0.5):end))

Файл results.txt я прикладываю для сверки результатов в MATLAB и С. results.txt

Помогите как сделать правильный расчет того же самое в чистом C.

Спасибо! С меня лайки))

Новые блоги и статьи Все статьи Все блоги /
Модель микоризы: классовый агентный подход 3 anaschu 06.01.2026 aa0a7f55b50dd51c5ec569d2d10c54f6/ O1rJuneU_ls https:/ / vkvideo. ru/ video-115721503_456239114	Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ФедосеевПавел 06.01.2026 Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .	Модель микоризы: классовый агентный подход 2 anaschu 06.01.2026 репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .	Расчёт токов в цепи постоянного тока igorrr37 05.01.2026 / * Дана цепь постоянного тока с сопротивлениями и напряжениями. Надо найти токи в ветвях. Программа составляет систему уравнений по 1 и 2 законам Кирхгофа и решает её. Последовательность действий:. . .
Новый CodeBlocs. Версия 25.03 palva 04.01.2026 Оказывается, недавно вышла новая версия CodeBlocks за номером 25. 03. Когда-то давно я возился с только что вышедшей тогда версией 20. 03. С тех пор я давно снёс всё с компьютера и забыл. Теперь. . .	Модель микоризы: классовый агентный подход anaschu 02.01.2026 Раньше это было два гриба и бактерия. Теперь три гриба, растение. И на уровне агентов добавится между грибами или бактериями взаимодействий. До того я пробовал подход через многомерные массивы,. . .	Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Programma_Boinc 28.12.2025 Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / ********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / ********/ gallery/ bKBkQFf Пост отсюда. . .	Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Нашел на реддите интересную статью под названием Anyone know where to get a free Desktop or Laptop? Ниже её машинный перевод. После долгих разбирательств я наконец-то вернула себе. . .