實(shí)驗(yàn)3-FFT變換信息隱藏算法(共6頁)

1、信息隱藏技術(shù) 實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)名稱實(shí)驗(yàn)3 FFT變換信息隱藏算法日期2013-05-23專業(yè)班級(jí)11信安1班實(shí)驗(yàn)人學(xué)號(hào)一、實(shí)驗(yàn)?zāi)康模?）了解信息隱藏算法的分類方式和分類依據(jù)（2）理解頻域變換信息隱藏算法的基本思想（3）理解DFT技術(shù)的基本理論（4）掌握DFT信息隱藏算法二、實(shí)驗(yàn)內(nèi)容載體圖像為8位bmp圖像Lena.bmp，嵌入的秘密信息為從屏幕上隨機(jī)輸入的英文字母，要求對(duì)載體圖像Lena.bmp進(jìn)行fft變換，將秘密信息轉(zhuǎn)換成二進(jìn)制流并嵌入到載體圖像的中頻區(qū)域，顯示原載體圖像、需要嵌入的秘密信息及其相應(yīng)的二進(jìn)制流、嵌入了秘密信息的偽裝載體。（編程語言不限）三、實(shí)驗(yàn)步驟和設(shè)計(jì)思想1用Matlab產(chǎn)生

2、正弦波,矩形波,以及白噪聲信號(hào)，并顯示各自時(shí)域波形圖2進(jìn)行FFT變換，顯示各自頻譜圖，其中采樣率，頻率、數(shù)據(jù)長度自選3做出上述三種信號(hào)的均方根圖譜,功率圖譜,以及對(duì)數(shù)均方根圖譜4用IFFT傅立葉反變換恢復(fù)信號(hào)，并顯示恢復(fù)的正弦信號(hào)時(shí)域波形圖四、程序清單fs=100;%設(shè)定采樣頻率N=128;n=0:N-1;t=n/fs;f0=10;%設(shè)定正弦信號(hào)頻率%生成正弦信號(hào)x=sin(2*pi*f0*t);figure(1);subplot(231);plot(t,x);%作正弦信號(hào)的時(shí)域波形xlabel('t');ylabel('y');title('正弦信號(hào)y

3、=2*pi*10t時(shí)域波形');grid;%進(jìn)行FFT變換并做頻譜圖y=fft(x,N);%進(jìn)行fft變換mag=abs(y);%求幅值f=(0:length(y)-1)'*fs/length(y);%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換figure(1);subplot(232);plot(f,mag);%做頻譜圖axis(0,100,0,80);xlabel('頻率(Hz)');ylabel('幅值');title('正弦信號(hào)y=2*pi*10t幅頻譜圖N=128');grid;%求均方根譜sq=abs(y);figure(1);subplot

4、(233);plot(f,sq);xlabel('頻率(Hz)');ylabel('均方根譜');title('正弦信號(hào)y=2*pi*10t均方根譜');grid;%求功率譜power=sq.2;figure(1);subplot(234);plot(f,power);xlabel('頻率(Hz)');ylabel('功率譜');title('正弦信號(hào)y=2*pi*10t功率譜');grid;%求對(duì)數(shù)譜ln=log(sq);figure(1);subplot(235);plot(f,ln);xlab

5、el('頻率(Hz)');ylabel('對(duì)數(shù)譜');title('正弦信號(hào)y=2*pi*10t對(duì)數(shù)譜');grid;%用IFFT恢復(fù)原始信號(hào)xifft=ifft(y);magx=real(xifft);ti=0:length(xifft)-1/fs;figure(1);subplot(236);plot(ti,magx);xlabel('t');ylabel('y');title('通過IFFT轉(zhuǎn)換的正弦信號(hào)波形');grid;%*2.矩形波*%fs=10;%設(shè)定采樣頻率t=-5:0.1:5;x=

6、rectpuls(t,2);x=x(1:99);figure(2);subplot(231);plot(t(1:99),x);%作矩形波的時(shí)域波形xlabel('t');ylabel('y');title('矩形波時(shí)域波形');grid;%進(jìn)行FFT變換并做頻譜圖y=fft(x);%進(jìn)行fft變換mag=abs(y);%求幅值f=(0:length(y)-1)'*fs/length(y);%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換figure(2);subplot(232);plot(f,mag);%做頻譜圖xlabel('頻率(Hz)');y

7、label('幅值');title('矩形波幅頻譜圖');grid;%求均方根譜sq=abs(y);figure(2);subplot(233);plot(f,sq);xlabel('頻率(Hz)');ylabel('均方根譜');title('矩形波均方根譜');grid;%求功率譜power=sq.2;figure(2);subplot(234);plot(f,power);xlabel('頻率(Hz)');ylabel('功率譜');title('矩形波功率譜'

8、);grid;%求對(duì)數(shù)譜ln=log(sq);figure(2);subplot(235);plot(f,ln);xlabel('頻率(Hz)');ylabel('對(duì)數(shù)譜');title('矩形波對(duì)數(shù)譜');grid;%用IFFT恢復(fù)原始信號(hào)xifft=ifft(y);magx=real(xifft);ti=0:length(xifft)-1/fs;figure(2);subplot(236);plot(ti,magx);xlabel('t');ylabel('y');title('通過IFFT轉(zhuǎn)換的矩形波

9、波形');grid;%*3.白噪聲*%fs=10;%設(shè)定采樣頻率t=-5:0.1:5;x=zeros(1,100);x(50)=100000;figure(3);subplot(231);plot(t(1:100),x);%作白噪聲的時(shí)域波形xlabel('t');ylabel('y');title('白噪聲時(shí)域波形');grid;%進(jìn)行FFT變換并做頻譜圖y=fft(x);%進(jìn)行fft變換mag=abs(y);%求幅值f=(0:length(y)-1)'*fs/length(y);%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換figure(3);subp

10、lot(232);plot(f,mag);%做頻譜圖xlabel('頻率(Hz)');ylabel('幅值');title('白噪聲幅頻譜圖');grid;%求均方根譜sq=abs(y);figure(3);subplot(233);plot(f,sq);xlabel('頻率(Hz)');ylabel('均方根譜');title('白噪聲均方根譜');grid;%求功率譜power=sq.2;figure(3);subplot(234);plot(f,power);xlabel('頻率(Hz

11、)');ylabel('功率譜');title('白噪聲功率譜');grid;%求對(duì)數(shù)譜ln=log(sq);figure(3);subplot(235);plot(f,ln);xlabel('頻率(Hz)');ylabel('對(duì)數(shù)譜');title('白噪聲對(duì)數(shù)譜');grid;%用IFFT恢復(fù)原始信號(hào)xifft=ifft(y);magx=real(xifft);ti=0:length(xifft)-1/fs;figure(3);subplot(236);plot(ti,magx);xlabel('t')