matlab求多项式的导函数,Matlab通过多项式拟合求贝塞尔函数的导函数

%产生拟合曲线，并求某点导数% hObject    handle to btn_ployder (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)x=0.01:0.01:14;steps=str2num(get(handles.edit_steps,'string')) ;if length(steps)==1%必须是单个曲线        switch get(handles.popupmenu_pick,'value')            case 1                           y=BESSELJ(steps(1),x);                   case 2                y=BESSELY(steps(1),x);                 end        a=polyfit(x,y,str2num(get(handles.edit_polyStep,'string')));%获取拟合多项式系数        a1=polyder(a);%多项式一阶导数        a2=polyder(a1);%多项式二阶导数        t=polyval(a,x);        hold on;        plot(x,y,'b');%绘制贝塞尔曲线        plot(x,t,'r:');%绘制拟合曲线        hold off;                xValue=str2num(get(handles.edit_xValue,'string')) ;        yValue=polyval(a,xValue);%函数值        dy1Value=polyval(a1,xValue);%一阶导数值        dy2Value=polyval(a2,xValue);%二阶导数值                set(handles.edit_yValue,'string',num2str(yValue));%显示        set(handles.edit_dy1Value,'string',num2str(dy1Value));        set(handles.edit_dy2Value,'string',num2str(dy2Value));end