проект D&E

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

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

новая тема

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

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 14-12-2008 22:48
здесь мы с Джином тихо в сторонке будем согласовывать работу над проектом. Дельные советы приветствуются.

магистр
Группа: Участники
Сообщений: 2148
Добавлено: 14-12-2008 22:50
общий проект это как?

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 14-12-2008 22:51
прогу делаем для рассчетов гаусса.

магистр
Группа: Участники
Сообщений: 2148
Добавлено: 14-12-2008 23:02
Тю ёпт. Я думал гаусс собираете. Ну если чё - я ж тож немного шарю в программировании. Что-нить когда-нить подскажу\посоветую.
Розскажи хоть основную соль программы...

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 14-12-2008 23:13
Ну там будет по сути мат. анализатор гаусс-процесса, т.е. без дескритизации времени как в фемме, а подстановкой значений в формулы Джина, которые чет та там в свою очередь выдают как результат другие формулы для конкретного гаусса. Короче это к Джину, я еще пока не доконца разобрался.

update
ах да, и там еще она быстрее фемма считать будет и точнее за счет гладкой мат. модели.

магистр
Группа: Участники
Сообщений: 648
Добавлено: 20-12-2008 14:46
Если делать "точнее", понадобится оочень много вычислений. А вот претендовать на создание очень быстро считающего прикидочного гаусскалькулятора мы можем (с точностью +-10%) :) Хотя впоследствии может и возможность точных вычислений подошьем.

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 21-12-2008 04:58
Хотя впоследствии может и возможность точных вычислений подошьем.

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

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

мастер
Группа: Участники
Сообщений: 109
Добавлено: 21-12-2008 06:23
этих прог по подсчетам многа и вообще я им не верю

http://www.coilgun.info/mark2/rlcsim.htm
http://www.gausspistol.com/Simulator.html

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 21-12-2008 07:45
коил ган симулятор, имхо, это так поржать просто: интерфейс гавеный и к тому же ограничение ввода кривое(десятичная точка не вводится), а постоитель графиков там строит график по переданому ему массиву значений, что не есть правильно. Да и много команд меню не пашет.

магистр
Группа: Участники
Сообщений: 2148
Добавлено: 24-01-2009 00:34
Ну что там?

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 24-01-2009 03:25
да хз Джин пропал куда. Интерфес и двигун готов, осталось заполнить библиотеку формулами.

магистр
Группа: Модераторы
Сообщений: 861
Добавлено: 24-01-2009 08:44
а жаль. так хорошо всё начиналось............

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 26-01-2009 01:26
дa не, придет он, щa сeсия вот он и нe зaxодит нaвeрно.

магистр
Группа: Участники
Сообщений: 648
Добавлено: 28-01-2009 13:29
1) анализ застопорился за отсутствием нормально работающих интернета и компа => не могу вести быстрый поиск нужной инфы и пользоваться феммом.
2)А может нуихнахер, вычисления? Может сделать прогу, выдающую оптимал в форме советов? Ведь для получения чисел вроде и ФЕММ'а хватает.
3) exe-dealer, а все-таки формулы нужно из опытов выводить, а не из глубокой теоретической жопы. А еще лучше делать как Чаплыгин - рыть тоннель с двух концов. Сдается мне, нужно снова скачать ФЕММ и сделать кое-какие теоретические расчеты.

Кстати, читал наше начало организации проекта, и вычитал... "испугать формулами"... Если надо, могу прислать то, что имею на данный момент - "испугаешься" наатличненько. Две основные формулы - это для силы тока и для магнитного поля. Это конечно не формулы, вычисляющие конечный результат, но формулы, являющиеся одними из основных.
Да, и чуть не забыл. Возможно B-H-кривая вообще будет не нужна - лишь только индукция насыщения... но это еще посмотрим.

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 28-01-2009 15:27
рaбoтaть пoкa нe могy ибo компa нeтy - я домa нa кaникyлax, a oн в oбщaгe ocтaлcя - сижy тут c мoбилкой дыркaюсь.
нaсчeт 2 нeпoнятнo.

магистр
Группа: Участники
Сообщений: 648
Добавлено: 28-01-2009 15:30
Я с мобилы вообще не пишу, хотя заглядываю на форум по несколько раз в день.
Насчет 2 - ну типа "провод для лучшего КПД желательно взять 0.9" или "Для достижения лучшей скорости необходимо 10 слоев по 21 витку в каждом. Лучший диаметр снаряда - 7мм" Что-то типа того. Короче, не будет конечных чисел (скорость, КПД, и т.д), а будут только те, при которых достигается луший КПД/лучшая скорость.

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 28-01-2009 15:43
ну мoжнo добaвить кaк eщe oдин тип анaлизa, тогдa придeтся с фeнoм комбинировaть.

магистр
Группа: Участники
Сообщений: 648
Добавлено: 28-01-2009 15:49
Ладно, пока у меня некоторые загвоздки, обсуждать проект бессмысленно. На данном этапе необходимо проанализировать кривые ускорения, скорости и перемещения от времени. Там видно будет.

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 15-02-2009 18:11
дай ка глянуть наработки (если не гемороно будет с пересылкой)

магистр
Группа: Участники
Сообщений: 648
Добавлено: 15-02-2009 20:20
Материться нет сил. Набивал формулы в ворде (причем на сегодня и анлим проплачен, мог бы с компа по почте переслать), никого не трогал, и тут - хуякс! и вырубается ганд... компьютер. Больше включаться не хочет. Думаю, в лучшем случае напишу тебе через неделю.

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 15-02-2009 20:28

как графики отображаться должны - по одному на координатную ось или несколько на 1 координатной оси?

магистр
Группа: Участники
Сообщений: 648
Добавлено: 15-02-2009 20:57
Может не оси а плоскости? Слушай, насчет графиков - очень понравилось как графики сделаны в гаусссимуляторе. Наводишь на кривую - и тут же значение функции в этой точке всплывает. А вообще вопрос ко всем: хорошо ли на ваш взгляд устроено построение графиков в гаусс симуляторе?

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 15-02-2009 21:03
Может не оси а плоскости?

ну ты меня понял. так как по одну графику или все вместе?

Наводишь на кривую - и тут же значение функции в этой точке всплывает.

хорошо, сделаю



магистр
Группа: Участники
Сообщений: 648
Добавлено: 15-02-2009 21:11
Эх, да вместе все, еще цветами их выделять, а сбоку где-нибудь легенду указывать, вроде красное - скорость, голубое - сила и тому подобное. Вопрос тогда один возникает - как размечать ось Oy?

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

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

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 15-02-2009 21:17
что нормально, как делать?

магистр
Группа: Участники
Сообщений: 648
Добавлено: 15-02-2009 21:25
Нормально в смысле нормально придумал - комбо бокс почему-то представляется над графиком. Тогда выделение цветом идет в топку вместе с легендой. При переключении между выводимыми графиками будет меняться разметка вертикальной оси.
Что-то всё-таки не так. Говорю же надо в гаусс ган симуляторе посмотреть как сделано.

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

там над осью у так и написано просто "ось y"

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


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

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

магистр
Группа: Участники
Сообщений: 648
Добавлено: 16-02-2009 19:00
Почему проца много жрёт? Что-то не понятно. Насчет точек - выводить точки - то есть писать значение функции в точке? Выражайся правильно! По выводу точек - подумаю. Если сильно грузит проц - значит надо сохранять вычисленную таблицу значений. Но тогда недалеко от ступенчатости уйдем... Так почему проц сильно грузится понял или нет? Если избежать его сильной загрузки, то можно было бы... интерполировать?! надо подумать. Хотя вообще ладно. Можно и без всплывающего значения функции.

магистр
Группа: Участники
Сообщений: 508
Добавлено: 16-02-2009 19:18
он и будет проц грузить, куда он денеться, у тебя сколько ходов по циклу получается? правильно дохуя+ вычисления в теле цикла+ сторонние прибамбасы + умножение вычисляется сложением набигает нормально. апроццессор он свои гигагерцы переплюнуть не может, или ты думал, что мега компы для инженерного моделирования процессов они так писюнами померяться?)))))))

магистр
Группа: Участники
Сообщений: 648
Добавлено: 16-02-2009 19:48
Ну ходов... один литой процесс, какие ходы...вычислений много из-за кратных интегралов. Гигагерцы не переплюнет - его в принципе и не заставляют. Насчет 'ты думал' - ДА ЙА ТАК И ДУМАЛЬ А ЧИТО НЕТ АТМО ЛЬБУЛЬБУАБЛЬСЯ ЫХШЫФЫШХ ШШШШ...
Если серьёзно - то вообще то я удивляюсь из-за того что у Эксе-Диллера даже нет моих наработок а нагрузить проц нашей прогой он уже умудрился

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 16-02-2009 19:51
Если сильно грузит проц - значит надо сохранять вычисленную таблицу значений. Но тогда недалеко от ступенчатости уйдем... Так почему проц сильно грузится понял или нет?

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

щас думаю займусь редактором ферромагнетиков. Какие параметры надо вводить для материала?
-плотность
-индукция насыщения
-BH пары
правильно?

магистр
Группа: Участники
Сообщений: 648
Добавлено: 16-02-2009 19:55
Да
Наводить мышь на график как снайпер? А знаешь, наверное всё-таки это лишнее. Зря я ляпнул.

магистр
Группа: Участники
Сообщений: 2148
Добавлено: 16-02-2009 19:59
А вы уверенны что прога будет считать правильно?

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 16-02-2009 20:02
никто ни в чем не уверен

магистр
Группа: Участники
Сообщений: 648
Добавлено: 16-02-2009 20:18
А вы уверенны что прога будет считать правильно?

Вы настоящие оптимисты и друзья - профессионально вдыхаете неуверенность в любого кто хочет что-то полезное сделать.

магистр
Группа: Участники
Сообщений: 2148
Добавлено: 16-02-2009 20:26
стараемся
ЗЫ я за графики в новой проге, как в гауссимуляторе

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

магистр
Группа: Участники
Сообщений: 648
Добавлено: 17-02-2009 07:21
Графики нужны для оптимизации системы. Так, видишь, что сила, тянущяя назад, слишком большая, значит надо укорачивать импульс - либо емкость кондёра поменьше, либо сопротивление катухи как-то поменьше. И так в принципе с графиком тока, напряжения, скорости, еще подумал интересный параметр - степень насыщенности снаряда, только я обломаюсь его считать. Вообщем, по графикам можно понять что не так, ПРОАНАЛИЗИРОВАВ ВРУЧНУЮ... Так! Опять мне думается можно замутить советную форму вроде 'Уменьшите колличество слоёв на ..' или 'Попробуйте постепенно уменьшать колличество слоёв в катушке, пока не получите оптимал'. А можно и без советной формы - все автоматически - программа сама будет вбивать то, что приблизит к оптималу, и заново считать. В принципе это замутить можно, главное - чтоб все расчетные формулы были. И будут

частый гость
Группа: Участники
Сообщений: 23
Добавлено: 18-02-2009 08:40
а матлаб здесь не к месту? там ещё есть и симулинк... кстати еси комуто нужен матлабец могу ключик выдать на 6.12

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 20-02-2009 00:42
Джин, как это ты насчитал что для хранения матрицы для интерполяции надо 100кб?

магистр
Группа: Участники
Сообщений: 648
Добавлено: 20-02-2009 07:22
В формулах напряженность H, максимум 6-го порядка, возводится в макс. степень (n-1). Если предположить, что пар вбивается не более 30, то максимальная степень получается 29. Итого степень сохраняемого числа 6*29=174 Учитывая, что 2^10~10^3 , получаем (10^3)^58~2^580 ; примерно 73 байта (я считал по 128). Всего чисел 30*31=930 На всю матрицу нужно 930*73=67890 байт.
Слухай, в программировании ты больше меня понимаешь, чета я напутал с форматом real - по-моему он не 128 байт , хыхыхы...

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 20-02-2009 07:51
чета я напутал с форматом real - по-моему он не 128 байт , хыхыхы...

ты хотел сказать бит? в dotnet 128 бит есть Decimal, его то я юзаю(тот что с 26 знаками после запятой).

если всего чисел 930 , то 930 * 128бит = 119040 бит = 14880 байт

магистр
Группа: Участники
Сообщений: 648
Добавлено: 20-02-2009 08:16
Нене, я именно байт имел ввиду, когда писал 'байт' :) Вообщем, нужны такие переменные, куда будет вмещаться до 10 в 180й степени. Кстати в супероптимизации интерполяция в нашем случае не нуждается - она будет выполняться всего один раз для одного материала; потом мы уже будем иметь формулу для связи B и H до индукции насыщения; после насыщения - график B(H) есть зависимость индукции от напряжённости в вакууме, "сдвинутая" вверх на индукцию насыщения материала.

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 20-02-2009 09:09
чтоб хранить 10E18 совсем не надо кучу байт, - 128 битный Decimal принимает значения от ±1.0*10E-28 до ±7.9*10E28.
http://ru.wikipedia.org/wiki/Числа_с_плавающей_точкой

магистр
Группа: Участники
Сообщений: 648
Добавлено: 20-02-2009 09:17
чтоб хранить 10E18 совсем не надо кучу байт, - 128 битный Decimal принимает значения от ±1.0*10E-28 до ±7.9*10E28.

ПЛЯЧУ - НЕ 10Е18 А 10Е180 !!!

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 20-02-2009 09:25
тогда 64 битный Double заюзаю, но точность снизиться до 15 знаков

магистр
Группа: Участники
Сообщений: 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, и гопов в лице чурок столько, что хоть из дома не выходи... и быдлы эти ебани повсюду и ржут все время полуумным смехом.

магистр
Группа: Участники
Сообщений: 648
Добавлено: 12-03-2009 13:04
Сама цепочка получается нехитровыежопездонистая, но вычисление индуктивностей и средних величин это ужас. Позже оформлю - перешлю. Черновой вариант ужо есть.

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 02-04-2009 01:40
какбэ чочо там с цепочкой то?

магистр
Группа: Участники
Сообщений: 648
Добавлено: 02-04-2009 09:33
Чочо упячка следит за нами, а цепочка... Ну нахер она нужна без основных звеньев? Нет энтузиазма её писать + времени нет бо планировать нихера не умею + за карточкой нетовской ищо идти надо + с самим соединением инетовским на уровне железа что-то не то. Проект замёрз :(((

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 02-04-2009 20:22
ааа, ну подождем.

магистр
Группа: Участники
Сообщений: 648
Добавлено: 18-04-2009 17:09
Только что глянул на редактор материалов. С твоей стороны все сделано отлично, а с моей... нет комментариев бо испацтала сложно комментировать

бывалый
Группа: Участники
Сообщений: 46
Добавлено: 22-04-2009 17:25
Народ, а можете хоть на черновой вариант проги (только б работала) ссылку выложить?

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 22-04-2009 23:58
она не считает ничего, только интерфейс написан, а плагинов для рассчетов нету... все еще надо?

какбе скрин

магистр
Группа: Участники
Сообщений: 1115
Добавлено: 23-04-2009 04:51
W@lk используй FEMM он уже проверен и легок в использовании. Ссылку можешь найти в спец. разделе.

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 23-04-2009 05:28
ога, и время приращения пониже ставь

бывалый
Группа: Участники
Сообщений: 46
Добавлено: 23-04-2009 07:49
Пасиб. А в FEMM'е тоже графики есть?

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

магистр
Группа: Участники
Сообщений: 648
Добавлено: 26-04-2009 14:27
Я вот чего подумал - фемм неплохо расчитывает гауссы. А сейчас я занимаюсь несколько иным - расчет магнитной пружины, о которой расскажу позже. Занимаюсь кстати со всех сторон - продвигаю одновременно и практику, и теорию.
И подумал: Exed, наверное мы впоследствии перепишем немного прогу - будет считать не простые гауссы, а... ну те, с активным снарядом...они более эффективные :) уж поверьте . Не буду торопить время, я потихоньку, но продвигаюсь в этом направлении - позже все будет, с картинками, с пояснениями и фотками.

магистр
Группа: Участники
Сообщений: 3209
Добавлено: 26-04-2009 14:35
с активным снарядом...они более эффективные :)

цитата с нового форума,в ответ на точно такую-же мою реплику:
А я наоборот, думаю, что никого толку от снаряда-катушки не будет...

http://gauss.getbb.ru/viewtopic.php?f=3&t=602&start=125

магистр
Группа: Участники
Сообщений: 648
Добавлено: 26-04-2009 15:22
цитата с нового форума,в ответ на точно такую-же мою реплику

Отписался. Беседы в кальянной продолжаются!

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 26-04-2009 15:30
наверное мы впоследствии перепишем немного прогу

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

магистр
Группа: Участники
Сообщений: 648
Добавлено: 26-04-2009 15:38
BH таблиц не будет! Равно как и железа!

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

магистр
Группа: Участники
Сообщений: 648
Добавлено: 26-04-2009 16:32
Хитро. Рад что когда-то обратился к тебе за помощью.

бывалый
Группа: Участники
Сообщений: 46
Добавлено: 28-04-2009 20:08
Ну дык как, будет прога?:-) Как начнёт работать выложите плиз ссылку

магистр
Группа: Участники
Сообщений: 648
Добавлено: 05-10-2009 20:08
Exed, не против продолжения работы?

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 05-10-2009 23:13
вообще некогда щяс(
от работы и так устаю шо ппц.
ну а если работы не будет, то думаю можно стереть пыль и вспомнить то что делал.

у тебя как? цепочка готова? или что щяс есть?

магистр
Группа: Участники
Сообщений: 648
Добавлено: 06-10-2009 16:56
Па суте есть все кроме итоговых формул, т.е. расчёт делается как в фемме - пошагово. Или еще есть метод вычисления при помощи сходимости ряда результатов. Проверял оба. Вродеб работают. Кстате перешёл на с++ борланд 6.0.

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 06-10-2009 21:55
расчёт делается как в фемме - пошагово.

а вот это уже не влезет в то что я делал(
изначально то задача стояла написать интерфейс для проведения быстрых рассчетов.
собсно все однопоточное( и чувствую если туда поять полезу, то с новыми знаниями что накопилось за эти полгода появится непреодолимое желание снести половину и написать по другому))

Кстате перешёл на с++ борланд 6.0.

брось каку, если под винду пишешь то на vs с++ лучше вроде.(только не под .net)
только для экспериментальных рассчетов не годится ни то ни другое. нужно думать о матиматике больше а не о том как перевести строку из юникода в ansi, правельно выделить память и т.д. для математики лучше юзать бэсик или паскаль.

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

магистр
Группа: Участники
Сообщений: 648
Добавлено: 07-10-2009 07:07

а вот это уже не влезет в то что я делал(
изначально то задача стояла написать интерфейс для проведения быстрых рассчетов.

Там все рассчёты в циклах пробегаются. Развеж все равно не влезет? А переписание программы с новыми знаниями - это вродеб тру?
Кстате перешёл на с++ борланд 6.0.

брось каку

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

Ога

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 07-10-2009 08:21
Там все рассчёты в циклах пробегаются. Развеж все равно не влезет?

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

А переписание программы с новыми знаниями - это вродеб тру?

тру то оно тру, вот только дата релиза всего этого заметно отдалится

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 27-10-2009 18:23
на прошлой недели родичи в астрахань ездили, а я не смог вместе с ними. печально кагбе( работаю блин, даж на форум зайти некогда щяс

магистр
Группа: Участники
Сообщений: 648
Добавлено: 27-10-2009 23:36
У меня как обычно круговоротный секс: чета матчасть есть, но прогу(сырье) пишу оочень медленно. Исчо надобно получше изучить сабж, а для этого его еще надо до ума довести. Это упирается как у всех в деньги. Но ничего - в последнее время меня так это заи.. достало, что решил плюнуть на часть принудиловки универа и продолжить развивать ацкие поделки.

магистр
Группа: Участники
Сообщений: 648
Добавлено: 03-03-2010 16:50
...продолжить развивать ацкие поделки.

Мое нормальное состояние
Exed, тут такое дело. Такая вот прога будет нужна - на конкурс в мае. Я постараюсь прислать формулы на этой неделе - замутишь прогу, накладно по времени не будет? Я конечно буду указывать тогда, что программер - такой-то человек из такого-то города, идея моя. Есть возможность выиграть бабло, коим клянусь поделиться в случае выигрыша по-братски.

магистр
Группа: Участники
Сообщений: 1586
Добавлено: 29-03-2010 15:48
мдя, поздно я зашел)) со временем вообще туго стало(

гроссмейстер
Группа: Участники
Сообщений: 206
Добавлено: 01-09-2010 16:00
проца жрет чет много
и медленно работает из-за того, что происходит непрерывный OLE-обмен с Экселем. Из-за этого и скорость... Примерно 99,99% ресурсов идёт на OLE-обмен между программами, а также на непрывную работу интерпретатора VBA, который, как любой скриптовый язык пасится во время выполнения. На компе с памятью в гигабайт и процессором 1,2 ГГц я запустил, пошёл курить, вернулся, не было выполнено и одного процента и выключил. Возможно, я что-то не так делал. :(

магистр
Группа: Участники
Сообщений: 648
Добавлено: 05-09-2010 17:54
И века не прошло как цепочка уже появилась! Но для лайнера. Сейчас я сам потихоньку возвращаюсь в программирование, потому вообщем-то сам смогу все написать и отладить. Тем более я еще сомневаюсь в верности базовых формул, но зато у меня есть самое главное - понимание хода вычислений! Уйти от пошагового моделирования так и не получилось


Группа: Участники
Сообщений: 8
Добавлено: 20-05-2011 14:51
Здравствуйте!
Интересуюсь той же темой, можно посмотреть Ваш результат?

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

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

KXK.RU