MATLAB —— 这篇介绍的是信号处理工具箱之波形产生函数 pulstran

mytomorrow

4 I4 G) c7 ?& u/ K; }/ I3 G前两篇博文和这篇博文有些许联系：
- z3 G: d! C' k1 c" a- o! r
MATLAB —— 认识一下信号处理工具箱之波形产生函数 rectpuls
) U) f! J" x. L4 @) V! v
) x1 m; l6 t" G9 k1 ?MATLAB —— 再介绍一下信号处理工具箱之波形产生函数 tripuls

' A* j/ e" t. P5 h+ _" wMATLAB帮助文档称pulstran函数为：脉冲串（Pulse train）产生函数。 pulstran从连续函数或采样原型脉冲生成脉冲序列。

" U" g1 w1 R0 \: B3 A* f: m& A语法格式：

3 m: }9 P1 y8 R
. m9 e) g. |4 X
7 X! G* q6 j2 \语法格式1：
' _( R5 _7 b# z7 ^8 l
0 D$ R: N7 o" \5 Dy = pulstran（t，d，'func'）基于连续函数'func'的样本生成脉冲序列，其中'func'是
* k+ n8 t5 h, m4 n* {+ z1 r
) r3 l* E6 x8 A- P. U'gauspuls'，用于产生高斯调制的正弦脉冲
+ h1 r& Y8 h8 f) B5 A'rectpuls'，用于生成采样的非周期矩形
'tripuls'，用于生成采样的非周期三角形
pulstran被评估长度（d）次并且返回评估的总和y = func（t-d（1））+ func（t-d（2））+ ...

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

; B+ u, \5 {& I  K+ n; J6 P通过将d指定为两列矩阵，可以将可选的增益因子应用于每个延迟的评估，其中第1列中定义的偏移和第2列中的相关增益。 请注意，行向量将被解释为仅指定延迟。

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

可以设定‘func’的相关参数，见语法格式2.
' R6 k- Y) w7 }/ h5 X3 b1 C. ?
语法格式2：
& ?! w6 o4 H+ W/ ]$ U9 a: ]4 p# V9 p, b
pulstran（t，d，'func'，p1，p2，...）允许根据需要将其他参数传递给'func'。 例如：
7 ~1 e* b% P0 u+ P6 x9 W% P* ^* m1 V
$ B( ?9 _* v! z; ^  W& lfunc（t-d（1），p1，p2，...）+ func（t-d（2），p1，p2，...）+ ...
+ T& s9 _' T$ w, q/ m
( I" i% D: a  h1 R( R针对语法格式2举个例子：
* z" s6 M% w6 U* E  {/ L
6 d7 t+ U- M- s

• % 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 all
• t = 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
  

   
" w8 l* a: T% u( Z8 s5 X1 A) q' H. z
* A) \9 h* e9 ^+ r! S. D

2 a' H6 q! ?( q; dy = pulstran(t,d,'tripuls',0.1,-1);

这个函数中的0.1以及-1都是‘tripuls’的参数，例子0.1是脉宽，-1表示脉冲最高幅值在最左边，如果将-1去掉，则默认为0.5.图像如下：
& F( y3 k0 ]0 w3 h4 c& [7 G
+ |# m4 z$ h9 q; {
+ v5 k% @2 z: q8 R" f
可见，这是一个正三角锯齿波脉冲串。

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

再举一个例子：
' E6 M2 Q' t) Y# @
要求都在程序开头的注释里，很简单，懒着翻译了。

• % 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 all
• t = 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'
  + F& Z) H) y3 |. P7 z

   

4 J5 U4 m) C/ N9 Z; D
/ U2 J! n8 c1 b2 }; g, {  f
语法格式3：

4 `9 q! n! f: @pulstran（t，d，p，fs）生成一个脉冲序列，它是矢量p中原型脉冲的多个延迟插值之和，以fs速率采样，其中p跨越时间间隔[0，（length（p） -1）/ fs]，其样本在此间隔之外相同。 默认情况下，线性插值用于生成延迟。
) }1 p3 j- d& O# p7 b7 _
4 V) J% f5 X8 q( y* }7 J6 z- D1 J不太明白，先放这里。

9 m6 J, l) r( N: C9 t6 u8 lpulstran（t，d，p）假设采样率fs等于1 Hz。
) C! g9 M6 c7 N6 p; o) u
也就是说如果省略了fs，则默认为1Hz。
% S! \2 q8 u9 v, ~; ^
pulstran（...，'func'）指定替代插值方法。

就这样吧！

* N3 o6 ~- Y' U# x
/ s5 u$ `7 v) c# a
9 d# t4 W4 |8 d$ R! U+ x

xixihahaheihei

谢谢分享