%此次代碼實現(xiàn)功能

? ? 1.生成一個單頻信號，并繪圖

????2. 產(chǎn)生信噪比從 -10dB:5dB:10dB 的噪聲并且繪圖表現(xiàn)出來

clear all

close all

clc

%% 產(chǎn)生聲波源信號

f0 = 1.e3;? ? ? ? ? % 單頻信號1頻率Hz

f1 = 2.e3;? ? ? ? ? % 單頻信號1頻率Hz

T_pulse = 5.e-2;? ? % 脈沖寬度s

fs = 7 * max(f0,f1);% 采樣率hz

n_number = T_pulse * fs;

t = 0: 1/fs :T_pulse;

s0 = cos(2*pi*f0*t);

s1 = 2*sin(2*pi*f1*t);

figure;plot(t,s0,'-k');hold on; plot(t,s1,'-r');

axis([-0.01 1.05*T_pulse -2.01 2.01]);? %控制坐標(biāo)軸的范圍

xlabel('時間/s');ylabel('幅度/v');grid on;

%% 產(chǎn)生信噪比從 -10dB:5dB:10dB 的噪聲并且繪圖表現(xiàn)出來

for? snr = -10:5:10

? ? y = awgn(s0,snr,'measured');? ? %在添加噪聲之前測量了x的能量

? ? figure;

? ? plot(t,y);xlim([-0.01 1.05*T_pulse]);

? ? xlabel('時間/s');ylabel('幅度/v');grid on;

end