%%%王丹妹-----電信112班------2011012489%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc;
%%%正弦信號(hào)的產(chǎn)生
fs=30000;%采樣頻率
t=（0:1/fs:0.005）;%圖形坐標(biāo)軸設(shè)計(jì)
f1=1000;f2=2000;f3=3000;%正弦信號(hào)頻率
s1=sin（2*pi*f1*t）;%產(chǎn)生正弦信號(hào)
s2=sin（2*pi*f2*t）;
s3=sin（2*pi*f3*t）;
figure（1）
subplot（321）;plot（ts1）;xlabel（‘單位：s‘）;ylabel（‘幅度‘）;title（‘正弦信號(hào)1‘）;%顯示正弦信號(hào)圖形
subplot（323）;plot（ts2）;xlabel（‘單位：s‘）;ylabel（‘幅度‘）;title（‘正弦信號(hào)2‘）;
subplot（325）;plot（ts3）;xlabel（‘單位：s‘）;ylabel（‘幅度‘）;title（‘正弦信號(hào)3‘）;
%%%頻域分析
N=1024;
Y1=fft（s1N）;Y1=fftshift（Y1）;%快速傅里葉變換得出頻譜函數(shù)
Y2=fft（s2N）;Y2=fftshift（Y2）;
Y3=fft（s3N）;Y3=fftshift（Y3）;???
f=（0:N-1）*fs/N-fs/2;
subplot（322）;plot（fabs（Y1））;xlabel（‘單位：HZ‘）;ylabel（‘幅度‘）;title（‘函數(shù)頻譜圖1‘）;
subplot（324）;plot（fabs（Y2））;xlabel（‘單位：HZ‘）;ylabel（‘幅度‘）;title（‘函數(shù)頻譜圖2‘）;
subplot（326）;plot（fabs（Y3））;xlabel（‘單位：HZ‘）;ylabel（‘幅度‘）;title（‘函數(shù)頻譜圖3‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%SSB調(diào)制
fc1=10000;fc2=20000;fc3=30000;fs1=100000;
sm1?=?modulate（s1fc1fs1‘a(chǎn)mssb‘）;???%對(duì)信號(hào)進(jìn)行SSB調(diào)制
sm2?=?modulate（s2fc2fs1‘a(chǎn)mssb‘）;???%對(duì)信號(hào)進(jìn)行SSB調(diào)制
sfm3?=?modulate（s3fc3fs1‘FM‘）;???%對(duì)信號(hào)進(jìn)行FM調(diào)制
figure（2）
subplot（321）;plot（tsm1）;xlabel（‘單位：s‘）;ylabel（‘幅度‘）;title（‘SSB1‘）;
subplot（323）;plot（tsm2）;xlabel（‘單位：s‘）;ylabel（‘幅度‘）;title（‘SSB2‘）;
subplot（325）;plot（tsfm3）;xlabel（‘單位：s‘）;ylabel（‘幅度‘）;title（‘FM3‘）;
%%%SSB調(diào)制的頻域分析
F1=fft（sm1N）;F1=fftshift（F1）;?
F2=fft（sm2N）;F2=fftshift（F2）;??????
FM3=fft（sfm3N）;FM3=fftshift（FM3）;?
f1=（0:N-1）*fs1/N-fs1/2;
subplot（322）;plot（f1abs（F1））;xlabel（‘單位：HZ‘）;ylabel（‘幅度‘）;title（‘SSB頻譜圖1‘）;
subplot（324）;plot（f1abs（F2））;xlabel（‘單位：HZ‘）;ylabel（‘幅度‘）;title（‘SSB頻譜圖2‘）;
subplot（326）;plot（f1abs（FM3））;xlabel（‘單位：HZ‘）;ylabel（‘幅度‘）;title（‘FM頻譜圖3‘）;
grid;?????
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%將三路信號(hào)疊加合成為一路在信道中傳輸%%%%%%%%%%%%%%%%%%
dj1=sm1+sm2+sfm3;%三路信號(hào)疊加
dj3=awgn（dj150）;%三路SSB+G的疊加							
figure（3）	
subplot（221）;plot（tdj1）;xlabel（‘單位：s‘）;ylabel（‘幅度‘）;title（‘三路信號(hào)疊加‘）;
subplot（223）;plot（tdj3）;xlabel（‘單位：s‘）;ylabel（‘幅度‘）;title（‘三路信號(hào)+G疊加‘）;
%%%三路信號(hào)疊加合成為一路信號(hào)的頻譜分析
DJ1=fft（dj1N）;DJ1=fftshift（DJ1）;??
DJ3=fft（dj3N）;DJ3=fftshift（DJ3）;??
subplot（222）;plot（f1abs（DJ1））;xlabel（‘單位：HZ‘）;ylabel（‘幅度‘）;title（‘三路信號(hào)疊加頻譜圖‘）;
subplot（224）;plot（f1abs（DJ3））;xlabel（‘單位：HZ‘）;ylabel（‘幅度‘）;title（‘三路信號(hào)+G疊加頻譜圖‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%帶通濾波器的設(shè)計(jì)
fs1=100000;
Rp=0.5;Rs=40;??%用切比雪夫2型設(shè)計(jì)帶通濾波器1；???????????????
Wp1=[9000?12000]/50000;???%數(shù)字頻率?fs1/2=50000?????????
Ws1=[8000?13000]/50000;
[n1Wn1]=cheb2ord（Wp1Ws1RpRs）;
[b1a1]=cheby2（n1RsWn1）;
[h1w1]=freqz（b1a1）;%顯示切比雪夫2型設(shè)計(jì)帶通濾波器頻率響應(yīng)
Wp2=[19000?23000]/50000;??%用切比雪夫2型設(shè)計(jì)帶通濾波器2；????
Ws2=[18000?24000]/50000;
[n2Wn2]=cheb2ord（Wp2Ws2RpRs）;
[b2a2]=cheby2（n2RsWn2）;
[h2w2]=freqz（b2a2）;
Wp3=[31000?35000]/50000;??%用切比

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-08-03?20:54??基于MATLAB的N路信號(hào)頻分復(fù)用系統(tǒng)的設(shè)計(jì)\
?????文件???????10577??2014-06-18?11:47??基于MATLAB的N路信號(hào)頻分復(fù)用系統(tǒng)的設(shè)計(jì)\123456.asv
?????文件???????10577??2014-06-18?11:50??基于MATLAB的N路信號(hào)頻分復(fù)用系統(tǒng)的設(shè)計(jì)\123456.m
?????文件??????260096??2015-08-03?19:49??基于MATLAB的N路信號(hào)頻分復(fù)用系統(tǒng)的設(shè)計(jì)\基于MATLAB的N路信號(hào)頻分復(fù)用系統(tǒng)的設(shè)計(jì).doc
?????文件??????395776??2015-08-03?19:48??基于MATLAB的N路信號(hào)頻分復(fù)用系統(tǒng)的設(shè)計(jì)\基于MATLAB的N路信號(hào)頻分復(fù)用系統(tǒng)的設(shè)計(jì).ppt
?????文件???????11638??2014-06-17?21:07??基于MATLAB的N路信號(hào)頻分復(fù)用系統(tǒng)的設(shè)計(jì)\完整版程序.asv
?????文件???????11634??2014-06-17?21:12??基于MATLAB的N路信號(hào)頻分復(fù)用系統(tǒng)的設(shè)計(jì)\完整版程序.m