S-функция, реализующая гистерезис

@Crazy-fill · Регистрация: 29.07.2013

Author24 — интернет-сервис помощи студентам

В общем задача такая. Необходимо сделать блок в Симулинке, моделирующий зависимость типа гистерезиса по методике Джона Чана.

Добавлено через 4 минуты

Matlab M

clc
Bs=0.49; %индукция насыщения
Br=0.17; %остаточная индукция
Hc=21;   %коэрцитивная сила 
Hmax =300 % максимальное значение намагничивающего поля
H = [-Hmax:1:Hmax]; % пределы изменения и шаг напряженности поля
Bup1 = Bs*(H+Hc); % пошаговое вычисление верхней ->
Bup2 = abs(H+Hc)+Hc*((Bs/Br)-1); % ->ветви петли гистерезиса
Bup = Bup1./Bup2; % верхняя ветвь петли гистерезиса 
Bdn1 = Bs*(H-Hc);  % пошаговое вычисление нижней-> 
Bdn2 = abs(H-Hc)+Hc*((Bs/Br)-1); % ->ветви петли гистерезиса
Bdn = Bdn1./Bdn2; % нижняя ветвь петли гистерезиса
Bmed = 0.5*(Bup+Bdn); % вычисление кривой начального намагничивания
Bupm = max(Bup); % максимальное значение индукции на верхней ветви предельного цикла 
Bmedm = max(Bmed); %координата точки поворота на кривой начального намагничивания
Delt = Bupm - Bmedm; % сдвиг верхней ветви 
Bup0 = Bup-Delt; % верхняя ветвь частного цикла 
Bdn0 = Bdn+Delt; % нижняя ветвь частного цикла

Добавлено через 7 минут
В общем задача такая. Необходимо сделать блок в Симулинке, моделирующий зависимость типа гистерезиса по методике Джона Чана. В моем понимании должна быть создана S-функция, на вход которой я задаю напряженность Н, а на выходе получаю индукцию В. На форуме нашел матлабовский код с реализацией гистерезиса по методу чана, однако как теперь по этому примеру сделать программу в s-функции - ума не приложу, так как мозгов не хватает. Может кто сталкивался с такой задачей или подскажете как это самому реализовать?

Добавлено через 4 минуты

Matlab M

clc
Bs=0.49; %индукция насыщения
Br=0.17; %остаточная индукция
Hc=21;   %коэрцитивная сила 
Hmax =300 % максимальное значение намагничивающего поля
H = [-Hmax:1:Hmax]; % пределы изменения и шаг напряженности поля
Bup1 = Bs*(H+Hc); % пошаговое вычисление верхней ->
Bup2 = abs(H+Hc)+Hc*((Bs/Br)-1); % ->ветви петли гистерезиса
Bup = Bup1./Bup2; % верхняя ветвь петли гистерезиса 
Bdn1 = Bs*(H-Hc);  % пошаговое вычисление нижней-> 
Bdn2 = abs(H-Hc)+Hc*((Bs/Br)-1); % ->ветви петли гистерезиса
Bdn = Bdn1./Bdn2; % нижняя ветвь петли гистерезиса
Bmed = 0.5*(Bup+Bdn); % вычисление кривой начального намагничивания
Bupm = max(Bup); % максимальное значение индукции на верхней ветви предельного цикла 
Bmedm = max(Bmed); %координата точки поворота на кривой начального намагничивания
Delt = Bupm - Bmedm; % сдвиг верхней ветви 
Bup0 = Bup-Delt; % верхняя ветвь частного цикла 
Bdn0 = Bdn+Delt; % нижняя ветвь частного цикла

	11.01.2015, 20:45