通信原理习题5.3

最新推荐文章于 2022-11-03 19:58:14 发布

chlorine27

最新推荐文章于 2022-11-03 19:58:14 发布

阅读量978

点赞数

文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/chlorine27/article/details/121067048

版权

解：（1）

a）由图得，

$s_{1}(t)=Rect(\frac{t}{T_{b}})$

对其做傅里叶变换得，

$S_{1}(f)=T_{b}Sa(\pi f T_{b})$

又s2(t)=-s1(t)，s1(t)表示0，s2(t)表示1，且0和1等概均为1/2，令该二进制序列为 $\left \{ a_{n} \right \}$ 即

$\left\{\begin{matrix} -1 & ,a_{n} \ equals\ to\ '1'\\ 1 & ,a_{n} \ equals\ to\ '0' \end{matrix}\right.$

计算均值 $m_{a}$ 与方差 $\sigma _{a^2}$ ，

$m_{a}=1\times 0.5+(-1)\times 0.5=0$

$\sigma_{a^2}=E[a_{n}^2]-m_{a}^2=1\times 0.5+1\times 0.5-0^2=1$

由功率谱密度公式，

$P_{s}(f)=\frac{\sigma_{a}^2}{T_{b}}\left |S_{1}(f) \right |^2+\frac{m_{a}^2}{T_{b}^2}\sum_{k}\left | S_{1}(\frac{k}{T_{b}}) \right|^2\delta (f-\frac{k}{T_{b}})=\frac{1}{T_{b}}\left |S_{1}(f) \right |^2=T_{b}Sa^2(\pi fT_{b})$

b）

①模型仿真图如下，

②用matlab进行仿真，得到基带波形图，双边功率谱密度图形如下。

波形分析：取 $T_{b}=1$ ，则得到带宽为0.5的门函数 $s_{1}(t)$ ，其幅值为1。令 $s_{2}(t)=-s_{1}(t)$ ，得到二进制符号“1”和“0”的两种表示方式，t从0开始，到截止时间T_stop为止共发送了T_stop/Tb个符号数，得到上图基带波形。再观察双边功率谱密度波形图Ps(f)，在f=0处取得最大值 $P_{s}(0)=T_{b}Sa^2(0)=1$ ，第一零点处f=1/Tb=1以此类推第n处f=n/Tb=n。

③matlab代码及注释

clear,close all
%---------------------------
%系统参数设置
%---------------------------
T_start=-5;%仿真开始时间
T_stop=10;%仿真结束时间
T=T_stop-T_start;%仿真持续时间
T_sample=1/1000;%时间间隔
f_sample=1/T_sample;%采样速率
N_sample=T/T_sample;%点数
T_b=1;%设比特时间间隔为1
M=T_stop/T_b;%发送符号数

%---------------------------
%变量t和f
%---------------------------
t=T_start:T_sample:T_stop-T_sample;%从-5开始以0.001为时间间隔采样10000个点，用于后续
                                   %近似模拟s1(t)
f=(-N_sample/2:N_sample/2-1)*f_sample/N_sample;%频率间隔=采样速率/采样点数=0.1，以
                                               %0.1为间隔，取10000个点，即长度为1000，
                                               %又最后为双边功率谱密度，关于y轴对称，
                                               %故f最终取(-500,500)区间
                                               
%---------------------------
%以下为第一题
%产生波形及傅里叶变化
%---------------------------
a_1=randsrc(1,M,[[0 1];[0.5 0.5]]);%生成M个随机数0，1（0与1等概）
s_t_1_1=rectpuls(t,T_b);%s1(t)
%以下for循环为生成基带信号s(t)
s_t_1=0;
for i=1:M%叠加M次g_T(t-i*Tb),求得s(t)
    s_t_flag_1=rectpuls(t-0.5-i+1,T_b);%矩形脉冲g_T(t-i*Tb)
    if a_1(i)==0
        s_t_1=s_t_1+s_t_flag_1;%an为0时对应s1(t)
    else
        s_t_1=s_t_1-s_t_flag_1;%an为1时对应s2(t)
    end
end
s_ff_1=abs(fft(s_t_1_1))*T/N_sample;%|S1(f)|，*T/N_sample用于幅度修正
%下两行对频谱位置进行修正
s_f_1(N_sample/2+1:N_sample)=s_ff_1(1:N_sample/2);
s_f_1(1:N_sample/2)=s_ff_1(N_sample/2+1:N_sample);

%----------------------------
%制图
%----------------------------
figure(1);
subplot(211);%2行1列位置1处制图
plot(t,s_t_1);%以t为横轴，s_t为纵轴制图
title('基带信号s(t)波形');%标题
xlabel('t');%x轴写t
ylabel('s(t)');%y轴写s(t)
axis([0,10,-1.5,1.5]);%修改xy轴呈现范围，便于观察
subplot(212);%2行1列位置2处制图
plot(f,T_b.*abs(s_f_1).*abs(s_f_1));%以f为横轴，Ps(f)=Tb|S1(f)|^2为纵轴制图
title('双边功率谱密度Ps(f)波形');%标题
xlabel('f');%x轴写f
ylabel('Ps(f)');%y轴写Ps(f)
axis([-5,5,0,1]);%修改xy轴呈现范围，便于观察

（2）

a）由图得，

$s_{1}(t)=Rect(\frac{t}{0.5T_{b}})$

对其做傅里叶变换得，

$S_{1}(f)=0.5T_{b}Sa(2\pi f T_{b})$

又s2(t)=0，s1(t)表示0，s2(t)表示1，且0和1等概均为1/2，令该二进制序列为 $\left \{ a_{n} \right \}$ 即

$\left\{\begin{matrix} 0 & ,a_{n} \ equals\ to\ '1'\\ 1 & ,a_{n} \ equals\ to\ '0' \end{matrix}\right.$

计算均值 $m_{a}$ 与方差 $\sigma _{a^2}$ ，

$m_{a}=1\times 0.5+0\times 0.5=0.5$

$\sigma_{a^2}=E[a_{n}^2]-m_{a}^2=1\times 0.5-0.5^2=0.25$

由功率谱密度公式，

$P_{s}(f)=\frac{\sigma_{a}^2}{T_{b}}\left |S_{1}(f) \right |^2+\frac{m_{a}^2}{T_{b}^2}\sum_{k}\left | S_{1}(\frac{k}{T_{b}}) \right|^2\delta (f-\frac{k}{T_{b}})$

$=\frac{T_{b}}{16}Sa^2(\frac{\pi fT_{b}}{2})+\frac{1}{4T_{b}^2}\sum_{m=-\infty }^{\infty}\frac{T_{b}^2}{4}Sa^2(\frac{m}{2})\delta (f-\frac{m}{T_{b}})$

b）

①模型仿真图如下，

②用matlab进行仿真，得到基带波形图，双边功率谱密度图形如下。

波形分析：取 $T_{b}=1$ ，则得到带宽为 $\frac{T_{b}}{4}$ 的门函数 $s_{1}(t)$ ，其幅值为1。令 $s_{2}(t)=0$ ，得到二进制符号“1”和“0”的两种表示方式，t从0开始，到截止时间T_stop为止共发送了T_stop/Tb个符号数，得到上图基带波形。再观察双边功率谱密度波形图Ps(f)，由冲激函数与Sa函数两个部分组成，由于叠加后画图二者会混叠在一起，不利于观察，故将两者分开绘制。Sa函数与冲激函数均在f=0处取得最大值，Sa函数最大值 $P_{sa\_s}(0)=\frac{\sigma_{a}^2T_{b}}{4}Sa^2(0)=\frac{1}{16}$ ，第一零点处f_sa=2/Tb=1以此类推第n处f=2n/Tb=n。冲激函数最大值 $P_{sa\_s}(0)=\frac{m_{a}^2}{4}Sa^2(0)=\frac{1}{16}$ ，当且仅当 $f= \frac {m}{T_{b}}$ ，其中m为整数是有冲激，且|m|为除0以外的偶数时Sa函数值为0，此时冲激也为0。

③matlab代码

%---------------------------
%以下为第二题，代码接上
%产生波形及傅里叶变化
%---------------------------
a_2=randsrc(1,M,[[0 1];[0.5 0.5]]);%生成M个随机数0，1（0与1等概）
s_t_2_1=rectpuls(t,T_b/2);%s1(t)
%以下for循环为生成基带信号s(t)
s_t_2=0;
for i=1:M%叠加M次g_T(t-i*Tb),求得s(t)
    s_t_flag_2=rectpuls(t-0.25-i+1,T_b/2);%矩形脉冲g_T(t-i*Tb)
    if a_2(i)==0
        s_t_2=s_t_2+s_t_flag_2;%an为0时对应s1(t)
    else
        s_t_2=s_t_2;%an为1时对应s2(t)
    end
end
s_ff_2=abs(fft(s_t_2_1))*T/N_sample;%|S1(f)|，*T/N_sample用于幅度修正
%下两行对频谱位置进行修正
s_f_2(N_sample/2+1:N_sample)=s_ff_2(1:N_sample/2);
s_f_2(1:N_sample/2)=s_ff_2(N_sample/2+1:N_sample);
%脉冲响应部分
delta=0;
for j=-M:M
    delta=delta+(sinc(j/2)).*(sinc(j/2)).*sign(dirac(f-j));%sign:f=j/T_b时产生冲激函数
end
P_s_1=0.25.*abs(s_f_2).*abs(s_f_2);
P_s_2=(1/16).*delta;

figure(2);
subplot(211);
plot(t,s_t_2);
title('基带信号s(t)波形');
xlabel('t');
ylabel('s(t)');
axis([0,10,-0.5,1.5]);
subplot(212);
plot(f,P_s_1);
hold on;
plot(f,P_s_2);
title('双边功率谱密度Ps(f)波形');
xlabel('f');
ylabel('Ps(f)');
axis([-10,10,0,1/16]);