信號(hào)分析與處理實(shí)驗(yàn)報(bào)告

實(shí)驗(yàn)一圖像信號(hào)頻譜分析及濾波一：實(shí)驗(yàn)原理FFT不是一種新的變化，而是DFT的快速算法?？焖俑道锶~變換能減少運(yùn)算量的根本因素在于它不停地把長(zhǎng)序列的離散傅里葉變換變?yōu)槎绦蛄械碾x散傅里葉變換，在運(yùn)用的對(duì)稱性和周期性使DFT運(yùn)算中的有些項(xiàng)加以合并，達(dá)成減少運(yùn)算工作量的效果。為了消除或削弱噪聲，提取有用信號(hào)，必須進(jìn)行濾波，能實(shí)現(xiàn)濾波功效的系統(tǒng)成為濾波器。按信號(hào)可分為模擬濾波器和數(shù)字濾波器兩大類。數(shù)字濾波器的核心是如何根據(jù)給定的技術(shù)指標(biāo)來(lái)得到能夠?qū)崿F(xiàn)的系統(tǒng)函數(shù)。從模擬到數(shù)字的轉(zhuǎn)換辦法諸多，慣用的有雙線性變換法和沖擊響應(yīng)不變法，本實(shí)驗(yàn)重要采用雙線性變換法。雙線性變換法是一種由s平面到z平面的映射過(guò)程，其變換式定義為：數(shù)字域頻率與模擬頻率之間的關(guān)系是非線性關(guān)系。雙線性變換的頻率標(biāo)度的非線性失真是能夠通過(guò)預(yù)畸變的辦法去賠償?shù)?。變換公式有Ωp=2/T*tan(wp/2)Ωs=2/T*tan(ws/2)二：實(shí)驗(yàn)內(nèi)容1．圖像信號(hào)的采集和顯示選擇一副不同彩色圖片，運(yùn)用Windows下的畫圖工具，設(shè)立成200*200像素格式。然后在Matlab軟件平臺(tái)下，運(yùn)用有關(guān)函數(shù)讀取數(shù)據(jù)和顯示圖像。規(guī)定顯示出原始灰度圖像、加入噪聲信號(hào)后的灰度圖像、濾波后的灰度圖像。2．圖像信號(hào)的頻譜分析規(guī)定分析和畫出原始灰度圖像、加入噪聲信號(hào)后灰度圖像、濾波后灰度圖像信號(hào)的頻譜特性。3.數(shù)字濾波器設(shè)計(jì)給出數(shù)字低通濾波器性能指標(biāo):通帶截止頻率fp＝10000Hz，阻帶截止頻率fs＝15000Hz，阻帶最小衰減Rs＝50dB，通帶最大衰減Rp＝3dB，采樣頻率40000Hz。三：實(shí)驗(yàn)程序clearallx=imread('D:\lan.jpg');%原始彩色圖像的數(shù)據(jù)讀取x1=rgb2gray(x);%彩色圖像值轉(zhuǎn)化為灰度圖像值[M,N]=size(x1);%數(shù)據(jù)x1的長(zhǎng)度,用來(lái)求矩陣的大小x2=im2double(x1);%unit8轉(zhuǎn)化為double型x3=numel(x2);%計(jì)算x2長(zhǎng)度f(wàn)igure(1);subplot(1,3,1);imshow(x2);title('原始灰度圖')z1=reshape(x2,1,x3);%將二維數(shù)據(jù)轉(zhuǎn)化成一維數(shù)據(jù)g=fft(z1);%對(duì)圖像進(jìn)行二維傅里葉變換mag=fftshift(abs(g));%fftshift是針對(duì)頻域的，將FFT的DC分量移到頻譜中心K=40000;Fs=40000;dt=1/Fs;n=0:K-1;f1=18000;z=0.1*sin(2*pi*f1*n*dt);x4=z1+z;%加入正弦噪聲f=n*Fs/K;y=fft(x4,K);z2=reshape(x4,M,N);%將一維圖轉(zhuǎn)換為二維圖subplot(1,3,2);imshow(z2);title('加入噪聲后')g1=fft(x4);mag1=fftshift(abs(g1));%設(shè)計(jì)濾波器ws=0.75*pi;wp=0.5*pi;fs=10000;wp1=2*fs*tan(wp/2);ws1=2*fs*tan(ws/2);rs=50;rp=3;%[n,wn]=buttord(wp/pi,ws/pi,rp,rs);%[bz,az]=butter(n,wn);[n,wn]=buttord(wp1,ws1,rp,rs,'s');[z,p,k]=buttap(n);[b,a]=zp2tf(z,p,k);[B,A]=lp2lp(b,a,wn);[bz,az]=bilinear(B,A,fs);[h,w]=freqz(bz,az,128,fs);L=numel(z2);z3=reshape(z2,1,L);x6=filter(bz,az,double(z3));x7=reshape(x6,M,N);subplot(1,3,3);imshow(x7);g2=fft(x6);mag2=fftshift(abs(g2));title('濾波后')%建立頻譜圖figure(2);subplot(1,3,1);plot(mag);title('原始Magnitude')subplot(1,3,2);plot(mag1);title('加噪聲Magnitude')subplot(1,3,3);plot(mag2);title('濾波后Magnitude')figure(3);subplot(1,2,1)plot(w,abs(h));xlabel('f');ylabel('h');title('濾波器幅譜');subplot(1,2,2);plot(w,angle(h));title('濾波器相譜');四：實(shí)驗(yàn)成果與分析圖一圖二分析：由圖二能夠懂得加入噪聲后的幅值譜和原始圖的幅值譜明顯多了兩條幅值