跪求————最大熵功率谱估计的matlab代码,..

来源：学生作业帮编辑：搜搜考试网作业帮分类：综合作业时间：2024/06/26 01:20:10

跪求————最大熵功率谱估计的matlab代码,哪位高人伸...
本熵代码
function [shang]=jss(xdate)
m=2;
n=length(xdate);
r=0.2*std(xdate);
cr=[];
gn=1;
gnmax=m;
while gn=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
k=1;
for i=1:n-m+1
for j=1:m
x2m(i,j)=xdate(i+j-1);
end
end
x2m;
for i=1:n-m+1
for j=1:n-m+1
if i~=j
d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
k=k+1;
end
end
k=1;
end
d;
for i=1:n-m+1
[k,l]=size(find(d(i,r));%将比R小的个数传送给L
cr1(1,i)=l;
end
cr1;
cr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1
if cr1(i)~=0
sum1=sum1+log(cr1(i));
end
end
cr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
gn=gn+1;
m=m+1;
end
cr;
shang=cr(1,1)-cr(1,2);
function [shang]=ybs(xdate)
m=2;
n=length(xdate);
r=0.2*std(xdate);
cr=[];
gn=1;
gnmax=m;
while gn=gnmax
x2m=zeros(n-m+1,m);%存放变换后的向量
d=zeros(n-m+1,n-m);% 存放距离结果的矩阵
cr1=zeros(1,n-m+1);%存放
k=1;
for i=1:n-m+1
for j=1:m
x2m(i,j)=xdate(i+j-1);
end
end
x2m;
for i=1:n-m+1
for j=1:n-m+1
if i~=j
d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离
k=k+1;
end
end
k=1;
end
d;
for i=1:n-m+1
[k,l]=size(find(d(i,r));%将比R小的个数传送给L
cr1(1,i)=l;
end
cr1;
cr1=(1/(n-m))*cr1;
sum1=0;
for i=1:n-m+1
sum1=sum1+cr1(i);
end
end
cr1=1/(n-m+1)*sum1;
cr(1,gn)=cr1;
gn=gn+1;
m=m+1;
end
cr;
shang=-log(cr(1,1)/cr(1,2));
春花含笑意爆竹增欢声喜气盈门

运用matlab进行随机信号的功率谱密度估计仿真【急求】 Matlab画图问题,自相关函数,功率谱密度,求具体代码. 求问matlab一段代码的意思! 求基于matlab的EMD代码, 求矩阵LU分解的matlab代码求MATLAB 拟合指数函数的代码~ 求函数matlab代码求矩阵秩的matlab代码求矩阵秩的matlab代码设总体X的概率密度为f(x)={(a+1)x^a,0}其中a>—1是未知参数……求a 的矩估计量和最大似然估计量（见下图求matlab代码解方程! 求大神解释一下matlab代码： matlab小段代码求解释