%%
%清除命令窗口，工作區
clc
clear?all?
%%
%導入原始振動數據和方位角數據
filename1=‘A.txt‘%振動數據
data1=textread（filename1‘%n‘）;
filename2=‘B.txt‘%方位角數據
data2=textread（filename2‘%n‘）;
fs=1024;%采樣頻率
f=6.77;%時標所在軸頻率
f_t=6.77;%感興趣時標所在軸頻率
n=length（data2）;
dic=zeros（n1）;
%%提取原時標對應振動信號的點數，存在dic矩陣中
for?i=1:n
????integer=fix（（data2（i1）*fs）/（13107200））;
????dicimal=（data2（i1）*fs）/（13107200）-integer;
????if?integer==（data2（i1）*fs）/（13107200）
???????dic（i1）=fix（（data2（i1）*fs）/（13107200））;
????else
????????dic（i1）=fix（（data2（i1）*fs）/（13107200））+1;
????end
end
%%等效時標轉換
k=length（dic）;%原時標個數
dicn=fix（（f_t/f）*（k-1））+1;%等效時標個數
dic_n=zeros（dicn1）;
for?j=1:dicn
????integer=fix（1+（j-1）*f/f_t）;%整數部分
????decimal=1+（j-1）*f/f_t-integer;%小數部分
????if?integer~=k
???????dic_n（j1）=dic（integer1）+round（decimal*（dic（integer+1）-dic（integer1）））;%round四舍五入，小數部分?
????else
????dic_n（j1）=dic（integer1）;
????end
end
dic_n=dic_n（2000:end1）;%我的數據里選擇從第dic_n（2000，1）個時標開始
????
??