|
楼主 |
发表于 2014-9-16 17:21
|
显示全部楼层
本帖最后由 whdwpy666 于 2014-9-16 17:22 编辑
rxy = zeros(ch,ch,n); thet = zeros(ch,ch);
for cch = 1:ch-1 % cch: colum of channel
for rch = cch+1:ch % rch: row of channel % 设置相干系数
r_20_2k = [0.1 0.8];
f_20_2k = log10([20 2000]);
r_xy = interp1(f_20_2k,r_20_2k,lgpf);
rxy(rch,cch,ran) = reshape(r_xy*cch*rch/(ch-1)/ch,1,1,length(pf)); % 设置相位角
thet(cch,rch) = pi/2;%第2点超前第1点相位角
% 设置互谱的值 R1(rch,cch,ran) = rxy(rch,cch,ran) .*sqrt(R1(rch,rch,ran)).*... sqrt(R1(cch,cch,ran))*exp(-1j*thet(cch,rch)); R1(cch,rch,ran) = conj(R1(rch,cch,ran));
end
end
R2=R1*9.8^2*df/2;%转为国际单位和双边谱
|
|