通信原理课程设计--抑制载波的双边带调制,河南工业大学.doc

课 程 设 计 课程设计名称： 常规双边带调幅信号的仿真与分析 专 业 班 级 ： 学 生 姓 名 ： 学 号 ： 指 导 教 师 ： 课程设计时间： 1 需求分析 对频率为1Hz、功率为1的余弦信源，进行常规双边带调幅,其中载波频率为10Hz，并对以调信号进行频谱分析。之后用相干解调将原始信号恢复，并与原信号进行对比分析，最后分析信号在加性高斯白信道下的特性 2 概要设计 原始信号：y1=A1*cos(2*pi*t); 载波信号：y2=A2*cos(20*pi*t) 调制原理：y3=y1*y2 解调原理：y3*y2（与y2同频同相即可） 由调制频谱线性搬移 产生9Hz 和 11Hz的边频，经过同频同相的信号解调后，由于再次搬移，在解调出原信号后，会产生19Hz 和21Hz的高频，根据频率分布，由低通滤波器将原信号滤出。 因为噪声是加性噪声，因此表达式上相加即可对其分析。 3 运行环境 PC机、Windows XP或 Windows 7操作系统 4 开发工具和编程语言 MATLAB工程软件 5 详细设计 clear all; clc; A1=sqrt(2); A2=2; w1=2*pi; w2=20*pi; fs=60; t=0:1/fs:4 y0=A1*cos(w1*t) %*************** 生成“高斯噪声”******************* k=randn(1,length(y0)); plot(length(y0),k); y1=A1*cos(w1*t)+k ; ;figure(1);subplot(1,1,1) ;plot(y0);title(通过加性高斯白系统原信号); y2=A2*cos(w2*t) %******************调制************************ y3=y1.*y2 figure(2);subplot(2,1,1) plot(t,y3) ;title(以调制信号时域图) %*****************进行FFT变换并做频谱图******************* N1=length(t) k=-N1/2:N1/2-1 f=k*fs/N1 y=fftshift(fft(y3)); %进行fft变换 mag1=abs(y); %求幅值 figure(2);subplot(2,1,2) plot(f,mag1.^2);%做频谱图 title(以调制信号频谱); %**************解调************************ y4=y3.*y2 N2=length(y4); n=0:N2-1; f=n*fs/N2-fs/2; F2=fftshift(fft(y4)); mag2=abs(F2);%求幅值 figure(3); subplot(3,1,1) plot(f,mag2);%做频谱图 title(解调滤波前频谱); % ********切比雪夫1型低通滤波器*************** ft=fs; fp=2;f1=4 wp=2*pi*fp/ft; ws=2*pi*f1/ft; omegp=tan(wp/2);omegs=tan(ws/2); omega=1;omegb=omegp/omegs; [n,wn]=cheb1ord(omega,omegb,1,40,s); [b,a] =cheby1(n,1, wn, s); [num,den]=bilinear(b,a,0.5); %*************以调信号滤波，滤出信号***************** h1=filter(num,den,y4); %解调信号频谱分析 N3=length(t); n=0:N3-1; f=n*ft/N3-ft/2; F3=fftshift(fft(h1)); mag3=abs(F3);%求幅值 figure(3);subplot(3,1,2) plot(f,mag3);%做频谱图 title(解调滤波后含加性高斯噪声信号频谱); %*****************原信号于虑出波形对比********************* N=length(y1); n=0:N-1; f=n*fs/N-fs/2; F1=fftshift(fft(y1)); mag=abs(F1);%求幅值 figure(3);subplot(3,1,3) plot(f,mag);%做频谱图 title(