仿真作业3：噪声通过DSB-SC解调器

  如下图所示，白噪声$n_i(t)$经过带通滤波器后与接收机本地载波相乘，随后进入低通滤波器。具体说明如下：

• $n_i(t)$为AWGN且单边功率谱密度为$n_0=10^{-6}$W/Hz；
• $h_1(t)$为带通滤波器，其主要作用是限制进入接收机的噪声功率，同时应使得有用信号$s(t)$失真尽可能小地通过。设其频域传递函数为$H_1(f)=\mathrm{R}\mathrm{e}\mathrm{c}\mathrm{t}(\frac{f-f_0}{B})+\mathrm{R}\mathrm{e}\mathrm{c}\mathrm{t}(\frac{f+f_0}{B})$，其中$f_0=100$HZ为其中心频率，$B$=20Hz为带宽。
• $n(t)=n_i(t)\ast h(t)$为经过带通滤波器$h_1(t)$之后的白噪声；
• $n_d(t)=n(t)\cos 2\pi f_{c}t$为乘法器输出信号，这里$f_c=100$Hz；
• $h_2(t)$为低通滤波器，$H_2(f)=\mathrm{R}\mathrm{e}\mathrm{c}\mathrm{t}(\frac{f}{B})$，B=20Hz，其输出$n_o(t)=n_d(t)\ast h_2(t)$。

试求：
1）分别写出$n_i(t)$、$n(t)$、$n_d(t)$、$n_o(t)$功率谱密度的表达式，并画出示意图。
2）分别计算$n_i(t)$、$n_i(t)$、$n_d(t)$、$n_o(t)$的平均功率。
3）用MATLAB仿真上述系统，观察记录$n_i(t)$、$n_i(t)$、$n_d(t)$、$n_o(t)$的功率谱密度，根据功率谱密度计算$n_i(t)$、$n_i(t)$、$n_d(t)$、$n_o(t)$各点功率大小。
【注】
（1）白噪声可用下面语句产生

 noise_i=wgn(1,N_sample,power_dB);

这里参数‘power_dB’为功率的分贝值（dBW）。
（2）功率谱密度建议直接法估计，便于观察强度大小的变化

PSD_Noise_i=abs(fft(noise_i)).^2*T_sample/T/f_sample;

（3）计算平均功率

P_noise_i=sum(PSD_Noise_i)/length(PSD_Noise_i)*f_sample;

注意仿真中产生的白噪声，频率范围为[-f_sample/2,f_sample/2]上，功率谱密度为$\frac{n_0}{2}$，因此平均功率为$\frac{n_0}{2}$f_sample。