matlab ellip,求大神帮我解释一下matlab最后几行是什么意思

该楼层疑似违规已被系统折叠隐藏此楼查看此楼

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); % 数据通过声卡转换为声音

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

% pause;

data1=x(:,1);

n=0:length(x)-1;

time=n/fs;

subplot(3,1,1) % 绘制2行1列的第1张子图

plot(time,data1) % 以时间为横轴，数据为纵轴作图

xlabel('Time (sec.)') % 标注横坐标

ylabel('Signal Level (Volts)') % 标注纵坐标

grid on % 添加网格

% pause;

% 对采集数据作滤波处理

blocksize =length(data1); % 计算窗函数长度

window = hanning(blocksize); % 计算汉宁窗函数(此函数为MATLAB自带)

data2=window.*data1; % 对数据先作加窗处理

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

wp=wp*2/fs;

wp=3000*2/fs;

[b,a]=ellip(4,0.1,20,wp); % 构造椭圆滤波器

data=filter(b,a,data2); % 求加窗处理后的数据经过滤波器的响应

subplot(3,1,2) % 绘制2行1列的第1张子图

plot(time,data) % 以时间为横轴，数据为纵轴作图

xlabel('Time (sec.)') % 标注横坐标

ylabel('Signal Level (Volts)') % 标注纵坐标

grid on % 添加网格

% 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'])

matlab ellip,求大神帮我解释一下matlab最后几行是什么意思相关推荐

matlab 最后一列,求大神帮我解释一下matlab最后几行是什么意思
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 clear; if nargin<1;action='initialized';end; [fname,pname]=uigetfile('*.wa ...
matlab fname pname,求大神帮我解释一下matlab最后几行是什么意思
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 clear; if nargin<1;action='initialized';end; [fname,pname]=uigetfile('*.wa ...
matlab中bits是什么意思,求大神帮我解释一下matlab最后几行是什么意思
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 clear; if nargin<1;action='initialized';end; [fname,pname]=uigetfile('*.wa ...
赫夫曼编码c语言排序部分,c语言构造哈夫曼树输出哈夫曼编码出错，跪求大神帮我找错...
c语言构造哈夫曼树输出哈夫曼编码出错 ,跪求大神帮我找错0 youxun0952016.09.07浏览120次分享举报 #include #include #include typedef struc ...
php输出后以密文显示,加密解密,_求大神帮看一下，如何根据下面的程序，把明文生成密文谢谢，加密解密 - phpStudy...
求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢 def parserLicense(self,encCode): data = base ...
c语言int32u的作用,求c语言大神帮我解释一下这段说的都是啥？
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 //滚动显示 } else if(Mode == 0xdd){ OffsetTemp = j + Offset; if(OffsetTemp<16) ...
Java修改全文数字,求大神帮我修改下java计算机代码，数字键只有0和8能用
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 package rj0211; import java.awt.BorderLayout; import java.awt.GridLayout; imp ...
java enum 报错_enum报错问题，求大神帮看下
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 public enum OrderStatusEnum { NOT_PAY("未付款", 0), NOT_YET_SHIPPEND(& ...
figure服务器无法显示,求大神帮帮忙，看一下为什么第二个figure出不来，只能显示第一个...
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 I=imread('C:\Users\Administrator\Desktop\123.jpg'); figure(1) subplot(3,2,1), ...

matlab ellip,求大神帮我解释一下matlab最后几行是什么意思

matlab ellip,求大神帮我解释一下matlab最后几行是什么意思相关推荐

最新文章

热门文章