iiv

во-первых, скалапак и лапак - разные вещи, если под линуксом, поставьте ACML и сам лапак с сорсов. Доступ к хедерам бласов и большей части лапака у вас будет из АСМЛя, а вот если надо будет совсем что-то извращенное (из 3 и более поздней версии лапака) то только ручками фортрановские вызовы придется дергать, без хедеров. Зато - это наиболее оптимальное по скорости решение. Атлас существенно тормозит. Если поставите бласы с лапаком с сорцов, скорость будет совсем никакая. Под МССтудио - однозначно АСМЛ, а если его не хватит, то только фортран ставить, или на цугвин переходить. Разумной альтернативы лапаку, к сожалению, до сих пор нет... ЗЫ: не пугайтесь, что АСМЛ - продукт АМДшника - он и на Интеле очень шустро работает, и с МатКернелЛибом (MKL) от интела по скорости даже на интеле почти не отличается. АСМЛ - софт - бесплатный, берите с сайта производителя.

Позвольте мне, пожалуйста, задать Вам еще один для меня пока еще очень актуальный вопрос! Если в дросселе максимально запасаемое количество энергии будет примерно равно энергии, запасаемой в каждом конденсаторе умножителя, например, при моих конденсаторах в 0.12мкФ на 1кВ, я возьму дроссель, в котором запасается 0.05Дж до насыщения, например, T200-2 от micrometals, и выведу это все на 50кГц, то смогу ли я работать и на не сильно высокой частоте, и на довольно большой мощности, чтобы и конденсаторы не горели и любой сильный всплеск с индуктора ими гасился? Заранее прошу прощение, если мой вопрос очень глуп, тогда, просто на меня за это не обижайтесь, пожалуйста! С уважением ИИВ

Я думал, что линейное предсказание, но, похоже, ошибался. Если представить $a_i$, $b_i$ как входные векторы, то результат в матричном виде и латехе можно записать примерно так: $$\left( \product_i (a_i I + b_i P) \right) \sum_i (a_i I + b_i P)^{-1} e,$$ где $e = (1, 0, ..., 0)^T$, а P матрица сдвига, как я сказал выше. Правда мне такое представление ничего не дало чтобы упростить или как-то серьезно перекаверкать этот алгоритм...

Да, действительно, каша в голове была, помоделировал разные ферриты, понял разницу, о которой спрашивал... короче для дросселя лучше ironpowder брать, а для транса - 3E5 или аналогичное... Собственно начал-то спрашивать из-за того, чтобы вычислить какую максимальную энергию можно запасти в феррите, и именно эта величина меня больше всего волновала, а в даташитах ее, к сожалению, не оказалось. Получилось, что для 3C90 при габаритах в 5см тора в нем можно запасти около 100мкДж, а в поделке от micrometals тех же габаритов - аж 0.1Дж. Огромное спасибо Тау, и всем-всем, кто наставлял на путь истинный! С уважением ИИВ

одна фирма - оффициальный CUDA Nvidia консультант год назад и с Альтерой и с Ксилинксом и с Латтисом этот вопрос обсуждали - в портфолио кудовской фирмы - под двадцать лет работы на массивно-параллельных платформах. Вопрос обсуждался на уровне деректоров по развитию. Решили забить, так как игра не стоит свеч - те кому надо - купят побыстрее проц и SSD, а остальные подождут.

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

a1 инициализирован только в первых 4 элементах. А дальше, в общем случае, тоже нули или модель просто не гоняли? Циклы в общем случае по i,j,k всегда до одного и того же N бегут? С виду очень походе на линейное предсказание с преобразованием исходных реальных векторов в комплексные, но и на сто процентов сказать не могу, если есть еще какая-то информация об алгоритме, говорите!

Основная идея почему я захотел уйти от полного моста в сторону буста вот в чем: у меня на входе в систему переменный ток 50Гц, который выпремляется до полупостоянного 100Гц. Так как мне надо заряжать конденсатор нагрузки не чаще 100 раз в секунду (обычно даже 20 раз в секунду достаточно), я хотел начинать зарядку слабым напряжением в момент когда в сети 0, в этом случае, система работает совсем в щадящем режиме, а потом, по наростанию напряжения на входе, поднимать и выходное напряжение. В этом случае есть надежда не пробить номиналы столбов и мосфета, а также не делать больших пульсаций на конденсаторах. Как я понимаю, с полномостовой схемой у меня такое не пройдет, так как + и питание p-мосфетов у меня должно будет гулять 100 раз в секунду, и это может как-то не хорошо сказаться на работе системы. По номиналам - что нарыл - добавлю, я последнее время пристрастился покупать через Арроу, у них пересылка только 24 бакса, в результате, тот же AVX с 120нФ 1кВ удалось достать по 23 цента + копейки за пересылку + растаможка

Именно вариант буст с дросселем + длинный умножитель (20-30 ступеней) я и хотел сравнить с вариантом полный мост с коротким умножителем (2-3 ступени).

Уважаемый Микроватт, очень Вам благодарен за простое и понятное разъяснение. Меня ввела в ступор фраза со схожего ресурса: me — относительная магнитная проницаемость сердечника (не путать с проницаемостью материала сердечника) http://ferrite.com.ua/site/page-Trancformatori_i_drocceli из-за которой я никак не мог сопоставить начальную проницаемость с относительной проницаемостью. Все-таки я что-то все еще не понимаю... помогите, пожалуйста, разобраться. Именно на том феррите, о котором я в начале говорил, я как-то собирал импульсник по полномостовой схеме, для повышения 12в в 120в, причем первички было 13 витков, вторички было 130, первичка была намотана 1мм проводом. Я тогда через него совершенно спокойно прокачивал 16А. Подставим I=16A и N=13, µ0=1.257*10-3 [мкГн/мм], le=89.7mm, и µe=2300 в B = 1000 * µ0 * µe * I * N / le получим 1000*1.257*10-3*2300*16*13/89.7=6704мТ, что примерно в 15 больше чем 430мТ. Скажите, пожалуйста, что я здесь снова не понял?

не, 1кВт - среднее потребление нагрузки, которое вначале, когда конденсатор нагрузки разряжен, может быть сколько схема сможет дать, а под конец зарядки - уже серьезно меньше 1кВт. На аналогичной схеме с трансформатором, правда с меньшей средней мощностью, у меня получалось, что в момент начала зарядки конденсатора система потребляла 200Вт, а к концу зарядки - всего-то 10Вт, средняя мощность потребления составляла 60Вт, а на выходе средняя мощность была около 20Вт. Сейчас дополнительно на источник наложено ограничение, что это бензиновый или дизельный генератор, с пиковой 5кВт мощностью, то есть потребление надо максимально сгладить. Понятно, что на том стенде, что на фотографии я больше 100Вт нагрузки цеплять побоюсь :)

Пусть у Вас есть матрица сдвига P (0 1 0) (0 0 1) (1 0 0) a и b - векторы, тогда матрица H= ( a, b ) ^T P ( a, b ) на диагонали будет содержать a0 и a2, а сумма внедиагональных членов будет равна a1. Сам такое много раз и в уравнениях Максвелла встречал, и в сигнал процессинге. Расскажите побольше про задачу, может удастся разобраться откуда уши ростут.

Спасибо конечно за совет про гугл. Я собственно сегодня пол дня гуглил, но, то и я не понимаю то, что в даташитах написано, то ли вот так на раз эти значения не вытаскиваются. Что смог нагуглить, так это фразу, что у большинства ферритов В достигается примерно при 300, но про µe надо гуглить совсем на этот феррит с его геометрией, а не огульно на 3c90. Вот я и вопрошаю, вдруг кто просто и понятно объяснит. Спасибо!

саму схему (симметриксом), или феррит на предмет запасания энергии на разных режимах (Максвеллом от Ансиса) или лучше и то и другое? То же думал над этим... похоже, только одно и осталось.

Уважаемые stells и gte, огромное Вам спасибо, что помогаете советами! То, что на фото - пока времянка, для того, чтобы убедиться, что вообще так можно, поэтому диоды только на 150мА и 2кВ. То, что хочется сделать на основе этого стенда, систему с 24кВ выходом около 1кВ мощности. Хотелось бы использовать те же конденсаторы (0.12мкФ, 1кВ), но другие диоды (тоже 50нс, но уже на 4А), которые вот-вот должны ко мне приехать. Нагрузка у меня как всегда - конденсатор (50нФ, 30кВ) который вначале заряжается, а потом ключом разряжается (хочется иметь систему до 24кВ, но реально использовать ее не выше 18кВ). Полевик больше 1500В не нашел, тот, что есть 4А (12А в пике), но можно привезди 8А (24А в пике). Как я понимаю, при качании 1кВатт, мой конденсатор в 0.12мкФ будет заряжаться за 60мкС, и для его контроля мне нужен АЦП раз в 10 быстрее... или ставить в начале какую-то очень не децкую керамику на 10мкФ 1кВ, тогда действительно, время зарядки будет растягиваться на милисекунду, и по нему можно будет все регулировать. Спасибо ИИВ

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

Пусть у вас система покоится (не ускоряется), положим, Вам надо получить обратную матрицу поворотов, запишем ее как H=(h_1, h_2, h_3), H \in \R^{3 \times 3}, h_i \in \R^3. Дополнительно потребуем, чтобы оси магнетометра и акселометра полностью совпадали. Положим, что Ваша система такова, что ось X направлена на север, Y - на восток, а Z - вверх. Тогда h_3=-m, h_1=a, где m - результат данных с магнетометра, a - результат данных с акселометра. Помня, что H - унитарная, и зная, что датчики безбожно врут, Вам надо будет перенормировать векторы h_3 и h_1, а, также возможно доортогонализовать их друг к другу, если есть ошибка измерения в датчиках (а она точно есть!!!) относительно положения осей магнетометра и акселометра. После того, как вы получите ортонормированные векторы h_3 и h_1, вычисление h_2 можно провести в виде векторного произведения h_3 и h_1. Сама же матрица поворота может быть вычислена из обратной H ее транспонированием. К сожалению, все это красиво в теории, но на практике, Вы столкнетесь с большими ошибками измерений, перенастройкой систем координат между датчиками... Для устойчивого решения этой задачи я обычно использую несколько датчиков, по-разному развернутых друг относительно друга, но там уже совсем другая математика :(

Всем привет, имеется кольцевой феррит, например, такого типа http://ru.farnell.com/jsp/search/productde...jsp?SKU=1784180 намотал на него несколько витков, пусть N, тогда его индуктивность будет составлять (Inductance Factor Al:2060nH) 2.06*N*N мкГн Теперь хочу посчитать его максимальный ток насыщения, например, исходя из формулы B = 1000 * µ0 * µe * I * N / le µ0=1.257*10-3 [мкГн/мм] le=89.7mm (Effective Magnetic Path Length:89.7mm) а где достать для 3C90 допустимые значения B и µe? Понятно, можно вытащить B/Н зависимость, вогнать геометрию феррита и провода в симулятор типа максвелла, долго и нудно считать, и получить то, что надо, но, может, есть какие-то табличные значения для ферритов из которых все это можно на вскидку получить? Посоветуйте, пожалуйста, кто знает! Спасибо ИИВ

Спасибо, Вам stells, я предполагал такую засаду... Интересно, если у меня нет вообще нагрузки, возможно можно выбрать такую высокую частоту, что в дроссель не будет заходить достаточно энергии, чтобы по получаемый с него бесконечно большой выброс не сглаживался бы первым конденсатором умножителя? Попробовал поднять частоту до 2МГц, (дальше побоялся, мог бы до 8МГц) но ничего хорошего не получил. Похоже мне потребуется очень быстрый АЦП, чтобы оценивать на сколько успело все зарядиться, или, при старте, включаться только когда входное напряжение после выпрямления меньше 100В, смотреть на поведение нагрузки, и, только если нагрузка потребляет хорошо, включаться дальнше на полную...

не, с нагрузкой все становится в порядке, но мой вопрос заключается в том, что если без нагрузки у меня на первой ступени получается 200В вместо 70В, которые я максимум хотел бы там видеть, то как только я воткну эту схему в 220VAC, то на тот же мосфет и первую ступень умножителя попадет 2-3кВ, и оно все сгорит, а мне этого бы не хотелось бы. Что я здесь не понимаю, посоветуйте, пожалуйста!

на сколько я понимаю, у меня хорошо умножитель работает только из-за очень больших номиналов конденсаторов (0.12мкФ на 1кВ, это 60мДж энергии в 2220 корпусе!!!) и реально шустных диодов (50нс). Причем все компоненты получились какими-то очень бюджетными - конденсаторы по 23цента, диоды по 10центов, правда цены там + доставка и разтаможка, но все равно, одна палка умножителя в 20 евро обошлась :) А вот с бустом я пока недоразобрался о чем и вопрошаю, кто знает, помогите, пожалуйста!

Чуток поэкспериментировал, смотрите картинку. Собрал управление на МК, пока без сильных и заумных способов управления по нагрузке. Ключевые номиналы, входное питание 35В, мосфет на 1500В от СТМа, дроссель на торе TX36/23/10-3C90 c 29 витками 1мм провода (1.7мГн), умножитель из 20 ступеней на 0.12мкФ 1кВ керамике и каких-то довольно шустрых диодах (нашлись в закромах Родины) со скоростью переключения около 50нс и пиковым током 150мА. На выходе пока нагрузки не было, если не считать делитель на 2ГОма. Что получилось: 1. можно работать в большом диапазоне частот от 1 до 1200кГц даже от МК и, при необходимости, измерять им напряжения в критических участках цепи, 2. на выходе с 20ой ступени получались напряжения от 4кВ (200-800кГц 50% ШИМ) до 10кВ (5кГц 50% ШИМ), ШИМы с меньшим или большим значениями только слегка меняли значения пикового напряжения, 3. при измерении напряжения на первой ступени и на 20ой, напряжение отличалось в 20 раз (в пределах измерения), что говорит о том, что умножитель на 20 ступеней может функционировать качественно и с минимальными потерями :) Что не понятно: без нагрузки не смог получить на выходе меньше 4кВ, при напряжении питания 35В как только не игрался с частотой и ШИМом... Это эквивалентно тому, что на первый конденсатор умножителя идет пульсирующее напряжение [0, 200]В. Если я на вход вместо 35В подам 320В, то получаемое пульсирующее непряжение в 1800В похерит все мои номиналы... Вдруг кто может что-то посоветовать, скажите, пожалуйста, в чем может быть причина таких высоких пульсаций, и как сделать так, чтобы пульсации были меньше (хочется при напряжении питания в 320В или 220VAC иметь пульсации около 900В). PS: шим на выходе с МК осциллографом проверял - все в порядке, левых всплесков или выбросов - нет, после МК стоит мосфет-драйвер TC4421 (два последовательно, так как купил, по глупости, инвертирующий). Спасибо всем откликавшимся, и заранее спасибо тем, кто откликнется! С уважением ИИВ

Спасибо большое за быстрый и исчерпывающий ответ, в Atmega328p стоит тот же конденсатор, похоже, что при потере не более 2 Ватта и скорости измерения 1кГц, я смогу измерять то, что хочу. При 70кВ и 2 Ваттах рассеяния мост будет с 2.5гОм верхним резистором и 125кОм нижним, то есть примерно при 1кГц оцифровке все будет ок. Конечно это на пределе разумного, но... что делать. Спасибо большое за простое и быстрое решение!!! С уважением ИИВ

Всем привет, собственно САБЖ, пожалуйста подскажите. Хочу воткнуть делитель в ВВ сеть и мерить напряжение Атмегой. Так как в сети имеется несколько десятков киловольт, то делитель приходится конструировать на номиналах 10ГОм-500кОм, что приводит к токам в делителе около 7мкА. Попробовал, похоже измеряемое напряжение занижается. Читал даташит, не нашел... Вдруг кто знает, посоветуйте, пожалуйста САБЖ. PS: 10ГОм собрал на последовательных 1ГОмах каждый на 10кВ расчитанный, все в масле, то есть глюков с пробоем быть не должно. Спасибо ИИВ

я боролся двумя методами на старой (50Ватт) установке. Первый - ставил очень большие номиналы на первую ступень каскада и использовал симметричный каскад, это в котором 3 конденсатора и 4 диода на ступень. Но вот в этой системе, где уже удалось получить 400Ватт, это дело не помогло, все успешно вылетело. Придумал простое и топорное решение - начальные несколько шагов ШИМа делать так, что как-будто мне на выходе половину напряжения хочется получить, а дальше, все как и раньше. Единственно, похоже что все эти численные извращения постепенно перерастают из атмеги, пробовал арм9 на управлении, но борда с ним очень неставильная, похоже только плиска и остается... Вы, как всегда, правы! Сильно задумался, надеюсь удастся соптимизировать