求助，相干函数的matlab程序问题

yaoyonghan

小弟现在在做振动方面的毕业设计，现在遇到一个问题就是，我自己用matlab编的程序求相干函数时，怎么也求不对啊，麻烦哪位高手给我看看阿,谢谢了！
    H=invG*Kt;
    Gqq=Gq*[1,exp(-j*omig*L/v);exp(-j*omig*L/v),1];
    OS=conj(H)*Gqq*(H)';
    OS55(i1)=conj([H(5,1),H(5,2)])*Gqq*[H(5,1),H(5,2)]';
    RS=Gqq*(H)';
    RS15(i1)=Gqq(1,:)*[H(5,1),H(5,2)]';
    RS25(i1)=Gqq(2,:)*[H(5,1),H(5,2)]';
    H51(i1)=H(5,1);
    H52(i1)=H(5,2);
    gama15(i1)=(abs(RS15(i1)))^2/Gqq(1,1)*OS55(i1);
    gama25(i1)=Gqq(2,2)*(abs(H52(i1)))^2/OS55(i1);
%plot(x,gama15);
plot(x,gama25);

其中：Gqq是系统输入的自谱，Gq是已知的，H是频响函数，Kt是已知的，OS是输出的自谱，OS55是第五个自由度的自谱，RS是互谱，RS15是1、5自由度的互谱，RS25类似，gama15是1、5自由度的相干函数，
gama25类似。

yaoyonghan

知道的请给我回复一下阿，本人急要，谢谢了！