- 1、本文档共11页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
MATLAB的语音信号频谱分析.
实验一
[y,fs,bit]=wavread(I do)%读取音乐片段,fs是采样率
size(y)%求矩阵的行数和列数
y1=y(:,1);%对信号进行分列处理
n1=length(y1);%取y的长度
t1=(0:n1-1)/fs;%设置波形图横坐标
figure
subplot(2,1,1);
plot(t1,y1); %画出时域波形图
ylabel(幅值);
xlabel(时间(s));
title(信号波形);
subplot(2,1,2);
Y1=fft(y1);
w1=2/n1*(0:n1-1);%设置角频率
plot(w1,abs(Y1));%画频谱图
title(信号频谱);
xlabel(数字角频率);
ylabel(幅度);
grid on;
sound(y,fs);
实验二
[y,fs,bit]=wavread(I do)
y1=y(:,1);
n1=length(y1);
D=2;%设置抽样间隔
y2=y1(1:D:n1);%减抽样
n2=length(y2);%减抽样后信号长度
t2=(0:n2-1)/fs;%设置横坐标
figure
subplot(2,2,1);
plot(t2,y2); %绘制减抽样信号波形图
ylabel(幅度);
xlabel(时间(s));
title(2:1减抽样信号波形);
Y2=fft(y2); %对y2进行n2点fft谱分析
w2=2/n2*[0:n2-1];
subplot(2,2,3);
plot(w2,abs(Y2));%绘制减抽样信号频谱图
title(2:1减抽样信号频谱);
xlabel(数字角频率w);
ylabel(幅度);
grid on;
sound(y2,fs/D);
D=5;%设置抽样间隔
y3=y1(1:D:n1);%减抽样
n3=length(y3);%减抽样后信号长度
t3=(0:n3-1)/fs;%设置横坐标
subplot(2,2,2);
plot(t3,y3); %绘制减抽样信号波形图
ylabel(幅度);
xlabel(时间(s));
title(5:1减抽样信号波形);
Y3=fft(y3); %对y2进行n2点fft谱分析
w3=2/n3*[0:n3-1];
subplot(2,2,4);
plot(w3,abs(Y3));%绘制减抽样信号频谱图
title(5:1减抽样信号频谱);
xlabel(数字角频率w);
ylabel(幅度);
grid on;
sound(y3,fs/D)
实验三
[y,fs,bit]=wavread(I do)
y1=y(:,1);
n1=length(y1);
n3=0:n1-1;
b1=cos(0.75*pi*n3);%设置调制信号
b2=cos(0.5*pi*n3);
c1=b1.*y1;%对原信号进行调制
c2=b2.*y1;
lc1=length(c1);
t=(0:lc1-1)/fs;
figure %用载波对信号进行调制,并对其做fft变换
subplot(2,2,1) %获取频谱,从图中可以观察到,调制后的
plot(t,c1); %信号频谱发生搬移
xlabel(时间(s));
ylabel(幅度);
title(调制后信号);
subplot(2,2,2) %获取频谱,从图中可以观察到,调制后的
plot(t,c2); %信号频谱发生搬移
xlabel(时间(s));
ylabel(幅度);
title(调制后信号);
w1=2/lc1*[0:lc1-1];%设置角频率W
C1=fft(c1);
C2=fft(c2);
subplot(2,2,3)
plot(w1,abs(C1));
xlabel(数字角频率w);
ylabel(幅度);
title(调制后信号的频谱(高频率调制));
grid on;
subplot(2,2,4)
plot(w1,abs(C2));
xlabel(数字角频率w);
ylabel(幅度);
title(调制后信号的频谱(低频率调制));
grid on;
sound(c1,fs);
实验四
clear all;cla;close all
[a,fs,bit]=wavread(I do);
y1=a(:,1);%去单列数据进行分析
f1=fft(y1);
n=length(f1);
tn=(0:n-1)/fs;
w=2/n*[0:n-1];
%sound(y1,fs);
figure(1)
subplot(2,3,1);plot(tn,y1);
grid on;
tit
您可能关注的文档
最近下载
- 2025年湖南科技职业学院单招英语模拟试题(附答案解析).pdf VIP
- 2025年高考数学重点题型归纳精讲精练5.4三角形四心和奔驰定理(新高考地区)(原卷版) .docx VIP
- 新生儿窒息复苏试卷附有答案.docx
- 三 气体压强.ppt VIP
- 图解---“健康中国2030”规划纲要-医学课件.pptx
- 历史:第6课《中古时代的欧洲》课件 (川教版九年级上).ppt
- DeepSeek-清华大学104页《DeepSeek:从入门到精通》.pptx
- FANUC 0i-MB操作说明书 完整版.pdf
- (八省联考)2025年新高考适应性演练 政治试卷合集(含答详解案).docx
- 赤泥坝施工方案.doc
文档评论(0)