基于Matlab的模拟调制与解调实验报告
基于matlab的-AM-FM调制与解调报告
基于matlab的-AM-FM调制与解调报告AM调制与解调100%% AMµ÷ÖÆfigure('Name','Ðźŵ÷Öƹý³ÌÖв¨Ðμ°ÆäƵÆ×','NumberTitle','off')a0=1;f0=2000;fc=20000;fs=1000000;t=[1:0.000001:4];am1=0*cos(2*pi*f0*t); %µ÷ÖÆÐźÅam=a0+am1;t1=cos(2*pi*fc*t); %Ôز¨s_am=am.*t1;AM1=fft(am1); T1=fft(t1); S_AM=fft(s_am);f=(0:3000000)*fs/3000001-fs/2;subplot(3,2,1); plot(t(1:1000),am1(1:1000)); title('ÐÅÏ¢ÐźŲ¨ÐÎ');subplot(3,2,2); plot(f,fftshift(abs(AM1))); title('ÐÅÏ¢ÐźÅƵÆ×');subplot(3,2,3); plot(t(1:1000),t1(1:1000)); title('Ôز¨ÐźÅ');subplot(3,2,4); plot(f,fftshift(abs(T1))); title('Ôز¨ÐźÅƵÆ×');subplot(3,2,5); plot(t(1:1000),s_am(1:1000)); title('Òѵ÷ÐźÅ');subplot(3,2,6); plot(f,fftshift(abs(S_AM))); title('Òѵ÷ÐźÅƵÆ×');%²úÉúÔëÉùfigure('Name','Ìí¼ÓÔëÉù¼°´øͨÂ˲¨¹ý³Ì²¨Ðμ°ÆäƵÆ×','NumberTitle','off ');snr=5;y=awgn(s_am,snr);fcuts=[16000 17500 22500 24000];mags=[0 1 0];devs=[0.05 0.01 0.05];[n,Wn,beta,ftype]=kaiserord(fcuts,mags,devs,fs);hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');[H,f]=freqz(hh,1,1024,fs);st_p=fftfilt(hh,y);Q=fft(st_p);Y=fft(y);f=(0:3000000)*fs/3000001-fs/2;subplot(2,2,1);plot(t(1500001:1501000),y(1500001:1501000));title('Ìí¼ÓÔëÉùºóÐźŲ¨ÐÎ');subplot(2,2,2);plot(f,fftshift(abs(Y)));title('Ìí¼ÓÔëÉùºóÐźÅƵÆ×');subplot(2,2,3);plot(t(1500001:1501000),st_p(1500001:1501000));title('´øͨÂ˲¨ºóÐźŲ¨ÐÎ');subplot(2,2,4); plot(f,fftshift(abs(Q)));title('´øͨÂ˲¨ºóÐźÅƵÆ×');%½âµ÷figure('Name','Ïà¸É½âµ÷ËùµÃ²¨Ðμ°ÆäƵÆ×','NumberTitle','off');ss_am=st_p.*t1;SS_AM=fft(ss_am)f=(0:3000000)*fs/3000001-fs/2;subplot(2,1,1);plot(t(1500001:1503000),ss_am(1500001:1503000));title( 'Ïà³ËÐźÅ');subplot(2,1,2); plot(f,fftshift(abs(SS_AM)));title('Ïà³ËÐźÅƵÆ×');fcuts1=[2500,30000];mags1=[1 0];devs1=[0.01 0.05];[n,Wn,beta,ftype]=kaiserord(fcuts1,mags1,devs1,fs);hh1=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');[H,f]=freqz(hh1,1,1024,fs);m0=fftfilt(hh1,ss_am);M0=fft(m0);f=(0:3000000)*fs/3000001-fs/2;figuresubplot(2,1,1);plot(t(1500001:1501000),m0(1500001:1501000));title('½âµ÷ÐźÅ');subplot(2,1,2); plot(f,fftshift(abs(M0)));title('½âµ÷ÐźÅƵÆ×');11.0005 1.001-101信息信号波形-505x 1050126信息信号频谱11.0005 1.001-101载波信号-505x 1050126载波信号频谱11.00051.001-202已调信号-55x 105012x 106已调信号频谱2.52.50052.501-4-2024添加噪声后信号波形-505x 105051015x 105添加噪声后信号频谱2.52.50052.501-2024带通滤波后信号波形-55x 10500.511.526带通滤波后信号频谱50%% AMµ÷ÖÆ2.5 2.5005 2.501 2.5015 2.502 2.5025 2.503 2.5035-1012相乘信号-5-4-3-2-1012345x 1050510155相乘信号频谱2.52.50012.50022.50032.50042.50052.50062.50072.50082.5009 2.501-0.500.51解调信号-5-4-3-2-1012345x 1050510155解调信号频谱figure('Name','Ðźŵ÷Öƹý³ÌÖв¨Ðμ°ÆäƵÆ×','NumberTitle','off')a0=2;f0=2000;fc=20000;fs=1000000;t=[1:0.000001:4];am1=0*cos(2*pi*f0*t); %µ÷ÖÆÐźÅam=a0+am1;t1=cos(2*pi*fc*t); %Ôز¨s_am=am.*t1;AM1=fft(am1); T1=fft(t1); S_AM=fft(s_am);f=(0:3000000)*fs/3000001-fs/2;subplot(3,2,1); plot(t(1:1000),am1(1:1000)); title('ÐÅÏ¢ÐźŲ¨ÐÎ');subplot(3,2,2); plot(f,fftshift(abs(AM1))); title('ÐÅÏ¢ÐźÅƵÆ×');subplot(3,2,3); plot(t(1:1000),t1(1:1000)); title('Ôز¨ÐźÅ');subplot(3,2,4); plot(f,fftshift(abs(T1))); title('Ôز¨ÐźÅƵÆ×');subplot(3,2,5); plot(t(1:1000),s_am(1:1000)); title('Òѵ÷ÐźÅ');subplot(3,2,6); plot(f,fftshift(abs(S_AM))); title('Òѵ÷ÐźÅƵÆ×');%²úÉúÔëÉùfigure('Name','Ìí¼ÓÔëÉù¼°´øͨÂ˲¨¹ý³Ì²¨Ðμ°ÆäƵÆ×','NumberTitle','off ');snr=5;y=awgn(s_am,snr);fcuts=[16000 17500 22500 24000];mags=[0 1 0];devs=[0.05 0.01 0.05];[n,Wn,beta,ftype]=kaiserord(fcuts,mags,devs,fs);hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');[H,f]=freqz(hh,1,1024,fs);st_p=fftfilt(hh,y);Q=fft(st_p);Y=fft(y);f=(0:3000000)*fs/3000001-fs/2;subplot(2,2,1);plot(t(1500001:1501000),y(1500001:1501000));title('Ìí¼ÓÔëÉùºóÐźŲ¨ÐÎ');subplot(2,2,2);plot(f,fftshift(abs(Y)));title('Ìí¼ÓÔëÉùºóÐźÅƵÆ×');subplot(2,2,3);plot(t(1500001:1501000),st_p(1500001:1501000));title('´øͨÂ˲¨ºóÐźŲ¨ÐÎ');subplot(2,2,4); plot(f,fftshift(abs(Q)));title('´øͨÂ˲¨ºóÐźÅƵÆ×');%½âµ÷figure('Name','Ïà¸É½âµ÷ËùµÃ²¨Ðμ°ÆäƵÆ×','NumberTitle','off');ss_am=st_p.*t1;SS_AM=fft(ss_am)f=(0:3000000)*fs/3000001-fs/2;subplot(2,1,1);plot(t(1500001:1503000),ss_am(1500001:1503000));title( 'Ïà³ËÐźÅ');subplot(2,1,2); plot(f,fftshift(abs(SS_AM)));title('Ïà³ËÐźÅƵÆ×'); fcuts1=[2500,30000];mags1=[1 0];devs1=[0.01 0.05];[n,Wn,beta,ftype]=kaiserord(fcuts1,mags1,devs1,fs);hh1=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');[H,f]=freqz(hh1,1,1024,fs);m0=fftfilt(hh1,ss_am);M0=fft(m0);f=(0:3000000)*fs/3000001-fs/2;figuresubplot(2,1,1);plot(t(1500001:1501000),m0(1500001:1501000));title('½âµ÷ÐźÅ');subplot(2,1,2); plot(f,fftshift(abs(M0)));title('½âµ÷ÐźÅƵÆ×');11.0005 1.001-101信息信号波形-505x 1050126信息信号频谱11.0005 1.001-101载波信号-505x 1050126载波信号频谱11.00051.001-505已调信号-55x 105024x 106已调信号频谱2.52.50052.501-505添加噪声后信号波形-505x 1050123x 106添加噪声后信号频谱2.52.50052.501-4-2024带通滤波后信号波形-55x 105012346带通滤波后信号频谱0%% AMµ÷ÖÆ2.5 2.5005 2.501 2.5015 2.502 2.5025 2.503 2.5035-10123相乘信号-5-4-3-2-1012345x 10501236相乘信号频谱2.52.50012.50022.50032.50042.50052.50062.50072.50082.5009 2.50100.511.5解调信号-5-4-3-2-1012345x 10501236解调信号频谱figure('Name','Ðźŵ÷Öƹý³ÌÖв¨Ðμ°ÆäƵÆ×','NumberTitle','off')a0=10^100;f0=2000;fc=20000;fs=1000000;t=[1:0.000001:4];am1=0*cos(2*pi*f0*t); %µ÷ÖÆÐźÅam=a0+am1;t1=cos(2*pi*fc*t); %Ôز¨s_am=am.*t1;AM1=fft(am1); T1=fft(t1); S_AM=fft(s_am);f=(0:3000000)*fs/3000001-fs/2;subplot(3,2,1); plot(t(1:1000),am1(1:1000)); title('ÐÅÏ¢ÐźŲ¨ÐÎ');subplot(3,2,2); plot(f,fftshift(abs(AM1))); title('ÐÅÏ¢ÐźÅƵÆ×');subplot(3,2,3); plot(t(1:1000),t1(1:1000)); title('Ôز¨ÐźÅ');subplot(3,2,4); plot(f,fftshift(abs(T1))); title('Ôز¨ÐźÅƵÆ×');subplot(3,2,5); plot(t(1:1000),s_am(1:1000)); title('Òѵ÷ÐźÅ');subplot(3,2,6); plot(f,fftshift(abs(S_AM))); title('Òѵ÷ÐźÅƵÆ×');%²úÉúÔëÉùfigure('Name','Ìí¼ÓÔëÉù¼°´øͨÂ˲¨¹ý³Ì²¨Ðμ°ÆäƵÆ×','NumberTitle','off ');snr=5;y=awgn(s_am,snr);fcuts=[16000 17500 22500 24000];mags=[0 1 0];devs=[0.05 0.01 0.05];[n,Wn,beta,ftype]=kaiserord(fcuts,mags,devs,fs);hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');[H,f]=freqz(hh,1,1024,fs);st_p=fftfilt(hh,y);Q=fft(st_p);Y=fft(y);f=(0:3000000)*fs/3000001-fs/2;subplot(2,2,1);plot(t(1500001:1501000),y(1500001:1501000));title('Ìí¼ÓÔëÉùºóÐźŲ¨ÐÎ');subplot(2,2,2);plot(f,fftshift(abs(Y)));title('Ìí¼ÓÔëÉùºóÐźÅƵÆ×');subplot(2,2,3);plot(t(1500001:1501000),st_p(1500001:1501000));title('´øͨÂ˲¨ºóÐźŲ¨ÐÎ');subplot(2,2,4); plot(f,fftshift(abs(Q)));title('´øͨÂ˲¨ºóÐźÅƵÆ×');%½âµ÷figure('Name','Ïà¸É½âµ÷ËùµÃ²¨Ðμ°ÆäƵÆ×','NumberTitle','off');ss_am=st_p.*t1;SS_AM=fft(ss_am)f=(0:3000000)*fs/3000001-fs/2;subplot(2,1,1);plot(t(1500001:1503000),ss_am(1500001:1503000));title( 'Ïà³ËÐźÅ');subplot(2,1,2); plot(f,fftshift(abs(SS_AM)));title('Ïà³ËÐźÅƵÆ×'); fcuts1=[2500,30000];mags1=[1 0];devs1=[0.01 0.05];[n,Wn,beta,ftype]=kaiserord(fcuts1,mags1,devs1,fs);hh1=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');[H,f]=freqz(hh1,1,1024,fs);m0=fftfilt(hh1,ss_am);M0=fft(m0);f=(0:3000000)*fs/3000001-fs/2;figuresubplot(2,1,1);plot(t(1500001:1501000),m0(1500001:1501000));title('½âµ÷ÐźÅ');subplot(2,1,2); plot(f,fftshift(abs(M0)));title('½âµ÷ÐźÅƵÆ×');11.0005 1.001-101信息信号波形-505x 1050126信息信号频谱1 1.0005 1.001-101载波信号-505x 1050126载波信号频谱11.00051.001-101x 10100已调信号-55x 105012x 10106已调信号频谱2.52.5005 2.501-1-0.500.51x 10100添加噪声后信号波形-505x 105051015x 10105添加噪声后信号频谱2.52.5005 2.501-2-1012100带通滤波后信号波形-55x 10501106带通滤波后信号频谱FM 调制与解调%%FMfigure('Name','FMµ÷ÖƲ¨ÐÎÓëƵÆ×')2.5 2.5005 2.501 2.5015 2.502 2.5025 2.503 2.5035-50510x 1099相乘信号-5-4-3-2-1012345x 105051015105相乘信号频谱2.52.50012.50022.50032.50042.50052.50062.50072.50082.5009 2.5013.93.913.923.93x 1099解调信号-1-0.500.51x 105123105解调信号频谱f0=2000; fc=20000; fs=1000000; Am=1; kf=0.5; Tc=8; Ta=0.001; dt=0.000001;t=[0:1/fs:3];f=(0:length(t)-1)*fs/(length(t))-fs/2;fm0=cos(2*pi*f0*t);mt=fm0;%»ý·ÖÆ÷Éè¼Æw1=0;w2=0;for m=1:length(t)w1=mt(m)+w2;w2=mt(m)+w1;fi(m)=w1/(2*fs);endfi=fi*2*pi/max(abs(fi));I=cos(kf*fi);Q=sin(kf*fi);y1=Am*cos(2*pi*fc*t).*I-Am*sin(2*pi*fc*t).*Q;subplot(2,1,1);plot(t,y1);title('²¨ÐÎ')axis([1e-3 4e-3 -2 2]);Y1=fft(y1);subplot(2,1,2);plot(f,fftshift(abs(Y1))/1e6); title('ƵÆ×') %%Ôز¨ÆµÆ×axis([-4e4 4e4 0 1]);figure('Name','FMµ÷Öƺó¼ÓÔëÉù²¨ÐÎÓë½âµ÷ºó²¨ÐÎÒÔ¼°Â˳ýÖ±Á÷·ÖÁ¿ºóµÄ²¨ÐÎ')y1o=awgn(y1,40);subplot(3,1,1);plot(t,y1o); title('¼ÓÔëÉùºó²¨ÐÎ') %%¼ÓÔëÉùºóµÄÐźÅaxis([1e-3 4e-3 -2 2]);%%´øͨÂ˲¨KSband=2*(3+1)*f0;fcutsb=[fc-KSband-2000 fc-KSband fc+KSbandfc+KSband+2000]; %%½ÓÊÕ»úÇ°¶Ë´øͨÂ˲¨magsb=[0 1 0];devsb=[0.05 0.01 0.05];[nb,Wnb,betab,ftypeb]=kaiserord(fcutsb,magsb,devsb,fs);hhb=fir1(nb,Wnb,ftypeb,kaiser(nb+1,betab),'noscale'); %´øͨÂ˲¨Æ÷£»st_pb=fftfilt(hhb,y1o);subplot(3,1,2);st_pb=st_pb/1e6;plot(t,st_pb); title('´øͨÂ˲¨Æ÷ºóµÄ²¨ÐÎ')axis([1e-3 4e-3 -2e-6 2e-6]);%΢·ÖÆ÷Éè¼Æfor i=1:length(t)-1 %½ÓÊÕÐźÅͨ¹ý΢·ÖÆ÷´¦Àídiff_st_pb(i)=(st_pb(i+1)-st_pb(i))/dt;endsfm=abs(hilbert(diff_st_pb));subplot(3,1,2);plot(t,[sfm*20 0]);axis([1e-3 4e-3 0 4]);%%¸ôÖ±% KSbandh=2*(3+1)*f0;fcutsh=[0.01 3000];magsh=[0 1];devsh=[0.01 0.05];[nh,Wnh,betah,ftypeh]=kaiserord(fcutsh,magsh,devsh,fs);hhh=fir1(nh,Wnh,ftypeh,kaiser(nh+1,betah),'noscale');sfm_out=fftfilt(hhh,sfm*20);subplot(3,1,3);plot(t,[sfm_out 0]);title('¸ôÖ±ºóµÄ²¨ÐÎ')axis([1e-3 4e-3 -2 2]);11.522.533.54x 10-3-2-1012波形-4-3-2-101234x 10400.51频谱11.522.533.54x 10-3-202加噪声后波形11.522.53 3.54x 10-302411.522.533.54x 10-3-202隔直后的波形。
基于MATLAB的模拟调制实验报告
基于MATLAB 的模拟调制实验报告一、实验目的1.进一步学习调制的知识，掌握调频与调角两种模拟调制技术。
2.进一步学习MATLAB 的编程，熟练使用MATLAB 进行作图。
二、实验原理1.调制的概念调制（modulation ）就是对信号源的信息进行处理加到载波上，使其变为适 合 于信道传输的形式的过程，是使载波随信号而改变的技术。
一般，用来传送消息的信号()t u c 叫作载波或受调信号，代表所欲传送消息的信号叫作调制信号，调制后的信号()t u 叫作已调信号。
用调制信号()t u Ω控制载波的某些参数，使之随()t u Ω而变化，就可实现调制。
2.调制的目的 ➢ 频谱变换当所要传送的信号的频率或者太低，或者频带很宽，对直接采用电磁波的形 式进行发送很不利，需要的天线尺寸很大，而且发射和接受短的天线与谐振回路的参数变化范围很大。
为了信息有效与可靠传输，往往需要将低频信号的基带频谱搬移到适当的或指定的频段。
这样可以提高传输性能，以较小的发送功率与较短的天线来辐射电磁波。
➢ 实现信道复用为了使多个用户的信号共同利用同一个有较大带宽的信道，可以采用各种复用技术。
如模拟电话长途传输是通过利用不同频率的载波进行调制。
将各用户话音每隔4 kHz 搬移到高频段进行传输。
➢ 提高抗干扰能力不同的调制方式，在提高传输的有效性和可靠性方面各有优势。
如调频广播系统，它采用的频率调制技术，付出多倍带宽的代价，由于抗干扰性能强，其音质比只占10 kHz 带宽的调幅广播要好得多。
扩频通信就是以大大扩展信号传输带宽，以达到有效抗拒外部干扰和短波信道多径衰落的特殊调制方式。
3.调制的种类根据()t u Ω和()t u c 的不同类型和完成调制功能的调制器传递函数不同，调制分为以下多种方式： （1）．按调制信号()t u Ω的类型分为：● 模拟调制：调制信号()t u Ω是连续变化的模拟量，如话音与图像信号。
● 数字调制：调制信号是数字化编码符号或脉冲编码波形。
matlab实验模拟调制解调器
实验四 模拟调制解调器一． 实验目的• 1.掌握模拟调制的方法原理• 2掌握模拟解调的方法原理二． 实验要求1熟悉matlab 中模拟调制解调函数的功能并能熟练调用；2信号为低频正弦信号3对其进行模拟调制，得到调制信号4假设信道为理想信道，噪声为加性、不相关白噪声 5对经信道的调制信号进行解调，恢复出原始信号 6对以上各个节点画出波形图三． 基本原理幅度是指对信号进行幅度调制。
一般做法就是先在原信号上叠加一个直流信号，以保证信号f(t)+A>0然后乘上一个高频的余弦信号，即得到g(t)=[f(t)+A]coswt 在频域上的效果就是将原信号的频谱移动到w 处，以适合信道传输的最佳频率范围。
g(t)的包络线即f(t)+A,用一个简单的包络检测电路就可以接收并还原信号了。
解调：将已调信号(t)V o 与载波(t)V c 相乘，一部分完整的信号频谱将会搬回低频上，通过一个低通滤波器就可以取出信号的频率成分。
四． 实验流程 1. 初始化各项参数。
2. 产生低频正弦信号 3. 进行幅度调制 4. 进行相位调制5. 对幅度调制信号和相位调制信号加入白噪声6. 分别对幅度调制信号和相位调制信号进行解调。
图4-1 幅度调制图4-2 相位调制五． 实验结果及分析-0.05-0.04-0.03-0.02-0.0100.010.020.030.040.05时间t原调制信号的时域图-500-400-300-200-100100200300400500频率f原调制信号的频谱图-0.05-0.04-0.03-0.02-0.0100.010.020.030.040.05时间t调幅信号的时域图-500-400-300-200-100100200300400500频率f调幅信号的频谱图调相信号的时域图-0.05-0.04-0.03-0.02-0.0100.010.020.030.040.05时间t调相信号的频谱图-500-400-300-200-1000100200300400500频率f调频信号的时域图-0.05-0.04-0.03-0.02-0.0100.010.020.030.040.05时间t调频信号的频谱图-500-400-300-200-1000100200300400500频率f已调信号的包络反映调制信号的变化。
基于MATLAB的模拟信号频率调制与解调分析
基于MATLAB的模拟信号频率调制与解调分析信号频率调制（FM）是一种将信息信号调制到载频波形上以便在传输过程中保持信号质量的技术。
本文将基于MATLAB对信号频率调制与解调进行分析与模拟。
首先，我们需要生成一个调制信号。
以正弦信号为例，通过改变该信号的频率来模拟调制信号。
我们可以使用MATLAB的信号处理工具箱中的`fmmod(`函数来实现这一点。
以下是一个示例代码：```matlabt = 0:1/fs:1; % 时间向量fc = 2000; % 载频频率fm = 100; % 调制信号频率m = sin(2*pi*fm*t); % 调制信号modulatedSignal = fmmod(m, fc, fs); % 使用fmmod进行调频调制subplot(2,1,1);plot(t, m);title('调制信号');xlabel('时间');ylabel('振幅');subplot(2,1,2);title('调制后信号');xlabel('时间');ylabel('振幅');```上述代码中，我们定义了采样频率、时间向量、载频频率和调制信号频率，并生成了调制信号。
然后，我们使用`fmmod(`函数将调制信号调制到载频波形上。
最后，我们用两个子图分别显示调制信号和调制后信号。
接下来，我们将对调制后的信号进行解调以还原原始信号。
我们可以使用MATLAB的信号处理工具箱中的`fmdemod(`函数。
以下是一个示例代码：```matlabdemodulatedSignal = fmdemod(modulatedSignal, fc, fs); % 使用fmdemod进行解调subplot(2,1,1);plot(t, modulatedSignal);title('调制后信号');xlabel('时间');ylabel('振幅');subplot(2,1,2);title('解调后信号');xlabel('时间');ylabel('振幅');```上述代码中，我们使用`fmdemod(`函数对调制后的信号进行解调。
MATLAB实现信号的调制与解调
实验8 信号调制与解调[实验目的]1．了解用MATLAB实现信号调制与解调的方法。
2．了解几种基本的调制方法。
[实验原理]由于从消息变换过来的原始信号具有频率较低的频谱分量，这种信号在许多信道中不适宜传输。
因此，在通信系统的发送端通常需要有调制过程，而在接收端则需要有反调制过程——解调过程。
所谓调制，就是按调制信号的变化规律去改变某些参数的过程。
调制的载波可以分为两类：用正弦信号作载波；用脉冲串或一组数字信号作为载波。
最常用和最重要的模拟调制方式是用正弦波作为载波的幅度调制和角度调制。
本实验中重点讨论幅度调制。
幅度调制是正弦型载波的幅度随调制信号变化的过程。
设正弦载波为式中——载波角频率——载波的初相位A——载波的幅度那么，幅度调制信号（已调信号）一般可表示为式中，m(t)为基带调制信号。
在MATLAB中，用函数y=modulate(x,fc,fs,’s’)来实现信号调制。
其中fc 为载波频率，fs为抽样频率，’s’省略或为’am-dsb-sc’时为抑制载波的双边带调幅，’am-dsb-tc’为不抑制载波的双边带调幅，’am-ss b’为单边带调幅，’pm’为调相，’fm’为调频。
[课上练习]产生AM FM PM signals[实验内容]0. 已知信号，当对该信号取样时，求能恢复原信号的最大取样周期。
设计MATALB 程序进行分析并给出结果。
1．有一正弦信号, n=[0:256]，分别以100000Hz的载波和1000000Hz的抽样频率进行调幅、调频、调相，观察图形。
2．对题1中各调制信号进行解调（采用demod函数），观察与原图形的区别3．已知线性调制信号表示式如下：⑴⑵式中，试分别画出它们的波形图和频谱图4．已知调制信号，载波为cos104t，进行单边带调制，试确定单边带信号的表示式，并画出频谱图。
[实验要求]1 自行编制完整的实验程序，实现对信号的模拟，并得出实验结果。
2 在实验报告中写出完整的自编程序，并给出实验结果和分析，学习demod 函数对调制信号进行解调的分析。
基于Matlab的模拟信号传输系统实验报告
实验七、基于Matlab 的模拟信号传输系统实验报告幅度的调制解调1，AM 调制原理2，AM 的解调原理Sam(t）Sp(t)Sd(t)c(t)=coswct相干解调的一般模型%041210308AM 调制解调部分代码close all ;clear;clc;……figure;yin=awgn(y1,3);subplot(2,1,1);plot(t,yin);title('加噪后信号');%%加噪声后的信号axis([1e-34e-3-55]);fcutsb=[16000175002250024000];%%接收机前端带通滤波magsb=[010];devsb=[0.050.010.05];[nb,Wnb,betab,ftypeb]=kaiserord (fcutsb,magsb,devsb,fs);hhb=fir1(nb,Wnb,ftypeb,kaiser(n b+1,betab),'noscale');%带通滤波器；[Hb,f]=freqz(hhb,1,1024,fs);st_pb=fftfilt(hhb,yin);subplot(2,1,2);plot(t,st_pb);title('加噪后频谱');axis([1e-34e-3-55]);figure;st_noise=st_pb.*xc;subplot(2,1,1);plot(t,st_noise);title('带通滤波后');axis([1e-34e-3-22]);fcutsl=[250030000];%%混频后低通滤波magsl=[10];devsl=[0.010.05];[nl,Wnl,betal,ftypel]=kaiserord(fcutsl,magsl,devsl,fs);hhl=fir1(nl,Wnl,ftypel,kaiser(n l+1,betal),'noscale');st_pl=fftfilt(hhl,st_noise);subplot(2,1,2);plot(t,st_pl);title('带通滤波后频谱');axis([1e-34e-301]);grid;LPF角度的调制解调FM的解调原理窄带的相干解调原理框图如下所示：Ｓfm(t)Si(t)C（t）窄带信号的相干解调FM的非相干解调模型%041210308FM调制解调部分代码……w1=0;w2=0;for m=1:length(t)w1=mt(m)+w2;w2=mt(m)+w1;fi(m)=w1/(2*fs);endfi=fi*2*pi/max(abs(fi));HI=cos(kfm*fi);HQ=sin(kfm*fi);yo=A*cos(2*pi*fc*t).*HI-A*sin(2 *pi*fc*t).*HQ;subplot(2,1,1);plot(t,yo);title('时域波形'); axis([1e-34e-3-22]);Yo=fft(yo);subplot(2,1,2);plot(f,fftshift(abs(Yo))/1e6); title('频谱图');%%载波频谱axis([-4e44e401]);figure;%%添加噪声yoo=awgn(yo,30);subplot(3,1,1);plot(t,yoo);title('加噪后波形');%%加噪声后的信号axis([1e-34e-3-22]);%%带通滤波KSband=2*(3+1)*f0;fcutsb=[fc-KSband-2000fc-KSband fc+KSband fc+KSband+2000];%%接收机前端带通滤波magsb=[010];devsb=[0.050.010.05];[nb,Wnb,betab,ftypeb]=kaiserord (fcutsb,magsb,devsb,fs);hhb=fir1(nb,Wnb,ftypeb,kaiser(n b+1,betab),'noscale');%带通滤波器；st_pb=fftfilt(hhb,yoo);%subplot(3,1,2);st_pb=st_pb/1e6;%plot(t,st_pb);%axis([1e-34e-3-2e-62e-6]);%figure;%Yst_pb=fft(st_pb);%subplot(2,1,1);%plot(f,fftshift(abs(Yst_pb))); %%载波频谱BPF LPF微分%axis([-4e44e401e6]);%figure('Name','FM解调波形')%微分器设计for i=1:length(t)-1%接收信号通过微分器处理diff_st_pb(i)=(st_pb(i+1)-st_pb (i))/dt;endst_noise=abs(hilbert(diff_st_pb) );subplot(3,1,2);plot(t,[st_noise*200]);title('解调后波形');axis([1e-34e-304]);%%隔直%KSbandh=2*(3+1)*f0;fcutsh=[0.013000];magsh=[01];devsh=[0.010.05];[nh,Wnh,betah,ftypeh]=kaiserord (fcutsh,magsh,devsh,fs);hhh=fir1(nh,Wnh,ftypeh,kaiser(n h+1,betah),'noscale');st_out=fftfilt(hhh,st_noise*20) ;subplot(3,1,3);plot(t,[st_out0]);title('解调后波形');axis([1e-34e-3-22]);。
基于MATLAB的GMSK调制与解调课设报告
基于Matlab的GMSK调制与解调1.课程设计目的（1）加深对GMSK基本理论知识的理解。
（2）培养独立开展科研的能力和编程能力。
（3）通过SIMULINK对BT=0.3的GMSK调制系统进行仿真。
2.课程设计要求（1）观察基带信号和解调信号波形。
（2）观察已调信号频谱图。
（3）分析调制性能和BT参数的关系。
3.相关知识3.1GMSK调制调制原理图如图2.2，图中滤波器是高斯低通滤波器，它的输出直接对VCO 进行调制，以保持已调包络恒定和相位连续。
图3.1 GMSK调制原理图为了使输出频谱密集，前段滤波器必须具有以下待性：1.窄带和尖锐的截止特性，以抑制FM调制器输入信号中的高频分量；2.脉冲响应过冲量小，以防止FM调制器瞬时频偏过大；3.保持滤波器输出脉冲响应曲线下的面积对应丁pi／2的相移。
以使调制指数为1／2。
前置滤波器以高斯型最能满足上述条件，这也是高斯滤波器最小移频键控(GMSK)的由来。
3.2GMSK 解调GMSK本是MSK 的一种，而MSK 又是是FSK 的一种，因此，GMSK 检波也可以采用FSK 检波器，即包络检波及同步检波。
而GMSK 还可以采用时延检波，但每种检波器的误码率不同。
GMSK 非相干解调原理图如图2.3，图中是采用FM 鉴频器（斜率鉴频器或相位鉴频器）再加判别电路，实现GMSK 数据的解调输出。
图3.2 GMSK 解调原理图4.课程设计分析 4.1信号发生模块因为GMSK 信号只需满足非归零数字信号即可，本设计中选用（Bernoulli Binary Generator ）来产生一个二进制序列作为输入信号。
图4.1 GMSK 信号产生器该模块的参数设计这只主要包括以下几个。
其中probability of a zero 设置为0.5表示产生的二进制序列中0出现的概率为0.5；Initial seed 为61表示随机数种子为61；sample time 为1/1000表示抽样时间即每个符号的持续时为0.001s 。
基于Matlab调制与解调的实现(DOC)
基于Matlab 调制与解调的实现一．实验目的1.熟悉Matlab 的使用2.掌握幅度调制、角度调制及FSK 调制的基本原理3.掌握解调的基本原理，并实现解调二．实验原理，仿真及结果分析 AM 调制与解调1.标准AM 波调制与解调的原理调制信号是只来来自信源的调制信号（基带信号），这些信号可以是模拟的，亦可以是数字的。
为首调制的高频振荡信号可称为载波，它可以是正弦波，亦可以是非正弦波(如周期性脉冲序列)。
载波由高频信号源直接产生即可，然后经过高频功率放大器进行放大，作为调幅波的载波，调制信号由低频信号源直接产生，二者经过乘法器后即可产生双边带的调幅波。
设载波信号的表达式为t c ωcos ，调制信号的表达式为t A t m m m ωcos )(= ，则调幅信号的表达式为t t m A t s c AM ωcos )]([)(0+=标准调幅波示意图从高频已调信号中恢复出调制信号的过程称为解调，又称为检波。
对于振幅调制信号，解调就是从它的幅度变化上提取调制信号的过程。
解调是调制的逆过程。
)(t m )(t可利用乘积型同步检波器实现振幅的解调，让已调信号与本地恢复载波信号相乘并通过低通滤波可获得解调信号。
2．matlab仿真% ======================载波信号===========================t=-1:0.00001:1;A0=10; %载波信号振幅f=6000; %载波信号频率w0=f*pi;Uc=A0*cos(w0*t); %载波信号figure(1);subplot(2,1,1);plot(t,Uc);title('载频信号波形');axis([0,0.01,-15,15]);subplot(2,1,2);Y1=fft(Uc); %对载波信号进行傅里叶变换plot(abs(Y1));title('载波信号频谱');axis([5800,6200,0,1000000]);% ======================调制信号============================== t=-1:0.00001:1;A1=5; %调制信号振幅f=6000; %载波信号频率w0=f*pi;mes=A1*cos(0.001*w0*t); %调制信号subplot(2,1,1);plot(t,mes);xlabel('t'),title('调制信号');subplot(2,1,2);Y2=fft(mes); % 对调制信号进行傅里叶变换plot(abs(Y2));title('调制信号频谱');axis([198000,202000,0,1000000]);% =======================AM已调信号=========================t=-1:0.00001:1;A0=10; %载波信号振幅A1=5; %调制信号振幅A2=3; %已调信号振幅f=3000; %载波信号频率w0=2*f*pi;m=0.15; %调制度mes=A1*cos(0.001*w0*t); %消调制信号Uam=A2*(1+m*mes).*cos((w0).*t); %AM 已调信号subplot(2,1,1);plot(t,Uam);grid on;title('AM调制信号波形');subplot(2,1,2);Y3=fft(Uam); % 对AM已调信号进行傅里叶变换plot(abs(Y3)),grid;title('AM调制信号频谱');axis([5950,6050,0,500000]);%=========================FIR低通滤波器=======================Ft=2000; %采样频率fpts=[100 120]; %通带边界频率fp=100Hz，阻带截止频率fs=120Hzmag=[1 0];dev=[0.01 0.05]; %通带波动1%，阻带波动5%[n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev,Ft);%kaiserord估计采用凯塞窗设计的FIR滤波器的参数b21=fir1(n21,wn21,Kaiser(n21+1,beta)); %由fir1设计滤波器[h,w]=freqz(b21,1); %得到频率响应plot(w/pi,abs(h));grid ontitle('FIR低通滤波器');%=========================AM信号解调=======================t=-1:0.00001:1;A0=10; %载波信号振幅A1=5; %调制信号振幅A2=3; %已调信号振幅f=3000; %载波信号频率w0=2*f*pi;m=0.15; %调制度k=0.5 ; %DSB 前面的系数mes=A1*cos(0.001*w0*t); %调制信号Uam=A2*(1+m*mes).*cos((w0).*t); %AM 已调信号Dam=Uam.*cos(w0*t); %对AM调制信号进行解调subplot(4,2,1);plot(t,Dam);title('滤波前AM解调信号波形');subplot(4,2,2);axis([187960,188040,0,200000]);Y5=fft(Dam); % 对AM解调信号进行傅里叶变换plot(abs(Y5)),grid;title('滤波前AM解调信号频谱');subplot(4,2,3);plot(t,z21);title('滤波后的AM解调信号波形');T5=fft(z21); %求AM信号的频谱subplot(4,2,4);plot(abs(T5));title('滤波后的AM解调信号频谱');axis([198000,202000,0,100000]);角度调制与解调角度调制是频率调制和相位调制的总称。
MATLAB仿真AM调制解调 无线通信实验报告
无线通信实验报告院系名称：信息科学与工程学院专业班级：电信班学生姓名：学号：授课教师：2014 年11 月 6 日实验一 高斯衰落信道建模一、基本原理QPSK 信号可以看成是对两个正交的载波进行多电平双边带调制后所得信号的叠加，因此可以用正交调制的方法得到QPSK 信号。
QPSK 信号的星座如图4.1.1所示：图1.1 QPSK 信号星座图从AWGN 信道中，在一个信号区间内接收到的带宽信号可以表示为()()()()()cos(2)()sin(2)m m c c s c r t u t n t u t n t f t n t f t ππ=+=+-这里()c n t 和()s n t 是加性噪声的两个正交分量。
可以将这个接收信号与1()()cos(2)T c t g t f t ψπ=，2()sin(2)T c g t f t ψπ=-给出的1()t ψ和2()t ψ作相关，两个相关器的输出产生受噪声污损的信号分量，它们可表示为22()m s s s m m r s n n n M Mππξξ=+=++ 式中c n 和s n 定义为 1()()2c T c n g t n t dt ∞-∞=⎰ 1()()2s T s n g t n t dt ∞-∞=⎰ 这两个正交噪声分量()c n t 和()s n t 是零均值，互不相关的高斯随机过程。
这样，()()0c s E n E n ==和()0c s E n n =。
c n 和s n 的方差是 220()()2c s N E n E n == 最佳检测器将接收信号向量r 投射到M 个可能的传输信号向量{m s }之一上去，并选取对应于最大投影的向量。
据此，得到相关准则为(,)m m C r s r s =•，m=0,1,…，M-1由于全部信号都具有相等的能量，因此，对数字相位调制一种等效的检测器标准是计算接收信号向量r=(c r ,s r )的相位为 arctan s r cr r θ= 并从信号集{m s }中选取其相位最接近r θ的信号。
matlab实验报告 模拟调制解调
实验报告姓名：李鹏博实验名称：模拟调制解调学号：2011300704 课程名称：数字信号处理班级：03041102 实验室名称：航海西楼303组号： 1 实验日期：2014.06.20一、实验目的、要求1．掌握掌握模拟调制以及对应解调方法的原理。
2．掌握模拟调制解调方法的计算机编程实现方法，即软件实现。
二、实验原理调制的目的是把要传输的模拟信号或数字信号变换成适合信道传输的信号，这就意味着把基带信号（信源）转变为一个相对基带频率而言频率非常高的代通信号。
三、实验环境PC机，Windows2000，office2000，Matlab6.5以上版本软件四、实验过程、数据记录、分析及结论实验过程1.根据参数产生消息信号s和载波信号。
2.编程实现调制过程。
调用函数y=ammod(s,Fc,Fs)完成幅度调制，y=fmmod(s,Fc,Fs,FREQDEV) 完成频率调制，y=pmmod(s,Fc,Fs, PHASEDEV) 完成相位调制。
3.编程实现信号的传输过程。
产生白噪声noise，并将其加到调制信号序列。
4.编程实现信号的解调。
调用函数x=amdemod(y,Fc,Fs)完成幅度调制信号的解调，x=fmdemod(y,Fc,Fs, FREQDEV) 完成频率调制信号的解调，x=pmdemod(y,Fc,Fs, PHASEDEV) 完成相位调制信号的解调。
数据记录消息信号s幅度调制和FFT 消息信号s频率调制和FFT消息信号s相位调制和FFT 幅度调制信号加入白噪声及FFT频率调制信号加入白噪声及FFT 相位调制信号加入白噪声及FFT三种调制方式后经过滤波如上图解调后信号及其FFT解调后信号及其FFT 解调后信号及其FFT五、讨论对于实验中的FREQDEV 、PHASEDEV 含义不清楚，具体数值不知道该如何确定。
通过这次的实验，我们对信息和通信系统有了更进一步的认识，尤其是在系统设计方面，尽管是非常基础的调制与解调的传输，也是经过若干设备协同工作，才能保证信号有效传输，而小到仅仅是一个参数，都有可能导致整个系统无法正常运行。
(代码)基于MATLAB的AM调制解调系统仿真报告
内部基于MATLAB的AM调制解调系统仿真报告XXXX-XXXX-XXXX-XXXXV1.0天津市智能信号与图像处理重点实验室2013年10月29日修订历史记录日期版本文档负责人修改内容2013-10-29 V1.0 刘亚洲创建文档编制姓名签字日期电话审查姓名签字日期电话审核姓名签字日期电话批准姓名签字日期电话文档评审负责人：参加评审人员：目录1引言 (5)1.1设计目的 (5)1.2术语定义 (5)1.3参考资料 (5)1.4文档组织 (5)2 AM调制解调 (6)2.1AM调制 (6)2.2AM解调 (7)3 基于MATLAB的AM仿真 (8)3.1仿真基本参数 (8)3.2生成调制信号 (8)3.3AM调制器 (8)3.4相干解调器 (9)4 仿真结果曲线 (10)4.1发送信号波形和频谱 (10)4.2载波信号波形和频谱 (12)4.3AM信号波形和频谱 (14)4.4相干解调波形和频谱 (16)4.5恢复信号波形和频谱 (18)5总结 (20)6程序附录 (20)1引言1.1设计目的本报告依照传统模拟调制的规范，给出了AM调制解调的具体流程，重点研究了系统中各阶段信号时域和频谱波形以及频谱的搬移变化，为AM调制解调系统信号波形的进一步深入研究做了基础。
1.2术语定义本文档使用以下关键术语和简略语。
英文缩写英文全称中文名称AM Amplitude Modulation 幅度调制AWGN Additive White Gaussian Noise 加性高斯白噪声1.3参考资料[1]通信原理（第六版）樊昌信曹丽娜编著国防工业出版社2007年1月1.4文档组织报告第二部分给出了AM调制解调的基本原理；第三部分给出了系统在MATLAB里面的程序调试及仿真；第四部分给出了各仿真模块输出时域和频域波形，并对比发射信号和接收信号的时域波形；第五部分对报告进行了总结。
2 AM 调制解调信源信号信宿信号AM 调制AM 解调信道加性噪声图1 AM 调制解调系统框图图1显示给出了用于AM 调制解调的系统框图。
实验四--基于matlab的FM调制与解调.wps
FM 解调模型
4.3.2 解调过程分析 输入调频信号为
解调模型
设相干载波为
uFM t Uc cos
0t k f
t 0
u
t
dt
ct cos2 pi fct
乘法器的作用是把调频信号变成有多种频率的波的混合，乘法器输出为
sp
t
1 2
sin
2ct
1 2
K
f
m
t
dt
1
cos
2ct
经低通滤波器后取出器低频分量为
diff_nsfm1(i)=(nsfm1(i+1)-nsfm1(i))./dt;
end
diff_nsfmn=abs(hilbert(diff_nsfm1)); %hilbert 变换，求绝对值得到瞬
时幅度（包络检波）
zero=(max(diff_nsfmn)-min(diff_nsfmn))/2;
subplot(2,1,2);
plot((1:length(diff_nsfmn3))./1000,diff_nsfmn3./400,'r'); xlabel('时间 t'); title('含高斯噪声条件下解调信号的时域图')；
五．实验结果：
调制结果：
高斯白噪声：
解调图形
实验总结：
通过这次的通信原理实验，我对系统的调制与解调有个详细而且深刻的认识，这使得我对通 信原理这门课程有了全新的理解。在实验中遇到了许多问题，通过网上查阅资料，在老师同 学的帮助下，完成了这次的实验。通过个 bask，pcm 等调制方式，我加深了对调制的认识。 并在 fm 解调方式中明白了解调的方式与方法。谢谢老师的悉心指导与同学的帮助。
基于MATLAB的模拟信号频率调制(FM)与解调分析.
课程设计任务书学生姓名：专业班级：电信指导教师：工作单位：武汉理工大学题目：信号分析处理课程设计－基于MATLAB的模拟信号频率调制（FM）与解调分析初始条件：1.Matlab6.5以上版本软件；2.先修课程：通信原理等；要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1、利用MATLAB中的simulink工具箱中的模块进行模拟频率（FM）调制与解调，观察波形变化2、画出程序设计框图，编写程序代码，上机运行调试程序，记录实验结果（含计算结果和图表等），并对实验结果进行分析和总结；3、课程设计说明书按学校统一规范来撰写，具体包括：⑴目录；⑵理论分析；⑶程序设计；⑷程序运行结果及图表分析和总结；⑸课程设计的心得体会（至少800字，必须手写。
）；⑹参考文献（不少于5篇）。
时间安排：周一、周二查阅资料，了解设计内容；周三、周四程序设计，上机调试程序；周五、整理实验结果，撰写课程设计说明书。
指导教师签名： 2013 年 7月 2 日系主任（或责任教师）签名： 2013年 7月 2日目录1 Simulink简介 (1)1.1 Matlab简介······················································错误！未定义书签。
模拟调制解调(matlab)-西工大
实验四 模拟调制解调一、实验目的让学生掌握掌握模拟调制以及对应解调方法的原理。
掌握模拟调制解调方法的计算机编程实现方法，即软件实现。
培养学生综合分析、解决问题的能力，加深对课堂内容的理解。
二、实验要求掌握模拟AM 、PM 和FM 的方法原理以及对应的解调原理；编制调制解调程序；完成对一个正弦信号的调制、传输、滤波、解调过程的仿真；实验后撰写实验报告。
三、实验内容、步骤实验内容已知消息信号为:频率为10f H z =的正弦信号；载波频率为8c f K H z =，采样频率为64KHz 。
编程实现一种调制、传输、滤波和解调过程。
实验步骤1． 根据参数产生消息信号s 和载波信号。
2． 编程实现调制过程。
调用函数y=ammod(s,Fc,Fs)完成幅度调制，y=fmmod(s,Fc,Fs,FREQDEV) 完成频率调制，y=pmmod(s,Fc,Fs, PHASEDEV) 完成相位调制。
3． 编程实现信号的传输过程。
产生白噪声noise ，并将其加到调制信号序列。
4． 编程实现信号的解调。
调用函数x=amdemod(y,Fc,Fs)完成幅度调制信号的解调，x=fmdemod(y,Fc,Fs, FREQDEV) 完成频率调制信号的解调，x=pmdemod(y,Fc,Fs, PHASEDEV) 完成相位调制信号的解调。
四、实验结果及分析（在matlab 下仿真的波形）仿真中，消息信号为:频率为10f H z =的正弦信号；载波频率为8c f K H z =，采样频率为64KHz 。
1. 消息信号和载波信号的波形和频谱如下图所示：00.51-1-0.500.51时间/s 幅度调制信号2.2 2.4 2.6 2.8x 10-3-1-0.500.51时间/s 幅度载波-500501234频率/Hz幅度调制信号频谱-4-2024x 104012344频率/Hz幅度载波频谱2.经过幅度调制、频率调制、相位调制后的时域波形如图所示：0.10.20.30.40.50.60.70.80.91-101时间/s 幅度调幅信号0.0050.010.0150.020.0250.030.035-11时间/s 幅度调频信号0.0050.010.0150.020.0250.030.0350.040.0450.05时间/s幅度调相信号3. 经过幅度调制、频率调制、相位调制后的频谱如图所示：-4-3-2-101234x 104012x 104频率/Hz 幅度调幅信号频谱-4-3-2-101234x 10405000频率/Hz 幅度调频信号频谱-4-3-2-101234x 104012x 104频率/Hz幅度调相信号频谱4.解调之后的时域波形如图所示：此时解调信号与原信号相比，恢复效果较好。
基于matlab的 AM,FM调制与解调报告
AM调制与解调100%% AMµ÷ÖÆfigure('Name','Ðźŵ÷Öƹý³ÌÖв¨Ðμ°ÆäƵÆ×','NumberTitle','off')a0=1;f0=2000;fc=20000;fs=1000000;t=[1:0.000001:4];am1=0*cos(2*pi*f0*t); %µ÷ÖÆÐźÅam=a0+am1;t1=cos(2*pi*fc*t); %Ôز¨s_am=am.*t1;AM1=fft(am1); T1=fft(t1); S_AM=fft(s_am);f=(0:3000000)*fs/3000001-fs/2;subplot(3,2,1); plot(t(1:1000),am1(1:1000)); title('ÐÅÏ¢ÐźŲ¨ÐÎ');subplot(3,2,2); plot(f,fftshift(abs(AM1))); title('ÐÅÏ¢ÐźÅƵÆ×');subplot(3,2,3); plot(t(1:1000),t1(1:1000)); title('Ôز¨ÐźÅ');subplot(3,2,4); plot(f,fftshift(abs(T1))); title('Ôز¨ÐźÅƵÆ×');subplot(3,2,5); plot(t(1:1000),s_am(1:1000)); title('Òѵ÷ÐźÅ');subplot(3,2,6); plot(f,fftshift(abs(S_AM))); title('Òѵ÷ÐźÅƵÆ×');%²úÉúÔëÉùfigure('Name','Ìí¼ÓÔëÉù¼°´øͨÂ˲¨¹ý³Ì²¨Ðμ°ÆäƵÆ×','NumberTitle','off ');snr=5;y=awgn(s_am,snr);fcuts=[16000 17500 22500 24000];mags=[0 1 0];devs=[0.05 0.01 0.05];[n,Wn,beta,ftype]=kaiserord(fcuts,mags,devs,fs);hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');[H,f]=freqz(hh,1,1024,fs);st_p=fftfilt(hh,y);Q=fft(st_p);Y=fft(y);f=(0:3000000)*fs/3000001-fs/2;subplot(2,2,1);plot(t(1500001:1501000),y(1500001:1501000));title('Ìí¼ÓÔëÉùºóÐźŲ¨ÐÎ');subplot(2,2,2);plot(f,fftshift(abs(Y)));title('Ìí¼ÓÔëÉùºóÐźÅƵÆ×');subplot(2,2,3);plot(t(1500001:1501000),st_p(1500001:1501000));title('´øͨÂ˲¨ºóÐźŲ¨ÐÎ');subplot(2,2,4); plot(f,fftshift(abs(Q)));title('´øͨÂ˲¨ºóÐźÅƵÆ×');%½âµ÷figure('Name','Ïà¸É½âµ÷ËùµÃ²¨Ðμ°ÆäƵÆ×','NumberTitle','off');ss_am=st_p.*t1;SS_AM=fft(ss_am)f=(0:3000000)*fs/3000001-fs/2;subplot(2,1,1);plot(t(1500001:1503000),ss_am(1500001:1503000));title( 'Ïà³ËÐźÅ');subplot(2,1,2); plot(f,fftshift(abs(SS_AM)));title('Ïà³ËÐźÅƵÆ×');fcuts1=[2500,30000];mags1=[1 0];devs1=[0.01 0.05];[n,Wn,beta,ftype]=kaiserord(fcuts1,mags1,devs1,fs);hh1=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');[H,f]=freqz(hh1,1,1024,fs);m0=fftfilt(hh1,ss_am);M0=fft(m0);f=(0:3000000)*fs/3000001-fs/2;figuresubplot(2,1,1);plot(t(1500001:1501000),m0(1500001:1501000));title('½âµ÷ÐźÅ');subplot(2,1,2); plot(f,fftshift(abs(M0)));title('½âµ÷ÐźÅƵÆ×');11.0005 1.001-101信息信号波形-505x 105012x 106信息信号频谱11.0005 1.001-101载波信号-505x 105012x 106载波信号频谱11.00051.001-202已调信号-55x 105012x 106已调信号频谱2.52.50052.501-4-2024添加噪声后信号波形-505x 105051015x 105添加噪声后信号频谱2.52.50052.501-2024带通滤波后信号波形-55x 10500.511.526带通滤波后信号频谱50%% AMµ÷ÖÆfigure('Name','Ðźŵ÷Öƹý³ÌÖв¨Ðμ°ÆäƵÆ×','NumberTitle','off')2.5 2.5005 2.501 2.5015 2.502 2.5025 2.503 2.5035-1012相乘信号-5-4-3-2-1012345x 1050510155相乘信号频谱2.52.50012.50022.50032.50042.50052.50062.50072.50082.5009 2.501-0.500.51解调信号-5-4-3-2-1012345x 1050510155解调信号频谱a0=2;f0=2000;fc=20000;fs=1000000;t=[1:0.000001:4];am1=0*cos(2*pi*f0*t); %µ÷ÖÆÐźÅam=a0+am1;t1=cos(2*pi*fc*t); %Ôز¨s_am=am.*t1;AM1=fft(am1); T1=fft(t1); S_AM=fft(s_am);f=(0:3000000)*fs/3000001-fs/2;subplot(3,2,1); plot(t(1:1000),am1(1:1000)); title('ÐÅÏ¢ÐźŲ¨ÐÎ');subplot(3,2,2); plot(f,fftshift(abs(AM1))); title('ÐÅÏ¢ÐźÅƵÆ×');subplot(3,2,3); plot(t(1:1000),t1(1:1000)); title('Ôز¨ÐźÅ');subplot(3,2,4); plot(f,fftshift(abs(T1))); title('Ôز¨ÐźÅƵÆ×');subplot(3,2,5); plot(t(1:1000),s_am(1:1000)); title('Òѵ÷ÐźÅ');subplot(3,2,6); plot(f,fftshift(abs(S_AM))); title('Òѵ÷ÐźÅƵÆ×');%²úÉúÔëÉùfigure('Name','Ìí¼ÓÔëÉù¼°´øͨÂ˲¨¹ý³Ì²¨Ðμ°ÆäƵÆ×','NumberTitle','off ');snr=5;y=awgn(s_am,snr);fcuts=[16000 17500 22500 24000];mags=[0 1 0];devs=[0.05 0.01 0.05];[n,Wn,beta,ftype]=kaiserord(fcuts,mags,devs,fs);hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');[H,f]=freqz(hh,1,1024,fs);st_p=fftfilt(hh,y);Q=fft(st_p);Y=fft(y);f=(0:3000000)*fs/3000001-fs/2;subplot(2,2,1);plot(t(1500001:1501000),y(1500001:1501000));title('Ìí¼ÓÔëÉùºóÐźŲ¨ÐÎ');subplot(2,2,2);plot(f,fftshift(abs(Y)));title('Ìí¼ÓÔëÉùºóÐźÅƵÆ×');subplot(2,2,3);plot(t(1500001:1501000),st_p(1500001:1501000));title('´øͨÂ˲¨ºóÐźŲ¨ÐÎ');subplot(2,2,4); plot(f,fftshift(abs(Q)));title('´øͨÂ˲¨ºóÐźÅƵÆ×');%½âµ÷figure('Name','Ïà¸É½âµ÷ËùµÃ²¨Ðμ°ÆäƵÆ×','NumberTitle','off');ss_am=st_p.*t1;SS_AM=fft(ss_am)f=(0:3000000)*fs/3000001-fs/2;subplot(2,1,1);plot(t(1500001:1503000),ss_am(1500001:1503000));title( 'Ïà³ËÐźÅ');subplot(2,1,2); plot(f,fftshift(abs(SS_AM)));title('Ïà³ËÐźÅƵÆ×'); fcuts1=[2500,30000];mags1=[1 0];devs1=[0.01 0.05];[n,Wn,beta,ftype]=kaiserord(fcuts1,mags1,devs1,fs);hh1=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');[H,f]=freqz(hh1,1,1024,fs);m0=fftfilt(hh1,ss_am);M0=fft(m0);f=(0:3000000)*fs/3000001-fs/2;figuresubplot(2,1,1);plot(t(1500001:1501000),m0(1500001:1501000));title('½âµ÷ÐźÅ');subplot(2,1,2); plot(f,fftshift(abs(M0)));title('½âµ÷ÐźÅƵÆ×');11.0005 1.001-101信息信号波形-505x 1050126信息信号频谱11.0005 1.001-101载波信号-505x 1050126载波信号频谱11.00051.001-505已调信号-55x 105024x 106已调信号频谱2.52.50052.501-505添加噪声后信号波形-505x 1050123x 106添加噪声后信号频谱2.52.50052.501-4-2024带通滤波后信号波形-55x 105012346带通滤波后信号频谱2.5 2.5005 2.501 2.5015 2.502 2.5025 2.503 2.5035-10123相乘信号-5-4-3-2-1012345x 10501236相乘信号频谱0%% AMµ÷ÖÆfigure('Name','Ðźŵ÷Öƹý³ÌÖв¨Ðμ°ÆäƵÆ×','NumberTitle','off')a0=10^100;f0=2000;fc=20000;fs=1000000;t=[1:0.000001:4];am1=0*cos(2*pi*f0*t); %µ÷ÖÆÐźÅam=a0+am1;t1=cos(2*pi*fc*t); %Ôز¨s_am=am.*t1;AM1=fft(am1); T1=fft(t1); S_AM=fft(s_am);f=(0:3000000)*fs/3000001-fs/2;subplot(3,2,1); plot(t(1:1000),am1(1:1000)); title('ÐÅÏ¢ÐźŲ¨ÐÎ');subplot(3,2,2); plot(f,fftshift(abs(AM1))); title('ÐÅÏ¢ÐźÅƵÆ×');2.52.50012.50022.50032.50042.50052.50062.50072.50082.5009 2.50100.511.5解调信号-5-4-3-2-1012345x 10501236解调信号频谱subplot(3,2,3); plot(t(1:1000),t1(1:1000)); title('Ôز¨ÐźÅ');subplot(3,2,4); plot(f,fftshift(abs(T1))); title('Ôز¨ÐźÅƵÆ×');subplot(3,2,5); plot(t(1:1000),s_am(1:1000)); title('Òѵ÷ÐźÅ');subplot(3,2,6); plot(f,fftshift(abs(S_AM))); title('Òѵ÷ÐźÅƵÆ×');%²úÉúÔëÉùfigure('Name','Ìí¼ÓÔëÉù¼°´øͨÂ˲¨¹ý³Ì²¨Ðμ°ÆäƵÆ×','NumberTitle','off ');snr=5;y=awgn(s_am,snr);fcuts=[16000 17500 22500 24000];mags=[0 1 0];devs=[0.05 0.01 0.05];[n,Wn,beta,ftype]=kaiserord(fcuts,mags,devs,fs);hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');[H,f]=freqz(hh,1,1024,fs);st_p=fftfilt(hh,y);Q=fft(st_p);Y=fft(y);f=(0:3000000)*fs/3000001-fs/2;subplot(2,2,1);plot(t(1500001:1501000),y(1500001:1501000));title('Ìí¼ÓÔëÉùºóÐźŲ¨ÐÎ');subplot(2,2,2);plot(f,fftshift(abs(Y)));title('Ìí¼ÓÔëÉùºóÐźÅƵÆ×');subplot(2,2,3);plot(t(1500001:1501000),st_p(1500001:1501000));title('´øͨÂ˲¨ºóÐźŲ¨ÐÎ');subplot(2,2,4); plot(f,fftshift(abs(Q)));title('´øͨÂ˲¨ºóÐźÅƵÆ×');%½âµ÷figure('Name','Ïà¸É½âµ÷ËùµÃ²¨Ðμ°ÆäƵÆ×','NumberTitle','off');ss_am=st_p.*t1;SS_AM=fft(ss_am)f=(0:3000000)*fs/3000001-fs/2;subplot(2,1,1);plot(t(1500001:1503000),ss_am(1500001:1503000));title( 'Ïà³ËÐźÅ');subplot(2,1,2); plot(f,fftshift(abs(SS_AM)));title('Ïà³ËÐźÅƵÆ×'); fcuts1=[2500,30000];mags1=[1 0];devs1=[0.01 0.05];[n,Wn,beta,ftype]=kaiserord(fcuts1,mags1,devs1,fs);hh1=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');[H,f]=freqz(hh1,1,1024,fs);m0=fftfilt(hh1,ss_am);M0=fft(m0);f=(0:3000000)*fs/3000001-fs/2;figuresubplot(2,1,1);plot(t(1500001:1501000),m0(1500001:1501000));title('½âµ÷ÐźÅ');subplot(2,1,2); plot(f,fftshift(abs(M0)));title('½âµ÷ÐźÅƵÆ×');11.0005 1.001-101信息信号波形-505x 1050126信息信号频谱1 1.0005 1.001-101载波信号-505x 1050126载波信号频谱11.00051.001-101x 10100已调信号-55x 105012x 10106已调信号频谱2.52.5005 2.501-1-0.500.51x 10100添加噪声后信号波形-505x 105051015x 10105添加噪声后信号频谱2.52.5005 2.501-2-1012x 10100带通滤波后信号波形-55x 1051x 10106带通滤波后信号频谱FM 调制与解调%%FMfigure('Name','FMµ÷ÖƲ¨ÐÎÓëƵÆ×')2.5 2.5005 2.501 2.5015 2.502 2.5025 2.503 2.5035-50510x 1099相乘信号-5-4-3-2-1012345x 105051015105相乘信号频谱2.52.50012.50022.50032.50042.50052.50062.50072.50082.5009 2.5013.93.913.923.93x 1099解调信号-1-0.500.51x 105123105解调信号频谱f0=2000; fc=20000; fs=1000000; Am=1; kf=0.5; Tc=8; Ta=0.001; dt=0.000001;t=[0:1/fs:3];f=(0:length(t)-1)*fs/(length(t))-fs/2;fm0=cos(2*pi*f0*t);mt=fm0;%»ý·ÖÆ÷Éè¼Æw1=0;w2=0;for m=1:length(t)w1=mt(m)+w2;w2=mt(m)+w1;fi(m)=w1/(2*fs);endfi=fi*2*pi/max(abs(fi));I=cos(kf*fi);Q=sin(kf*fi);y1=Am*cos(2*pi*fc*t).*I-Am*sin(2*pi*fc*t).*Q;subplot(2,1,1);plot(t,y1);title('²¨ÐÎ')axis([1e-3 4e-3 -2 2]);Y1=fft(y1);subplot(2,1,2);plot(f,fftshift(abs(Y1))/1e6); title('ƵÆ×') %%Ôز¨ÆµÆ×axis([-4e4 4e4 0 1]);figure('Name','FMµ÷Öƺó¼ÓÔëÉù²¨ÐÎÓë½âµ÷ºó²¨ÐÎÒÔ¼°Â˳ýÖ±Á÷·ÖÁ¿ºóµÄ²¨ÐÎ')y1o=awgn(y1,40);subplot(3,1,1);plot(t,y1o); title('¼ÓÔëÉùºó²¨ÐÎ') %%¼ÓÔëÉùºóµÄÐźÅaxis([1e-3 4e-3 -2 2]);%%´øͨÂ˲¨KSband=2*(3+1)*f0;fcutsb=[fc-KSband-2000 fc-KSbandfc+KSbandfc+KSband+2000]; %%½ÓÊÕ»úÇ°¶Ë´øͨÂ˲¨magsb=[0 1 0];devsb=[0.05 0.01 0.05];[nb,Wnb,betab,ftypeb]=kaiserord(fcutsb,magsb,devsb,fs);hhb=fir1(nb,Wnb,ftypeb,kaiser(nb+1,betab),'noscale'); %´øͨÂ˲¨Æ÷£»st_pb=fftfilt(hhb,y1o);subplot(3,1,2);st_pb=st_pb/1e6;plot(t,st_pb); title('´øͨÂ˲¨Æ÷ºóµÄ²¨ÐÎ')axis([1e-3 4e-3 -2e-6 2e-6]);%΢·ÖÆ÷Éè¼Æfor i=1:length(t)-1 %½ÓÊÕÐźÅͨ¹ý΢·ÖÆ÷´¦Àídiff_st_pb(i)=(st_pb(i+1)-st_pb(i))/dt;endsfm=abs(hilbert(diff_st_pb));subplot(3,1,2);plot(t,[sfm*20 0]);axis([1e-3 4e-3 0 4]);%%¸ôÖ±% KSbandh=2*(3+1)*f0;fcutsh=[0.01 3000];magsh=[0 1];devsh=[0.01 0.05];[nh,Wnh,betah,ftypeh]=kaiserord(fcutsh,magsh,devsh,fs);hhh=fir1(nh,Wnh,ftypeh,kaiser(nh+1,betah),'noscale');sfm_out=fftfilt(hhh,sfm*20);subplot(3,1,3);plot(t,[sfm_out 0]);title('¸ôÖ±ºóµÄ²¨ÐÎ')axis([1e-3 4e-3 -2 2]);11.522.533.54x 10-3-2-1012波形-4-3-2-101234x 10400.51频谱11.522.533.54x 10-3-202加噪声后波形11.522.53 3.54x 10-302411.522.533.54x 10-3-202隔直后的波形。
基于Matlab的模拟信号传输系统实验(幅度调制)
1）、实验原理及框图调制：DSB信号实质上就是基带信号与载波直接相乘。
设均值为零的模拟基带信号为m(t)，双边带调幅（DSB）信号为:s(t)=m(t)cosωt其原理框图如下：解调：由于m(t)均值为0，因此调制后的信号不含离散的载波分量，若接收端能恢复出载波分量，则可采用相干解调再用低通滤波器滤去高频分量,即可无失真地恢复出原始信号。
其原理框图如下：2）、实验结果：附：程序源代码fs=1e6;%采样率N=8192;t=0:1/fs:(N-1)/fs;%时域采样点f=(0:(N-1))*fs/N-fs/2;%频域采样点fm=2000;%基带频率fss=20000;%载波频率m=cos(fm*2*pi*t);%基带信号subplot(4,1,1);plot(t,m);title('基带信号','FontWeight','bold');xlabel('t/s','FontSize',12);axis([0.001,0.004,-2,2]);m_w=fft(m);%基带频谱subplot(4,1,2);plot(f,fftshift(abs(m_w)));title('基带频谱','FontWeight','bold');xlabel('f/Hz','FontSize',12);axis([-10000,10000,0,3500]);s=cos(fss*2*pi*t);%载波信号subplot(4,1,3);plot(t,s);title('载波信号','FontWeight','bold');xlabel('t/s','FontSize',12);axis([0.001,0.004,-2,2]);s_w=fft(s);%载波频谱subplot(4,1,4);plot(f,fftshift(abs(s_w)));title('载波频谱','FontWeight','bold');xlabel('f/Hz','FontSize',12);axis([-30000,30000,0,3500]);sm=s.*m;%调制信号figure;subplot(4,1,1);plot(t,sm);title('调制信号','FontWeight','bold');xlabel('t/s','FontSize',12);axis([0.001,0.004,-2,2]);sm_w=fft(sm);%调制信号频谱subplot(4,1,2);plot(f,fftshift(abs(sm_w)));title('调制信号频谱','FontWeight','bold');xlabel('f/Hz','FontSize',12);axis([-30000,30000,0,2000]);sm1=awgn(sm,5);subplot(4,1,3);plot(t,sm1);title('加噪声后的调制信号','FontWeight','bold');xlabel('t/s','FontSize',12);axis([0.001,0.004,-2,2]);sm1_w=fft(sm1);subplot(4,1,4);plot(f,fftshift(abs(sm1_w)));title('加噪声后的调制信号频谱','FontWeight','bold');xlabel('f/Hz','FontSize',12);axis([-30000,30000,0,2000]);fsamp =1e6;fcuts =[16000175002250024000];mags =[010];devs =[0.050.010.05];[n,Wn,beta,ftype]=kaiserord(fcuts,mags,devs,fsamp);hh =fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');sm2=fftfilt(hh,sm1);figure;subplot(4,1,1);plot(t,sm2);title('滤掉部分噪声后的调制信号','FontWeight','bold');xlabel('t/s','FontSize',12);axis([0.001,0.004,-2,2]);sm2_w=fft(sm2);subplot(4,1,2);plot(f,fftshift(abs(sm2_w)));title('滤掉部分噪声后的调制信号频谱','FontWeight','bold');xlabel('f/Hz','FontSize',12);axis([-30000,30000,0,2000]);sp=2*sm2.*s;subplot(4,1,3);plot(t,sp);title('与本地载波相乘后的信号','FontWeight','bold');xlabel('t/s','FontSize',12);axis([0.001,0.004,-2,2]);sp_w=fft(sp);subplot(4,1,4);plot(f,fftshift(abs(sp_w)));title('与本地载波相乘后的信号频谱','FontWeight','bold');xlabel('f/Hz','FontSize',12);axis([-45000,45000,0,1200]);fsamp =1e6;fcuts =[300020000];mags =[10];devs =[0.010.05];[n,Wn,beta,ftype]=kaiserord(fcuts,mags,devs,fsamp);hh1=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');sd=fftfilt(hh1,sp);figure;subplot(4,1,1);plot(t,sd);title('低通滤波后的信号','FontWeight','bold');xlabel('t/s','FontSize',12);axis([0.001,0.004,-2,2]);sd_f=fft(sd);subplot(4,1,2);plot(f,fftshift(abs(sd_f)));title('低通滤波后的信号频谱','FontWeight','bold');xlabel('f/Hz','FontSize',12);axis([-10000,10000,0,3500]);。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
基于Matlab的模拟调制与解调（开放实验）一、实验目的（一）了解AM、DSB和SSB 三种模拟调制与解调的基本原理（二）掌握使用Matlab进行AM调制解调的方法1、学会运用MATLAB对基带信号进行AM调制2、学会运用MATLAB对AM调制信号进行相干解调3、学会运用MATLAB对AM调制信号进行非相干解调（包络检波）（三）掌握使用Matlab进行DSB调制解调的方法1、学会运用MATLAB对基带信号进行DSB调制2、学会运用MATLAB对DSB调制信号进行相干解调（四）掌握使用Matlab进行SSB调制解调的方法1、学会运用MATLAB对基带信号进行上边带和下边带调制2、学会运用MATLAB对SSB调制信号进行相干解调二、实验环境MatlabR2020a三、实验原理（一）滤波法幅度调制（线性调制）（二）常规调幅（AM）1、AM表达式2、AM波形和频谱3、调幅系数m（三）抑制载波双边带调制（DSB-SC）1、DSB表达式2、DSB波形和频谱（四）单边带调制（SSB）（五）相关解调与包络检波四、实验过程（一）熟悉相关内容原理 （二）完成作业已知基带信号()()()sin 10sin 30m t t t ππ=+，载波为()()cos 2000c t t π= 1、对该基带信号进行AM 调制解调（1）写出AM 信号表达式，编写Matlab 代码实现对基带进行进行AM 调制，并分别作出3种调幅系数（1,1,1m m m >=<）下的AM 信号的时域波形和幅度频谱图。
代码 基带信号fs = 10000; % 采样频率 Ts = 1/fs; % 采样时间间隔t = 0:Ts:1-Ts; % 时间向量m = sin(10*pi*t) + sin(30*pi*t); % 基带信号载波信号fc = 1000; % 载波频率c = cos(2*pi*fc*t); % 载波信号AM调制Ka = [1, 0.5, 2]; % 调制系数m_AM = zeros(length(Ka), length(t)); % 存储AM调制信号相干解调信号r = zeros(length(Ka), length(t));绘制AM调制信号的时域波形和幅度频谱图figure;for i = 1:length(Ka)m_AM(i, :) = (1 + Ka(i)*m).*c; % AM调制信号subplot(3, 2, i);plot(t, m_AM(i, :));title(['AM调制信号(Ka = ' num2str(Ka(i)) ')']);xlabel('时间');ylabel('幅度');ylim([-2, 2]);subplot(3, 2, i+3);f = (-fs/2):fs/length(m_AM(i, :)):(fs/2)-fs/length(m_AM(i, :));M_AM = fftshift(abs(fft(m_AM(i, :))));plot(f, M_AM);title(['AM调制信号的幅度频谱图(Ka = ' num2str(Ka(i)) ')']);xlabel('频率');ylabel('幅度');r(i, :) = m_AM(i, :) .* c; % 相干解调信号end绘制相干解调信号的时域波形和幅度频谱图figure;for i = 1:length(Ka)subplot(length(Ka), 1, i);plot(t, r(i, :));title(['相干解调信号(Ka = ' num2str(Ka(i)) ')']);xlabel('时间');ylabel('幅度');end图像（2）编写Matlab代码实现对AM调制信号的相干解调，并作出图形。
（3*）编写Matlab代码实现对AM调制信号的包络检波，并作出图形。
代码定义时间和信号t = 0:0.0001:0.1;m = sin(10*pi*t) + sin(30*pi*t);c = cos(2000*pi*t);计算AM调制信号s = (1+m).*c;包络检波env = envelope(s);绘图figure;subplot(3,1,1);plot(t, m);title('基带信号');xlabel('时间（秒）');ylabel('幅度');subplot(3,1,2);plot(t, c);title('载波信号');xlabel('时间（秒）');ylabel('幅度');subplot(3,1,3);plot(t, s);hold on;plot(t, env, 'r');title('AM调制信号及其包络'); xlabel('时间（秒）');ylabel('幅度');legend('AM调制信号', '包络'); 图像2、对该基带信号进行DSB调制解调（1）写出DSB信号表达式，编写Matlab代码实现对基带进行DSB调制，并分别作出DSB信号的时域波形和幅度频谱图。
代码基带信号fs = 10000; % 采样频率Ts = 1/fs; % 采样时间间隔t = 0:Ts:1-Ts; % 时间向量m = sin(10*pi*t) + sin(30*pi*t); % 基带信号载波信号fc = 1000; % 载波频率c = cos(2*pi*fc*t); % 载波信号DSB调制m_DSB = m .* c; % DSB调制信号% 绘制DSB调制信号的时域波形和幅度频谱图figure;subplot(2, 1, 1);plot(t, m_DSB);title('DSB调制信号');xlabel('时间');ylabel('幅度');subplot(2, 1, 2);f = (-fs/2):fs/length(m_DSB):(fs/2)-fs/length(m_DSB); M_DSB = fftshift(abs(fft(m_DSB)));plot(f, M_DSB);title('DSB调制信号的幅度频谱图');xlabel('频率');ylabel('幅度');相干解调r = m_DSB .* c; % 相干解调信号绘制相干解调信号的时域波形和幅度频谱图figure;subplot(2, 1, 1);plot(t, r);title('相干解调信号');xlabel('时间');ylabel('幅度');subplot(2, 1, 2);R = fftshift(abs(fft(r)));plot(f, R);title('相干解调信号的幅度频谱图');xlabel('频率');ylabel('幅度');图像（2）编写Matlab代码实现对DSB调制信号的相干解调，并作出图形。
3、对该基带信号进行SSB调制解调（1）写出SSB信号表达式，编写Matlab代码实现对基带进行SSB调制，并分别作出上边带调制和下边带调制信号的时域波形和幅度频谱图。
基带信号fs = 10000; % 采样频率Ts = 1/fs; % 采样时间间隔t = 0:Ts:1-Ts; % 时间向量m = sin(10*pi*t) + sin(30*pi*t); % 基带信号载波信号fc = 1000; % 载波频率c = cos(2*pi*fc*t); % 载波信号SSB调制(上边带调制)m_H = hilbert(m); % m(t) 的希尔伯特变换m_SSB_upper = m .* cos(2*pi*fc*t) - real(m_H) .* sin(2*pi*fc*t);绘制上边带调制信号的时域波形和幅度频谱图figure;subplot(2, 1, 1);plot(t, m_SSB_upper);title('上边带调制信号');xlabel('时间');ylabel('幅度');subplot(2, 1, 2);f = (-fs/2):fs/length(m_SSB_upper):(fs/2)-fs/length(m_SSB_upper); M_SSB_upper = fftshift(abs(fft(m_SSB_upper)));plot(f, M_SSB_upper);title('上边带调制信号的幅度频谱图');xlabel('频率');ylabel('幅度');SSB调制(下边带调制)m_SSB_lower = m .* cos(2*pi*fc*t) + imag(m_H) .* sin(2*pi*fc*t);绘制下边带调制信号的时域波形和幅度频谱图figure;subplot(2, 1, 1);plot(t, m_SSB_lower);title('下边带调制信号');xlabel('时间');ylabel('幅度');subplot(2, 1, 2);M_SSB_lower = fftshift(abs(fft(m_SSB_lower)));plot(f, M_SSB_lower);title('下边带调制信号的幅度频谱图');xlabel('频率');ylabel('幅度');相干解调(上边带调制)r_upper = m_SSB_upper .* cos(2*pi*fc*t);绘制上边带相干解调信号的时域波形和幅度频谱图figure;subplot(2, 1, 1);plot(t, r_upper);title('上边带相干解调信号');xlabel('时间');ylabel('幅度');subplot(2, 1, 2);R_upper = fftshift(abs(fft(r_upper)));plot(f, R_upper);title('上边带相干解调信号的幅度频谱图');xlabel('频率');ylabel('幅度');相干解调(下边带调制)r_lower = m_SSB_lower .* cos(2*pi*fc*t);绘制下边带相干解调信号的时域波形和幅度频谱图figure;subplot(2, 1, 1);plot(t, r_lower);title('下边带相干解调信号');xlabel('时间');ylabel('幅度');subplot(2, 1, 2);R_lower = fftshift(abs(fft(r_lower)));plot(f, R_lower);title('下边带相干解调信号的幅度频谱图'); xlabel('频率');ylabel('幅度');图像（2）编写Matlab代码实现对SSB调制信号的相干解调，并作出图形。