clc;close?all;clear?all;
%sig_test=csvread（‘H:\Qjt.dat‘）;?%讀入數據生成（N1）矩陣

data=fopen（‘H:\Qjt.dat‘‘r‘）;
sig_test=fread（data）;
fclose（data）;

nfft=2048;???????????????????????%設置分段fft長度
overlap=0;??????????????????????%設置疊合長度
wn=bartlett（nfft）;????????????????%設置選用的窗函數
fs=50;???????????????????????????%設置抽樣頻率
[PxF]=pwelch（sig_testwnoverlapnfftfs）;???%進行Welch譜估計

log_Px=10*log10（Px）;
plot（Flog_Px）;
title（‘Welch法估計‘）;grid?on;ylabel（‘Px‘）;xlabel（‘F‘）;

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

?????文件????????514??2014-11-11?03:11??welch\pwelch_test.m

?????文件???????1335??2014-11-11?04:06??welch\P_welch.m

?????目錄??????????0??2014-11-14?22:51??welch

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

?????????????????1849????????????????????3