%讀入信號
%對實測管道數據計算
%讀取文件數據
FS=500*1024;%實際信號的采樣頻率
TS=1/FS;

filename=‘泄漏定位4-3\73091217_07_06.dat‘;
%%%filename=‘泄漏定位8-3\73091203_18_03.dat‘;
filename=[‘D:\20100303PipeData\‘?filename];
[dataSampleF]?=ReadData（filename）;%讀取數據
temp=data/SampleF.Gain;%消除增益
N=3264;%N是每次處理的數據段長度
lendata=size（temp1）;%取出總數據長度
t=0:TS:（lendata-1）*TS;

figure（1）;
plot（t（1:N）temp（1:N））;
ylim（[-0.050.05]）;
xlabel（‘時間/s‘）;
ylabel（‘幅度/V‘）;
title（‘管道泄漏原始信號‘）;


dim=8;
tau=18;
startPosToZero=5;
%對時間序列加窗或者分段實行降噪處理，總數據長度261120
N=3264;%N是每次處理的數據段長度
xs=zeros（2611201）;%存放降完噪的純凈信號
for?j=1:5?%迭代奇異值分解的迭代次數
for?iD=1:1%---------------------由于有26112個數據，id的范圍應該選擇為1:80---------
%對時間序列加窗或者分段實行降噪處理，總數據長度261120
?x=temp（1+N*（iD-1）:N*iD）;%每段取N個數據存放在列向量中
????%newx?=?getDenoiseBySVD（xdimstartPosToZero）;
????%SVD降噪
????????????%?參數?dim表示相空間重構的維數
????????????%?參數?startPosToZero表示從第幾個奇異值開始之后就全為零
????????????%X?=?PhaSpaRecon（xdim10）;%相空間重構
????????????[XD]=reconstructDistanceM（x‘dimtau）;%相空間重構
????????????[USV]?=?svd（X）;
????????????sv?=?diag（S）;???%獲取奇異值向量?sigular?vector
%?%?????????????%?改變該奇異值向量到?newSv;
????????????newSv?=?sv;
????????????newSv（startPosToZero:end）?=?0;
%?%?????????????%?復原
????????????[mn]?=?size（S）;
????????????newS?=?[diag（newSv）zeros（mn-m）];
????????????newX?=?U*newS*V‘;
????????????%newx?=?PhaSpaRecov（newXdim1）;
????????????[newxnum]=RecoverToTimeseries（dimtaunewX）;
????????????clear?U?S?V
%?%?????
???xs（1+N*（iD-1）:N*iD）=newx‘;
end
???tempnew（1+N*（iD-1）:N*iD）=xs（1+N*（iD-1）:N*iD）;
end;
DN=tempnew-xs;
save（‘D:\！！MY20100308畢業論文仿真\8-3SVDDN.txt‘‘DN‘‘-ASCII‘）;
%?%?
figure（3）;
plot（t（1:N）xs（1:N））;
ylim（[-0.050.05]）;
xlabel（‘時間/s‘）;
ylabel（‘幅度/V‘）;
title（‘經過svd降噪后的管道泄漏信號‘）;

%?%?D=128;%控制fft變換點數
%?%?ymax=fft（xsD*1024）;?
%?%?py=ymax.*conj（ymax）/（D*1024）;%功率譜的計算
%?%?ff=FS*（0:D*1024/2）/（D*1024）;%算頻率坐標
%?%?plot（ffpy（1:（D*1024/2+1）））;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????2942112??2011-03-23?09:09??奇異值分解降噪\B007_0.mat

?????文件????2910768??2011-03-23?09:09??奇異值分解降噪\IR007_0.mat

?????文件???????2203??2011-03-23?09:07??奇異值分解降噪\MainSVDDenoise.m

?????文件????3903344??2011-03-23?09:09??奇異值分解降噪\Normal0.mat

?????文件????2928192??2011-03-23?09:09??奇異值分解降噪\OR007@6_0.mat

?????文件???????1009??2011-03-23?09:07??奇異值分解降噪\qiyizhi.m

?????文件???????2842??2011-03-23?09:07??奇異值分解降噪\ReadData.m

?????文件????1318969??2011-03-23?09:07??奇異值分解降噪\reconstructDistanceM.dll

?????文件???????1642??2011-03-23?09:07??奇異值分解降噪\RecoverToTimeseries.m

?????文件???????1228??2011-03-23?09:12??奇異值分解降噪\Untitled.m

?????文件??????26624??2011-03-23?09:07??奇異值分解降噪\說明文檔.doc

?????目錄??????????0??2013-08-23?21:38??奇異值分解降噪

-----------?---------??----------?-----??----

?????????????14039151????????????????????13