|
Форум посвящен вопросам релейной защиты и автоматики (РЗА). Обмену опытом и общению релейщиков. |
Вы не вошли. Пожалуйста, войдите или зарегистрируйтесь.
Советы бывалого релейщика → Учимся делать расчёты → Моделирование в РЗА
|
|
Страницы Назад 1 2 3 4 5 6 7 … 15 Далее
Чтобы отправить ответ, вы должны войти или зарегистрироваться
Спасибо, что откликнулись и уделили внимание. Кому то будет интересно. А по поводу этого
Да и MATLAB для этого использовать – как микроскопом гвозди забивать.
Иногда хочется и микроскопом, отложив кувалду в сторонку :-)
Маленький вопрос по Matlab
подскажите, кто владеет данным продуктом, как взять производную функции в символьном виде, а затем поставить в нее численные значения переменных и вычислить?
чтобы было понятно, простой пример:
y=x^2
берем производную:
y'=2x
x=2
y'(x)=2*2=4
PS Если у кого-нибудь есть поделитесь, пожалуйста, чтобы не изобретать: программа для решения систем нелинейных уравнений методом Ньютона-Рафсона в Matlab.
http://rzia.ru/extensions/hcs_image_uploader/uploads/10000/7500/17521/thumb/p16f4vdmb6a2k14ie1hmng011b4u1.JPG
взять производную: d/dx, к ней приписать дифференцируемую функцию: d(x^2)/dx -> 2x. (стрелка из закладки Symbolic - символьные решения)
Задать х и получим ответ.
взять производную: d/dx, к ней приписать фиддеренцируемую функцию: d(x^2)/dx -> 2x. (стрелка из закладки Symbolic - символьные решения)
извините, а Вы случайно не про Mathcad говорите? просто там именно так это делается?
и по форме записи видно это
Или использовать символьное преобразование. Поставить курсок в дифференцируемом уравнеии к диффреенцируемой ыеличине и в Symbilic/Variable/Differenciate.
ЗЫ: Сорри. Спутал с Маткадом. 
:(
| Присоединяйтесь!!! Мы в социальных сетях и на Ютуб. |
Матлаба на работе нет. Вечером дома гляну как делать в Матлабе.
Хотя в Гуугле должно быть.
Спасибо за отклик. В mathcad я умею данные расчеты производить. А вот с Matlab не приходилось пока работать. вот видимо настал момент изучать.
Взять в производную в символьном виде могу:
к примеру:
gg2=sym('gg2');
y=sym('y');
O=sym('O');
y=(O^2)/(O^4-1);
gg2=diff(y);
disp(gg2);
результат: y' = (2*O)/(O^4 - 1) - (4*O^5)/(O^4 - 1)^2
непонятно вот как теперь вместо "символьной" "О" подставить значение и вычислить?
Deleted
Для Вашего случая нужно будет записать (подставляем O = 5):
Благодарю)) это именно то, что нужно))
clc;
gg2=sym('gg2');
y=sym('y');
P=sym('P');
y=(P^2)/(P^4-1);
gg2=diff(y);
disp(gg2);
z=subs(gg2, P, 6);
disp(z);
Результат:
(2*P)/(P^4 - 1) - (4*P^5)/(P^4 - 1)^2
-0.0093
И маленький совет - в программировании использование переменной O не поощряется, очень легко спутать с нулем.
совет принят, подкорректировал
Для экспорта данных в COMTRADE используется процедура write_comtrade.m (прилагается). Результат проверялся на WaveWin и TransView (Omicron).
При написании процедуры write_comtrade.m была выявлена проблема с ранее выложенной процедурой read_comtrade.m (сложности с пониманием символа возврата каретки и новой строки CR/LF). Поэтому также выкладываю исправленную функцию read_comtrade.mНу и главное – собственно осциллограмма 1p_fault.zip
Извините, если не сложно, то можно ли пошагаво описать каким образом использовать данную функцию, начиная с запуска модели *.mdl.
Маленькая модель показывающая как используя статистические наблюдения получить функциональную зависимость. Например, формулу зависимости числа отказов РЗА какого-либо производителя, в зависимости от срока эксплуатации(условно конечно сказано). в общем тут показываю, что такое аппроксимация и с чем ее едят.
к Mathcad надо прикладывать фонт, иначе русский не разобрать 
:) (мой 15.10 не понимает)
Deleted
Несколько сомнений по поводу адекватности математических моделей трансформаторов тока (ТТ), в том числе и в "секретных" моделях...
Мысль следующая.
1. Любая модель использует характеристику намагничивания ТТ полученную на частоте 50Гц, а она может существенно отличаться при отклонении частоты либо при наличии высших гармонических составляющих в первичном токе. А моделируют многие в предположении, что и на других частотах (либо в первичном сигнале 50Гц с другими частотами 100Гц, 150Гц и т.п.) эта характеристика намагничивания такая же...
2. Гистерезис кривой намагничивания (помучавшись, конечно удастся учесть);
3. Остаточное намагничивание (учитывают ли?);
Deleted
Олег, а что думаете о Jiles-Atherton и о пункте 1 в №94?
Как правильно пользоваться файлом CT_langevin.mdl? Необходимо ввести исходные данные, например, для Rloss?
К тому же ошибка в Lookup Table (Output values)
Deleted
Deleted
Ошибка в Lookup Table, скорее всего, из-за версии Simulink. У меня стоит версия 7.4 (MATLAB R2009b). Специально проверил - никаких библиотечных линков в модели нет, должна работать в stand-alone режиме на любом компе.
Открыл в Matlab 7.12 (R2011a), все работает на ура!
Страницы Назад 1 2 3 4 5 6 7 … 15 Далее
Чтобы отправить ответ, вы должны войти или зарегистрироваться
Советы бывалого релейщика → Учимся делать расчёты → Моделирование в РЗА
