1信号间的相互关系

1 信号间的相互关系

1.2 新的信号间运算类

1. classdef CSigCalc<handle
2. % relation between two signals
3. %
4. properties
5. m_sig1;
6. m_sig2;
7. end
9. methods
10. function obj = CSigCalc(sig1, sig2)
11. dFs1 = sig1.m_dFs;
12. dFs2 = sig2.m_dFs;
13. % 判断两信号的采样率一致
14. if abs(dFs1 - dFs2) < eps
15. obj.m_sig1 = sig1;
16. obj.m_sig2 = sig2;
17. else
18. disp('Sample frequency of the two signals are not equal.')
19. disp('Signal object is not constructed.')
20. return;
21. end
22. end
24. function relat = Relation(obj)
25. % 计算相关系数
26. vecSig1 = obj.m_sig1.m_vecSig;
27. [~, colNum] = size(vecSig1);
28. if colNum == 1 % to be sure the row vector
29. vecSig1 = vecSig1';
30. end
31. vecSig2 = obj.m_sig2.m_vecSig;
32. [~, colNum] = size(vecSig2);
33. if colNum == 1 % to be sure the row vector
34. vecSig2 = vecSig2';
35. end
37. if length(vecSig1) > length(vecSig2)
38. vecSig1 = vecSig1(1:length(vecSig2));
39. elseif length(vecSig1) < length(vecSig2)
40. vecSig2 = vecSig2(1:length(vecSig1));
41. else
43. end
45. relat = (vecSig1*vecSig2')/...
46. sqrt(vecSig1*vecSig1')/sqrt(vecSig2*vecSig2');
47. end
49. function Disp2SigCurve(obj)
50. Fs = obj.m_sig1.m_dFs;
51. vecSig1 = obj.m_sig1.m_vecSig;
52. vecT1 = (0:length(vecSig1)-1)/Fs;
53. vecSig2 = obj.m_sig2.m_vecSig;
54. vecT2 = (0:length(vecSig2)-1)/Fs;
55. figure;
56. plot(vecT1,vecSig1,vecT2,vecSig2,'linewidth',2)
57. xlabel('时间 S')
58. ylabel('幅值')
59. legend('1','2')
60. end
62. function [vecRelFunc, vecLags] = RelatFunc(obj)
63. strArrScaleOpt = {'biased','unbiased','coeff','none'};
64. ind = 3;
65. dFs = obj.m_sig1.m_dFs;
66. [vecRelFunc, vecLags] = xcorr(obj.m_sig1.m_vecSig, obj.m_sig2.m_vecSig,strArrScaleOpt{ind});
67. vecLags = vecLags/dFs;
68. figure;
69. plot(vecLags, vecRelFunc)
70. title(strArrScaleOpt{ind})
71. grid;
73. [~,indexMax] = max(vecRelFunc);
74. disp(['Time lag = ' num2str(vecLags(indexMax))])
75. end
76. end
77. end

1.2 脚本程序

1. % Exer_SigProcess_Script.m
2. close all;
4. strCurrentPath = cd;
5. addpath([strCurrentPath '\classdef']);
7. % 构造信号1，sig1
8. vecAmp=[1 3 5 7];
9. vecFreq = [0.25 2 3 4];
10. vecPhi = [0 pi/3 pi/4 pi/2];
12. Fs = 2^10;
13. t = 0:1/Fs:2*pi;
15. ind = 1;
16. xt = vecAmp(ind)*cos(2*pi*vecFreq(ind)*t + vecPhi(ind));
17. sig1 = CSignal(Fs, xt);
19. xt = vecAmp(ind)*cos(2*pi*vecFreq(ind)*t + vecPhi(4));
20. sig2 = CSignal(Fs, xt);
22. sigCalc = CSigCalc(sig1,sig2);
23. relat = sigCalc.Relation;
24. disp(['relat = ' num2str(relat)])
25. Disp2SigCurve(sigCalc);
28. rmpath([strCurrentPath '\classdef']);