simomo 发表于 2006-5-21 10:21

[求助]经过FFT后得到的响应和激励谱线求自功率谱和互功率谱的问题

就是<BR>对加窗后的响应信号进行FFT，得响应谱线 Ｘi;<BR>对加窗后的力信号进行FFT，得力谱线Fi;<BR><BR>得到ＧＦＦi=Fi*Fi'(Fi的共扼)；ＧＦＸi=Xi*Fi'<BR><BR>下一步要将通过Ｍ次测得的的ＧＦＦ和ＧＦＸ求和<BR>就是下面的式子怎么用ＭＡＴＬＡＢ实现？<BR>我用<BR>clc<BR>N=1000;<BR>n=0:N-1;<BR>Xk=fft(xn);　　　　　％xn是采集到的响应信号,fn是采集到的激励信号<BR>Xf=fft(fn);<BR>Gff=Xf*Xf'; %求自功率谱<BR>Gfx=Xk*Xf'; %求互功率谱<BR>m=100; %m是采集点的个数<BR>问题就在下面,要怎么才能实现我上面说的那些呢?请各位指点一下,谢谢了(我是本科生,做柴油机故障诊断的毕业设计其中要用到ANYSY,PRO/E和MATLAB,导师给了个研究生的课题,要学的东西实在是太多了,自己精力有限,只要来这里求助了)<BR>sym Gfx m;<BR>sym Gfx;<BR>SFX=simple(symsum(Gfx,1,m)) <BR><a href="http://bbs.matwav.com/upload/2006/05/18/13613970.jpg" target="_blank" ><IMG src="http://bbs.matwav.com/upload/2006/05/18/13613970.jpg" align=absMiddle border=0></A> <BR>

yangzj 发表于 2006-5-21 10:26

<P>经过累加就可得到呀</P>

simomo 发表于 2006-5-21 10:28

<P>上面的问题可能就出在这里,就是如果我采集了很多个信号,用x1,x2,x3.........xm来表示<br>怎样才能实现在程序中用xn (n=1,2,3....m)来表示x1,x2,x3.........xm<br><br></P>
[此贴子已经被作者于2006-5-21 10:31:06编辑过]

simomo 发表于 2006-5-21 10:33

回复：（yangzj）经过累加就可得到呀

<DIV class=quote><B>以下是引用<I>yangzj</I>在2006-5-21 10:26:44的发言：</B><BR>
<P>经过累加就可得到呀</P></DIV>
<br>请问累加使用什么语句?我哪个对吗?<BR>要是不对能不能帮我指点一下那里错了,怎么改,谢谢

yangzj 发表于 2006-5-21 10:34

用一个l*m的矩阵表示信号即可.一列代表一个信号
页: [1]
查看完整版本: [求助]经过FFT后得到的响应和激励谱线求自功率谱和互功率谱的问题