function?[]=vsb（）
??
%?殘留邊帶調制
%?輸入信號為模擬信號
%?碼元寬度?0.000999s
%?整個仿真總計1000個采樣點
t0=0.000999;??
ts=0.000001;
t=（0:ts:t0）;??%時間變量

%噪聲參數調節(jié)
snr=100;?%信噪比，以dB來表示
snr_lin?=10^（snr/10）;?%dB轉換
fs=1/ts;
df=0.3;??%頻率分辨率

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

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

%DSB載波調制
%相乘器
dsb=c.*m;

%將雙邊帶信號由時域轉換到頻域
[DSBdsbdf1]=fftseq（dsbtsdf）;
DSB=DSB/fs;
f=（0:df1:df1*（length（dsb）-1））-fs/2;

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

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

%頻率正半軸互補對稱頻率分量
%變量名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

%頻率負半軸互補對稱頻率分量
%變量名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

%合成調制濾波器
H=H1+H2+H3;
clear?H1?H2?H3