基于MATLAB的N路信号频分复用系统的设计

通信课程设计作业，有演示PPT、设计报告和程序。

%%%王丹妹-----电信112班------2011012489%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc;
%%%正弦信号的产生
fs=30000;%采样频率
t=（0:1/fs:0.005）;%图形坐标轴设计
f1=1000;f2=2000;f3=3000;%正弦信号频率
s1=sin（2*pi*f1*t）;%产生正弦信号
s2=sin（2*pi*f2*t）;
s3=sin（2*pi*f3*t）;
figure（1）
subplot（321）;plot（ts1）;xlabel（‘单位：s‘）;ylabel（‘幅度‘）;title（‘正弦信号1‘）;%显示正弦信号图形
subplot（323）;plot（ts2）;xlabel（‘单位：s‘）;ylabel（‘幅度‘）;title（‘正弦信号2‘）;
subplot（325）;plot（ts3）;xlabel（‘单位：s‘）;ylabel（‘幅度‘）;title（‘正弦信号3‘）;
%%%频域分析
N=1024;
Y1=fft（s1N）;Y1=fftshift（Y1）;%快速傅里叶变换得出频谱函数
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（‘函数频谱图1‘）;
subplot（324）;plot（fabs（Y2））;xlabel（‘单位：HZ‘）;ylabel（‘幅度‘）;title（‘函数频谱图2‘）;
subplot（326）;plot（fabs（Y3））;xlabel（‘单位：HZ‘）;ylabel（‘幅度‘）;title（‘函数频谱图3‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%SSB调制
fc1=10000;fc2=20000;fc3=30000;fs1=100000;
sm1 = modulate（s1fc1fs1‘amssb‘）;   %对信号进行SSB调制
sm2 = modulate（s2fc2fs1‘amssb‘）;   %对信号进行SSB调制
sfm3 = modulate（s3fc3fs1‘FM‘）;   %对信号进行FM调制
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调制的频域分析
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;     
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%将三路信号叠加合成为一路在信道中传输%%%%%%%%%%%%%%%%%%
dj1=sm1+sm2+sfm3;%三路信号叠加
dj3=awgn（dj150）;%三路SSB+G的叠加							
figure（3）	
subplot（221）;plot（tdj1）;xlabel（‘单位：s‘）;ylabel（‘幅度‘）;title（‘三路信号叠加‘）;
subplot（223）;plot（tdj3）;xlabel（‘单位：s‘）;ylabel（‘幅度‘）;title（‘三路信号+G叠加‘）;
%%%三路信号叠加合成为一路信号的频谱分析
DJ1=fft（dj1N）;DJ1=fftshift（DJ1）;  
DJ3=fft（dj3N）;DJ3=fftshift（DJ3）;  
subplot（222）;plot（f1abs（DJ1））;xlabel（‘单位：HZ‘）;ylabel（‘幅度‘）;title（‘三路信号叠加频谱图‘）;
subplot（224）;plot（f1abs（DJ3））;xlabel（‘单位：HZ‘）;ylabel（‘幅度‘）;title（‘三路信号+G叠加频谱图‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%带通滤波器的设计
fs1=100000;
Rp=0.5;Rs=40;  %用切比雪夫2型设计带通滤波器1；               
Wp1=[9000 12000]/50000;   %数字频率 fs1/2=50000         
Ws1=[8000 13000]/50000;
[n1Wn1]=cheb2ord（Wp1Ws1RpRs）;
[b1a1]=cheby2（n1RsWn1）;
[h1w1]=freqz（b1a1）;%显示切比雪夫2型设计带通滤波器频率响应
Wp2=[19000 23000]/50000;  %用切比雪夫2型设计带通滤波器2；    
Ws2=[18000 24000]/50000;
[n2Wn2]=cheb2ord（Wp2Ws2RpRs）;
[b2a2]=cheby2（n2RsWn2）;
[h2w2]=freqz（b2a2）;
Wp3=[31000 35000]/50000;  %用切比

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-08-03 20:54  基于MATLAB的N路信号频分复用系统的设计\
     文件       10577  2014-06-18 11:47  基于MATLAB的N路信号频分复用系统的设计\123456.asv
     文件       10577  2014-06-18 11:50  基于MATLAB的N路信号频分复用系统的设计\123456.m
     文件      260096  2015-08-03 19:49  基于MATLAB的N路信号频分复用系统的设计\基于MATLAB的N路信号频分复用系统的设计.doc
     文件      395776  2015-08-03 19:48  基于MATLAB的N路信号频分复用系统的设计\基于MATLAB的N路信号频分复用系统的设计.ppt
     文件       11638  2014-06-17 21:07  基于MATLAB的N路信号频分复用系统的设计\完整版程序.asv
     文件       11634  2014-06-17 21:12  基于MATLAB的N路信号频分复用系统的设计\完整版程序.m