| ← Октябрь 2011 → | ||||||
|
1
|
2
|
|||||
|---|---|---|---|---|---|---|
|
3
|
4
|
5
|
6
|
7
|
8
|
9
|
|
10
|
11
|
12
|
13
|
14
|
15
|
16
|
|
17
|
18
|
20
|
21
|
22
|
23
|
|
|
24
|
25
|
26
|
27
|
28
|
30
|
|
|
31
|
||||||
За последние 60 дней ни разу не выходила
Сайт рассылки:
http://rfpro.ru/
Открыта:
11-11-2009
Статистика
0 за неделю
RFpro.ru: Алгоритмы и теория программирования
Хостинг портала RFpro.ru: РАССЫЛКИ ПОРТАЛА RFPRO.RU
Лучшие эксперты по данной тематике
/ КОМПЬЮТЕРЫ И СОФТ / Программирование / Алгоритмы и теория программирования
Консультация # 184310: Уважаемые эксперты! Пожалуйста, ответьте на вопрос: Есть набор экспериментально полученных точек(200 они тут). По ним строиться петля гистерезиса магнитного материала. У построенной кривой нужно посчитать коэрцитивную силу, то есть найти серединые точки петли петли, и посчитать в этом м... Уважаемые эксперты! Пожалуйста, ответьте на вопрос:
Дата отправки: 24.10.2011, 12:58 Консультирует Лысков Игорь Витальевич (Старший модератор): Здравствуйте, Посетитель - 383763! Код : #include <iostream>
#include <malloc.h>
#include <stdio.h>
#include <string.h>
#include <math.h>
using namespace std;
typedef struct POINT
{
double x;
double y;
}POINT;
//Загрузка всех точек в массив
//возвращает заполненный массив и количество точек
int ReadPoints(char* fName, POINT** points)
{
char line[256];
char *pline;
POINT one;
int max = 1000;
int count;
*points = (POINT*)malloc((sizeof POINT)*max);
FILE * pFile = fopen( fName, "r" );
if ( !pFile )
return 0;
for (count=0 ;0!=fgets( line, 256, pFile );)
{
one.x = atof(line);
pline = strchr(line, ' ');
one.y = atof(pline);
(*points)[count++] = one;
if (count >= max)
{
max += 1000;
*points = (POINT*)realloc(*points, (sizeof POINT)*max);
}
}
fclose( pFile );
return count;
}
//расчет коэрцитивной силы
//параметры: число точек и массив с точками
double Calc(int count, POINT *points)
{
double ymin, ymax, ymid, ycurr;
double x1, x2;
int i;
//найдем максимальный и минимальный Y
ymax = -10000000; //max
ymin = 10000000; //min
for (i=0; i<count; i++)
{
if (points[i].y > ymax)
ymax = points[i].y;
if (points[i].y < ymin)
ymin = points[i].y;
}
//теоритическое среднее
ymid = (ymax - ymin)/2;
//ищем первую точку
//ycurr - текущее, ближайшее к ymid Y
for (ycurr=100000,i=0; i<count; i++)
{
//идем по точкам, пока идет уменьшение
if (fabs(ymid - points[i].y) < fabs(ymid - ycurr))
{
//запоминаем
ycurr = points[i].y;
x1 = points[i].x;
} //учтем возможную "неправильность" данных на краях
else if (fabs(ymid - points[i].y) < 100)
//выходим из поиска, если получаем вблизи точки увеличение расстояния
break;
}
//"разбегаемся" :), находим наиболее удаленную точку
for (; i<count; i++)
{
if (fabs(ymid - points[i].y) > fabs(ymid - ycurr))
ycurr = points[i].y;
else
break;
}
//и "с разбега", аналогично предыдущему, находим вторую точку
for (; i<count; i++)
{
if (fabs(ymid - points[i].y) < fabs(ymid - ycurr))
{
ycurr = points[i].y;
x2 = points[i].x;
}
else if (fabs(ymid - points[i].y) < 100)
break;
}
//осталось только посчитать модуль разности
return (fabs(x2 - x1));
}
int main(int argc, char *argv[ ])
{
system("chcp 1251 >> nul");
int count;
POINT * points = 0;
if (argc > 1)
{
if (0 != (count = ReadPoints(argv[1], &points)))
{
cout << "Коэрцитивная сила = " << Calc(count, points) << endl;
free(points);
}
else
cout << "Ошибка чтения файла" << endl;
}
else
cout << "Введите: CalcF "<имя файла>" << endl;
system("pause");
return 0;
}
Оценить выпуск | Задать вопрос экспертам
главная страница
|
стать участником
|
получить консультацию
© 2001-2011, Портал RFPRO.RU, Россия
Авторское право: ООО "Мастер-Эксперт Про" Калашников О.А. | Гладенюк А.Г. Хостинг: Компания "Московский хостер" Версия системы: 2011.6.33 от 10.09.2011 |
| В избранное | ||
