[求助]fibonacci数列的置乱程序

游鱼儿 · 发表于 2006-3-22 16:15

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

求助]fibonacci数列的置乱程序[em10]

happy · 发表于 2006-3-22 21:33

fibonacci数列是不是满足 a[k+2]=a[k]+a[k+1]，k=1,2,.. ；且a[1]=a[2]=1 的数列？ 另外置乱是什么意思？

游鱼儿 · 发表于 2006-3-23 18:33

定义1 设F0=1,F1=1,Fn=Fn-1+ Fn-2,n\3,则数列{Fn}称为Fibonacci数列" 定义2[9] 设Fn,Fn+1是Fibonacci数列中的两 个相邻数,称以下变换: Sk=(kFn+r)modFn+1 为Fibonacci变换"这里r可作算法中密钥的一部 分:r可取0 k=0,1,2,,,Fn+1-1 由定义可以看出,该变换可将数列{Q}=(0,1,2, ,,Fn+1-1)变换成另一新数列{S}=(S0,S1,S2, ,,SFn+1-1),且可以证明数列{S}是数列{Q}的一个 新的伪随机置换 用Fibonacci变换将水印图像的像素位置打 乱;将水印图像的行列分别作变换,选取 Fn+1为载体图像的宽Fn+1和Fn+1高,置乱算法可表述 如下:Xi=(i*Fn)modFn+1 Yi=（i*Fn）modFn+1 这样做的目的是将水印图像的像素扩散到载体 图像大小的空间域上,使得水印图像像素有更广泛 的分散空间,且置乱后的水印图像与载体图像尺寸 相同"例如数列{0,1,2,3}可用计算式Xi=(i@ 89)mod144分散到144的空间上为{0，89，34，123} 还有我没讲明白的你再问，拜托，拜托，帮帮忙，我做毕业设计，不会用MATLAB[em04]

happy · 发表于 2006-3-24 08:04

还是不明白你的变换具体是什么回事，见笑了 不过FIBONACCI数列是搞清楚了，可以有下面两种方式产生 function out=fibo1(n) r1=(1+sqrt(5))/2; r2=(1-sqrt(5))/2; out=(r1^(n-1)-r2^(n-1))/sqrt(5); function out=fibo2(n) if n==1 out=0; return; elseif n==2 out=1; return; else out=fibo(n-1)+fibo(n-2); end

游鱼儿 · 发表于 2006-3-24 14:07

变换就是：产生一个随机的64个数的数列，然后用公式 K*89MOD144做变换（K取1到64）来确定K位置上的数到被换到的新的位置，结果就是把这64个数伪随机的散列到了144的空间上，例如，K=1时，1*89MOD144=89，既第一个位置上的数要到89位置上去

happy · 发表于 2006-3-25 08:57

你要的应该就是这个伪随机序列吧，看红色标注 function x = rndgen(n,type,options) %RNDGEN Uniformly distributed random numbers generator. % TYPE = 'Mul' Linear congruential random number generators % OPTIONS = [a,c,m,x0] % X is a random number vector get from % X = (a*X+c) mod m % X0 is the inital value of X which is also called SEED. % For example IMSL uses it with m=2^31-1, a=16807 and c=0; % a=84589, c=45989 and m=217728 will give a random seriers % which will overflow at 2^35. % % TYPE = 'Add' Fibonacci random number generators % OPTIONS = [,m,x0,x1] % X is a random number vector get from % X(i) = (X(i-1)+X(i-2)) mod m % X0 and X1 are inital values of X which are also called SEED. % % See also RAND, RNDGENM %--------------------------------------------------------------------------- % xywu % Fudan university. % % $Revision: 1.0 $ % $Date: 03/08/2002 $ % %--------------------------------------------------------------------------- if nargin>3 error('Too many arguments! '); end if nargin<1 error('Requires at least one input argument.'); end if nargin==1 type = 'Mul'; options = [16807,0,2^31-1,fix(rand(1)*now*62)]; end switch type case 'Mul' a = options(1); c = options(2); m = options(3); xt = options(4); x = []; for i = 1:n xt = mod(a*xt+c,m); x = [x,xt/m]; end x = x(end-n+1:end); case 'Add' m = options(1); xt1 = options(2); xt2 = options(3); x = []; for i = 1:n xt1 = mod(xt1+xt2,m); xt2 = mod(xt1+xt2,m); x = [x,xt1/m,xt2/m]; end otherwise x = []; error('Illegal type!'); end

游鱼儿 · 发表于 2006-3-26 20:30

这么复杂啊，我研究研究，谢谢了[em07]

游鱼儿 · 发表于 2006-3-27 13:55

a=imread('c:\baboon64.bmp'); b=imresize(a,[64,64]); c=round(double(b)./256); subplot(331) imshow(c) new=c; for i=1:64 for j=1:64 k=j*89; b=rem(k,144); new(i,b)=c(i,j); end end for j=1:64 i=1:64 m=i*89; n=rem(m,144); new(n,j)=c(i,j); end end subplot(333) imshow(new) 看看这段程序。我的大概意思就是这样，就是要把图C 的各点分散了，成为一个新的图NEW 但是循环语句那有错误，麻烦看一下，错在哪了

happy · 发表于 2006-3-27 18:18

可以运行，请说明问题所在

游鱼儿 · 发表于 2006-4-2 17:36

上个问题解决了，十分感谢 请教个新问题 利用融合技术公式S=T*W+（1-T）*G T取0。1 W为水印图象，G为一段音频，求S的程序原代码帮忙写一下，谢了

happy · 发表于 2006-4-3 16:27

<DIV class=quote>以下是引用游鱼儿在2006-4-2 17:36:39的发言： 上个问题解决了，十分感谢 请教个新问题 利用融合技术公式S=T*W+（1-T）*G T取0。1 W为水印图象，G为一段音频，求S的程序原代码帮忙写一下，谢了</DIV>
这个问题应该只是一个公式运算啊，需要找源代码吗？

游鱼儿 · 发表于 2006-4-5 10:55

可是我按照公式写出的程序运行结果没什么变化，显示的音频和原来的一样 for i=1:20736 m=0.1*yiwei(i); n=0.9*wav(i); newwav(i)=m+n; end 
 这样的循环语句就可以了吗？

happy · 发表于 2006-4-5 14:39

<DIV class=quote>以下是引用游鱼儿在2006-4-5 10:55:01的发言： 
可是我按照公式写出的程序运行结果没什么变化，显示的音频和原来的一样 for i=1:20736 m=0.1*yiwei(i); n=0.9*wav(i); newwav(i)=m+n; end 
 这样的循环语句就可以了吗？</DIV>
 你的音频和图像数据都是向量？不太可能吧，图片一般都是二维或者三维的

账号		自动登录	找回密码
密码			我要加入

[求助]fibonacci数列的置乱程序

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（游鱼儿）[求助]fibonacci数列的置乱程序

FIBONACCI变换

回复：（游鱼儿）[求助]fibonacci数列的置乱程序

回复：（游鱼儿）[求助]fibonacci数列的置乱程序

回复：（游鱼儿）[求助]fibonacci数列的置乱程序

融合技术

回复：（游鱼儿）融合技术

回复：（happy）回复：（游鱼儿）融合技术

回复：（游鱼儿）回复：（happy）回复：（游鱼儿）融合...