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）;????????????