EDA365电子论坛网

标题: MATLAB —— 这篇介绍的是信号处理工具箱之波形产生函数 pulstran [打印本页]

---

作者: mytomorrow    时间: 2019-10-21 10:42
标题: MATLAB —— 这篇介绍的是信号处理工具箱之波形产生函数 pulstran

前两篇博文和这篇博文有些许联系：
# ]! P# [, C2 w. Y9 Q" i% I! k0 F
* Q, X) |' K& l9 z; d* yMATLAB —— 认识一下信号处理工具箱之波形产生函数 rectpuls

MATLAB —— 再介绍一下信号处理工具箱之波形产生函数 tripuls
( k  L+ X! |) X! m! j) B
: @, m0 w0 i+ n1 u2 LMATLAB帮助文档称pulstran函数为：脉冲串（Pulse train）产生函数。 pulstran从连续函数或采样原型脉冲生成脉冲序列。
8 z2 b; N9 L4 D. m" _' E
语法格式：
; j1 H3 \  }) G: c1 a
% B" k7 C0 `1 a- [: ~
$ N; R. v: F/ L5 m6 e
3 W+ O) m1 u, n; ]8 ?5 a1 Z语法格式1：

' ?; |( m/ p" Sy = pulstran（t，d，'func'）基于连续函数'func'的样本生成脉冲序列，其中'func'是

'gauspuls'，用于产生高斯调制的正弦脉冲
& `7 O1 P8 ]9 @# z! T) f'rectpuls'，用于生成采样的非周期矩形
'tripuls'，用于生成采样的非周期三角形
) a: F* w6 M  x6 n9 I# o/ i- tpulstran被评估长度（d）次并且返回评估的总和y = func（t-d（1））+ func（t-d（2））+ ...
! e# V7 o8 g9 h: R
在移除从向量d获取的标量参数偏移之后，在数组t中指定的参数值范围内计算该函数。 请注意，func必须是一个矢量化函数，可以将数组t作为参数。
- D/ |, {! C9 C1 S' j+ d
2 \% ^: M; c, i% |通过将d指定为两列矩阵，可以将可选的增益因子应用于每个延迟的评估，其中第1列中定义的偏移和第2列中的相关增益。 请注意，行向量将被解释为仅指定延迟。
" t+ [3 L+ T* c  M' u* y
使用上面这种格式，'func'的相关参数都是默认的。

2 M! M. H/ G3 F7 I可以设定‘func’的相关参数，见语法格式2.

语法格式2：
( y+ E( I3 v5 o3 K3 e" U2 A
8 y# v% J! h, wpulstran（t，d，'func'，p1，p2，...）允许根据需要将其他参数传递给'func'。 例如：

1 ^, [% ^$ H9 H* O3 V* Nfunc（t-d（1），p1，p2，...）+ func（t-d（2），p1，p2，...）+ ...
% u1 j7 Z, [! J) z0 R
针对语法格式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 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
  

   
8 B# o6 k- O" `/ w

4 j, `" m- N: i
3 N" ^/ L8 x8 V1 G! N; C% }y = pulstran(t,d,'tripuls',0.1,-1);
7 ?6 h4 E2 e& w6 L. W, r
& Z7 s: I3 x  Y2 \* r5 t5 C这个函数中的0.1以及-1都是‘tripuls’的参数，例子0.1是脉宽，-1表示脉冲最高幅值在最左边，如果将-1去掉，则默认为0.5.图像如下：
7 ~6 b: n8 {5 I" \

3 n; K3 i* e8 P& B9 X" L7 q
4 P$ Q* t1 z: B* \: Z可见，这是一个正三角锯齿波脉冲串。
6 M& _) K5 N3 }5 |) _# ~
8 |! a9 g6 s4 p如果将‘func’的所有参数都使用默认值的话，那么就成了语法格式1了。

3 H8 \* P( z# A( m( ~再举一个例子：

0 k. a+ C' U4 ]要求都在程序开头的注释里，很简单，懒着翻译了。

• % 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'
  : j5 K, H; Z: l. d7 C/ C

   

- n2 H+ {4 I1 o' w. x3 [
+ q9 {. o# Y1 u% k$ p6 N
( u+ ^6 r* `9 D3 D4 ~语法格式3：

pulstran（t，d，p，fs）生成一个脉冲序列，它是矢量p中原型脉冲的多个延迟插值之和，以fs速率采样，其中p跨越时间间隔[0，（length（p） -1）/ fs]，其样本在此间隔之外相同。 默认情况下，线性插值用于生成延迟。
$ T: e& [. p! J% F
5 {( K9 K0 y5 @不太明白，先放这里。

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

+ T# T6 r4 g, l4 c$ q! ?* u也就是说如果省略了fs，则默认为1Hz。
) x  z, y  P: T8 g8 t. Y, g- ^. {1 t
pulstran（...，'func'）指定替代插值方法。

就这样吧！


6 p' T/ m4 b: q/ N' {5 H
# t4 {; s6 E) s. _+ ]2 \

---

作者: xixihahaheihei    时间: 2019-10-23 10:11
谢谢分享

---

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)

Powered by Discuz! X3.2