why does the following code give error?
1 次查看(过去 30 天)
显示 更早的评论
u=[3 4 30 50];% Desired Vector
b=u;
[R,C]=size(b);
P=C/2;
M=2*C;
% calculate observed vector
xo=zeros(1,M);
for k=1:M
for i=1:P
xo(1,k)=xo(1,k)+1*exp(1i*((-pi/2)*sin(u(P+i))*(k-1)+(pi/(16*u(i)))*cos^2(u(P+i))*(k-1)^2)));
end
end
xe=zeros(1,M);
for k=1:M
for i=1:P
xe(1,k)=xe(1,k)+1*exp(1i*((k-1)*(-pi/2)*sin(b(P+i))+(pi/(16*b(i)))*cos^2(b(P+i))*(k-1)^2));
end
end
abc=0.0;
for m1=1:M
abc=abc+(abs(xo(1,m1)-xe(1,m1))).^2;
end
abc=abc/M;
e=abc
0 个评论
采纳的回答
Torsten
2023-12-13
Use
xo(1,k)=xo(1,k)+1*exp(1i*((-pi/2)*sin(u(P+i))*(k-1)+(pi/(16*u(i)))*cos(u(P+i))^2*(k-1)^2));
xe(1,k)=xe(1,k)+1*exp(1i*((k-1)*(-pi/2)*sin(b(P+i))+(pi/(16*b(i)))*cos(b(P+i))^2*(k-1)^2));
instead of
xo(1,k)=xo(1,k)+1*exp(1i*((-pi/2)*sin(u(P+i))*(k-1)+(pi/(16*u(i)))*cos^2(u(P+i))*(k-1)^2)));
xe(1,k)=xe(1,k)+1*exp(1i*((k-1)*(-pi/2)*sin(b(P+i))+(pi/(16*b(i)))*cos^2(b(P+i))*(k-1)^2));
5 个评论
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Phase-Locked Loops 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!