|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
小弟最近在研究SSI模态参数识别的程序,由于本人在Matlab方面是菜鸟,故有问题请教。
在求互功率谱密度是有一段程序:
fni=input('NEXT:','s');
fid=fopen(fni,'r');
sf=fscanf(fid,'%f',1);
np=fscanf(fid,'%d',1);
fno=fscanf(fid,'%s',1);
%第一行为参照点,第二行为响应数据点
b=fscanf(fid,'%f',[2,inf]);
status=fclose(fid);
x=b(1,:);
y=b(2,:);
t=0:1/sf:(np-1)/sf;
nfft=2^nextpow2(2*np);
p=csd(x,y,nfft); %求互功率谱密度
p(nfft/2+1)=real(p(nfft/2));
p(nfft/2+2:nfft)=conj(p(nfft/2:-1:2));
g=ifft(p); %进行逆傅里叶变化,求互相关函数
r=real(g(1:np));
plot(t,r);
xlabel('时间(s)');
ylabel('幅值');
grid on
其中,
p(nfft/2+1)=real(p(nfft/2));
p(nfft/2+2:nfft)=conj(p(nfft/2:-1:2));
这两个语句没看懂,好像是求共轭什么的,请问一下各位,谢谢。 |
|