[ На главную ] -- [ Список участников ] -- [ Правила форума ] -- [ Зарегистрироваться ]

проект D&E

  Вход на форум   логин       пароль   Забыли пароль? Регистрация
On-line:  

Раздел: 
Gauss2k - gauss gun у тебя дома / Gauss tecnology / проект D&E

Страницы: << Prev 1 2 3 Next>> новая тема

Автор Сообщение

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 20-02-2009 09:30
А ты не мог бы подставить какую-нибудь таблицу ВН и проверить насколько адекватно подбирается формула?

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 20-02-2009 09:34
щас попробую.
Как с компом у тебя?

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 20-02-2009 09:36
Видюха ебнутая. Потому никак.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 21-02-2009 21:08
жепь какая то выходит опять - график вообще вниз уходит круто. юзал те пары что ты прислал с vb-проектом

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 21-02-2009 22:10
Я вспомнил! Последние мои работы в этом направлении успехом так и не увенчались. Что любопытно: для небольшого колличества пар, менее 5 чтоли, все работает.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 21-02-2009 22:22
да и для 5 не пашет как надо

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 21-02-2009 22:24
Значит у тебя устаревшая версия :D

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 21-02-2009 22:33
ну вот то что ты на бэсике мне прислал я перекодил на c#.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 22-02-2009 09:13
подрубил к проге готовый код для интерполяции сплайнами который лежал по той ссылке которые ты мне давал для ознакомления, и все заработало как надо - на графике очень пологая кривая и проходит там где надо.

update
ссылку не знаю ибо вся переписка удалилась, а код тот я додумался сохранить на компе.

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 22-02-2009 10:06
[Доволен]
Я вывел формулу для тока среднего, максимального, для времени, за которое ток достигает наибольшего значения. В силу моей увлеченности средними значениями, здесь также не без них - индуктивность берётся средняя за весь процесс. Как вычислять её - ещё будет видно. Вообщем-то самое сложное я уже понял - как вычислять индуктивность контура с сердечником. Позже перешлю всё это дело как нибудь.
А ссылку я тебе могу переслать, только геморно искать нужное письмо через сотовый.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 22-02-2009 10:17
А ссылку я тебе могу переслать

да не не стоит, все что нужно мне было на том сайте я сразу сохранил

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 22-02-2009 12:42
Кстати, свой метод интерполяции я начал мутить не потому, что посчитал себя шибко умным, а потому что с интерполяцией кусками в нашем случае работать очень сложно.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 22-02-2009 14:27
Кстати, свой метод интерполяции я начал мутить не потому, что посчитал себя шибко умным, а потому что с интерполяцией кусками в нашем случае работать очень сложно.

а че с ним работать? Функцию выдает, считает быстро, и нормально так ведь?

вротмненоги, что ж мне так не везет. нечайно удалил проект - благо бэкап сделал вчера, щас сидел дописывал.

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 22-02-2009 15:20
а че с ним работать? Функцию выдает, считает быстро, и нормально так ведь?

Не разбирался с тамошними исходниками, но если выдает функцию, тогда норм. Кстати, что по твоему значит 'выдает функцию'? Строит график? Сам график без промежуточных циферок нинахер не нужен. А интерполяция необходима для вычислений, для получения зависимости B=B(H)

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 22-02-2009 15:38
там короче как у тебя все - сначала я 1 раз получаю вектор неизвестных а потом в другой процедуре юзаю его для интерполяции. Только там как то сложно все я нихера не понимаю че там написано, и хлама много всякого - сразу видно математик писал.

Кстати, что по твоему значит 'выдает функцию'?

в результате я имею функцию Function B(H as Double) as Double (если по бэсику), а уже от нее и график и анализ идти будет.

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 22-02-2009 17:18
вектор неизвестных

А можно ссылочку где так написано? Образования ради. Да и спасибо за то что у тебя хватает терпения работать со мной :)

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 22-02-2009 19:52
препод так говорил о матрице X когда решал слау методом гаусса.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 23-02-2009 10:40
фу она плохая эта интерполяция сплайнами - я щас ее на других точках потестил, она их просто прямыми отрезками соединяет и все. У тебя красиво было гладко но не то выдает почемуто

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 23-02-2009 13:22
Я поищу ошибку, может вкралась куда? КСТАТЕ!!!11 ПЕРВАЯ ПАРА ОБЯЗАТЕЛЬНО (1 ;0.01 )!!!

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 23-02-2009 17:26
для чистого железа вообще круть

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 23-02-2009 17:50
Ты как откладывал? По оси абсцисс должно быть от нуля до ста тысяч например! Да и вообще я вижу что-то совсем лажа - вниз уходит, да еще и 0.5 Тл при 2.4 А/м! Эх, без компьютера как без рук.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 23-02-2009 18:31
По оси абсцисс должно быть от нуля до ста тысяч например!

я ничего не откладывал, все отрисовщик сам делает, я ему просто функцию передаю и все, а он считает только то что видно на экране, и пересчитывает все при сдвиге графика. Кстати первую пару ставил как ты написал и все равно все графики примерно такие. Скажи как эта интерполяция твоя называется? Мож сам чего нагуглю почитаю.

Эх, без компьютера как без рук

ога не хорошо
/ *вспоминает себя на каникулах*/

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 24-02-2009 00:41
при 2.4 А/м!

почему А/м? В/м ведь вроде

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 24-02-2009 06:26
Не понимаю о чем вы
http://ru.wikipedia.org/wiki/Напряжённость_магнитного_поля
Аппроксимация - так вот что это такое! По сути это и есть подбор эмпирической формулы. У меня есть странное ощущение, что ВН-кривая очень похожа на график какой-то общеизвестной функции...
Скажи как эта интерполяция твоя называется? Мож сам чего нагуглю почитаю.

Я ж тебе говорил, что это собсно сварганенный метод. В интернете почитать - разве только я выложу :)

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 24-02-2009 07:12
А понятно, я думал то напряженность электрического поля.
Я ж тебе говорил, что это собсно сварганенный метод
Апроксимация выдает фунцию которая не через все исходные точки проходит, или как там она, щас надо почитать.

update
А разве она подойдет? Интерполяцию всеж надо, поэтому наверно она у тебя и не работает что не через все точки проходит. Хотя сам думай в математике я тебе не советник.

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 24-02-2009 07:29
Апроксимация выдает фунцию которая не через все исходные точки проходит, или как там она, щас надо почитать.

Да и похер! Ведь близко к реальности!
Вот накопал из глубин шкафа журнал "Reader's Digest" за октябрь 1997 года. Там статья - 'Насколько точен приблизительный ответ' Может где в интернете есть электронная версия?
Добавлено позже: а нет

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 24-02-2009 07:33
а мож тогда интерполяцию сплайнами оставить? это на крутых графиках заметна ее ломанность а на BH кривой(если она вообще кривая а не прямая) ничего не заметно.

магистр
Группа: Модераторы
Сообщений: 861
 Добавлено: 24-02-2009 11:31
http://ru.wikipedia.org/wiki/Гамма_распределение
Там внизу
"Примеры абсолютно непрерывных распределений"
очень похоже на график распределения какого-нить.

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 24-02-2009 17:00
Да, за ссылку спасибо. Уфф, начинаю новое погружение в неизвестность...

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 27-02-2009 06:25
Кстати, проблемы с интерполяцией могут быть из-за того, что в ячейках А сначала хранятся очень большие значения, а потом очень маленькие. Отсюда возможны и внушительные погрешности.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 05-03-2009 13:13
Во, для значений
x y
0 0
1 1
2 4
порабола чоткая получается, но если больше точек добавить уже не то.


а полиномиальная интерполяция не подойдет? я для нее готовый код нашел, только щас не разобрался еще как его подрубить.

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 05-03-2009 13:29
а полиномиальная интерполяция не подойдет?

Не знаю её особенностей. Попробуй.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 05-03-2009 19:32
уффф....
А вот те крестики в формулах это векторное произведение или постое?
И как полное сопротивление RLC цепи считать, все эти внутренние ебкости, esr?
А формула L(x) будет, или она там не нужна?

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 05-03-2009 20:09
Крестики - простое произведение. Там нигде векторов нет. Полное сопротивление цепи = сопр. обмотки + 0.42ом (также сделано и в скриптах фемма). L(x) скорее всего не будет.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 06-03-2009 13:27
1. В формуле I(t) параметр U это напряжение кондера в начале процесса или функция U(t)?
2. Как вычислять Lmin и Lmax, и какое среднее значение брать?
3. Для вычисления определенного инеграла юзать формулу трапеций? и с какой точностью вычислять?
4. Как прога будет находит производную и неопределенный интеграл?

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 06-03-2009 14:37
1) напряжение в начале процесса
2)А разве там не написано? Индуктивность - это полный магнитный поток деленый на силу тока. Полный магнитный поток - это сумма потоков через все витки (т.е. через их площади). Вот риальне проблема с сердечником - ведь он насыщается постепенно, и определить индукцию в каждой точке для случая с сердечником я чёйто не могу. Среднее значение... а попробуй арифметическое для начала, потом будет видно, куда нужно сдвигать среднее
3)Мне даже метода прямоугольников хватало. А точности сейчас хватит даже 1%
4) Это нам, людям, всё выводить надо (интегралы, пятна и т.д.), а компьютеру достаточно способности быстро считать. Привожу два примера с их решением.
пр.1. Нужно найти I = int (t1 to t2) (int (x1 to x2) f(x)dx) dt , двойной интеграл. Что можно сделать:
I=0
I2=0
For t=t1 to t2 step 0.01
I2=0
For x=x1 to x2 step 0.01
I2=I2 + f(x)*0.01 ' в интегральную сумму добавляется площадь прямоугольника
Next x
I=I+I2*0.01
Next t
В I соответственно имеем ответ.
пр. 2. Если появляется ацкое желание найти производную функции y=f(x) в точке x0, нужно найти предел приращения функции к приращению аргумента в этой точке. Что- то типа такого:
y'(x0)=(y(x0+0.001) - y(x0))/0.001

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 06-03-2009 19:55
чет ни хуище не пойму че куда подставлять.

В формуле магнитного потока под неопределенным интегралом параметр B, а интеграл по идее от функции берется - тож непонятно

I=0
I2=0
For t=t1 to t2 step 0.01
I2=0
For x=x1 to x2 step 0.01
I2=I2 + f(x)*0.01 ' в интегральную сумму добавляется площадь прямоугольника
Next x
I=I+I2*0.01
Next t

ВО! вот так со мной разговаривай чтоб понятно было
это вычисление определенного интеграла по формуле прямоугольников, а как неопределенный?

Можешь формулы в цепочку разложить? Чтоб были формулы от исходных параметров до конечных результатов(а где еще нет звеньев то просто какую нибудь заглушку), и указать список функций которые нужно отобразить в графиках - как бэ чтоб алгоритмы можно было составить.

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 06-03-2009 22:47
:ржу нах с предыдущего поста:
В формуле магнитного потока под неопределенным интегралом параметр B

Как бы сказать... расчет состоит не только из функций, но и некоторых значений, вычисляемых по сложным формулам. Если их описывать как функции от нескольких переменных, то этож пальцы отсохнут, да и нахер это не надо. В магнитом потоке там вообще B(H), если тебе так удобней. Да, по поводу неопределённого интеграла - я думал ты сможешь разобраться. Вообщем, Ф=int(0 to R1)(int(0 to 2pi)(B(H)dfip))dR, т.е. dS=d(fip)dR. Понимаешь?
а как неопределенный?

Неопред. интеграл - это совок-ть функций. Его нельзя вычислить. А записано как неопределённый потому, что интеграл не по переменной, а по площади, имеющей пределы ;)
Можешь формулы в цепочку разложить? Чтоб были формулы от исходных параметров до конечных результатов(а где еще нет звеньев то просто какую нибудь заглушку), и указать список функций которые нужно отобразить в графиках - как бэ чтоб алгоритмы можно было составить.

Бомба - этим я и займусь позже - упорядочиванием известного.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 06-03-2009 23:41
давай знаешь как сделаем - я все рассчеты вот эти буду в отдельной dll-шке писать на vb.net а потом как плагин подключать (я уже добавил такую возможность), чтоб мы на одном языке говорили, а инетрфейс/двиг на С# останется.

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 07-03-2009 06:28
Да в принципе с программной точки зрения делай как считаешь нужным.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 08-03-2009 09:16
Помоги с формулой разобраться, она нужна для полиномиальной интерполяции

Я так понял под знаком произведения стоит k неравно j, и что это значит?

Ф=int(0 to R1)(int(0 to 2pi)(B(H)dfip))dR, т.е. dS=d(fip)dR. Понимаешь?

нет. R1 это внешний радиус катушки? d(fip)dR это что? А Н откуда брать? Не, чувствую дело с места не сдвинется пока цепочки не будет ибо гавеный я математик.

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

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 08-03-2009 16:49
Слухай, мне проблематично цитировать, поэтому по-другому буду отвечать.
1)это значит, что k от 0 до n, перепрыгивая через значение j
2)Радиус не катушки, а слоя
3)d(fip)dR - это произведение дифференциалов
Короче, жди цепочки.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 08-03-2009 19:35
UPDATED

Вообщем ниче так показали себя интерполяционый полином Лагранжа и интерполяция по алгоритму Невиля, но проблеема в том что если большой разрыв между точками, то график между ними впукливается или выпукливается сильно. Барицентрическая интерполяция гавно.


На графиках интерполяция Невилля
<p x="1" y="1E-08" />
тут просидает ниже оси ОХ
<p x="510" y="0,100015744" />
<p x="750" y="0,199993788" />
<p x="950" y="0,499953055" />
<p x="1200" y="0,900003463" />
<p x="1400" y="1,099997252" />
<p x="1500" y="1,19996273" />
<p x="1750" y="1,300053872" />
<p x="2250" y="1,399956518" />
<p x="3250" y="1,499984828" />
<p x="6000" y="1,599950307" />
тут вообще разрывы и скачки, ибо смотри какой пробел между иксами
<p x="4130000" y="13,67095459" />




<p x="1" y="1E-08" />
тут просидает ниже оси ОХ
<p x="13,8984" y="0,227065" />
<p x="27,7967" y="0,45413" />
<p x="42,3974" y="0,681195" />
<p x="61,4157" y="0,90826" />
<p x="82,3824" y="1,13533" />
вот тут выпукливается вверх
<p x="144,669" y="1,36239" />
опять же разрывы
<p x="897,76" y="1,58935" />
<p x="4581,74" y="1,81236" />
<p x="17736,2" y="2,01004" />
<p x="41339,3" y="2,13316" />
<p x="68321,8" y="2,19999" />
<p x="95685,5" y="2,25479" />
<p x="123355" y="2,29993" />
<p x="151083" y="2,34251" />
<p x="178954" y="2,37876" />
<p x="206825" y="2,41501" />
<p x="234696" y="2,45126" />
<p x="262568" y="2,4875" />
<p x="290439" y="2,52375" />
<p x="318310" y="2,56" />
<p x="1000000" y="3,15" />
<p x="10000000" y="11,2" />


ну думаю можно будет заполнить пробелы точками, взятыми из интерполяции сплайнами.


Теория:

Полиномиальная интерполяция является наиболее известным из методов одномерной интерполяции. Её достоинствами являются простота реализации и хорошее качество получаемых интерполянтов. Впрочем, она не лишена недостатков (некоторые из них будут рассмотрены ниже), и в последнее время испытывает сильное давление со стороны альтернативных методов интерполяции: сплайнов и рациональных функций. Но, несмотря на это, полиномиальная интерполяция по-прежнему остается одним из главных инструментов численного анализа.
_______________

Методы полиномиальной интерполяции

Простейшей формой интерполяционного полинома, приводящейся в начале каждого учебника по численному анализу, является интерполяционный полином Лагранжа:



Эта форма интерполяционного полинома интуитивно понятна, но не очень подходит для практического применения. В практике численного анализа для полиномиальной интерполяции обычно используется алгоритм Невилля. Этот алгоритм представляет полином N-ой степени P 0,1,..,N-1,N , проходящий через N точек (с 0-ой по N-ую), как функцию двух полиномов N-1-ой степени по формуле



К полученным полиномам рекурсивно применяется та же формула, до тех пор, пока мы не дойдем до полиномов вида Pi , которые вычисляются по формуле Pi = yi . Этот алгоритм позволяет вычислить значение интерполяционного полинома N-ой степени за O(N 2) действий. Достоинством алгоритма Невилля является простота реализации и тот факт, что алгоритм легко может быть расширен для вычисления не только значения интерполяционного полинома, но и его производных. Недостатком - сравнительно невысокое быстродействие.

Следует отметить, что алгоритм Невилля не является последним словом в области полиномиальной интерполяции. Из новых достижений в этой области следует отметить методы, основанные на т.н. барицентрическом представлении интерполяционного полинома:



В общем случае приведенная выше формула задает рациональную функцию, проходящую через точки (xi , yi ) - независимо от того, какие весовые коэффициенты wj мы выбираем. Таким образом, эти коэффициенты влияют только на то, как интерполянт приближает интерполируемую функцию между точками. Конкретно для нашей задачи важны два факта. Во-первых, то, что при выборе



рациональная функция r(x) становится полиномом N-ой степени. А во-вторых, то, что в некоторых важных частных случаях (если интерполяция осуществляется на равномерной сетке или на Чебышевских узлах) существует алгоритм, который может вычислить все wj за O(N) операций. Таким образом, задача вычисления значения интерполирующего полинома решается в N раз быстрее, чем алгоритмом Невилля.

Замечание #1
Следует отметить, что прямолинейное применение описанной выше формулы барицентрической интерполяции может привести к неприятным последствиям. Если xj = 0, а x находится слишком близко к xj (но не равно ему), то возможно переполнение из-за деления на слишком малое число. В практической работе такие случаи возникают очень редко, но все же по какой-то необъяснимой причине это происходит намного чаще, чем это могло бы следовать из теории вероятностей. Но можете не беспокоиться: всё это было напечатано лишь для того, чтобы вы знали - алгоритм, реализованный в этом модуле, корректно обрабатывает такие случаи, переключаясь при необходимости на устойчивый к погрешностям переполнения вариант формулы.

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 08-03-2009 21:30
exe-dealer, а можно я буду звать тебя Exed? Как бы дружеское сокращение. Цепочку сейчас начал создавать, но не думаю, что сегодня её и напишу. В случае чего потом как-нибудь на мыло закину.

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 08-03-2009 22:05
exe-dealer, а можно я буду звать тебя Exed?

да называй

Цепочку сейчас начал создавать, но не думаю, что сегодня её и напишу.

ога, буду ждать. Я во вторник в таганрог уезжаю аттестат забрать из института, так что меня дня 3 со вторника не будет.

магистр
Группа: Участники
Сообщений: 648 		Добавлено: 09-03-2009 15:02
аттестат забрать из института

???

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 09-03-2009 17:10
ну я в мухосранск свой вернулся жеш, втопку вуз

магистр
Группа: Участники
Сообщений: 882 		Добавлено: 10-03-2009 09:58
Ну зачем свой родной город называть так пренебрежительно (я бы даже сказал обзывать не хорошими словами). Ты что его так не навидиь? Или это всё влияние Маасквы и маасквичей, которые все остальные города (если город не столица) называют мухосрансками, а людей от туда считают второсортными?

магистр
Группа: Участники
Сообщений: 3209
 Добавлено: 10-03-2009 12:04
ыыыТЬ!Ё не все,только самые беспринципные,коих 90+%

магистр
Группа: Участники
Сообщений: 1586
 Добавлено: 10-03-2009 15:29
Ну зачем свой родной город называть так пренебрежительно (я бы даже сказал обзывать не хорошими словами). Ты что его так не навидиь? Или это всё влияние Маасквы и маасквичей, которые все остальные города (если город не столица) называют мухосрансками, а людей от туда считают второсортными?

Да что родной не родной, один хер, мне везде мухосранск ибо справедливости нет нигде, а кругом одно лицемерие и конформизм.
И детали тут хер найдешь, и компьтеры дорогие, и инетернет дорогой, а зарплата средняя здесь тыщ 5, и гопов в лице чурок столько, что хоть из дома не выходи... и быдлы эти ебани повсюду и ржут все время полуумным смехом.

Страницы: << Prev 1 2 3 Next>> новая тема
Раздел: 
Gauss2k - gauss gun у тебя дома / Gauss tecnology / проект D&E

Отвечать на темы могут только зарегистрированные пользователи

KXK.RU