clear all
close all
N=100;
fc=1.0e9;
fs=20*fc;
w=2*pi*fc/fs;
n=0:N-1;
s=sqrt(1.0)*exp(i*w*n);
wn=sqrt(0.1/2)*(randn(1,N)+i*randn(1,N));
x=s+wn;
w_est=0;
M=N;
for m=1:M-1
w_est=w_est+(1/(M-1))*arg((x(m))'*x(m+1));
end
error_1=w-w_est;
w_est=0;
M=N;
for m=1:M-1
p(m)=6*(m)*(M-m)/(M*(M^2-1));
w_est=w_est+p(m)*arg((x(m))'*x(m+1));
end
error_2=w-w_est;

 채택된 답변

DGM
DGM 2021년 5월 4일

0 개 추천

Use angle() instead of arg() to get the argument of complex numbers.

댓글 수: 1

Sadiq Akbar
Sadiq Akbar 2021년 5월 4일
Thank you very much dear DGM for your prompt respose.

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

질문:

2021년 5월 4일

댓글:

2021년 5월 4일

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by