作者q33339955 (羊弟)
看板MATLAB
标题[绘图] 有关beamforming mic的
时间Thu Apr 4 14:24:46 2013
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