力と位置 (クーロンの法則) のプロットを作成し、定数 e0 を推定したいと考えています。私は料金、q1 = 1、q2 = 1の値を持っています。私は例えば
position=[0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1];
force=[0.08,0.015,0.013,0.0062,0.0016,0.00519,-0.00159,0.00118,...
0.0061,0.00155,0.00143];
クーロンは F= (1/4*pi*e0) * q1*q2/r^2 です。したがって、次の形式になります。
y=ax^-m 、ここで a= (q1*q2/4*pi*e0)
私がやっている:
p=polyfit(-log10(position),log10(force),1); % I am not sure about '1' and minus
m=p(1);
a=10^(p(2)); % I am not sure about a
xp=0.1:0.1:1.1;
yp=a*xp.^(-m);
plot(position,force,'o',xp,yp)
e0=q1*q2/4*pi*a
e0 の正しい値が見つかりません。何か間違っていますか? m値は-2である必要がありますが、私は取っています:
メートル=
1.6287 - 0.2008i