实验3-FFT变换信息隐藏算法

信息隐藏技术实验报告实验名称实验三快速傅立叶变换信息隐藏算法日期，2013-05-23专业课11一级校长实验者学生身份证一、实验目的(1)了解信息隐藏算法的分类方法和基础(2)了解频域变换信息隐藏算法的基本思想(3)了解离散傅立叶变换技术的基本理论(4)掌握离散傅立叶变换信息隐藏算法二、实验内容载体图像是8位的bmp图像，嵌入的秘密信息是从屏幕上随机输入的英文字母。需要对载体图像Lena.bmp进行fft变换，将秘密信息转换成二进制流并嵌入载体图像的中频区域，以显示原始载体图像、待嵌入的秘密信息及其对应的二进制流，以及嵌入有秘密信息的伪装载体。(编程语言不受限制)三、实验步骤和设计思路1.用Matlab生成正弦波、矩形波和白噪声信号，并显示它们的时域波形图2.进行快速傅立叶变换以显示各自的频谱图，其中选择采样率、频率和数据长度。3.制作上述三种信号的均方根谱、功率谱和对数均方根谱4.通过IFFT傅里叶逆变换恢复信号，并显示恢复的正弦信号的时域波形四.程序清单fs=100%设置采样频率N=128n=0:N-1；t=n/fs。f0=10%设置正弦信号频率%生成正弦信号x=sin(2 * pi * F0 * t)；图(1)；子情节(231)；图(t，x)；%作为正弦信号的时域波形XL Abel(t)；y label(y)；标题(正弦信号y=2*pi*10t时域波形)；网格；%快速傅里叶变换和频谱图y=fft(x，N)；% fft变换mag=ABS(y)；%寻道幅度f=(0: length(y)-1)* fs/长度(y)；%用于相应的频率转换图(1)；子情节(232)；绘图(f，mag)；% do光谱图轴(0，100，0，80)；Xlabel(频率(赫兹)；伊拉贝尔(振幅)；标题(正弦信号y=2 * 10t频谱N=128)；网格；%均方根谱sq=ABS(y)；图(1)；子情节(233)；地块(f，sq)；Xlabel(频率(赫兹)；均方根谱；标题(正弦信号y=2*pi*10t均方根频谱)；网格；%功率谱功率=sq。2；图(1)；子情节(234)；图(f，功率)；Xlabel(频率(赫兹)；功率谱；标题(正弦信号的功率谱y=2 * pi * 10t)；网格；%对数光谱ln=对数(sq)；图(1)；子情节(235)；图(f，ln)；Xlabel(频率(赫兹)；Ylabel(对数谱)；标题(正弦信号的对数频谱y=2 * pi * 10t)；网格；%使用IFFT恢复原始信号xif ft=IFFT(y)；magx=实数(xifft)；ti=0: length(xif ft)-1/fs；图(1)；子情节(236)；图(ti，magx)；XL Abel(t)；y label(y)；标题(IFFT转换的正弦信号波形)；网格；% * * * * * * * * * * * * * * 2。矩形波* * * * * * * * * * * * *%fs=10%设置采样频率t=-5:0.1:5；x=矩形脉冲(t，2)；x=x(1:99)；图(2)；子情节(231)；图(t(1:99)，x)；%作为矩形波的时域波形XL Abel(t)；y label(y)；标题(矩形波时域波形)；网格；%快速傅里叶变换和频谱图y=FFT(x)；% fft变换mag=ABS(y)；%寻道幅度f=(0: length(y)-1)* fs/长度(y)；%用于相应的频率转换图(2)；子情节(232)；绘图(f，mag)；% do光谱图Xlabel(频率(赫兹)；伊拉贝尔(振幅)；标题(矩形振幅谱)；网格；%均方根谱sq=ABS(y)；图(2)；子情节(233)；地块(f，sq)；Xlabel(频率(赫兹)；均方根谱；标题(矩形波的均方根谱)；网格；%功率谱功率=sq。2；图(2)；子情节(234)；图(f，功率)；Xlabel(频率(赫兹)；功率谱；标题(矩形波的功率谱)；网格；%对数光谱ln=对数(sq)；图(2)；子情节(235)；图(f，ln)；Xlabel(频率(赫兹)；Ylabel(对数谱)；标题(矩形波的对数谱)；网格；%使用IFFT恢复原始信号xif ft=IFFT(y)；magx=实数(xifft)；ti=0: length(xif ft)-1/fs；图(2)；子情节(236)；图(ti，magx)；XL Abel(t)；y label(y)；标题(IFFT转换的矩形波波形)；网格；% * * * * * * * * * * * * * * 3。白噪声* * * * * * * * * * * * * * *%fs=10%设置采样频率t=-5:0.1:5；x=零(1，100)；x(50)=100000；图(3)；子情节(231)；图(t(1:100)，x)；%作为白噪声时域波形XL Abel(t)；y label(y)；标题(白噪声时域波形)；网格；%快速傅里叶变换和频谱图y=FFT(x)；% fft变换mag=ABS(y)；%寻道幅度f=(0: length(y)-1)* fs/长度(y)；对应频率c的百分比Xlabel(频率(赫兹)；伊拉贝尔(振幅)；标题(白噪声幅度谱)；网格；%均方根谱sq=ABS(y)；图(3)；子情节(233)；地块(f，sq)；Xlabel(频率(赫兹)；均方根谱；标题(白噪声的均方根谱)；网格；%功率谱功率=sq。2；图(3)；子情节(234)；图(f，功率)；Xlabel(频率(赫兹)；功率谱；标题(白噪声功率谱)；网格；%对数光谱ln=对数(sq)；图(3)；子情节(235)；图(f，ln)；Xlabel(频率(赫兹)；Ylabel(对数谱)；标题(白噪声对数谱)；网格；%使用IFFT恢复原始信号xif ft=IFFT(y)；magx=实数(xifft)；ti=0: length(xif ft)-1/fs；图(3)；子情节(236)；图(ti，magx)；XL Abel(t)；y label(y)；标题