Matlab的相关练习题

Question

 

Answer 1

x = linspace(0,pi,200);
fx = sqrt(x.^2+1);

%% 1
gx = diff(fx)./diff(x);  % 前向差分求导数
figure
plot(x(2:end),gx)

%% 2
[~,ind] = min( abs(x-pi/6) );
gx1 = gx(ind)           

%% 3
fp = polyfit(x,fx,4);   % 多项式拟合
gp = polyder(fp);       % 多项式求导
gx2 = polyval(gp,pi/6)

结果：

gx1 =

    0.4675


gx2 =

    0.4599