求具有滚降系数 a 的余弦滚降特性 H(f) 相应的冲激响应 h(t)

该博客详细介绍了具有滚降系数α的余弦滚降特性H(f)的冲激响应h(t)的求解过程。通过傅里叶变换的对称性、分段函数的积分计算,结合三角函数的和差化积公式,逐步推导出h(t)的最终表达式为 Sa(πt/Tb)(1−4α^2t^2/Tb^2)cos(παt/Tb),其中Sa(x)是sinc函数。
摘要由CSDN通过智能技术生成

问题的提出

求 具 有 滚 降 系 数 α 的 余 弦 滚 降 特 性 H ( f ) 相 应 的 冲 激 响 应 h ( t ) 。 H ( f ) = { T b   , 0 ≤ ∣ f ∣ ≤ ( ( 1 − α ) / 2 T b ) T b 2 { 1 + cos ⁡ [ π T b α ( ∣ f ∣ − 1 − α 2 T b ) ] }   , 1 − α 2 T b ≤ ∣ f ∣ ≤ 1 + α 2 T b 0   , ∣ f ∣ ≤ ( ( 1 + α ) / 2 T b ) 其 中   ,   0 ≤ α ≤ 1 。 求具有滚降系数 \alpha的余弦滚降特性H( f ) 相应的冲激响应 { h(t)}。\\ \\ \\ H(f)= \begin{cases} T_b\,,0 \leq |f| \leq ((1-\alpha)/2T_b) \\ \large {T_b \over 2}\{1+\cos[{\pi T_b \over \alpha}(|f|-{1-\alpha \over 2T_b})]\} \,,{1-\alpha \over 2T_b} \leq |f| \leq {1+\alpha \over 2T_b} \\ 0 \,, |f| \leq ((1+\alpha)/2T_b) \end{cases} \\ \\ 其中\,,\, 0 \le \alpha \le 1。 αH(f)h(t)H(f)=Tb,0f((1α)/2Tb)2Tb{ 1+cos[απTb(f2Tb1α)]},2Tb1αf2Tb1+α0,f((1+α)/2Tb),0α1

求解分段函数的第1部分对应的h1(t)

利用傅里叶变换的对称性:

{ f ( t )   ⇔ F ( ω ) F ( t )   ⇔ 2 π f ( − ω ) \begin{cases} f(t)\, \Leftrightarrow F(\omega) \\ F(t)\, \Leftrightarrow 2\pi f(-\omega) \end{cases} { f(t)F(ω)F(t)2πf(ω)

故由:

f ( t ) = { E   , ∣ t ∣ ≤ τ 2 0   , 其 他 ⇔ F ( ω ) = 2 E sin ⁡ ω τ 2 ω f(t)=\begin{cases} E\,,|t| \leq {\large \tau \over 2} \\ 0\,,其他 \end{cases} \Leftrightarrow F(\omega)=2E{\sin{\large \omega \tau \over 2} \over \omega} f(t)={ E,t2τ0,F(ω)=2Eωsin2ωτ

推出:(1)

H 1 ( ω ) = { T b   , 0 ≤ ∣ ω ∣ ≤ π ( 1 − α ) T b 0   , 其 他 ⇔ h 1 ( t ) = 2 T b sin ⁡ π ( 1 − α ) t T b 2 π t = T b sin ⁡ π ( 1 − α ) t T b π t ⇒ H 1 ( f ) = { T b   , 0 ≤ ∣ f ∣ ≤ ( 1 − α ) 2 T b 0   , 其 他 ⇔ h 1 ( t ) = 2 T b sin ⁡ π ( 1 − α ) t T b 2 π t = T b sin ⁡ π ( 1 − α ) t T b π t \begin{aligned} H_1(\omega) &= \begin{cases} T_b\,,0\leq |\omega| \leq \large{ \pi (1-\alpha) \over T_b} \\ 0\,,其他 \end{cases} \Leftrightarrow \large h_1(t)=2T_b{ {\sin \large{ \pi (1-\alpha)t} \over T_b } \over 2\pi t} = T_b{ {\sin \large{ \pi (1-\alpha)t \over T_b }} \over \pi t} \\ \Rightarrow H_1(f) &= \begin{cases} T_b\,,0\leq |f| \leq \large{ (1-\alpha) \over 2T_b} \\ 0\,,其他 \end{cases} \Leftrightarrow \large h_1(t)=2T_b{ {\sin \large{ \pi (1-\alpha)t }\over T_b } \over 2\pi t} = T_b{ {\sin \large{ \pi (1-\alpha)t} \over T_b } \over \pi t} \end{aligned} H1(ω)H1(f)={ Tb,0ωTbπ(1α)0,h1(t)=2Tb2πtsinTbπ(1α)t=TbπtsinTbπ(1α)t={ Tb,0f2Tb(1α)0,h1(t)=2Tb2πtsinTbπ(1α)t=TbπtsinTbπ(1α)t

求解分段函数的第2部分对应的h2(t)

又由:

H 2 ( f ) = T b 2 { 1 + cos ⁡ [ π T b α ( ∣ f ∣ − 1 − α 2 T b ) ] }   , 1 − α 2 T b ≤ ∣ f ∣ ≤ 1 + α 2 T b ⇒ H 2 ( ω ) = T b 2 { 1 + cos ⁡ [ T b 2 α ( ∣ ω ∣ − π ( 1 − α ) T b ) ] }   , π ( 1 − α ) T b ≤ ∣ ω ∣ ≤ π ( 1 + α ) T b \large \begin{aligned} H_2(f) &= {T_b \over 2}\{1+\cos[{\pi T_b \over \alpha}(|f|-{1-\alpha \over 2T_b})]\} \,,{1-\alpha \over 2T_b} \leq |f| \leq {1+\alpha \over 2T_b} \\ \Rightarrow H_2(\omega) &= {T_b \over 2}\{1+\cos[{ T_b \over 2\alpha}(|\omega|-{\pi(1-\alpha) \over T_b})]\} \,,{\pi(1-\alpha) \over T_b} \leq |\omega| \leq {\pi(1+\alpha) \over T_b} \end{aligned} H2(f)H2(ω)=2Tb{ 1+cos[απTb(f2Tb1α)]},2Tb1αf2Tb1+α=2Tb{ 1+cos[2αTb(ωTbπ(1α))]},Tbπ(1α)ωTbπ(1+α)

推出:(2)

对 应 的 h ( t ) 如 下 : 2 π h 2 ( t ) = ∫ − ∞ ∞ T b 2 { 1 + cos ⁡ [ T b 2 α ( ∣ ω ∣ − π ( 1 − α ) T b ) ] } e j ω t d ω ⇒ 4 π T b h 2 ( t ) = ∫ − ∞ ∞ { 1 + cos ⁡ [ T b 2 α ( ∣ ω ∣ − π ( 1 − α ) T b ) ] } e j ω t d ω = ∫ − 1 + α T b π − 1 − α T b π { 1 + cos ⁡ [ T b 2 α ( − ω − π ( 1 − α ) T b ) ] } e j ω t d ω   + ∫ 1 − α T b π 1 + α T b π { 1 + cos ⁡ [ T b 2 α ( ω − π ( 1 − α ) T b ) ] } e j ω t d ω = e j ω t j t ∣ − 1 + α T b π − 1 − α T b π   + e j ω t j t ∣ 1 − α T b π 1 + α T b π   + ∫ − 1 + α T b π − 1 − α T b π cos ⁡ [ T b 2 α ( − ω − π ( 1 − α ) T b ) ] e j ω t d ω   + ∫ 1 − α T b π 1 + α T b π cos ⁡ [ T b 2 α ( ω − π ( 1 − α ) T b ) ] e j ω t d ω = 2 [ sin ⁡ π ( 1 + α ) T b t t − sin ⁡ π ( 1 − α ) T b t t ]   + ∫

在MATLAB中,绘制升余弦滚降滤波器的单位冲激响应通常涉及到信号处理中的数字滤波器设计。升余弦滚降滤波器是一种线性相位滤波器,它结合了低通滤波器的滚降特性,用于减小过渡带内的幅度失真。 以下是简单的步骤: 1. **导入必要的库**:首先需要加载MATLAB内置的信号处理工具箱,如果尚未加载,可以使用`load(' Communications Toolbox')`。 2. **设定滤波器参数**:确定滤波器的截止频率、滚降因子等。例如,假设截止频率为0.5π rad/s,滚降因子为0.6。 ```matlab fs = 1; % Sampling frequency (Hz) fc = 0.5 * pi; % Cutoff frequency (rad/s) alpha = 0.6; % Roll-off factor ``` 3. **计算滤波器系数**:使用`fir1`函数设计一个FIR(有限 impulse response)滤波器,它会产生一个近似升余弦滚降特性的滤波器。 ```matlab [num, den] = fir1(4096, alpha*fc/fs, 'low'); ``` 这里,`4096`是滤波器长度(即采样点数),可以根据需要调整。 4. **绘制单位冲激响应**:创建一个单位冲激序列,然后通过滤波器得到响应,最后使用`impz`函数绘制幅频响应曲线,并显示滤波器的单位冲激响应。 ```matlab % Create unit impulse and apply the filter u = impz(num, den); % Unit impulse response % Plot the impulse response and magnitude spectrum figure; stem(u); xlabel('Sample index'); ylabel('Amplitude'); title('Unit Impulse Response'); figure; freqz(num, den, [0 fs/2]); % Frequency response grid on; xlabel('Frequency (Hz)'); ylabel('Magnitude'); title('Filter Magnitude Response'); ``` 完成以上步骤后,你会看到两个图形,一个展示了滤波器的单位冲激响应,另一个是其频率响应。
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值