matlabを使ってヒステリシスをプロットする方法は？

Question

私は、Matlabを使って指数モードでその式を使ってmemristorをモデル化したいと思います。私の目的は、Matlabでヒステリシスプロット（i-vプロット）を取得することです。式はf(x)は窓関数であり、h(v)が多項式関数である、Aが一定である

よう

i=x^n*b*sinh(av)+m(exp(gv)-1) 
x'=Af(x)h(v) 

です。

f(x)=1-(2x-1)^2 
h(v)=cv+dv^3; 

cとdはc<0 & d>0 a=2;b=0.01;g=4;n=4;A=25;m=0.001;

x一定であるMathWorks社のMATLABでプロットを取得することが可能device..isの内部状態ですか？ http://webee.technion.ac.il/people/skva/Memristor%20Models/MATLAB/memristor.m

あなたは（非線形イオンドリフトモデルを使いたいんモデルの種類が選択できます。私はそれがエラー

Answer 1

私は非常に簡単に検索した後、Googleで見つけたこのスクリプトを示しているMATLABで試してみました あなたのケースでは）、ウィンドウ関数のタイプ（あなたはJogelkarの私が想定している）と非線形の電圧 - 電流関係。

少し自分のスクリプトを変更した後、あなたがこのような何かを取得する必要があります：

b=0.01; g=4; n=4; a=2;A=25; m=0.001; 
c=-1;d=4;%your constants 
numOfPoints = 10000; 
t = linspace(-1, 1,numOfPoints); 
dt = t(2) - t(1); 
volt = .003*sin(2*pi*t); 
x = zeros(size(t)); 
curr = x; 
for i=2:numOfPoints 
    x(i) = x(i-1) + A*(c*volt(i)+d*(volt(i))^3)*(1-(2*x(i-1)-1)^2)*dt; 
    curr(i)=x(i)^n*b*sinh(a*volt(i))+m*(exp(g*volt(i))-1); 
end; 
fig = figure(1); 
plot(volt,curr); 
xlabel('Voltage');ylabel('Current'); 
figure(fig); 

幸運を！