function?[]=vsb（）
??
%?殘留邊帶調(diào)制
%?輸入信號(hào)為模擬信號(hào)
%?碼元寬度?0.000999s
%?整個(gè)仿真總計(jì)1000個(gè)采樣點(diǎn)
t0=0.000999;??
ts=0.000001;
t=（0:ts:t0）;??%時(shí)間變量

%噪聲參數(shù)調(diào)節(jié)
snr=100;?%信噪比，以dB來(lái)表示
snr_lin?=10^（snr/10）;?%dB轉(zhuǎn)換
fs=1/ts;
df=0.3;??%頻率分辨率

fc=50000;????%載波頻率
fc1=50000;???%本地提取載波頻率
%載波信號(hào)生成
c=cos（2*pi*fc.*t）;?????

%生成消息信號(hào)（正弦信號(hào)）
fc_message=5000;???????????????%調(diào)制信號(hào)頻率
m=cos（2*pi*fc_message.*t）;?????%調(diào)制信號(hào)為單頻正弦信號(hào)

%DSB載波調(diào)制
%相乘器
dsb=c.*m;

%將雙邊帶信號(hào)由時(shí)域轉(zhuǎn)換到頻域
[DSBdsbdf1]=fftseq（dsbtsdf）;
DSB=DSB/fs;
f=（0:df1:df1*（length（dsb）-1））-fs/2;

%繪制雙邊帶信號(hào)時(shí)域波形和頻譜
%subplot（211）
%plot（tdsb（1:length（t）））
%subplot（212）
%plot（ffftshift（abs（DSB）））
%殘留邊帶調(diào)制濾波器設(shè)計(jì)
f_band=10000;?????????????%確定互補(bǔ)部分的帶寬
f_cutoff1=fc-f_band;??????????%下截止頻率??中心頻率為載波頻率
f_cutoff2=fc+f_band;??????????????%上截止頻率??中心頻率為載波頻率
n_cutoff1=floor（f_cutoff1/df1）;
%調(diào)制濾波器，變量名H
H=zeros（size（f））;????????

%常數(shù)增益部分，變量名H1
H1（1:n_cutoff1）=ones（1n_cutoff1）;
H1（length（f）-n_cutoff1+1:length（f））=ones（1n_cutoff1）;

%頻率正半軸互補(bǔ)對(duì)稱頻率分量
%變量名H2
H2=zeros（size（f））;?
H2=-1/（（5*10^5-f_cutoff1）-（5*10^5-f_cutoff2））.*（f+（5*10^5-f_cutoff2））;
for?i=1:1:floor（f_cutoff1/df1）
????H2（i）=0;
end

for?i=floor（f_cutoff2/df1）:1:length（f）
????H2（i）=0;
end

%頻率負(fù)半軸互補(bǔ)對(duì)稱頻率分量
%變量名H3
H3=zeros（size（f））;?
H3=1/（（5*10^5-f_cutoff1）-（5*10^5-f_cutoff2））.*（f-（5*10^5-f_cutoff2））;
for?i=（length（f）-n_cutoff1+1）:1:length（f）
????H3（i）=0;
end

for?i=1:1:（length（f）-floor（f_cutoff2/df1））
????H3（i）=0;
end

%合成調(diào)制濾波器
H=H1+H2+H3;
clear?H1?H2?H3