Matlab学习笔记(4)——微积分及概率论相关函数讲解
练习一—绘制如下图形
参考实现:
a=[20 -7 5 10];%创建多项式系数矩阵a
b=[0 4 12 -3]; %创建多项式系数矩阵b
c=conv(a,b); %对a、b求卷积,得到多项式系数相乘的矩阵
d=polyder(c); %对多项式求微分
x=-2:0.01:1;
y=polyval(c,x); %多项式计算
f=polyval(d,x); %多项式计算
plot(x,f,'r',x,y,'--','LineWidth',2);
xlabel('x');
ylabel('f(x)');
set(gca,'FontSize',14);
相关知识点:
1.conv():卷积和多项式乘法
w = conv(u,v) ;%返回向量 u 和 v 的卷积。如果 u 和 v 是多项式系数的向量,对其卷
%积与将这两个多项式相乘等效
w = conv(u,v,shape);%返回如 shape 指定的卷积的分段。例如,conv(u,v,'same')
%仅返回与 u 等大小的卷积的中心部分,而 conv(u,v,'valid') 仅返回计算的没有补零边
%缘的卷积部分
2.polyder():多项式微分
k = polyder(p);%返回 p 中的系数表示的多项式的导数
k = polyder(a,b);%返回多项式 a 和 b 的乘积的导数
3.polyval():多项式计算
y = polyval(p,x);%计算多项式 p 在 x 的每个点处的值。参数 p 是长度为 n+1 的向
%量,其元素是 n 次多项式的系数(降幂排序
执行效果:
练习二—做出两点的差
程序实现:
x=[1 2];
y=[5 7];
delta=diff(y)/diff(x);
相关知识点:
1.diff():
Y = diff(X);%计算沿大小不等于 1 的第一个数组维度的 X 相邻元素之间的差分:
%如果 X 是长度为 m 的向量,则 Y = diff(X) 返回长度为 m-1 的向量。Y 的元素是 X %相邻元素之间的差分。
%Y = [X(2)-X(1) X(3)-X(2) ... X(m)-X(m-1)]
%如果 X 是不为空的非向量 p×m 矩阵,则 Y = diff(X) 返回大小为 (p-1)×m 的矩阵,
%其元素是 X 的行之间的差分。
%Y = [X(2,:)-X(1,:); X(3,:)-X(2,:); ... X(p,:)-X(