1

私は単純なオード関数を解決するために使用しています。コードを以下に示します。

function window_sine

% Some constant values to be used later in the code.............

epsilon0=8.85*10^-12;
d_mos=0.65*10^-9;
epsilon_mos_min=5*epsilon0;
epsilon_mos_max=20*epsilon0;
d_g=10*10^-9;
epsilon_g=30*epsilon0;
epsilon_mos=5*epsilon0;
vt=0;
e=1.6*10^-19;  
n=[];
i=1;
u=30; % cm^2/v*S
h=1.05*10^-34;  % j*s
cond_d=e^2/h;

%tIME Step...........

step=0.1;
t = 0:step:10;

%Some other constant values..simple..............
c_g=(epsilon_g/d_g);
c_mos=(epsilon_mos/d_mos);
c_t=1/((1/c_g)+(1/c_mos));
% This below equation involves time "t" which is defined above.
vs=(1-(c_t/c_g))*t;

%ode function calling............

p0=0.01;
[t,p] = ode45(@state,t,p0);

%plotting the result..................

figure
plot(t,p)
title('State Variable')

%Now the Ode45 function.................

function dpdt = state( t,p)

nc=6.8*10^12;
smooth=1;   
vg=t;
k=10^-1;

window1=1-((2*p)-1).^(2*smooth);

% How to incorporate the time dependent "n" value here.......... 

dpdt=k*(n-nc)*window1;   

end

end

問題は関数の「n」です。その時間依存とその値は、次のコードを使用して個別に計算できます。

while i<length(t)
n(i)=((c_g*(t(i)-vs(i)-vt))/e)*(10^-4);% this 10^-4 is for cm^2 unit
i=i+1;
end

SHORT では、「n」の値は「t」に比例して増加します。「t」はどこでも同じであり、導関数で使用される時間であることに注意してください。

問題は、この時間依存の "n" 値を、"n" と記述した ode 関数に組み込む方法です。? コードのこの部分、関数ブロックまたはメイン コードのどこに書き込むか。ありがとう

4

1 に答える 1

0

ode-function の最初の引数は時間であることに注意してください。これは、まさに時間依存変数やレオノミック制約のようなものです。したがって、定数を関数に入力し、ode-function 内で時間依存の計算を行うだけです。

例えば：

function dydt = myodefun(t,y,c)
n=1+t*c; % calculate your time-dependent value n
dydt = [y(2)*c; (1-y(1)^2)*y(2)-y(1)];

たとえば、次のように実行します。

c=1; % a constant, or a structure of constants
myodefunn=@(t,y)myodefun(t,y,c)
[tout,yout]=ode45(myodefunn,[0 20],[2 0]);
于 2015-10-02T12:35:35.233 に答える