那位大哥大姐能帮我检查密集谱校正的程序

S0704082

程序是采用“利用FFT校正两个密集信号的频率和相位  作者：方体莲 洪一”的方法做得，关键问题出现在方程求解的传递过程中，请各位赐教，谢谢！！

songzy41

原帖由 S0704082 于 2009-2-24 15:47 发表
程序是采用“利用FFT校正两个密集信号的频率和相位  作者：方体莲 洪一”的方法做得，关键问题出现在方程求解的传递过程中，请各位赐教，谢谢！！

1，程序中
a2=Q11*Q2*yn-Q1*I22*xn-I11*Q2+I1*I22*yn+I1*Q22*xn+I2*I11;
是否应改为
a2=Q11*Q2*yn-Q1*I22*xn-I11*Q2+I1*I22*yn+I1*Q22*xn+I2*Q11;
2，不明白“在方程求解的传递过程中”指的是什么。
solve('x1^2+y1^2=1','a0+a1*x1+a2*y1+a4=0')
给出的解便是x1,y1,x2,y2的字符表达式，把a0,a1,a2,a4代入便给出了x1,y1,x2,y2的其体数值，这样就能进一步求出dk1和dk2。
3，在已知x1,y1,x2,y2后，可用反正切直接求出dk1和dk2。