データのベクトル (ydot という名前) と時間のベクトルがあります。ydot を時間に対してプロットすると、周期的な正弦関数のような図が得られます。時間平均 ydot を計算するにはどうすればよいですか?
次のコードは、ydot と time の正確な値を計算します。
T=(2*pi)/(160e6)
tspan=linspace(0,2*T,1500)
current=linspace(0,1e-6,40);
for k=1:length(current)
f = @(y, t) (current(k)/3.2911e-016)-(2.6151e+009)*sin(y)+(4.8448e+008)*sin(y+0.5697)+(5.2266e+008)*sin((160e6)*t)*cos(y);
[t{k}, y{k}] = ode45(f,tspan,2e22);
end
y1=cell2mat(y);
t1=cell2mat(t);
for k=1:length(tspan)
for j=1:length(current)
ydot(k,j)=(current(j)/3.2911e-016)-(2.6151e+009)*sin(y1(k,j))+(4.8448e+008)*sin(y1(k,j)+0.5697)+(5.2266e+008)*sin((160e6)*t1(k,j))*cos(y1(k,j));
end
end
これにより、40 の異なる電流の ydot が得られます。次のコードは、特定の電流 (k) (k=1:40) の ydot/時間をプロットします。
plot(t1(:,k),ydot(:,k))