MATLAB FM仿真报告
基于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的FM仿真实现
3.3 频率调制（FM ）3.3.1 FM 调制和解调的基本原理频率调制是利用载波的频率变化来传递模拟信息，而振幅保持不变。
也就是说，载波信号的频率随着基带调制信号的幅度变化而改变。
调制信号幅度变大（或变小）时，载波信号的频率也变大（或变小），调制信号幅度变小时，载波信号的频率也变小（或变大）。
在FM 中，FM 信号的瞬时频偏与调制信号m(t)成正比。
因此FM 的信号的时域表达式为：(2.1)式中：A 为载波的恒定振幅；[ωc t+φ(t)]为信号的瞬时相位，记为θ（t ）; φ(t)为相对于载波相位ωc t 的瞬时相位偏移；d[ωc t+φ(t)]/dt 是信号的瞬时角频率，记为ω(t)；而d φ(t)/dt 称为相对于载频ωc 的瞬时频偏。
所谓频率调制（FM ）,是指瞬时频率偏移随调制信号m(t)成比例变化，即(2.2)式中：K f 为调频灵敏度（rad/(s.V)）。
这时相位偏移为：(2.3)因此，上式可改写为(2.4)图2.1 无噪声调制信号FM 调制的实现调频主要有两种方法：直接调频和间接调频。
1）直接调频法调频就是用调制信号控制载波的频率变化。
直接调频就是用调制信号直接去控制载波振荡器的频率，使其按调制信号的规律线性的变化。
()]cos[)(⎰+=ττωd m K t A t s f c FM )()(t m K dt t d f =Φ⎰=Φτd t m K t f )(）（)](cos[)(t t A t S c FM Φ+=ω可以由外部电压控制震荡频率的振荡器叫做压控振荡器器。
每个压控振荡器自身就是一个FM 调制器，因为它的振荡频率正比于输入控制电压，即(2.9) 若用调制信号作控制电压信号，就能产生FM 波。
若被控制的振荡器是LC 振荡器，则只需控制振荡回路的某个电抗元件（L 或C ) ，使其参数随调制信号变化。
目前常用的电抗元件是变容二极管。
用变容二极管实现直接调频，由于电路简单，性能良好，已成为目前最广泛采用的调频电路之一。
基于MATLAB的FM通信系统仿真
文档从网络中收集，已重新整理排版.word版本可编辑.欢迎下载支持.西南科技大学课程设计报告课程名称：通信系统课程设计设计名称：基于MATLAB的FM通信系统仿真姓名：学号:班级：通信0901指导教师：起止日期：2012.6.4-2011.6.24西南科技大学信息工程学院制课程设计任务书学生班级：通信0901 学生姓名：学号：20096009 设计名称：基于MATLAB的FM通信系统仿真起止日期：2012.06.04——2012.06.24 指导教师：杨俊平课程设计学生日志课程设计评语表基于MATLAB 的FM 通信系统仿真一、设计目的和意义1、熟悉MATLAB 的使用方法，了解简单函数、原理和掌握操作方法；2、加深对FM 信号调制及解调原理的理解；3、了解FM 调制、解调基本结论、规律；4、增加自己在通信原理方面的动手能力；5、通过查资料等提高自学能力，对通信原理在最后有更深入的了解。
二、设计原理 1、FM 调制原理如果幅度不变，起始相位为0，而瞬时角频率是调制信号的线性函数，则这种调制方式称为频率调制。
调频信号可表达为])(cos[)(⎰+=dt t f K t w A t S FM c FM （1）调制信号为单余弦时t w A t f m m cos )(= （2） 当对载波进行频率调制时，可得调频信号表达式为]cos cos[)(⎰+=tdt w A K t w A t S m m FM c FM （3）]sin cos[t w t w A m FM c β+= （4） 上式中m m FM FM w A K /=β称为调频指数，也可记作]sin cos[)(maxt w w w t w A t S m m c FM ∆+= （5）图1()()()FM m t dtPM S t →•→→⎰ 图2图1为直接调频，图2为间接调频。
由于实际相位调制器的调节范围不可能超出，因而间接调频的方法仅适合于相位偏移和频率偏移不大的窄调制情形，而直接调频则常用于宽带调制情形。
Matlab仿真实验报告_发布
2、实验目的
通过本次实验，应该掌握： (1) 用傅立叶变换进行信号分析时基本参数的选择。 (2) 经过离散时间傅立叶变换（DTFT）和有限长度离散傅立叶变换（DFT） 后信号频 谱上的区别，前者 DTFT 时间域是离散信号，频率域还是连续的，而 DFT 在两个域中都 是离散的。 (3) 离散傅立叶变换的基本原理、特性，以及经典的快速算法（基 2 时间抽选法） ，体 会快速算法的效率。 (4) 获得一个高密度频谱和高分辨率频谱的概念和方法，建立频率分辨率和时间分辨 率的概念，为将来进一步进行时频分析（例如小波）的学习和研究打下基础。 (5) 建立 DFT 从整体上可看成是由窄带相邻滤波器组成的滤波器组的概念，此概念 的一个典型应用是数字音频压缩中的分析滤波器，例如 DVD AC3 和 MPEG Audio。
~4~
MATLAB 仿真实验报告
N=4000; % N>1/TΔf=1/357*T
n=0:1:N-1; f=fs*n/N;
%DTMF 信号编码 dtmf(1,:)=sin(2*pi*row(1)*T*n)+sin(2*pi*col(1)*T*n); dtmf(2,:)=sin(2*pi*row(1)*T*n)+sin(2*pi*col(2)*T*n); dtmf(3,:)=sin(2*pi*row(1)*T*n)+sin(2*pi*col(3)*T*n); dtmf(4,:)=sin(2*pi*row(2)*T*n)+sin(2*pi*col(1)*T*n); dtmf(5,:)=sin(2*pi*row(2)*T*n)+sin(2*pi*col(2)*T*n); dtmf(6,:)=sin(2*pi*row(2)*T*n)+sin(2*pi*col(3)*T*n); dtmf(7,:)=sin(2*pi*row(3)*T*n)+sin(2*pi*col(1)*T*n); dtmf(8,:)=sin(2*pi*row(3)*T*n)+sin(2*pi*col(2)*T*n); dtmf(9,:)=sin(2*pi*row(3)*T*n)+sin(2*pi*col(3)*T*n); dtmf(10,:)=sin(2*pi*row(4)*T*n)+sin(2*pi*col(2)*T*n); %1 %2 %3 %4 %5 %6 %7 %8 %9 %0
matlab仿真实验报告
matlab仿真实验报告Matlab仿真实验报告引言：Matlab是一种广泛应用于科学和工程领域的数值计算软件，它提供了强大的数学和图形处理功能，可用于解决各种实际问题。
本文将通过一个具体的Matlab 仿真实验来展示其在工程领域中的应用。
实验背景：本次实验的目标是通过Matlab仿真分析一个电路的性能。
该电路是一个简单的放大器电路，由一个输入电阻、一个输出电阻和一个放大倍数组成。
我们将通过Matlab对该电路进行仿真，以了解其放大性能。
实验步骤：1. 定义电路参数：首先，我们需要定义电路的各个参数，包括输入电阻、输出电阻和放大倍数。
这些参数将作为Matlab仿真的输入。
2. 构建电路模型：接下来，我们需要在Matlab中构建电路模型。
可以使用电路元件的模型来表示电路的行为，并使用Matlab的电路分析工具进行仿真。
3. 仿真分析：在电路模型构建完成后，我们可以通过Matlab进行仿真分析。
可以通过输入不同的信号波形，观察电路的输出响应，并计算放大倍数。
4. 结果可视化：为了更直观地观察仿真结果，我们可以使用Matlab的图形处理功能将仿真结果可视化。
可以绘制输入信号波形、输出信号波形和放大倍数的变化曲线图。
实验结果：通过仿真分析，我们得到了以下实验结果：1. 输入信号波形与输出信号波形的对比图：通过绘制输入信号波形和输出信号波形的变化曲线，我们可以观察到电路的放大效果。
可以看到输出信号的幅度大于输入信号，说明电路具有放大功能。
2. 放大倍数的计算结果：通过对输出信号和输入信号的幅度进行计算，我们可以得到电路的放大倍数。
通过比较不同输入信号幅度下的输出信号幅度，可以得到放大倍数的变化情况。
讨论与分析：通过对实验结果的讨论和分析，我们可以得出以下结论：1. 电路的放大性能：根据实验结果，我们可以评估电路的放大性能。
通过观察输出信号的幅度和输入信号的幅度之间的比值，可以判断电路的放大效果是否符合设计要求。
MATLAB FM仿真报告
实验二FM仿真一实验题目假设基带信号，载波频率为20kHz，FM的频率偏移常数为5kHz/V。
仿真产生FM信号，观察已调信号的波形和频谱。
二基本原理单音频信号经FM调制后的表达式为其中调制指数。
同实验一中相仿，定义必要的仿真参数，在此基础上可得到载波信号和调制信号。
根据可得到频偏，由此可写出最终的FM信号的表达式进行仿真计算。
对FM信号进行傅里叶变换可得频谱特性，变换依旧使用实验一中给出的t2f.m函数。
三仿真方案四仿真源代码cleart1=0.1; %调制信号的时域范围fs=600000; %抽样频率ts=1/fs; %采样率t=-t1:ts:t1;fc=20000; %设定载波频率fc=20kHzs=cos(2*pi*fc*t); %生成载波figureplot(t,s);xlabel('时间');ylabel('幅度');title('载波波形');grid onaxis([-0.0001 0.0001 -2 2]);fm=250; %设定调制信号频率kf=5000;x=sin(2*pi*4*fm*t)+2*cos(2*pi*2*fm*t)+4*sin(2*pi*fm*t+pi/3);%生成调制信号figure;plot(t,x);grid onxlabel('t');ylabel('x');axis([0 0.01 -10 10])title('调制信号波形')figurea=fftshift(fft(x));f=linspace(-fs/2,fs/2,length(t)); plot(f,abs(a));grid onxlabel('ffrequence(Hz)');ylabel('powerspectrum(x)');axis([-1500 1500 0 100]);title('调制信号频谱图');figurephi=2*pi*kf*cumsum(x)*ts;y=cos(2*pi*fc*t+phi);plot(t,y);grid onxlabel('t(s)');ylabel('y');axis([0 0.01 -5 5]);b=fft(y,1024);f=(0:length(b)-1)*fs/length(b)-fs/2; title('已调信号波形');figureplot(f,abs(b));grid onxlabel('frequence(Hz)');ylabel('powerspectrum(x)');axis([-300000 300000 0 200]); title('已调信号频谱图');grid on五实验结果及分析1.调制信号波形极其频谱图根据调制信号表达式可知信号x由三个正弦波相加而成，其中三个正弦波频率符合频率频谱图所示。
基于matlab的FM通信系统仿真设计与实现 课程设计报告
成绩课程设计报告通信系统仿真设计与实现matlab 地FM 题 目： 基于学生姓名： 杨丽君 1008030317学生学号： 电气信息工程学院系 别：电子信息工程专 业： 届 别： 14届马立宪指导教师：电气信息工程学院制 月52013年通信系统仿真设计与实现FM 基于matlab 地 学生：杨丽君 指导教师：马立宪 电气信息工程学院 电子信息工程 课程设计地任务与要求1 课程设计地任务1.1.MATLAB 文件中M 文件地使用方法，包括函数、原理和方法地应用(1)熟悉 .FM 信号调制原理地理解(2)加深对 .MATLAB 画出基于地FM 通信系统仿真设计与实现设计地原理图(3) 课程设计地要求1.2.软件地安装MATLAB 学会(1)(2)在做完FM 调制仿真之后，在今后遇到类似地问题，学会对所面对地问题进行系统地分析，并能从多个层面进行比较..TLAB 软件地使用，学会输入程序并且加以运行(3) 熟练并且掌握对MA 课程设计地研究基础1.3通信地目地是传输信号.通信系统地作用是将信息从信息源发送到一个或者多个目地地.模拟信号是时间和幅值上都连续地信号.调制是用原始信号即调制信号去控制高频载波信号地某一参数，是指随着原始信号幅度地变化而变化.而FM频率调制是高频载波信号地频率随着原始信号幅度变化而变化.解调是将已调制地信号恢复成原始信号即基带调制信号.以下是通信系统地一般模型：1通信系统地一般模型图信息源1)（根据消息地种类不同，信息源可以分信息源（简称信源）地作用是把各种消息转换成电信号.为模拟信号源和数字信号源，模拟信号源输出模拟信号；数字信号源输出数字信号（本次课程设计.是模拟信号源））发送设备（2发送设备地作用是产生适合在信道中传输地信号即使发送信号地特性与信道特性相匹配，具因此，发送设备涵盖地内容很.有抗信道干扰能力，并且具有足够地功率以满足远距离传输地需要.多，包含变换、放大、滤波、编码、调制地过程）信道（3信道分为无线信道和有线信信道是一种物理媒介，用于将来自发送设备地信号发送到接收端.信道地固有特性及...在有线信道中可以是明线、电缆、光纤在无线信道中，信道可以是自由空间道.引入地干扰与噪声直接关系到通信地质量噪声源4（）噪声源是信道中及分布在系统中地其他各处噪声集中表示，噪声是随机地、形式是多样地，.它地出现直接干扰信号地传输）接收设备5（目地是从受到减损地信号中正确.接收设备地功能是将信号放大及反变换（如译码、解调等）．恢复出原始信号，减少在传输过程中噪声与干扰所带来地影响.通信系统方案制定2 FM方案提出2.12模拟通信系统模型图模拟通信系统模型中对于调制信号通过调制器产生地调频波有两种方法其一本次课程设计FM解调器对应地解调方法也有两种其一为相干.为直接调制即宽带调制，其二为间接调制及窄带调制所以据此我提出了四种方案：.解调另外一种为非相干解调方案一：模型中调制器中地调制方法为直接调制即宽带调制，解调器对应地解调方法为相干解.调方案二：模型中调制器中地调制方法为直接调制即宽带调制，解调器对应地解调方法为非相干.解调方案三：模型中调制器中地调制方法为间接调制即窄带调制，解调器对应地解调方法为非相干解调方案四：模型中调制器中地调制方法为间接调制即宽带调制，解调器对应地解调方法为相干解.调方案论证2.2此外，它地最大频率. 窄带调频地应用更广泛与宽带调频，我们对此也更为熟悉，技术也更为成熟.所以本次课程设计选择窄带调频偏移较小，占据带宽较窄、抗干扰性能更好等.所以在以上选择地前提下，本次课程设计地解调方一般情况下，相干解调法较适用于窄带调频.法选择相干解调.针对以上地分析选择地最佳方案为方案四通信系统方案设计3 FM通信系统模型设计3.1 FM．)n(t)s/n(oo)n(t噪声解调器信宿)n(t o)s/n(i模拟通信系统模型设计图3通信系统各部分地功能3.2 FM.滤除调制器输出地无用信号使信号与信道相匹配, 有利于信号在信道中传输.发滤波器: 调制器:.为窄带白噪声N(t)为高斯白噪声,则Ni(t)收滤波器: 滤除信号频带以外地噪声.一般设以利于信号在信道上地传输，调制是将用原始信在通信系统中一般需要将信号进行相应调制,调制可分为线性调制和号去控制高频振荡信号地某一参数，使之随原始信号地变化而成规律变化.调制通信系FM、PM等，这里主要讨论FM非线性调制.线性调制有AM、DSB等，非线性调制有统通信系统参数地计算及原理3.3 FM调制原理FM1）（角调制不是线性调制，角调制中已调信号和调制信号频谱之间不是线性关系而是产生出新地与调制中瞬时角频率是关.FM频谱搬移不同地新地频率分量，呈现非线性特性，故又成为非线性调制于调制信号地线性函数，)(tw?w?KFMf KFM)tw??KFMf(, 瞬时角频率偏移量则，瞬时角频率为:为频偏指数。
matlab仿真实验报告,Matlab仿真及其应用实验报告.doc
matlab仿真实验报告,Matlab仿真及其应⽤实验报告.doc Matlab仿真及其应⽤ 实验报告温州⼤学物理与电⼦信息⼯程学院Matlab仿真及其应⽤ 实验报告课程名称：Matlab仿真及其应⽤班 级：10电信姓名：吴** 学号：1011000****实验地点：5B305⽇期：12.25实验⼆ Matlab 基本编程基础[实验⽬的和要求]熟悉MATLAB环境与⼯作空间熟悉变量与矩阵的输⼊、矩阵的运算熟悉M⽂件与M函数的编写与应⽤熟悉MATLAB控制语句与逻辑运算掌握if语句、switch语句、try语句的使⽤。
掌握利⽤for语句、while语句实现循环结构的⽅法。
[实验内容]1⾏100列的Fibonacc 数组a，a(1)=a(2)=1,a(i)=a(i-1)+a(i-2),⽤for循环指令来寻求该数组中第⼀个⼤于10000的元素，并之处其位置i。
编写M函数表⽰曲线以及它的包络线，并从命令窗⼝输⼊命令语句绘制曲线。
t的取值范围是[0，4π]。
设，编写⼀个M函数⽂件，使得调⽤f(x)时，x可⽤矩阵代⼊，得出的f(x)为同阶矩阵。
根据，求时的最⼤n值；与(1)的n值对应的y值。
已知求中，最⼤值、最⼩值、各数之和，以及正数、零、负数的个数。
输⼊⼀个百分制成绩，要求输出成绩等级A,B,C,D,E。
其中，90~100分为A，80~89分为B，70~79分为C，60~69分为D，60分以下为E。
求分段函数的值。
⽤if语句实现输出x=-5.0, -3.0, 1.0, 2.0, 2.5, 3.0, 5.0时的y值。
编写⼀M函数，实现近似计算指数，其中x为函数参数输⼊，当n+1步与n步的结果误差⼩于0.00001时停⽌。
编写⼀M函数，a和x作为M函数参数输⼊，函数⾥⾯分别⽤if结构实现函数表⽰实验结果及分析：1.a=ones(1,100); %定义数组for i=3:100a(i)=a(i-1)+a(i-2);if(a(i)>10000)a(i),break;endend ,i2.function y=ff(t)y1=exp(-t/3);y2=exp(-t/3).*sin(3*t); y=[y1;y2]3.function y=f(x);a=input('输⼊a值：');x=input('输⼊x值：');if(x<=-a)y=-1;elseif(x-a)y=x/a;elsey=1;endend4.for n=1:100f(n)=1./(2*n-1);y=sum(f)if y>=3my=y-f(n)breakendendmy5.f(1)=1,f(2)=0,f(3)=1; for n=4:100f(n)=f(n-1)-2*f(n-2)+f(n-3);enda=sum(f);b=max(f);c=min(f);p=f==0,d=sum(p);%p等于f为0的个数p1=f>0,e=sum(p1);p2=f<0,f=sum(p2);a,b,c,d,e,f6.clear;n=input('输⼊成绩：');m=floor(n/10);%取整switch mcase num2cell(9:10)disp('A'); %显⽰在控制框case 8disp('B');case 7disp('C');case 6disp('D');case num2cell(0:5)disp('E');otherwisedisp('error')end7.function y=ex3_4(x)for i=1:length(x)if (x(i)<0)&(x(i)~=-3)y(i)=x(i)^2+x(i)-6elseif (x(i)>=0)&(x(i)<5)&(x(i)~=2)&(x(i)~=3) y(i)=x(i)^2-5*x(i)+6else y(i)=x(i)^2-x(i)-1 endendy8.function t=ex3_4(x) n=0;t=1;y=1;x=input(‘’);while y>=0.00001n=n+1;y=x^n/factorial(n);t=t+y;endn9.function y=f(x);a=input('输⼊a值：'); x=input('输⼊x值：'); if。
基于matlab的FM通信系统仿真设计与实现_课程设计报告
成绩课程设计报告题目：基于matlab的FM通信系统仿真设计与实现学生姓名：杨丽君学生学号：1008030317系别：电气信息工程学院专业：电子信息工程届别：14届指导教师：马立宪电气信息工程学院制2013年5月基于matlab的FM通信系统仿真设计与实现学生：杨丽君指导教师：马立宪电气信息工程学院电子信息工程1课程设计的任务与要求1.1 课程设计的任务(1)熟悉MATLAB文件中M文件的使用方法，包括函数、原理和方法的应用。
(2)加深对FM信号调制原理的理解。
(3)画出基于MATLAB的FM通信系统仿真设计与实现设计的原理图。
1.2 课程设计的要求(1) 学会MATLAB软件的安装。
(2)在做完FM调制仿真之后，在今后遇到类似的问题，学会对所面对的问题进行系统的分析，并能从多个层面进行比较。
(3) 熟练并且掌握对MATLAB软件的使用，学会输入程序并且加以运行。
1.3 课程设计的研究基础通信的目的是传输信号。
通信系统的作用是将信息从信息源发送到一个或者多个目的地。
模拟信号是时间和幅值上都连续的信号。
调制是用原始信号即调制信号去控制高频载波信号的某一参数，是指随着原始信号幅度的变化而变化。
而FM频率调制是高频载波信号的频率随着原始信号幅度变化而变化。
解调是将已调制的信号恢复成原始信号即基带调制信号。
以下是通信系统的一般模型：图1通信系统的一般模型（1)信息源信息源（简称信源）的作用是把各种消息转换成电信号。
根据消息的种类不同，信息源可以分为模拟信号源和数字信号源，模拟信号源输出模拟信号；数字信号源输出数字信号（本次课程设计是模拟信号源）。
（2）发送设备发送设备的作用是产生适合在信道中传输的信号即使发送信号的特性与信道特性相匹配，具有抗信道干扰能力，并且具有足够的功率以满足远距离传输的需要。
因此，发送设备涵盖的内容很多，包含变换、放大、滤波、编码、调制的过程。
（3）信道信道是一种物理媒介，用于将来自发送设备的信号发送到接收端。
基于MATLAB的AM与FM的调制系统仿真
总732期第三十四期2020年12月河南科技Henan Science and Technology基于MATLAB的AM与FM的调制系统仿真张博文林君（延边大学，吉林延吉133002）摘要：本文通过MATLAB进行AM信号与FM信号的调制系统仿真实验，分析了AM和FM调制的原理，观察了AM、FM信号调制过程中各环节的时域和频域波形。
最后，在仿真的基础上分析并比较了这两种调制方式的性能。
关键字：MATLAB；AM调制；FM调制；调制系统中图分类号：TN911.3文献标识码：A文章编号：1003-5168（2020）34-0037-03 Modulation System Simulation Based on MATLAB AM and FMZHANG Bowen LIN Jun（Yanbian University，Yanji Jilin133002）Abstract:In this paper,the Modulation system simulation experiment of AM signal and FM signal was carried out by MATLAB,the principle of AM and FM modulation was analyzed,and the time-domain and frequency-domain wave⁃form of each link in the modulation process of AM and FM signal was observed.Finally,the performance of the two modes was analyzed and compared on the basis of simulation.Keywords:MATLAB；AM modulation；FM modulation；modulation system通信系统包含两种调制方式，即模拟调制方式和数字调制方式［1］。
matlab FM调制仿真
Matｌab FＭ调制仿真ﻬ目录引言ﻩ错误!未定义书签。
一、课程设计得目得与要求 ............. 错误!未定义书签。
1、1课程设计得目得 .............. 错误!未定义书签。
1、2课程设计得要求ﻩ错误!未定义书签。
二、FM调制解调系统设计.............. 错误!未定义书签。
2、1FM调制模型得建立......... 错误!未定义书签。
2、2调制过程分析.................... 错误!未定义书签。
2、3FM解调模型得建立......... 错误!未定义书签。
2、４解调过程分析ﻩ错误!未定义书签。
２、５高斯白噪声信道特性ﻩ错误!未定义书签。
2、6调频系统得抗噪声性能分析ﻩ错误!未定义书签。
三、仿真实现ﻩ错误!未定义书签。
3、１ＭATLAB源代码............ 错误!未定义书签。
3、2仿真结果 .......................... 错误!未定义书签。
四、心得体会ﻩ错误!未定义书签。
五、参考文献..................................... 错误!未定义书签。
引言本课程设计用于实现ＤSＢ信号得调制解调过程。
信号得调制与解调在通信系统中具有重要得作用。
调制过程就就是一个频谱搬移得过程,它就就是将低频信号得频谱搬移到载频位置。
解调就就是调制得逆过程,即就就是将已调制得信号还原成原始基带信号得过程。
信号得接收端就就就是通过解调来还原已调制信号从而读取发送端发送得信息。
因此信号得解调对系统得传输有效性与传输可靠性有着很大得影响。
调制与解调方式往往决定了一个通信系统得性能。
双边带DSB信号得解调采用相干解调法,这种方式被广泛应用在载波通信与短波无线电话通信中。
一、课程设计得目得与要求1、1 课程设计得目得通过《ＦM调制解调系统设计与防真》得课程设计,掌握通信原理中模拟信号得调制与解调、数字基带信号得传输、数字信号得调制与解调,模拟信号得抽样、量化与编码与信号得最佳接收等原理。
基于Matlab的模拟(AM、FM、PM)调制系统仿真
通信系统模拟调制系统仿真一 课题内容 AM FM PM 调制 二 设计要求1.掌握AM FM PM 调制和解调原理。
2.学会Matlab 仿真软件在AM FM PM 调制和解调中的应用。
3.分析波形及频谱1.AM 调制解调系统设计1.振幅调制产生原理所谓调制，就是在传送信号的一方将所要传送的信号附加在高频振荡上，再由天线发射出去。
这里高频振荡波就是携带信号的运载工具，也叫载波。
振幅调制，就是由调制信号去控制高频载波的振幅，直至随调制信号做线性变化。
在线性调制系列中，最先应用的一种幅度调制是全调幅或常规调幅，简称为调幅（AM ）。
在频域中已调波频谱是基带调制信号频谱的线性位移；在时域中，已调波包络与调制信号波形呈线性关系。
设正弦载波为)cos()(0ϕω+=t A t c c式中，A 为载波幅度；c ω为载波角频率；0ϕ为载波初始相位(通常假设0ϕ=0).调制信号（基带信号）为)(t m 。
根据调制的定义，振幅调制信号（已调信号）一般可以表示为)cos()()(t t Am t s c m ω=设调制信号)(t m 的频谱为)(ωM ，则已调信号)(t s m 的频谱)(ωm S ：)]()([2)(c c m M M AS ωωωωω-++=2.调幅电路方案分析标准调幅波（AM ）产生原理调制信号是只来来自信源的调制信号（基带信号），这些信号可以是模拟的，亦可以是数字的。
为首调制的高频振荡信号可称为载波，它可以是正弦波，亦可以是非正弦波(如周期性脉冲序列)。
载波由高频信号源直接产生即可，然后经过高频功率放大器进行放大，作为调幅波的载波，调制信号由低频信号源直接产生，二者经过乘法器后即可产生双边带的调幅波。
设载波信号的表达式为t c ωcos ，调制信号的表达式为t A t m m m ωcos )(= ，则调幅信号的表达式为t t m A t s c AM ωcos )]([)(0+=图5.1 标准调幅波示意图 3.信号解调思路从高频已调信号中恢复出调制信号的过程称为解调(demodulation )，又称为检波(detection )。
FM信号仿真实验
《实验四FM信号仿真实验》实验（上机）报告实验名称FM信号仿真实验验证型实验目的及要求：1．掌握FM信号系统工作原理；2.掌握FM信号的的Matlab/Simulink仿真方法。
3.FM系统的Simulink模型；4.参数设置说明；5. 仿真结果及说明；6. 意见和建议；实验内容：1.进一步熟悉并掌握Matlab/Simulink基本库、通信库和DSP库中较为重要的一些功能模块的作用以及相应功能参数的物理意义与设置方法。
2.搭建FM信号仿真模型如图1所示。
设置系统参数并调试，同时观测并记录各点的时域波形及频域图。
3.改变Random Number模块所产生的噪声大小，重测图所示模型中各点的时域波形以及频域图。
请写出实验内容2的操作步骤：2.搭建FM信号仿真模型如图1所示。
设置系统参数并调试，同时观测并记录各点的时域波形及频域图。
3.改变Random Number模块所产生的噪声大小，重测图所示模型中各点的时域波形以及频域图。
Signal generatror设置为：FM Modulator Passband设置为：FM Demodulator Passband设置为：Zero-Order Hold设置为：Specrtum Scope的设置为：Simulation time的设置为：操作步骤（续）：Scope中所得波形：Scope1中所得波形：从FM Modulator Passband中输出的频谱为：信号源的频谱为：从FM Demodulator Passband中输出的频谱为：使用的FM仿真模型图为：实验中遇到的问题和总结：在实验过程中参数不知道怎么设置，以至于看到的结果不是很明显或是看不到调试结果。
所以通过实验了解到了各个参数设置，得到仿真结果。
教师批语与成绩评定：评阅教师：年月日。
基于MATLAB的模拟信号频率调制(FM)与解调分析
课程设计任务书学生姓名：杨刚专业班级：电信1302指导教师: 工作单位：武汉理工大学题目信号分析处理课程设计—基于MATLAB的模拟信号频率调制（FM）与解调分析初始条件:1.Matlab6.5以上版本软件；2.先修课程：通信原理等；要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1、利用MATLAB中的simulink工具箱中的模块进行模拟频率（FM）调制与解调，观察波形变化2、画出程序设计框图，编写程序代码，上机运行调试程序，记录实验结果（含计算结果和图表等），并对实验结果进行分析和总结；3、课程设计说明书按学校统一规范来撰写，具体包括：⑴目录；⑵理论分析；⑶ 程序设计；⑷ 程序运行结果及图表分析和总结；⑸课程设计的心得体会（至少800字，必须手写。
）；⑹参考文献（不少于5篇）。
时间安排：周一、周二查阅资料，了解设计内容；周三、周四程序设计，上机调试程序；周五、整理实验结果，撰写课程设计说明书2013系主任（或责任教师）签名: 2013 年7月2日指导教师签名:目录1 Simulink 简介 (1)1.1 Matlab简介.................................... 错误！未定义书签。
1.2 Simulink介绍 .................................. 错误！未定义书签。
2原理分析........................................... 错误！未定义书签。
2.1通信系统....................................... 错误！未定义书签。
2.1.1通信系统的一般模型........................ 错误！未定义书签。
2.1.2模拟通信系统 (3)2.2 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 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的FM仿真实现
基于Matlab的FM仿真实现摘要本次设计主要是以Matlab为基础平台，对FM信号进行仿真。
介绍了FM信号，及其调制和解调的基本原理，并设计M文件，分析在混入噪声环境下的波形失真，以及分析FM的抗噪声性能。
本设计的主要目的是对Matlab的熟悉和对模拟通信理论的更深化理解。
关键词：Matlab；FM；噪声前言通信按照传统的理解就是信息的传输。
在当今高度信息化的社会，信息和通信已成为现代社会的命脉。
信息作为一种资源，只有通过广泛传播与交流，才能产生利用价值，促进社会成员之间的合作，推动社会生产力的发展，创造出巨大的经济效益。
而通信作为传输信息的手段或方式，与传感技术、计算机技术相融合，已成为21世纪国际社会和世界经济发展的强大动力。
可以预见，未来的通信对人们的生活方式和社会的发展将会产生更加重大和意义深远的影响。
在通信系统中，从消息变换过来的原始信号所占的有效频带往往具有频率较低的频谱分量(例如语音信号)，如果将这种信号直接在信道中进行传输，则会严重影响信息传送的有效性和可靠性，因此这种信号在许多信道中均是不适宜直接进行传输的。
在通信系统的发射端通常需要有调制过程，将调制信号的频谱搬移到所希望的位置上，使之转换成适于信道传输或便于信道多路复用的已调信号；而在接收端则需要有解调过程，以恢复原来有用的信号。
调制解调方式常常决定了一个通信系统的性能。
随着数字化波形测量技术和计算机技术的发展，可以使用数字化方法实现调制与解调过程。
调制在通信系统中具有重要的作用。
通过调制，不仅可以进行频谱搬移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合于信道传输或便于信道多路复用的已调信号，而且它对系统的传输有效性和传输可靠性有着很大的影响。
调制方式往往决定了一个通信系统的性能。
调制技术是指把基带信号变换成传输信号的技术。
基带信号是原始的电信号，一般是指基本的信号波形，在数字通信中则指相应的电脉冲。
在无线遥测遥控系统和无线电技术中调制就是用基带信号控制高频载波的参数（振幅、频率和相位），使这些参数随基带信号变化。
基于matlab的FM通信系统仿真设计与实现_课程设计报告
题目：基于matlab地FM通信系统仿真设计与实现学生姓名：杨丽君学生学号： 1008030317系别：电气信息工程学院专业：电子信息工程届别： 14届指导教师：马立宪电气信息工程学院制2013年5月基于matlab地FM通信系统仿真设计与实现学生：杨丽君指导教师：马立宪电气信息工程学院电子信息工程1课程设计地任务与要求1.1 课程设计地任务(1)熟悉MATLAB文件中M文件地使用方法，包括函数、原理和方法地应用.(2)加深对FM信号调制原理地理解.(3)画出基于MATLAB地FM通信系统仿真设计与实现设计地原理图.1.2 课程设计地要求(1) 学会MATLAB软件地安装.(2)在做完FM调制仿真之后，在今后遇到类似地问题，学会对所面对地问题进行系统地分析，并能从多个层面进行比较.(3) 熟练并且掌握对MA TLAB软件地使用，学会输入程序并且加以运行.1.3 课程设计地研究基础通信地目地是传输信号.通信系统地作用是将信息从信息源发送到一个或者多个目地地.模拟信号是时间和幅值上都连续地信号.调制是用原始信号即调制信号去控制高频载波信号地某一参数，是指随着原始信号幅度地变化而变化.而FM频率调制是高频载波信号地频率随着原始信号幅度变化而变化.解调是将已调制地信号恢复成原始信号即基带调制信号.以下是通信系统地一般模型：图1通信系统地一般模型（1)信息源信息源（简称信源）地作用是把各种消息转换成电信号.根据消息地种类不同，信息源可以分为模拟信号源和数字信号源，模拟信号源输出模拟信号；数字信号源输出数字信号（本次课程设计是模拟信号源）.（2）发送设备发送设备地作用是产生适合在信道中传输地信号即使发送信号地特性与信道特性相匹配，具有抗信道干扰能力，并且具有足够地功率以满足远距离传输地需要.因此，发送设备涵盖地内容很多，包含变换、放大、滤波、编码、调制地过程.（3）信道信道是一种物理媒介，用于将来自发送设备地信号发送到接收端.信道分为无线信道和有线信道.在有线信道中可以是明线、电缆、光纤.在无线信道中，信道可以是自由空间.信道地固有特性及引入地干扰与噪声直接关系到通信地质量.（4）噪声源噪声源是信道中及分布在系统中地其他各处噪声集中表示，噪声是随机地、形式是多样地，它地出现直接干扰信号地传输.（5）接收设备接收设备地功能是将信号放大及反变换（如译码、解调等）.目地是从受到减损地信号中正确恢复出原始信号，减少在传输过程中噪声与干扰所带来地影响.2 FM通信系统方案制定2.1 方案提出下图为大体模拟通信系统模型：图2模拟通信系统模型本次课程设计FM模拟通信系统模型中对于调制信号通过调制器产生地调频波有两种方法其一为直接调制即宽带调制，其二为间接调制及窄带调制.解调器对应地解调方法也有两种其一为相干解调另外一种为非相干解调.所以据此我提出了四种方案：方案一：模型中调制器中地调制方法为直接调制即宽带调制，解调器对应地解调方法为相干解调.方案二：模型中调制器中地调制方法为直接调制即宽带调制，解调器对应地解调方法为非相干解调.方案三：模型中调制器中地调制方法为间接调制即窄带调制，解调器对应地解调方法为非相干解调方案四：模型中调制器中地调制方法为间接调制即宽带调制，解调器对应地解调方法为相干解调.2.2 方案论证窄带调频地应用更广泛与宽带调频，我们对此也更为熟悉，技术也更为成熟.此外，它地最大频率偏移较小，占据带宽较窄、抗干扰性能更好等.所以本次课程设计选择窄带调频.一般情况下，相干解调法较适用于窄带调频.所以在以上选择地前提下，本次课程设计地解调方法选择相干解调针对以上地分析选择地最佳方案为方案四.3 FM通信系统方案设计3.1 FM通信系统模型设计图3 模拟通信系统模型设计3.2 FM 通信系统各部分地功能调制器: 使信号与信道相匹配, 有利于信号在信道中传输.发滤波器: 滤除调制器输出地无用信号.收滤波器: 滤除信号频带以外地噪声.一般设N(t)为高斯白噪声,则Ni(t)为窄带白噪声.在通信系统中一般需要将信号进行相应调制,以利于信号在信道上地传输，调制是将用原始信号去控制高频振荡信号地某一参数，使之随原始信号地变化而成规律变化.调制可分为线性调制和非线性调制.线性调制有AM 、DSB 等，非线性调制有FM 、PM 等，这里主要讨论FM 调制通信系统 3.3 FM 通信系统参数地计算及原理 （1）FM 调制原理角调制不是线性调制，角调制中已调信号和调制信号频谱之间不是线性关系而是产生出新地与频谱搬移不同地新地频率分量，呈现非线性特性，故又成为非线性调制.FM 调制中瞬时角频率是关于调制信号地线性函数， 瞬时角频率偏移量)(t KFMf w =∆ 则， 瞬时角频率为:)(t KFMf w w c +=。
实验一 模拟通信的MATLAB仿真..
实验一 模拟通信的MATLAB 仿真姓名：左立刚 学号：031040522简要说明：实验报告注意包括AM ，DSB ，SSB ，VSB ，FM 五种调制与解调方式的实验原理，程序流程图，程序运行波形图，simulink 仿真模型及波形，心得体会，最后在附录中给出了m 语言的源程序代码。
一．实验原理1．幅度调制（AM ）幅度调制（AM ）是指用调制信号去控制高频载波的幅度，使其随调制信号呈线性变化的过程。
AM 信号的数学模型如图3-1所示。
图2-1 AM 信号的数学模型为了分析问题的方便，令δ=0，1.1 AM 信号的时域和频域表达式()t S AM=[A 0+m ()t ]cos t cω （2-1）()t S AM =A 0π[()()ωωωωδC C ++-]+()()[]ωωωωc c M M ++-21（2-2）AM 信号的带宽2=BAMfH（2-3）式中，fH为调制信号的最高频率。
2.1.3 AM 信号的功率P AM 与调制效率ηAMP AM=()2222t m A +=PP mc + （2-4）式中，P C=2A为不携带信息的载波功率；()22t m P m=为携带信息的边带功率。
()()t t m A m PP AMCAM222+==η（2-5） AM 调制的优点是可用包络检波法解调，不需要本地同步载波信号，设备简单。
AM 调制的最大缺点是调制效率低。
2.2、双边带调制（DSB ）如果将在AM 信号中载波抑制，只需在图3-1中将直流 A 0去掉，即可输出抑制载波双边带信号。
2.2.1 DSB 信号的时域和频域表达式()()t t m t cDSB S ωcos= （2-6）()()()[]ωωωωωC C DSBM M S ++-=21 （2-7） DSB 信号的带宽fB BHAM DSB2== （2-8）DSB 信号的功率及调制效率由于不再包含载波成分，因此，DSB 信号的功率就等于边带功率，是调制信号功率的一半，即()()t t m PS P CDSB DSB 2221=== （2-9） 显然，DSB 信号的调制效率为100%。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验二FM仿真
一实验题目
假设基带信号，载波频率为20kHz，FM的频率偏移常数为5kHz/V。
仿真产生FM信号，观察已调信号的波形和频谱。
二基本原理
单音频信号
经FM调制后的表达式为
其中
调制指数。
同实验一中相仿，定义必要的仿真参数，在此基础上可得到载波信号和调制
信号。
根据可得到频偏，由此可写出最终的FM信号的
表达式进行仿真计算。
对FM信号进行傅里叶变换可得频谱特性，变换依旧使
用实验一中给出的t2f.m函数。
三仿真方案
四仿真源代码
clear
t1=0.1; %调制信号的时域范围
fs=600000; %抽样频率
ts=1/fs; %采样率
t=-t1:ts:t1;
fc=20000; %设定载波频率fc=20kHz
s=cos(2*pi*fc*t); %生成载波
figure
plot(t,s);
xlabel('时间');
ylabel('幅度');
title('载波波形');
grid on
axis([-0.0001 0.0001 -2 2]);
fm=250; %设定调制信号频率
kf=5000;
x=sin(2*pi*4*fm*t)+2*cos(2*pi*2*fm*t)+4*sin(2*pi*fm*t+pi/3);%生成调制信号figure;
plot(t,x);
grid on
xlabel('t');ylabel('x');
axis([0 0.01 -10 10])
title('调制信号波形')
figure
a=fftshift(fft(x));
f=linspace(-fs/2,fs/2,length(t)); plot(f,abs(a));
grid on
xlabel('ffrequence(Hz)');
ylabel('powerspectrum(x)');
axis([-1500 1500 0 100]);
title('调制信号频谱图');
figure
phi=2*pi*kf*cumsum(x)*ts;
y=cos(2*pi*fc*t+phi);
plot(t,y);
grid on
xlabel('t(s)');ylabel('y');
axis([0 0.01 -5 5]);
b=fft(y,1024);
f=(0:length(b)-1)*fs/length(b)-fs/2; title('已调信号波形');
figure
plot(f,abs(b));
grid on
xlabel('frequence(Hz)');
ylabel('powerspectrum(x)');
axis([-300000 300000 0 200]); title('已调信号频谱图');
grid on
五实验结果及分析
1.调制信号波形极其频谱图
根据调制信号表达式可知信号x由三个正弦波相加而成，其中三个正弦波频率符合频率频谱图所示。
2.载波波形
载波波形符合程序设计的20kHz的正弦波。
3.仿真FM信号波形，其形状为疏密波，最大频偏5kHz/V。
由图可读出带宽为70kHz，符合公式计算B≈2（β+1）fm=2（6.1×5+4）=69kHz
4.已调信号频谱图
六实验心得
通过这次通信原理实验，我发现自己的基础知识和应用能力都很差。
刚拿到题目时，在如何实现FM调制的地方就卡住了。
最后通过各种网上的
资料解决了问题。
通过了这次实验，我加深了课堂上所学的理论知识，提
高了知识的应用能力。
熟悉了MATLAB的运行环境，掌握了一些MATLAB语
言和函数的调用，很大程度上提高了自己的能力。