基于matlab的4FSK系統(tǒng)設(shè)計仿真

1、 課程設(shè)計任務(wù)書課程名稱： 專業(yè)課程設(shè)計 題 目： 基于matlab的4FSK系統(tǒng)設(shè)計仿真 學(xué) 院： 電子系 系： 信息工程學(xué)院 專業(yè)班級： 學(xué) 號： 學(xué)生姓名： 指導(dǎo)教師：職稱： 學(xué)院審核（簽名）： 審核日期： 設(shè)計基本原理與系統(tǒng)框圖以前學(xué)過2FSK信號的產(chǎn)生，知道它有兩種方法：調(diào)頻法和開關(guān)法，前者是用二進制基帶矩形脈沖信號去調(diào)制一個調(diào)頻器，使其能輸出兩個不同頻率的碼元；后者是用一個受基帶脈沖控制的開關(guān)電路去選擇兩個獨立頻率源的振蕩作為輸出。2FSK鍵控法調(diào)頻原理圖如下：2fsk鍵控原理圖振蕩器f1選通開關(guān)反相器振蕩器f2選通開關(guān)相加器這里我們要研究的是4FSK信號，是采用第二種方法得到的，

2、即用基帶四進制信號去鍵控四個頻率不同的載波，就可以得到四進制頻移鍵控信號，其中4FSK是采用四個不同的頻率分別表示四進制的的四個碼元00、01、10、11，每個碼元都含有2bit的信息，其波形如圖1-1所示，這時仍和2FSK時的條件相同，即要求每個載頻之間的距離足夠大，使不同頻率的碼元頻譜能夠用濾波器分離開，或者說使不同頻率的碼元相互正交。4FSK調(diào)制原理如下： 傳“0”信號(或00)時，發(fā)送頻率為f1的載波； 傳“1”信號(或10)時，發(fā)送頻率為f2的載波； 傳“2”信號(或11)時，發(fā)送頻率為f3的載波； 傳“3”信號(或01)時，發(fā)送頻率為f4的載波。圖1-1 4FSK信號波形系統(tǒng)方框圖

3、如圖1-2所示4FSK晶振4FSK調(diào)制串/并轉(zhuǎn)換NRZ圖1-2系統(tǒng)方框圖222濾波器濾波器濾波器各單元電路圖設(shè)計本次系統(tǒng)設(shè)計大致可分為四大模塊：晶體振蕩器與信源共用，位于信源單元；多級分頻電路；4FSK調(diào)制中的邏輯電路單元；二進制基帶信號的串/并轉(zhuǎn)換模塊。信源單元電路本模塊是整個課程設(shè)計系統(tǒng)的發(fā)終端，模塊內(nèi)部只使用+5V電壓，其原理方框圖如圖1-3所示本單元用來產(chǎn)生晶振信號和NRZ信號，圖1-4為信源單元電路圖，它上面的元器件與圖1-3上各單元對應(yīng)關(guān)系如下：晶振 CRY：晶體； U1：反相器7404分頻器 U2：計數(shù)器74161；U3：計數(shù)器74193；U4：計數(shù)器40160并行碼產(chǎn)生器 K1

4、、K2、K3：8位手動開關(guān)，從左至右依次與幀同步碼、數(shù)據(jù)1、數(shù)據(jù)2相對應(yīng)；發(fā)光二極管：左起分別與一幀中的24位代碼相對應(yīng)八選一 U5、U6、U7：8位數(shù)據(jù)選擇器4512三選一 U8：8位數(shù)據(jù)選擇器4512倒相器 U20：非門74LS04抽樣 U9：D觸發(fā)器74LS74CLKBS-OUTNRZ-OUTNRZFSBSS5S4S3S2S1晶振并行碼產(chǎn)生器八選一八選一八選一三選一倒相器抽 樣分 頻 器圖1-3信源單元方框圖從晶振產(chǎn)生一個4.096KHz的信號，一路做時鐘信號CLK，一路送到74LS193，74193完成2、4、8運算，輸出BS、S1、S2、S3等4個信號。BS為位同步信號，頻率為204

5、8MHz，S1、S2、S3為3個選通信號，頻率分別為BS信號頻率的1/2、1/4和1/8。74193是一個4位二進制加/減計數(shù)器，當CD= PL =1、MR=0時，可在QA、QB、QC及QD端分別輸出上述4個信號。、多級分頻電路它是由晶振電路和一個具有分頻功能的74LS193芯片構(gòu)成的。系統(tǒng)要求產(chǎn)生四個頻率不同的脈沖序列，所以要通過分頻電路將主時鐘輸出的信號分別進行2分頻、4分頻、8分頻。根據(jù)頻率的不同，采用三個分頻器進行分頻。CLK信號由一個晶振電路實現(xiàn)，信號分別是由三個串接的二分頻器而得到的，它們的頻率分別為CLK頻率的1/2、1/4、1/8，其電路圖如圖2-1所示。圖2-1 分頻電路圖，

6、、4FSK調(diào)制中的邏輯電路單元類似于2FSK調(diào)制模塊，4FSK調(diào)制模塊是采用一個四選一數(shù)據(jù)選擇器，經(jīng)過多次分頻產(chǎn)生的四個不同頻率的信號作為數(shù)據(jù)選擇器的四個輸入端，經(jīng)過串并轉(zhuǎn)換的兩路并行信號作為數(shù)據(jù)選擇器的兩個地址端，對應(yīng)的每一種地址選中不同頻率的輸入信號，即“00”時選中頻率為的信號；“01”時選中頻率為的信號；“10”時選中頻率為的信號；“11”時選中頻率為的信號。圖3-1 邏輯電路圖3-1中的邏輯電路即為雙四選一數(shù)據(jù)選擇器，接口分別接來自多次分頻的頻率不同的四種信號附加控制端接地，分別接串/并轉(zhuǎn)換電路的兩個輸出端，輸出端即為4FSK調(diào)制信號。、二進制基帶信號的串/并轉(zhuǎn)換模塊串/并轉(zhuǎn)換器：移

7、位寄存器除了具有存儲代碼的功能以外，還具有移位功能。所謂移位功能，是指寄存器里存儲的代碼能在移位脈沖的作用下依次左移或右移，因此移位寄存器還可以用來實現(xiàn)數(shù)據(jù)的串并轉(zhuǎn)換。 圖5.1所示電路是由對稱式多諧振蕩器和二位移位寄存器（邊沿觸發(fā)器結(jié)構(gòu)的D觸發(fā)器組成）構(gòu)成的串并轉(zhuǎn)換模塊電路，其中對稱式多諧振蕩器是用來產(chǎn)生串行輸入信號和移位脈沖的，第一個觸發(fā)器的輸入端接收串行輸入信號，其余的每個寄存器輸入端均與前邊一個觸發(fā)器的Q端相連。圖5-1串并轉(zhuǎn)換電路圖 因為從CP上升沿到達開始到輸出端新狀態(tài)的建立需要經(jīng)過一段傳輸延遲時間，所以當CP的上升沿同時作用于所有的觸發(fā)器時，它們輸入端的狀態(tài)還沒有改變。于是按原來

8、的狀態(tài)反轉(zhuǎn)。 例如，在4個時鐘周期內(nèi)輸入代碼依次為1011，而移位寄存器的初始狀態(tài)為=00，那么在移位寄存器的作用下，移位寄存器里的代碼移動情況如表5-2所示，圖5-3給出了各觸發(fā)器在移位過程中的電壓波形圖。 可以看到，經(jīng)過四個CP信號以后，串行輸入的4個代碼全部移入了移位寄存器中，同時在4個觸發(fā)器的輸出端得到了并行輸出的代碼。表5-2 CP的順序輸入的 01234 0 1 0 1 1001 1 01 1圖5.3 電路的電壓波形1101CPtcTttt總電路圖如圖5.4所示圖5.4所示 總電路圖三、4FSk的解調(diào)原理4FSK信號的相干解調(diào)法原理框圖如圖1.2所示。其原理是：4FSK信號先經(jīng)過帶

9、通濾波器去除調(diào)制信號頻帶以外的在信道中混入的噪聲，此后該信號分為四路，每路信號與相應(yīng)載波相乘，再經(jīng)過低通濾波器去除高頻成分，得到包含基帶信號的低頻信號，將其送入抽樣判決器中進行抽樣判決，抽樣判決器的輸出分別得到兩路原基帶信號表示四進制得到原始碼元。4FSK已調(diào)信號帶通濾波器相乘器低通濾波器抽樣判決器解調(diào)信號1載波f1載波f2帶通濾波器帶通濾波器帶通濾波器相乘器相乘器相乘器低通濾波器低通濾波器低通濾波器載波f2載波f2加法器抽樣判決器帶通濾波器加法器圖1.2 4FSK鍵控法解調(diào)原理框圖四、 4FSK 調(diào)制算法分析(1)、將輸入的二進制序列按奇位、偶位進行串并轉(zhuǎn)換。(2)、根據(jù)DMR標準中的符號和

10、比特的對應(yīng)關(guān)系表14，將二進制的0、1序列映射為相應(yīng)的四電平符號流。(3)、將這些符號流每符號插入8個數(shù)值點，并輸入平方根升余弦濾波器進行平滑處理，則可得到輸入調(diào)制信號m(n)。濾波器為平方根升余弦濾波器4，奈奎斯特升余弦濾波器的一部分用于抑制鄰道干擾，另一部分用于接收機抑制噪聲。抑制鄰道干擾濾波器的輸入包含一系列脈沖，這些脈沖之間的間隔為208，33ms(1/4800s)。通過定義根升余弦濾波器的頻率響應(yīng)為奈奎斯特升余弦濾波器的平方根，來定義奈奎斯特升余弦濾波器的分割。濾波器的群延遲在帶通范圍|f|2880Hz內(nèi)是平滑的。濾波器的的幅頻響應(yīng)由下面公式近似給出4：F( f ) =1 當 |f|

11、 1920 HzF( f ) = cos( f / 1920) 當 1920 Hz2880 Hz其中F( f ) 代表平方根升余弦濾波器的幅頻響應(yīng)。則該濾波器的傳輸頻率df 為2880-1920=960Hz，濾波器的等效截止頻率F0 為2400Hz。(4)、將m(n)輸入頻率調(diào)制器進行FM調(diào)制。則可得到4FSK調(diào)制輸出信號。調(diào)頻信號： 式中: kf 為調(diào)頻指數(shù)。將其離散化,在nTs的時間內(nèi)對信號m(n)累加求和,得:由(2)式得該4FSK系統(tǒng)調(diào)制實現(xiàn)方框圖如下:五、4FSK 解調(diào)算法分析解調(diào)過程跟調(diào)制過程恰好相反，將經(jīng)過信道傳輸?shù)浇邮斩说男盘柾ㄟ^頻率解調(diào)器進行解調(diào)。對于同一種數(shù)字調(diào)制信號，采用

12、相干解調(diào)方式的誤碼率低于采用非相干解調(diào)方式的誤碼率，所以這里采用相干解調(diào)方式。將解調(diào)后的信號經(jīng)濾波器后，再經(jīng)抽樣判決則可輸出四電平符號，將其按表1反映射，即可輸出二進制比特。解調(diào)框圖如下所示：解調(diào)中F(f)濾波器的幅頻響應(yīng)與調(diào)制中相同。將(2)式cos部分按三角公式展開得：將其經(jīng)過低通濾波器，濾掉高頻分量，則(8)式剩下1/ 2A* I (nTs)分量，(10)式只剩下1/ 2A*Q(nTs)分量。將(13)式經(jīng)相位校正后，再由(14)即可解調(diào)出m(n)。由于在數(shù)字域內(nèi)，頻率和相位的關(guān)系是簡單的一階差分關(guān)系，如公式(13)所示。在實際中，要想準確實現(xiàn)調(diào)頻信號的解調(diào)，差分鑒頻必須滿足以下要求67

13、：(a)、相鄰的兩個相位差應(yīng)限制在 , 內(nèi)，否則出現(xiàn)相位跳變。(b)、當相位差大于2 時，應(yīng)進行模2 處理。4FSK 相干解調(diào)的實現(xiàn)方框圖如下：最后將解調(diào)后的信號m(n)經(jīng)過匹配濾波器后，再經(jīng)過抽樣判決，可獲得四電平符號，由標準中的表1 所示對應(yīng)關(guān)系，可將四電平符號再反映射為二進制比特流。六、 軟件實現(xiàn)方案按照上節(jié)的算法分析及調(diào)制、解調(diào)方框圖，做出軟件實現(xiàn)方案圖如圖7 所示。 FSK調(diào)制解調(diào)MATLAB源代碼關(guān)鍵詞：FSK 高斯白噪聲 調(diào)制 眼圖 信噪比function FSKFc=10; %載頻Fs=40; %系統(tǒng)采樣頻率Fd=1; %碼速率N=Fs/Fd;df=10;numSymb=20;

14、%進行仿真的信息代碼個數(shù)M=4; %進制數(shù)SNRpBit=80;%信噪比SNR=SNRpBit/log2(M);seed=12345 54321;numPlot=20;x=randsrc(numSymb,1,0:M-1);%產(chǎn)生20個二進制隨機碼figure(1)stem(0:numPlot-1,x(1:numPlot),bx);title(輸入四進制隨機序列)xlabel(Time);ylabel(Amplitude);%調(diào)制y=dmod(x,Fc,Fd,Fs,fsk,M,df);%4FSK調(diào)制函數(shù)numModPlot=numPlot*Fs;t=0:numModPlot-1./Fs;figu

15、re(2)plot(t,y(1:length(t),b-);axis(min(t) max(t) -1.5 1.5);title(調(diào)制信號輸出)xlabel(Time);ylabel(Amplitude);%在已調(diào)信號中加入高斯白噪聲randn(state,seed(2);y=awgn(y,SNR-10*log10(0.5)-10*log10(N),measured,dB);%在已調(diào)信號中加入高斯白噪聲figure(3)plot(t,y(1:length(t),b-);%畫出經(jīng)過信道的實際信號axis(min(t) max(t) -1.5 1.5);title(加入高斯白噪聲后的已調(diào)信號)xl

16、abel(Time);ylabel(Amplitude);%相干解調(diào)figure(4)z1=ddemod(y,Fc,Fd,Fs,fsk/eye,M,df);title(相干解調(diào)后的信號的眼圖)%帶輸出波形的相干M元頻移鍵控解調(diào)figure(5)stem(0:numPlot-1,x(1:numPlot),bx);hold on;stem(0:numPlot-1,z1(1:numPlot),ro);hold off;axis(0 numPlot -0.5 1.5);title(相干解調(diào)輸出信號與原序列的比較)legend(原輸入二進制隨機序列,相干解調(diào)后的信號)xlabel(Time);ylabe

17、l(Amplitude);%非相干解調(diào)figure(6)z2=ddemod(y,Fc,Fd,Fs,fsk/eye/noncoh,M,df);title(非相干解調(diào)后的信號的眼圖)%帶輸出波形的非相干M元頻移鍵控解調(diào)figure(7)stem(0:numPlot-1,x(1:numPlot),bx);hold on;stem(0:numPlot-1,z2(1:numPlot),ro);hold off;axis(0 numPlot -0.5 1.5);title(非相干解調(diào)輸出信號與原序列的比較)legend(原輸入二進制隨機序列,非相干解調(diào)后的信號)xlabel(Time);ylabel(Am

18、plitude);%誤碼率統(tǒng)計errorSym ratioSym=symerr(x,z1);figure(8)simbasebandex(0:1:5);title(相干解調(diào)后誤碼率統(tǒng)計)errorSym ratioSym=symerr(x,z2); figure(9)simbasebandex(0:1:5);title(非相干解調(diào)后誤碼率統(tǒng)計)%濾除高斯白噪聲Delay=3;R=0.5; %滯后3syf,tf=rcosine(Fd,Fs,fir,R,Delay); %升余弦函數(shù)yo2,to2=rcosflt(y,Fd,Fs,filter,yf);%加入高斯白噪聲后的已調(diào)信號和經(jīng)過升余弦濾波器后

19、的已調(diào)信號t=0:numModPlot-1./Fs;figure(10)plot(t,y(1:length(t),r-);hold on;plot(to2,yo2,b-);hold off;axis(0 20 -1.5 1.5);xlabel(Time);ylabel(Amplitude);legend(加入高斯白噪聲后的已調(diào)信號,經(jīng)過升余弦濾波器后的已調(diào)信號)title(升余弦濾波前后波形比較)eyediagram(yo2,N);%眼圖title(加入高斯白噪聲后的已調(diào)信號的眼圖)七、 軟件實現(xiàn)結(jié)果本文給出了MATLAB軟件實現(xiàn)4FSK 調(diào)制解調(diào)過程中的詳細結(jié)果。7.1 調(diào)制過程仿真實現(xiàn)結(jié)果

20、先將截取的二進制比特流按奇、偶位分開，再按照前面表1 所示的映射關(guān)系，將二進制序列映射為四電平符號，得到的四電平符號局部如下圖8 所示：圖 8將二進制比特流映射為四電平符號對四電平符號進行插值濾波，并進行平滑濾波后得到輸入調(diào)制信號m(n)，濾波器即為前面介紹的平方根升余弦濾波器。插值濾波后生成的局部的與上述顯示四電平符號對應(yīng)的輸入調(diào)制信號m(n)顯示如下圖9： 圖 9 輸入調(diào)制信號m(n)做輸入調(diào)制信號m(n)的頻譜(由于重在觀察調(diào)頻前后頻譜的搬移情況，為便于觀察，取頻譜的絕對值，并將零頻時圖形轉(zhuǎn)換至圖的正中來顯示)，如下圖10： 圖 10 輸入調(diào)制信號m(n)的頻譜做出調(diào)制輸出信號的頻譜如下

21、圖11 所示(為便于觀察頻譜的搬移情況，取了頻譜的絕對值，并將零頻時的圖形轉(zhuǎn)換到了圖的正中間顯示，若未轉(zhuǎn)換，頻譜邊帶在載頻+fc 與-fc周圍分布)，顯然可以看到，調(diào)制過程實現(xiàn)了頻譜的搬移，將輸入調(diào)制信號的頻譜由低頻搬移到了中頻fc 處。圖 11調(diào)制輸出信號頻譜7.2 解調(diào)過程仿真實現(xiàn)結(jié)果將調(diào)制輸出信號輸入到傳輸信道進行傳輸。在信道中由于存在干擾，則到達接收端的在調(diào)制輸出信號基礎(chǔ)上附加了高斯白噪聲。下面開始解調(diào)過程。我們采用相干解調(diào)方式對其進行解調(diào)。將接收到信號分別乘上本地載波 與 后通過低通濾波器，濾掉高頻分量，然后對其進行抽樣，可得到Ik2 與Qk2。用Ik2 與Qk2 計算相位，應(yīng)用(1

22、4)式，就可以得到解調(diào)后的序列。顯示如下圖12。 圖 12 解調(diào)后得到的信號做解調(diào)后得到的信號的頻譜(為便于觀察， 取了頻譜的絕對值，并將零頻轉(zhuǎn)換至了圖的正中顯示)，顯示如下圖13： 圖 13 解調(diào)后得到信號的頻譜對產(chǎn)生的解調(diào)后信號每8 點取一個抽樣點，得一新序列，對按下面三個門限+2，0，-2依次進行抽樣判決，即可得到解調(diào)后的四電平符號。解調(diào)出的四電平符號局部顯示如下圖14： 圖 14 解調(diào)后得出的四電平符號按照與表1 相反的過程，將該四電平符號反映射為二進制序列，即可完成解調(diào)過程。 MATLAB仿真眼圖6.3 誤碼率計算輸入二進制比特流與輸出二進制比特流均為二進制序列，則二者對應(yīng)不同的位即

23、為誤碼位，統(tǒng)計出誤碼位的個數(shù)，則可求出此過程的誤碼率來。使信噪比發(fā)生改變，則誤碼率也相應(yīng)發(fā)生變化，可以做出誤碼率與信噪比之間的關(guān)系圖。誤碼率與信噪比關(guān)系如下圖15 所示： 圖 17 信噪比與誤碼率關(guān)系圖由上面誤碼率與信噪比的關(guān)系圖，我們可以看到，在這種相干解調(diào)方式中，隨著信噪比的增大，誤碼率逐漸降低，與實際信道是相吻合的。八、心得體會首先，通過練習(xí)熟悉和掌握通信仿真軟件System View的使用，在練習(xí)的同時我又重新復(fù)習(xí)了通信原理很多知識，進一步加深了對通信系統(tǒng)知識的掌握，把所學(xué)的知識理論與實踐相結(jié)合。其次，在通信系統(tǒng)仿真的過程中，我遇到了很多的問題，我通過自己查資料，請教老師和同學(xué)一一的解

24、決遇到的各種問題。我不僅學(xué)會獨立思考問題、分析問題、解決問題，而且加深了對理論知識的掌握，還增強了自己的動手能力。再者，對于信息工程專業(yè)的學(xué)生來說，通信系統(tǒng)的設(shè)計是我們的很重要專業(yè)課，我很認真的完成各項設(shè)計，這次實習(xí)中學(xué)到的知識對我們增強專業(yè)知識的幫助很大。 附錄資料：不需要的可以自行刪除各類濾波器的MATLAB程序理想低通濾波器IA=imread(lena.bmp);f1,f2=freqspace(size(IA),meshgrid);Hd=ones(size(IA);r=sqrt(f1.2+f2.2);Hd(r0.2)=0;Y=fft2(double(IA);Y=fftshift(Y);Y

25、a=Y.*Hd;Ya=ifftshift(Ya);Ia=ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA);subplot(2,2,2),imshow(uint8(Ia);figuresurf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong); 二、理想高通濾波器IA=imread(lena.bmp);f1,f2=freqspace(size(IA),meshgrid);Hd=ones(size(IA);r=sqrt(f1.2+f2.2);Hd(r0.2)=0;Y=fft2(double(I

26、A);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA);subplot(2,2,2),imshow(uint8(Ia);figuresurf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong); Butterworth低通濾波器IA=imread(lena.bmp);f1,f2=freqspace(size(IA),meshgrid);D=0.3;r=f1.2+f2.2;n=4;for i=1:s

27、ize(IA,1) for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=1/(tn+1); endendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA);subplot(2,2,2),imshow(uint8(Ia);figuresurf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong); Butterworth高通濾波器IA

28、=imread(lena.bmp);f1,f2=freqspace(size(IA),meshgrid);D=0.3;r=f1.2+f2.2;n=4;for i=1:size(IA,1) for j=1:size(IA,2) t=(D*D)/r(i,j); Hd(i,j)=1/(tn+1); endendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA);subplot(2,2,2),imshow(uint8(Ia);f

29、iguresurf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong); 高斯低通濾波器IA=imread(lena.bmp);IB=imread(babarra.bmp);f1,f2=freqspace(size(IA),meshgrid);D=100/size(IA,1);r=f1.2+f2.2;Hd=ones(size(IA);for i=1:size(IA,1) for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=exp(-t); endendY=fft2(double(IA);Y=fftshi

30、ft(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA);subplot(2,2,2),imshow(uint8(Ia);figuresurf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong); 高斯高通濾波器IA=imread(lena.bmp);IB=imread(babarra.bmp);f1,f2=freqspace(size(IA),meshgrid);%D=100/size(IA,1);D=0.3;r=f1

31、.2+f2.2;for i=1:size(IA,1) for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=1-exp(-t); endendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA);subplot(2,2,2),imshow(uint8(Ia);figuresurf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong);

32、梯形低通濾波器IA=imread(lena.bmp);IB=imread(babarra.bmp);f1,f2=freqspace(size(IA),meshgrid);%D=100/size(IA,1);D0=0.1;D1=0.4;r=sqrt(f1.2+f2.2);Hd=zeros(size(IA);Hd(r=D0 & r(i,j)=D1 Hd(i,j)=(D1-r(i,j)/(D1-D0); end endendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot

33、(2,2,1),imshow(uint8(IA);subplot(2,2,2),imshow(uint8(Ia);figuresurf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong); 梯形高通濾波器IA=imread(lena.bmp);IB=imread(babarra.bmp);f1,f2=freqspace(size(IA),meshgrid);%D=100/size(IA,1);D0=0.1;D1=0.4;r=sqrt(f1.2+f2.2);Hd=ones(size(IA);Hd(r=D0 & r(i,j)=D1 Hd(i

34、,j)=(D0-r(i,j)/(D0-D1); end endendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA);subplot(2,2,2),imshow(uint8(Ia);figuresurf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong); 用其他方法編寫的理想低通、理想高通、Butterworth低通、同態(tài)濾波程序理想低通i1=imread(lena.bmp);i2=imnoise(i1,salt & pepper,0.1);f=double(i2);k=fft2(f);g=fftshift(k);N1,N2=size(g);d0=50;u0=floor(N1/2)+1;v0=floor(N2/2)+1;for i=1:N1 for j=1:N2 d=sqrt(i-u0)2+(j-v0)2); if d=d0 h=1; else h=0; end y(i,j)=g(i,j)*h; endendy=ifftshift(y);E1=ifft2(y);E2=real(E1);figuresubplot(2,2