【 MATLAB 】信号处理工具箱之波形产生函数 pulstran

前两篇博文和这篇博文有些许联系：

【 MATLAB 】信号处理工具箱之波形产生函数 rectpuls

【 MATLAB 】信号处理工具箱之波形产生函数 tripuls

MATLAB帮助文档称pulstran函数为：脉冲串（Pulse train）产生函数。 pulstran从连续函数或采样原型脉冲生成脉冲序列。

语法格式：

语法格式1：

y = pulstran（t，d，'func'）基于连续函数'func'的样本生成脉冲序列，其中'func'是

'gauspuls'，用于产生高斯调制的正弦脉冲
'rectpuls'，用于生成采样的非周期矩形
'tripuls'，用于生成采样的非周期三角形
pulstran被评估长度（d）次并且返回评估的总和y = func（t-d（1））+ func（t-d（2））+ ...

在移除从向量d获取的标量参数偏移之后，在数组t中指定的参数值范围内计算该函数。请注意，func必须是一个矢量化函数，可以将数组t作为参数。

通过将d指定为两列矩阵，可以将可选的增益因子应用于每个延迟的评估，其中第1列中定义的偏移和第2列中的相关增益。请注意，行向量将被解释为仅指定延迟。

使用上面这种格式，'func'的相关参数都是默认的。

可以设定‘func’的相关参数，见语法格式2.

语法格式2：

pulstran（t，d，'func'，p1，p2，...）允许根据需要将其他参数传递给'func'。例如：

func（t-d（1），p1，p2，...）+ func（t-d（2），p1，p2，...）+ ...

针对语法格式2举个例子：

% This example generates an asymmetric sawtooth waveform with a repetition frequency of 3 Hz and a sawtooth width of 0.1 s.
% The signal length is 1 s and the sample rate is 1 kHz.clc
clear
close allt = 0 : 1/1e3 : 1;         % 1 kHz sample freq for 1 s
d = 0 : 1/3 : 1;           % 3 Hz repetition frequency
y = pulstran(t,d,'tripuls',0.1,-1);plot(t,y)
xlabel 'Time (s)', ylabel Waveform

y = pulstran(t,d,'tripuls',0.1,-1);

这个函数中的0.1以及-1都是‘tripuls’的参数，例子0.1是脉宽，-1表示脉冲最高幅值在最左边，如果将-1去掉，则默认为0.5.图像如下：

可见，这是一个正三角锯齿波脉冲串。

如果将‘func’的所有参数都使用默认值的话，那么就成了语法格式1了。

再举一个例子：

要求都在程序开头的注释里，很简单，懒着翻译了。

% This example generates a periodic Gaussian pulse signal at 10 kHz with 50% bandwidth.
% The pulse repetition frequency is 1 kHz, the sample rate is 50 kHz,
% and the pulse train length is 10 ms. Each pulse has 80% of the amplitude of the preceding pulse.clc
clear
close allt = 0 : 1/50e3 : 10e-3;
d = [0 : 1/1e3 : 10e-3 ; 0.8.^(0:10)]';
y = pulstran(t,d,'gauspuls',10e3,0.5); plot(t,y)
xlabel 'Time (s)', ylabel 'Periodic Gaussian pulse'

语法格式3：

pulstran（t，d，p，fs）生成一个脉冲序列，它是矢量p中原型脉冲的多个延迟插值之和，以fs速率采样，其中p跨越时间间隔[0，（length（p） -1）/ fs]，其样本在此间隔之外相同。默认情况下，线性插值用于生成延迟。

不太明白，先放这里。

pulstran（t，d，p）假设采样率fs等于1 Hz。

也就是说如果省略了fs，则默认为1Hz。

pulstran（...，'func'）指定替代插值方法。

就这样吧！