%隨機減量法

clear
close?all
format?long

fprintf（‘Samples?In?Test?:\n‘）;
path_base?=?sprintf（‘test‘）;
dir_sample?=?dir（path_base）;
[nsampledummy]?=?size（dir_sample）;
nsample?=?nsample?-?2;
for?i?=1:nsample
????fprintf（strcat（num2str（i）‘:?????‘））;
????fprintf（dir_sample（i+2）.name）;
????fprintf（‘\n‘）;
end
SampleNum?=?input（‘Please?select?the?sample?number:\n‘）;

path_sample?=?strcat（path_base‘\‘dir_sample（SampleNum+2）.name）;?%strcat：連接字符串的函數
name_sample?=?dir_sample（SampleNum+2）.name;

%提取測試數據
data（:1）=load（path_sample）;?
%所選的時域樣本內數據個數
ntdata?=?length（data（:1））;?
Cvtdata=?zeros（ntdata1）;


%【1】正常信號處理:
Cvtdata?=?data?;
%【1】正常信號處理end;


%?%【2】Narada信號轉換:
%?for?j=1:length（data（:1））
%?????data（j1）?=?data（j1）?*?5.0?/?65535.0;??%Convert?the?data?轉化為電壓信號
%?end???????
%?
%?%增益系數
%?gain?=?10;
%?for?j=1:ntdata??????
%?????Cvtdata（j1）=?data（j1）/（gain*200/1000）;
%?end
%?Avgdata?=?zeros（11）;
%?Avgdata?=?mean（?Cvtdata（:1）??）;?%求數據平均值
%?%【2】Narada信號轉換end;

sf=800;
np?=?100;
x?=?Cvtdata‘;

t=0:1/sf:（np-1）/sf;????%輸出時間向量
nt=length（x）;????%提取輸入數據的長度
s=1.5*std（x）;????%設置截取振幅A的大小
m=0;
y=zeros（1np）;
for?k=2:nt-np
????a=abs（x（k-1）-s）;
????b=abs（x（k）-s）;
????c=abs（x（k+1）-s）;
????if?b????????y（1:np）=y（1:np）+x（k:k+np-1）;
????????m=m+1;
????end
end

y=y./m;???%對疊加結果進行平均

%繪圖
plot（ty）;
xlabel（‘時間（s）‘）;
ylabel（‘幅值‘）;
grid?on;
%?fid=fopen（fno‘w‘）;
%?for?k=1:np
%????fprintf（fid‘%f\n‘y（k））;
%?end
%?status=fclose（fid）;

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

?????文件???????1928??2015-09-16?13:05??RD法\outRDT.mat

?????文件???????1657??2015-10-21?13:47??RD法\RDT.m

?????文件??????25600??2015-10-21?12:03??RD法\test\result1_2.txt

?????文件??????25600??2015-10-21?12:03??RD法\test\result2_2.txt

?????文件?????350023??2015-09-16?13:11??RD法\x值.xlsx

?????文件??????80192??2007-11-18?21:58??RD法\y.mat

?????目錄??????????0??2015-10-21?13:44??RD法\test

?????目錄??????????0??2015-10-21?13:44??RD法

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

???????????????485000????????????????????8