2022年MATLAB变声程序代码 .pdf
( 4.5 )《2022年MATLAB变声程序代码 .pdf》由会员分享，可在线阅读，更多相关《2022年MATLAB变声程序代码 .pdf（4页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、%ct1 clear all,close all,clc;%定义常数FL=80;%帧长WL=240;%窗长P=10;%预测系数个数s,fs=wavread(sunday_2.wav);%载入语音s s=s/max(s);%归一化L=length(s);%读入语音长度FN=floor(L/FL)-2;%计算帧数%预测和重建滤波器exc=zeros(L,1);%激励信号（预测误差）zi_pre=zeros(P,1);%预测滤波器的状态s_rec=zeros(L,1);%重建语音zi_rec=zeros(P,1);%合成滤波器exc_syn=zeros(L,1);%合成的激励信号（脉冲串）s_syn
2、=zeros(L,1);%合成语音last_syn=0;%存储上一个（或多个）段的最后一个脉冲的下标zi_syn=zeros(P,1);%合成滤波器的状态%变调不变速滤波器exc_syn_t=zeros(L,1);%合成的激励信号（脉冲串）s_syn_t=zeros(L,1);%合成语音last_syn_t=0;%存储上一个（或多个）段的最后一个脉冲的下标zi_syn_t=zeros(P,1);%合成滤波器的状态%变速不变调滤波器（假设速度减慢一倍）v=.5;exc_syn_v=zeros(vL,1);%合成的激励信号（脉冲串）s_syn_v=zeros(vL,1);%合成语音last_syn
3、_v=0;%存储上一个（或多个）段的最后一个脉冲的下标zi_syn_v=zeros(P,1);%合成滤波器的状态hw=hamming(WL);%汉明窗%依次处理每帧语音for n=3:FN%计算预测系数（不需要掌握）s_w=s(n*FL-WL+1:n*FL).*hw;%汉明窗加权后的语音 A E=lpc(s_w,P);%用线性预测法计算P 个预测系数%A 是预测系数，E 会被用来计算合成激励的能量名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页，共 4 页 -if n=27%(3)观察预测系统的零极点图 zplane(1,A);end s_f=s(n-1)*FL+1:n*FL);%本帧
4、语音，下面就要对它做处理%(4)用 filter 函数 s_f 计算激励，注意保持滤波器状态 exc1,zi_pre=filter(A,1,s_f,zi_pre);exc(n-1)*FL+1:n*FL)=exc1;%计算得到的激励%(5)用 filter 函数和 exc 重建语音，注意保持滤波器状态 s_rec1,zi_rec=filter(1,A,exc1,zi_rec);s_rec(n-1)*FL+1:n*FL)=s_rec1;%计算得到的重建语音%注意下面只有在得到exc 后才会计算正确 s_Pitch=exc(n*FL-222:n*FL);PT=findpitch(s_Pitch);%
5、计算基音周期PT（不要求掌握）G=sqrt(E*PT);%计算合成激励的能量G（不要求掌握）%方法 3：本段激励只能修改本段长度 tempn_syn=1:n*FL-last_syn;exc_syn1=zeros(length(tempn_syn),1);exc_syn1(mod(tempn_syn,PT)=0)=G;%某一段算出的脉冲 exc_syn1=exc_syn1(n-1)*FL-last_syn+1:n*FL-last_syn);s_syn1,zi_syn=filter(1,A,exc_syn1,zi_syn);exc_syn(n-1)*FL+1:n*FL)=exc_syn1;%计算得
6、到的合成激励 s_syn(n-1)*FL+1:n*FL)=s_syn1;%计算得到的合成语音 last_syn=last_syn+PT*floor(n*FL-last_syn)/PT);%(11)不改变基音周期和预测系数，将合成激励的长度增加一倍，再作为filter%的输入得到新的合成语音，听一听是不是速度变慢了，但音调没有变。FL_v=floor(FL/v);tempn_syn_v=1:n*FL_v-last_syn_v;exc_syn1_v=zeros(length(tempn_syn_v),1);exc_syn1_v(mod(tempn_syn_v,PT)=0)=G;%某一段算出的脉冲
7、exc_syn1_v=exc_syn1_v(n-1)*FL_v-last_syn_v+1:n*FL_v-last_syn_v);s_syn1_v,zi_syn_v=filter(1,A,exc_syn1_v,zi_syn_v);last_syn_v=last_syn_v+PT*floor(n*FL_v-last_syn_v)/PT);exc_syn_v(n-1)*FL_v+1:n*FL_v)=exc_syn1_v;%计算得到的加长合成激励 s_syn_v(n-1)*FL_v+1:n*FL_v)=s_syn1_v;%计算得到的加长合成语音%(13)将基音周期减小一半，将共振峰频率增加150Hz，
- 配套讲稿：
如PPT文件的首页显示word图标，表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制：
部分文档作品中含有的国旗、国徽等图片，仅作为作品整体效果示例展示，禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词：
- 2022年MATLAB变声程序代码 2022 MATLAB 变声 程序代码
淘文阁 - 分享文档赚钱的网站所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。
限制150内