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
您好,我是有问必答小助手,您的问题已经有小伙伴解答了,您看下是否解决,可以追评进行沟通哦~
如果有您比较满意的答案 / 帮您提供解决思路的答案,可以点击【采纳】按钮,给回答的小伙伴一些鼓励哦~~
ps:问答VIP仅需29元,即可享受5次/月 有问必答服务,了解详情>>>https://vip.csdn.net/askvip?utm_source=1146287632