各位大侠 我看不懂下面的程序 能帮我解释一下吗 谢谢啦

jiaxiangning

function p = gfp(beta,deta,n,e)
%Function to generate fractal profile using IFT.
%refer to Gallant,Moore,Hutchinson and Gessler,1994
%Input:
%       beta: beta=2H+1
%       deta: the sample interval of time series
%       n: number of time series
%       e:stadard devation
%Writen by Zhang Jinling
%Institute of geology and geophysics china academy of sciences
pi=3.1415926;
f=(1:n)/deta/n;
seta=2*pi*rand(1,n);
amp=e^2*f.^(-1*beta/2);
fp=amp.*exp(sqrt(-1)*seta);
p=real(ifft(fp,n));
figure;plot(p);
end

jiaxiangning

回复 1 # jiaxiangning 的帖子

运行了一下  总是出错

jiaxiangning

回复 2 # jiaxiangning 的帖子

f=(1:n)/deta/n;   就看不懂了

ChaChing

...总是出错...

求助完整格式：出错代码和出错提示

还有LZ这样问不理想, 建议看看

5.建议提问的网友分清 编程问题 和 专业问题 http://forum.vibunion.com/thread-36746-1-1.html
6.提问的智慧！！！！(发帖前请认真阅读) http://forum.vibunion.com/thread-21991-1-1.html

个人水平有限, LZ没相关说明, 仅猜测好像在处理信号!?:@)

f=(1:n)/deta/n;

求相对应频率??

321forever

回复 2 # jiaxiangning 的帖子

不是很懂，但，fp=amp.*exp(sqrt(-1)*seta); 是不是对-1求根，还是对 -1*seta求根