24小时热门版块排行榜    

查看: 765  |  回复: 3

mymarik

铁杆木虫 (小有名气)

[求助] 用频谱法求分形维数,这里的频谱要怎么求呢(图片里有具体方法说明,但不是很理解)

这里的P(w)是要怎么求呢 ,下面这样求对吗?

testD=Zcunew(:,20); %Zcunew(:20)是一个分形序列点

Pw=[];
w=[];
ZcuPw=[];
absZcuPw=[];

ZcuPw=fft(testD);
absZcuPw=abs(ZcuPw);

for i=1:100
    Pw(i)=absZcuPw(i);
    w(i)=i;
end

    lgPw=log10(Pw);
    lgw=log10(w);
    lgPwT=lgPw';
    lgwT=lgw';
    plot(lgw,lgPw,'-*');
    hold on
    [p1,p2]=polyfit(lgw,lgPw,1);
   
    for i=1:100
        plot(lgw(i),p1(1)*lgw(i)+p1(2),'b-o')
    end
   
    kk(ii)=p1(1);
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mymarik

铁杆木虫 (小有名气)

2楼2017-05-16 22:17:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mymarik

铁杆木虫 (小有名气)

送红花一朵
图片在二楼

发自小木虫IOS客户端
3楼2017-05-16 22:18:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mymarik

铁杆木虫 (小有名气)

另外图片里 推倒的 D=2.5-a 是不是错的 ,我算出来是(5+a)/2

发自小木虫IOS客户端
4楼2017-05-17 10:10:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mymarik 的主题更新
信息提示
请填处理意见