1.信号的尺度变换、翻转、时移
这些运算实际上是函数的自变量运算。略略略。
2.离散序列的差分
调用diff函数:y=diff(f);
3.离散序列的求和
调用sum函数:y=sum(f(k1:k2));
4.连续信号的微分
也可以用diff函数近似计算,求 y = x 2 y=x^2 y=x2的导数
>> h=.01;x=0:h:2;
y=diff(x.^2)*1/h;
plot(x(1:length(x)-1),y) %前向差分,向量x和y长度必须相同
结果如下:
5.连续信号的积分
利用quadl函数或integral函数。例如,求三角波 f ( t ) f(t) f(t)的积分
%tripl.m文件
function y=tripl(t)
y=tripuls(t,4,0.5);
end
命令行窗口执行以下语句
%利用quadl函数
t=-3:0.1:3;
for x=1:length(t)
y2(x)=quadl(@tripl,-3,t(x));
end
plot(t,y2)
title('integral of f(t)')
%利用integral函数
t=-3:0.1:3;
for x=1:length(t) %要画图,所以要这么写;
y2(x)=integral(@tripl,-3,t(x));
end
plot(t,y2)
title('integral of f(t)')
结果如下:
ps:单纯计算结果 Int=integral(@tripl,-3,3);
或 Int=quadl(@tripl,-3,3);
那么,quadl函数和integral函数有什么区别呢???
收藏里有!!!
6.总结
(1)了解了构建函数f(t)的几种方法。如:利用内联函数inline();编写.m文件;匿名函数法;
(2)在diff函数中体会到了向量长度必须相同的要求;
(3)求积分的几种方法以及它们之间的区别。