Question: Hysteresis function and plot

Has someone an idea how to create a hysteresis plot that models a system with hysteresis?

A simple plot can be created by plotting two arctan functions (g1,g2)

But, I'd rather need a continous function as in the second example (f1..3).

 

 

 

restart;

g1:= x -> 10*arctan(x-6);
g2:= x -> 10*arctan(x+6);
Plot_g1:=plot(g1(x),x=-30..+30);
Plot_g2:=plot(g2(x),x=-30..+30);
plots[display](Plot_g1,Plot_g2);

 

 

 

 

 

f1:= x -> 1.5+sin(x);
Plot_f1:=plot(f1(x),x=0..2*Pi,color=green);
f2:= x -> 1.5+sin(x+(0.5*Pi));
Plot_f2:=plot(f2(x),x=0..2*Pi,color=blue);
plots[display](Plot_f1,Plot_f2);

f3:= x -> f1(x)-f2(x);
Plot_f3:=plot([f3(t),f1(t),t=0..2*Pi]);
plots[display](Plot_f3);

 

 

 

 

 

 

 

 

 

 

 

This yields an ellipse but not a hysteresis plot.

Is there some literature available how to modify trigonometric functions to get the desired result?

 

Please Wait...