您现在的位置:首页 >> 技术文章 >> MATLAB技术 >> 内容

MATLAB代做-python代做-FPGA代做-频率估计

时间:2019-6-26 1:17:56 点击:

  核心提示:MATLAB代做-python代做-FPGA代做-频率估计 ...

%频率估计--esprit

%----------------------------------------------

%初始条件

%----------------------------------------------

clear; 

clc;


for L=1;


n=1:128;               

xn=2*cos(2*pi*0.0005*n);

p=6;                            %阶数;

m=1.5*p;                          

%-------------------------------------------------------------------------

%1、由xn求r(0)—r(m)

%-------------------------------------------------------------------------

N=length(xn);

rxx1=xcorr(xn,'biased');          %求自相关函数;

rxx=rxx1(N:m+N);              %取r(0)—r(m);

%-------------------------------------------------------------------------

%2、构造自相关矩阵Rxx和Rxy

%-------------------------------------------------------------------------

Rxxx=toeplitz(rxx);              %m+1阶自相关矩阵

Rxx=Rxxx(1:m,1:m);            %m阶  Rxx

Rxy=Rxxx(1:m,2:m+1);          %m阶  Rxy

%-------------------------------------------------------------------------

%3、对Rxx特征分解,求出最小特征值t

%-------------------------------------------------------------------------

ev=eig(Rxx);                    %矩阵分解;

[s i]=min(ev);

[V D]=eig(Rxx);

t=D(i,i);

%-------------------------------------------------------------------------

% 4、计算Rx=(Cxx,Cxy)=(Rxx-t*I,Rxy-t*Z)

%-------------------------------------------------------------------------

%4.1求Z

II=eye(m+1);

I=eye(m);

Z=II(1:m,2:m+1);

%4.2求(Cxx,Cxy)

Cxx=Rxx-t*I;

Cxy=Rxy-t*Z;

syms x;

Rx=Cxx-x*Cxy;

%-------------------------------------------------------------------------

%5、求广义特征根w 与频率分量f

%-------------------------------------------------------------------------

y=det(Rx);

yy=sym2poly(y);             %提取系数

w=roots(yy); 


for i=1:m;

f(L,i)=angle(w(i))/(2*pi);               

end;


end



联系:highspeedlogic

QQ :1224848052

微信:HuangL1121

邮箱:1224848052@qq.com

网站:http://www.mat7lab.com/

网站:http://www.hslogic.com/

微信扫一扫:

作者:频率估计 来源:频率估计
  • 您是如何找到本站的?
  • 百度搜索
  • Google搜索
  • 查阅资料过程中
  • 论坛发现
  • 百度贴吧发现
  • 朋友介绍
本站最新成功开发工程项目案例
相关文章
  • 没有相关文章
相关评论
发表我的评论
  • 大名:
  • 内容:
  • matlab代做|matlab专业代做|matlab淘宝代做(www.hslogic.com) © 2019 版权所有 All Rights Reserved.
  • Email:highspeed_logic@163.com 站长QQ: 1224848052

    专业代做/代写/承接、MATLAB、SIMULINK、FPGA项目、博士/硕士/本科毕业设计、课题设计、论文,毕业论文,Coursework、Eassy、Assignment