# MATLAB代做|FPGA代做|Python代做--f形谱宽ｄｅｌｔａa 谱差

### 时间：2019-1-13 2:59:04 点击：

核心提示：MATLAB代做|FPGA代做|Python代做--f形谱宽ｄｅｌｔａa 谱差...

% clc
clear
a=rgb2gray(a);
x=double(a);
%   x1=x(65:192,65:192)
x1=x/sum(sum(x));
for q=-60:2:60
w=(q+60)/2+1;

for k=1:8
xe(k)=log(2^k/128);
fun=inline('sum(sum(Y))');
b=blkproc(x1,[2^k,2^k],fun);
a1=b.^q;
b1=sum(sum(a1));
x12=b;
[row vic]=find(x12==0);
ss=size(row);
for i=1:ss
x12(row(i),vic(i))=1;
end
a2=log(x12);

a3=a1.*a2;
b3=sum(sum(a3));
rq(k)=b3/b1;
fq(k)=b3/b1*q-log(b1);
end
Sxx=sum(xe.^2)-(sum(xe))^2/8;

Sxr=xe*rq'-sum(xe)*sum(rq)/8;
Sxf=xe*fq'-sum(xe)*sum(fq)/8;

r(w)=Sxr/Sxx;
f(w)=Sxf/Sxx;

end
q1=-60:2:60;
% figure,plot(r,f);
q=-60:2:60;
% figure,plot(q,r);%%%%%%%%%%%%%ｑ和阿尔法a
% figure,plot(q,f);%%%%%%%%%%%%%ｑ和f(a)

%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf
%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf
%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf
rmax=max(r);
rmin=min(r);
deltar=(rmax-rmin)

% fmax=max(f);
% fmin=min(f);
% (fmax-fmin)

ffmin=find(r==rmin);
fmin=f(ffmin);

ffmax=find(r==rmax);
fmax=f(ffmax);

deltaf=fmin-fmax%%%%%%这个才是真正的谱差啊啊啊啊啊啊啊

%      amin=find(isnan(hurstH));
%      hurstH(hursti)=[];

%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf
%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf
%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf
figure,plot(r,f);
xlabel('奇异指数a'),ylabel('多重分行谱f(a)')
% q=-60:2:60;
% figure,plot(q,r);%%%%%%%%%%%%%ｑ和阿尔法a
% xlabel('权重因子q'),ylabel('奇异指数a')
% figure,plot(q,f);%%%%%%%%%%%%%ｑ和f(a)
% xlabel('权重因子q'),ylabel('多重分行谱f(a)')

% ff=find(isnan(f));
% f(ff)=[]

QQ ：1224848052

• 百度搜索
• 查阅资料过程中
• 论坛发现
• 百度贴吧发现
• 朋友介绍
• 没有相关文章
• 大名：
• 内容：