各位大侠 我看不懂下面的程序 能帮我解释一下吗 谢谢啦
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 回复 1 # jiaxiangning 的帖子
运行了一下总是出错
回复 2 # jiaxiangning 的帖子
f=(1:n)/deta/n; 就看不懂了
本帖最后由 ChaChing 于 2011-12-29 22:36 编辑
...总是出错...求助完整格式：出错代码和出错提示
还有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;
求相对应频率??
回复 2 # jiaxiangning 的帖子
不是很懂，但，fp=amp.*exp(sqrt(-1)*seta); 是不是对-1求根，还是对 -1*seta求根
页:
[1]