Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

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

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

求助]fibonacci数列的置乱程序

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

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

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

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

FIBONACCI变换

定义1 设F0=1,F1=1,Fn=Fn-1+ Fn-2,n\3,则数列{Fn}称为Fibonacci数列" 定义2 设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

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

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

还是不明白你的变换具体是什么回事，见笑了 不过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

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

你要的应该就是这个伪随机序列吧，看红色标注 function x = rndgen(n,type,options) %RNDGEN Uniformly distributed random numbers generator. % TYPE = 'Mul' Linear congruential random number generators % OPTIONS = % 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 = ; 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 = ; 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 = ; end otherwise x = []; error('Illegal type!'); end

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

这么复杂啊，我研究研究，谢谢了

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

a=imread('c:\baboon64.bmp'); b=imresize(a,); 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

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

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

游鱼儿 发表于 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

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

可是我按照公式写出的程序运行结果没什么变化，显示的音频和原来的一样 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

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

<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>
 你的音频和图像数据都是向量？不太可能吧，图片一般都是二维或者三维的

页: [1]

声振论坛's Archiver

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

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

FIBONACCI变换

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

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

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

融合技术

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

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

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