数学建模MATLAB 台大课程10 数值微积分

目录

多项式微积分

用到的一些函数

多项式的值polyval（）

多项式微分polyder()

多项式积分polyint()

多项式相乘conv（）

练习1

其他公式的微分

数值微分法

diff()

微分法

练习2

练习3

高阶微分

数值积分法

面积积分法

 矩形

梯形

曲边梯形

使用函数积分计算

获取函数handle

定积分integral()

---

视频：

MATLAB教學 - 10数值微积分_哔哩哔哩_bilibili

多项式微积分

用到的一些函数

多项式的值polyval（）

y = polyval(p,x) 计算多项式 p 在 x 的每个点处的值。参数 p 是长度为 n+1 的向量，其元素是 n 次多项式的系数（降幂排序）。

多项式微分polyder()

k = polyder(p) 返回 p 中的系数表示的多项式的导数；
k = polyder(a,b) 返回多项式 a 和 b 的乘积的导数；
[q,d] = polyder(a,b) 返回多项式 a 和 b 的商的导数(q为分子，d为分母)。

多项式积分polyint()

q = polyint(p,k) 返回对多项式积分后的多项式系数，常数项为k；
q = polyint(p) 常数项 k = 0。

多项式相乘conv（）

conv（）函数

练习1

画f(x)=(5x^3-7x^2+5x+10)(4x^2+12x-3)  -2<=x<=1 和他的导数图 

%f(x)=(5x^3-7x^2+5x+10)(4x^2+12x-3)  -2<=x<=1 
clc
g = colormap(lines);hold on
a = [5,-7,5,10];
b = [4,12,-3];
x = linspace(-2,1);
y = conv(a,b);
f = polyval(y,x);
h = polyval(polyder(y),x);
plot(x,f,'--','LineWidth',2);
plot(x,h,'LineWidth',2);box on;
xlabel('x');%ylabel('f(x)');
set(gca,'FontSize',14)
set(gca,'YLim',[-800,800])
H = legend('f(x)','f\prime(x)');