数值微积分
1.polyval()---建立多项式
>> a = [1,0,0];----------f = x^2,用列表a来表示多项式的系数
>> x = -100:0.01:100;----用x来表示自变量
>> f = polyval(a, x);----创建多项式f(x)
>> plot(x, f);------画图
2.polyder(x)--对多项式求导,x表示多项式的系数
>> polyder(a)------f(x) = x^2,求导后f‘(x)=2x
ans =
2 0
>> polyval(polyder(a),2)-----求出f‘(2)
ans =
4
3.conv()做多项式乘法,deconv()做多项式除法
f(x) = (20x^3 - 7x^2 + 5x + 10)(4x^2 + 12x -3) -2<=x<=1
>> a = [20 -7 5 10];
b = [4 12 -3];
c = conv(a,b);
x = -2:0.1:1;
f1 = polyval(c,x);
f2 = polyval(polyder(c),x);
plot(x,f1,'-b',x,f2,'r');
legend('f(x)','f`(x)')
4.polyint(x,k)--对多项式求积分,x表示多项式的系数,k表示积分后的常数项
>> a = [5 0 -2 0 1];
>> polyint(a,3);
>> polyint(a,3)
ans =
1.0000 0 -0.6667 0 1.0000 3.0000