实验一正弦波加高斯白噪声并通过窄带系统

实验目的
1、通过实验深入理解窄带系统；
2、掌握信号加高斯白噪声的方法；
3、练习根据理论分析自行设计实验方法的能力。
实验内容
1、设定正弦信号的频率为10HZ，抽样频率为100HZ；
2、设定高斯白噪声参数为N(0,0.25)，及噪声功率为0.25W；
3、将噪声叠加到正弦信号上，观察其三者时域波形和频谱；
4、信号通过中心频率为ω0的窄带系统，观察通过后的时域波形和频谱。

%加性高斯白噪声信道-用awgn(x,snr)实现
fs=100; %采样频率
T_N=1;%总时间
t=1/fs:1/fs:T_N;%产生时间矢量t
L=T_N*fs;%采样数量
x=sin(20*pi*t);%生成正弦信号x，频率为10Hz
snr=0.25;  %设定加性高斯白噪声的功率snr

zs=wgn(L,1,snr);%生成高斯白噪声
y=awgn(x,snr);  %在正弦信号上叠加高斯白噪声 

figure(1)
subplot(321);%定义第一个图像
plot(t,x);%画出原始信号
title('正弦信号x（时域）');%对图像主题的说明
xlabel('时间/s');%对横轴进行说明
ylabel('幅度/V');%对纵轴进行说明
subplot(322);%定义第一个图像
fft_x=fft(x);%快速傅里叶变换之后的噪声
P_x = abs(fft_x/L);%取幅频特性，除以L
P_x = P_x(1:L/2+1);%截取前半段
P_x(2:end-1)=2*P_x(2:end-1);%单侧频谱非直流分量记得乘以2
f = fs*(0:(L/2))/L;%频率，最多到一半（奈奎斯特采样定理）
plot(f,P_x);%画出原始信号的频域
title('正弦信号x（频域）');%对图像主题的说明
xlabel('频率/Hz');%对横轴进行说明
ylabel('幅度/V');%对纵轴进行说明

subplot(323);%定义第二个图像
plot(t,zs);%画出高斯白噪声信号
title('高斯白噪声zs（时域）');%对图像主题的说明
xlabel('时间/s');%对横轴进行说明
ylabel('幅度/V');%对纵轴进行说明               
subplot(324);%定义第二个图像
fft_zs=fft(zs);%快速傅里叶变换之后的噪声
P_zs = abs(fft_zs/L);%取幅频特性，除以L
P_zs = P_zs(1:L/2+1);%截取前半段
P_zs(2:end-1)=2*P_zs(2:end-1);%单侧频谱非直流分量记得乘以2
f = fs*(0:(L/2))/L;%频率，最多到一半（奈奎斯特采样定理）
plot(f,P_zs);%画出原始信号的频域
title('高斯白噪声zs（频域）');%对图像主题的说明
xlabel('频率/Hz');%对横轴进行说明
ylabel('幅度/V');%对纵轴进行说明 

subplot(325);%定义第三个图像 
plot(t,y);%画出叠加后的信号
title('叠加后的信号y（时域）');%对图像主题的说明
xlabel('时间/s');%对横轴进行说明
ylabel('幅度/V');%对纵轴进行说明
subplot(326);%定义第三个图像 
fft_y=fft(y);%快速傅里叶变换之后的噪声
P_y = abs(fft_y/L);%取幅频特性，除以L
P_y = P_y(1:L/2+1);%截取前半段
P_y(2:end-1)=2*P_y(2:end-1);%单侧频谱非直流分量记得乘以2
f = fs*(0:(L/2))/L;%频率，最多到一半（奈奎斯特采样定理）
plot(f,P_y);%画出原始信号的频域
title('叠加后的信号y（频域）');%对图像主题的说明
xlabel('频率/Hz');%对横轴进行说明
ylabel('幅度/V');%对纵轴进行说明

figure(2)
[b,a]=butter(8,[8/(fs/2),12/(fs/2) ]);
lvbo_y=filter(b,a,y);
subplot(211)
plot((lvbo_y))
xlabel("时间/s")
ylabel("幅度/v")
title("窄带高斯噪声（时域）")
fft_lvbo_y=fft(lvbo_y);%傅里叶变换
P_yf = abs(fft_lvbo_y/L);%取幅频特性，除以L
P_yf = P_yf(1:L/2+1);%截取前半段
P_yf(2:end-1)=2*P_yf(2:end-1);%单侧频谱非直流分量记得乘以2
subplot(212)
plot(f,P_yf)
xlabel("频率/Hz")
ylabel("幅度/v")
title("窄带高斯噪声（频域）")