首先,MATLAB很容易求出导数
syms x y
y=x^2;
差异(y ) )。
执行结果:
ans=
2*x
求出一点导数也很简单:
subs(diff ),1 ) ) ) ) ) ) ) ) ) ) )。
执行结果:
ans=
2
今天,求出正态分布反函数的导数,如上求出,发现提示错误。 因为正态分布反函数的导数只有[ 0,1 ]有值。 而且,没办法,只能用最愚蠢的定义法。 例如,以下求出0.5下的导数值。
y1=norminv(0.50000000001 )-norminv ) 0.5 );
但是,有一个奇怪的地方。 区间越小,求导越应该准确,但区间太小,反而会让人不满。
但是,最后的验证是相等的,是快乐的
y2=1/(norm pdf (norm inv (0.5,0,1 ),0,1 ) );
DISP(y1;
DISP(y2;