[繪圖] 有關beamforming mic的
clear all
close all
clc
b=0;
M=8;
N=50;
n=0:N-1;
f = 1000;
c = 33000;
lumbda = c/f;
d = lumbda/2;
fs = 32 * f;
del_s=d*fs/c;
Om=2*pi*f/fs;
mv=0:M-1;
a=ones(M,1)/M;
S1=sin((-90:90)*pi/180);
for id = 1:181
qb=del_s*S1(id);
% X=sin(Om*(n(:,ones(1,M)) + qb*mv(ones(N,1),:) ));
X=sin(Om*(qb*mv(ones(N,1),:) ));
for m =1:M
qmb=eye(1,(m-1)*b+M);
Xqmb(:,m)=filter(fliplr(qmb),qmb,X(:,m));
end;
Y=Xqmb*a;
P(id)=sum(Y.^2);
end;
P_dB=10*log10(P);
plot(id,P_dB,'-o');
xlabel('DOA (Degree)');
ylabel('Power (dB)');
grid;
請問一下要如何把每一度的點都各別顯示呀,不要全部的圖都疊在181度的地方!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 112.104.151.58