语音的调制与解调

clear;
[sfs]=wavread （‘C:\Users\Administrator\Desktop\1.wav‘）;   %打开已录好的语音文件
wavplay（sfs）;          %回放语音
fs=44100; 
duration=10;
t=0:duration*fs-1;    
figure（1）;     
subplot（211）;
plot（ts）;title（‘原始语音信号时域波形‘）;
xlabel（‘t‘）;ylabel（‘幅度‘）;
grid; 
S=fft（s）;     %快速傅里叶变换
subplot（212）;stem（t（abs（S））‘.‘）;title（‘原始语音信号频域波形‘）; 
xlabel（‘f‘）;ylabel（‘幅度‘）;
grid;
fc=4000;
Uam=s‘.*cos（2*pi*fc*t/fs）;  %进行DSB调制
C=abs（fft（Uam））;
figure（2）; 
subplot（211）;
plot（tUam）;title（‘已调信号时域波形‘）;
grid; 
subplot（212）;
stem（tC‘.‘）;title（‘已调信号频域波形‘）; grid;

ns1=awgn（Uam20）;       %给已调信号添加噪声SNR=20；
ns2=awgn（Uam40）;       %给已调信号添加噪声SNR=40；
ns3=awgn（Uam80）;       %给已调信号添加噪声SNR=80；
NS1=abs（fft（ns1））;
NS2=abs（fft（ns2））;
NS3=abs（fft（ns3））;

figure（3）; 
subplot（311）;
plot（ns1）;title（‘添加噪声SNR=20后信号时域波形‘）;
subplot（312）;
plot（ns2）;title（‘SNR=40‘）;
grid; 
subplot（313）;
plot（ns3）;title（‘SNR=80‘）;

figure（4）
subplot（311）;
stem（tNS1‘.‘）;title（‘添加噪声SNR=20后信号频域波形‘）; grid;
subplot（312）;
stem（tNS2‘.‘）;title（‘SNR=40‘）; grid;
subplot（313）;
stem（tNS3‘.‘）;title（‘SNR=80‘）; grid;

Rp=0.4;Rs=45;               %设计带通滤波器                            
Wp1=[40008000]/（fs/2）;                       
Ws1=[38008500]/（fs/2）;
[n1wc1]=cheb2ord（Wp1Ws1RpRs）;
[b1a1]=cheby2（n1Rswc1）;
[h1w1]=freqz（b1a1）;
H1=20*log10（abs（h1））;

y1=filter（b1a1ns1）;
y01=y1.*cos（2*pi*fc*t/fs）;