基于MATLAB調(diào)制解調(diào)仿真設(shè)計(jì)

1、 PAGE11 / NUMPAGES11數(shù)字信號(hào)處理 實(shí)驗(yàn)八調(diào)制解調(diào)系統(tǒng)的實(shí)現(xiàn)一、實(shí)驗(yàn)?zāi)康? 深刻理解濾波器的設(shè)計(jì)指標(biāo)與根據(jù)指標(biāo)進(jìn)行數(shù)字濾波器設(shè)計(jì)的過程了解濾波器在通信系統(tǒng)中的應(yīng)用二、實(shí)驗(yàn)步驟:1.通過SYSTEMVIEW軟件設(shè)計(jì)與仿真工具，設(shè)計(jì)一個(gè)FIR數(shù)字帶通濾波器，預(yù)先給定截止頻率和在截止頻率上的幅度值， 通過軟件設(shè)計(jì)完后，確認(rèn)濾波器的階數(shù)和系統(tǒng)函數(shù)，畫出該濾波器的頻率響應(yīng)曲線，進(jìn)行技術(shù)指標(biāo)的驗(yàn)證。通過仿真驗(yàn)證，原理圖如下：輸入方波與鋸齒波，都為10HZ，載波100hz與300hz正弦波，仿真的結(jié)果如下：還是可以比較好的恢復(fù)信號(hào)。基帶信號(hào)1建立一個(gè)兩載波幅度調(diào)制與解調(diào)的通信系統(tǒng)，將該濾波

2、器作為兩個(gè)載波分別解調(diào)的關(guān)鍵部件，驗(yàn)證其帶通的頻率特性的有效性。系統(tǒng)框圖如下：sin1基帶信號(hào)1+基帶信號(hào)2帶通濾波器1中心頻率1帶通濾波器2中心頻率2基帶信號(hào)2sin2sin1sin2低通濾波低通濾波規(guī)劃整個(gè)系統(tǒng)，確定系統(tǒng)的采樣頻率、觀測(cè)時(shí)間、細(xì)化并設(shè)計(jì)整個(gè)系統(tǒng)，仿真調(diào)整并不斷改進(jìn)達(dá)到正確調(diào)制、正確濾波、正確解調(diào)的目的。（參考文件zhan3.svu）設(shè)計(jì)的思路是：基帶信號(hào)乘上一個(gè)高頻信號(hào)，稱為調(diào)制，實(shí)現(xiàn)頻譜搬移，與另一調(diào)制信號(hào)疊加，再分別通過以W為中心頻率的帶通FIR數(shù)字濾波器，再乘以原來(lái)的高頻信號(hào)，實(shí)現(xiàn)再頻譜搬移，最后通過IIR低通濾波器得到解調(diào)信號(hào)。本實(shí)驗(yàn)是通過編程的方式完成的。1、首先

3、，產(chǎn)生信號(hào)：n=1;f1=100;f2=300;fs=1000;%采樣頻率t=0:1/fs:n;fre=10; y1=square(2*fre*pi*t)/2+1.1;y2=sawtooth(fre*2*pi*t)/2+1.1;觀察圖形與頻譜：2、基帶信號(hào)乘以一個(gè)高頻載波：yy1=y1.*z1;yy2=y2.*z2;觀察頻譜：3、兩調(diào)制信號(hào)相加：yy3=yy1+yy2;4、設(shè)計(jì)帶通濾波器：設(shè)計(jì)100-200,330-430hz的FIR數(shù)字濾波器，用漢明窗實(shí)現(xiàn)。fp1=100;fp2=200;%FIR濾波器100-200hzfs1=50;fs2=250;As=15;Ws1=(fp1+fs1)/f

4、s;Ws2=(fp2+fs2)/fs;w=(fp1-fs1)/fs;M=ceil(As-7.95)/(14.36*w);hamming=Hamming(M+1);b=fir1(M,Ws1,Ws2,hamming);%figure(2);%freqz(b,1,fs,fs);t=0:1/fs:n;yyy1=filter(b,2,yy3);zz1=filter(b,2,z1);fp1=330;fp2=430;%FIR濾波器330-430hzfs1=200;fs2=490;As=15;Ws1=(fp1+fs1)/fs;Ws2=(fp2+fs2)/fs;w=(fp1-fs1)/fs;M=ceil(As-

5、7.95)/(14.36*w);hamming=Hamming(M+1);b=fir1(M,Ws1,Ws2,hamming);%figure(4);%freqz(b,1,fs,fs);t=0:1/fs:n;yyy2=filter(b,2,yy3);zz2=filter(b,2,z2);5、濾波后的信號(hào)再乘載波信號(hào)：k1=yyy1.*zz1;k2=yyy2.*zz2;6、設(shè)計(jì)低通濾波器為100hz的巴特沃斯低通濾波器。N=8;%8階巴特沃斯低通濾波器上限頻率100hzWn=100/(fs/2);b,a = butter(N,Wn,low);kk1=filter(b,a,k1);%figure(5

6、);%H,W=freqz(b,a); %返回頻率響應(yīng)%subplot(1,2,1);plot(W*fs/(2*pi),abs(H); xlabel(頻率HZ); ylabel(幅值);grid on;%subplot(1,2,2);plot(W*fs/(2*pi),20*log10(abs(H); xlabel(頻率HZ);ylabel(幅值dB);grid on;7、調(diào)制信號(hào)經(jīng)過低通濾波器后的到的信號(hào)就是解調(diào)信號(hào)。8、觀察調(diào)制信號(hào)的頻譜：結(jié)論與體會(huì)：調(diào)制解調(diào)實(shí)際就是頻譜的搬移，達(dá)到高頻傳播，本次實(shí)驗(yàn)的關(guān)鍵就是濾波器的設(shè)計(jì)，濾波器設(shè)計(jì)不好解調(diào)不了波形導(dǎo)致失真，通過寫代碼使我的編程能力加強(qiáng)了，學(xué)

7、會(huì)了移植代碼，matlab是個(gè)很不錯(cuò)的仿真軟件。附錄：n=1;f1=100;f2=300;fs=1000;%采樣頻率t=0:1/fs:n;fre=10; y1=square(2*fre*pi*t)/2+1.1;dt=1/fs; %定義時(shí)間步長(zhǎng)。n1=length(t); %樣點(diǎn)個(gè)數(shù)%y1=cos(2*pi*fre*t); %余弦信號(hào)f_end=1/dt; %頻率軸的顯示圍f=(0:n1-1)*f_end/n1-f_end/2; %頻率自變量Xf=dt*fftshift(fft(y1); %頻譜figure(1);subplot(211);plot(t,y1);xlabel(t);title(時(shí)

8、間波形);%時(shí)間波形subplot(212);plot(f,abs(Xf);xlabel(f);title(方波幅度頻譜);%頻譜波形y2=sawtooth(fre*2*pi*t)/2+1.1;Xf1=dt*fftshift(fft(y2); %頻譜figure(2);subplot(211);plot(t,y2);xlabel(t);title(時(shí)間波形);%時(shí)間波形subplot(212);plot(f,abs(Xf1);xlabel(f);title(三角波幅度頻譜);%頻譜波形z1=10*sin(2*pi*f1*t);z2=10*sin(2*pi*f2*t);yy1=y1.*z1;yy

9、2=y2.*z2;yy3=yy1+yy2;Xf=dt*fftshift(fft(yy1); %頻譜figure(3);subplot(211);plot(t,yy1);xlabel(t);title(時(shí)間波形);%時(shí)間波形subplot(212);plot(f,abs(Xf);xlabel(f);title(調(diào)制信號(hào)頻譜);%頻譜波形Xf=dt*fftshift(fft(yy2); %頻譜figure(7);subplot(211);plot(t,yy2);xlabel(t);title(時(shí)間波形);%時(shí)間波形subplot(212);plot(f,abs(Xf);xlabel(f);titl

10、e(調(diào)制信號(hào)頻譜);%頻譜波形fp1=100;fp2=200;%FIR濾波器100-200hzfs1=50;fs2=250;As=15;Ws1=(fp1+fs1)/fs;Ws2=(fp2+fs2)/fs;w=(fp1-fs1)/fs;M=ceil(As-7.95)/(14.36*w);hamming=Hamming(M+1);b=fir1(M,Ws1,Ws2,hamming);%figure(2);%freqz(b,1,fs,fs);t=0:1/fs:n;yyy1=filter(b,2,yy3);zz1=filter(b,2,z1);fp1=330;fp2=430;%FIR濾波器330-430

11、hzfs1=200;fs2=490;As=15;Ws1=(fp1+fs1)/fs;Ws2=(fp2+fs2)/fs;w=(fp1-fs1)/fs;M=ceil(As-7.95)/(14.36*w);hamming=Hamming(M+1);b=fir1(M,Ws1,Ws2,hamming);%figure(4);%freqz(b,1,fs,fs);t=0:1/fs:n;yyy2=filter(b,2,yy3);zz2=filter(b,2,z2);k1=yyy1.*zz1;k2=yyy2.*zz2;N=8;%8階巴特沃斯低通濾波器上限頻率100hzWn=100/(fs/2);b,a = but

12、ter(N,Wn,low);kk1=filter(b,a,k1);figure(10);H,W=freqz(b,a); %返回頻率響應(yīng)subplot(1,2,1);plot(W*fs/(2*pi),abs(H); xlabel(頻率HZ); ylabel(幅值);grid on;subplot(1,2,2);plot(W*fs/(2*pi),20*log10(abs(H); xlabel(頻率HZ);ylabel(幅值dB);grid on;figure(4);subplot(211);plot(t,kk1);grid on;title(解調(diào)信號(hào)kk1);kk2=filter(b,a,k2);subplot(212);plot(t,kk2);grid on;title(解調(diào)信號(hào)kk2);Xf2=dt*fftshift(fft(kk1); %頻譜figure(5);subplot(211);plot(t,kk1);xlabel(t);title(時(shí)間波形);%時(shí)間波形subplot(212);pl