該樓層疑似違規(guī)已被系統(tǒng)折疊?隱藏此樓查看此樓

clear;

if nargin<1;action='initialized';end;

[fname,pname]=uigetfile('*.wav','Open Wave File');

file=[pname,fname];

[x,fs,bits]=wavread(file); % 讀入聲音文件(*.wav)

sound(x,fs,bits); % 數(shù)據(jù)通過(guò)聲卡轉(zhuǎn)換為聲音

%===========================================================

% pause;

data1=x(:,1);

n=0:length(x)-1;

time=n/fs;

subplot(3,1,1) % 繪制2行1列的第1張子圖

plot(time,data1) % 以時(shí)間為橫軸，數(shù)據(jù)為縱軸作圖

xlabel('Time (sec.)') % 標(biāo)注橫坐標(biāo)

ylabel('Signal Level (Volts)') % 標(biāo)注縱坐標(biāo)

grid on % 添加網(wǎng)格

% pause;

% 對(duì)采集數(shù)據(jù)作濾波處理

blocksize =length(data1); % 計(jì)算窗函數(shù)長(zhǎng)度

window = hanning(blocksize); % 計(jì)算漢寧窗函數(shù)(此函數(shù)為MATLAB自帶)

data2=window.*data1; % 對(duì)數(shù)據(jù)先作加窗處理

wp=[70,400]; %100Hz--400Hz

wp=wp*2/fs;

wp=3000*2/fs;

[b,a]=ellip(4,0.1,20,wp); % 構(gòu)造橢圓濾波器

data=filter(b,a,data2); % 求加窗處理后的數(shù)據(jù)經(jīng)過(guò)濾波器的響應(yīng)

subplot(3,1,2) % 繪制2行1列的第1張子圖

plot(time,data) % 以時(shí)間為橫軸，數(shù)據(jù)為縱軸作圖

xlabel('Time (sec.)') % 標(biāo)注橫坐標(biāo)

ylabel('Signal Level (Volts)') % 標(biāo)注縱坐標(biāo)

grid on % 添加網(wǎng)格

% pause;

[xmax,index]=max(data1);

timewin=floor(0.015*fs);

xwin=data1(index-timewin:index+timewin);

[y,lags]=xcov(xwin);

subplot(3,1,3)

plot(lags,y)

grid on

ylen=length(y);

halflen=(ylen+1)/2 +30;

yy=y(halflen: ylen);

[ymax,maxindex] = max(yy);

fmax=fs/(maxindex+30);

disp(['Maximum occurred at ', num2str(fmax), ' Hz'])

總結(jié)

以上是生活随笔為你收集整理的matlab 最后一列,求大神帮我解释一下matlab最后几行是什么意思的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。