残夜痕 残夜痕
关注数: 2 粉丝数: 39 发帖数: 2,055 关注贴吧数: 12
减谱法程序求帮忙改下,或者给我一个简单减谱法的程序 %基本减谱法 clear; %[xx,fs]=wavread('E:\mywhisper\shu.wav'); %[xx,fs]=wavread('E:\speech\x\w1xun_01.wav'); %[xx,fs]=wavread('E:\speech\耳语音切割\b\w1ba_5'); [xx,fs]=wavread('D:/工具/11.wav');% 读取音频文件yuan.wav,并返回采样数据给变量xx及采样率Fs [team,row]=size(xx);%将数组xx的行数赋给team,列数赋给row if row==2 x=(xx(:,1)+xx(:,2))/2; yy=x; %如果语音信号xx为2列,即信号为双声道,则将其转换成单声道信号,即取两列的平均值赋给x,并将x的值赋给yy else x=xx; yy=x; %若语音信号xx为单声道,则将xx的值赋给x,并将x的值赋给yy end x=x-mean(x)+0.1*rand(length(x),1); N=length(x);%将语音信号长度赋给变量N n=220;%对语音信号进行分帧,帧长为220 n1=160;%帧移为160 frame=floor((N-n)/(n-n1));%将分帧数赋给变量frame %frame=floor(N/n); for i=1:frame y1=x((i-1)*(n-n1)+1:(i-1)*(n-n1)+n).*hamming(n); %对每段分帧进行加窗处理 fy=fft(y1,n); nen(i,:)=abs(fy).^2;% 将频域信号功率赋给矩阵变量nen ang(i,:)=angle(fy);%将频域信号的相位角赋给矩阵变量ang end yuzhi=sum(sum(nen(2:5,:)))/(4*n); for i=1:frame nen(i,:)=nen(i,:)-yuzhi; nen(i,find(nen(i,:)<0))=0; %chuli=nen(i,1:n/2); %chuli=chuli-yuzhi; %chuli(find(chuli<0))=0; %nen(i,:)=[chuli,fliplr(chuli)]; % nen(i,:)=filter(1,[0.5 0.5],nen(i,:)); % nen(i,find(nen(i,:)<0))=0; end for i=1:frame nen(i,:)=sqrt(nen(i,:));%将纯语音功率谱开根,得到频域值 jie=nen(i,:).*exp(j*ang(i,:)); out(i,:)=real(ifft(jie))/hamming(n)'; %对纯语音频谱进行逆傅里叶变换,并取其实部,并进行去窗处理 end zong=out(1,:)';%将第一帧中未重叠部分记入数组zong jiewei=n; for i=2:frame zong(jiewei-n1+1:jiewei)=(zong(jiewei-n1+1:jiewei)+out(i,1:n1)')/2; jiewei=jiewei+n-n1;% 使指针jiewei依次指向下一帧的帧尾 zong=[zong;out(i,n1+1:end)'];%将从第二帧开始的每一帧中未重叠部分记入数组zong end %zong=out(1,:)'; %for i=2:frame %zong=[zong;out(i,:)']; %end % %for i=1:frame % zong=[zong,nen(i,:)']; % zong(i*(n-n1)+1:(i-1)*(n-n1)+n) =zong(i*(n-n1)+1:(i-1)*(n-n1)+n)/2; figure(1); %创建图1 subplot(211); %把图形窗口分成2*1个小窗口,取第1个小窗口 plot(x);%以数组x绘图基本二维曲线 axis([1,(n-n1)*frame+n,min(x),max(x)]);% 对当前二维图形对象的X轴和Y轴进行标定,x轴的范围为1到最后一个分帧结尾,y轴的范围为带噪语音时域最小值到最大值 subplot(212); %把图形窗口分成2*1个小窗口,取第2个小窗口 specgram(x,fs,1024,n,n1);%画出语音的语谱图 figure(2); %创建图2 subplot(211); %把图形窗口分成2*1个小窗口,取第1个小窗口 plot(zong); %以数组zong绘图基本二维曲线 axis([1,(n-n1)*frame+n,min(zong),max(zong)]); % 对当前二维图形对象的X轴和Y轴进行标定,x轴的范围为1到最后一个分帧结尾,y轴的范围为纯语音时域最小值到最大值 subplot(212); %把图形窗口分成2*1个小窗口,取第2个小窗口 specgram(zong,fs,1024,n,n1); %画出语音的语谱图 wavplay(x,fs);%播放单声道带噪语音音频 wavplay(zong,fs);%播放单声道纯净语音音频 报错Index exceeds matrix dimensions
凭心而说,真不希望看你们黑大中华 凭心而说,真不希望看你们黑大中华。你们可以黑共肠蛋,可以黑小三,可以黑李刚,可以黑等等等等等,但不应该黑大中华。每个国家都有自己很不好的一面,但我们却没怎么发现外国的坏?为什么,因为他们都不说,他们的国家,他们的淫都只知道自己国家有多好多好。但我们是亲自挖出不好的一次次展现在外国面前,然后告诉他们我们多不民主,我们多不好。 二战后,很多有野心的国家都预见到战争短期内是不可能再爆发了,于是采取的是各种同化政策。尤其对我大中华,文化入侵,科技入侵,还有即将就发生的市场和经济的入侵。他们成功了,DS们向往美利坚的民主,美利坚的发达,喜欢美利坚的产品,虽然说美国确实是有其好的地方,但是DS们对其的崇拜达到了盲目。国货始终不敢买,一定要到国外过个品牌才肯买。很多大品牌的的代加工都是中国,我们国货的质量不是那么差的。 还有,很多事情不在其位根本不知其苦。比如整天黑的城管,我自己亲眼目睹了一幕,你们敢黑吗?“、 城管来收了,因为摆在了街道两旁,使街道宽度少了3分1,影响了出行,摊主是个老人,因为没来得及跑,就在原地装傻充愣,城管说“老人家,你的摊子违反规定了,要没收”老人家继续不说话,估计是怕老人,城管又说,“摊子我们不没收了,您别再摆了好不好。。。。。”大概意思就说有人举报他也难办,然后老人家就推着车走了,城管一走,她马上又回来。后来城管估计是收队了,拉着一堆东西又路过那,看到老人家又在那。就上前去说不是不让摆了,老人又当没听到,另一个城管估计脾气比较差,就说跟他说那么多干神马,直接收了算了,就要去收摊,手干抓在篮上老人就开始拒绝,一直推那城管,城管估计也不敢还手,就任老人打,结果刚把篮子弄上车,老人就倒在地上,眼泪一直流,城管一看也荒了,就把老人拉上车送医院。。。然后我就不知道怎么解决了。 这种情况很多次了,我们每次一看那些惊心动魄的图就指责城管,但事实上呢,城管是工作,他如果不管估计我们街道都不能通车,他要管,而有些摊子为了谋生又必须摆,然后起了争执,最后动手,相信城管被打也不是少数。 这只是一种情况,还有很多职业,甚至可以说是没有一个职业没被黑过,但是当我们自己亲自去做的时候,又有多少的无可奈何? 大中华是有很多问题,但不是无可救药,而且不管再怎么样都是祖国不是吗?每次我看到大美利坚,天朝这种嘲讽的词我都有一种很绝望的感觉,甚至也想出国,可见言论毒害的深刻。 所以我觉得我们DS应该去发掘黑暗,但不该放弃光明,不要再建设一个绝望的中华了,你们敢说10后不是继承了一些我们腐朽的思想并且发扬的?DS应该调正心态,救中华,而不该恨中华,除非以后全部出国,改国籍大美利坚吧 来喷吧,能解释的我一一解释,不能解释就当我吃翔了吧
1 下一页