这个程序怎么改啊？

chy19840208

　　本设计步骤如下：
　　1.录制一段自己的语音信号，并对录制的信号进行采样;画出采样后的语音信号的频谱图;
　　2. 插入噪音;画出含有噪音的频谱图;
　　3.给定滤波器的性能指标，采用窗函数法，划出滤波器的频域响应
　　4.用该滤波器对含噪音的信号进行滤波，画出滤波后信号的频谱;
　　5.回放语音信号。
　　程序如下：
　　fs=8000;
　　x1=wavread('E:\我的论文\论文\luyin.wav');
　　y1=fft(x1);
　　figure(1);
　　plot(abs(y1));
　　title('原始信号频谱');
　　xlabel('频率');
　　ylabel('频谱幅度');
　　t=0:1/8000:(length(x1)-1)/8000;
　　d=sin(2*pi*2000*t);
　　x2=x1+d;
　　sound(x2,fs);
　　y2=fft(x2);
　　figure(2);
　　plot(abs(y2));
　　title(‘含噪音的信号频谱’);
　　xlabel(‘频率’);
　　ylabel(‘频谱幅度’);
　　wp=0.25*pi;
　　ws=0.3*pi;
　　wdelta=ws-wp;
　　N=ceil(6.6*pi/wdelta); %取整
　　wn=(0.2+0.3)*pi/2;
　　b=fir1(N,wn/pi,hamming(N+1)); %选择窗函数，并归一化截止频率
　　f2=filter(bz,az,x2)
　　sound(f2,8000); %播放滤波后的语音信号
　　y3=fft(f2,);
　　figure(3)
　　y2=fft(x2);
　　subplot(3,1,1);
　　plot(abs(y2));
　　title('滤波前的频谱')
　　xlabel('频率');
　　ylabel('频谱幅度');
　　subplot(3,1,2)
　　f2=plot(abs(y3));
　　title('滤波后的频谱')
　　xlabel('频率');
　　ylabel('频谱幅度');
　　但运行到x2=x1+d是就说+错误之类
　　谁能该啊

pengzk

很可能x1 和d一个是列向量，一个是纵向量，如果是这样，加上一句 d= d'，试试！

chy19840208

　　这个程序怎么改啊?
　　本设计步骤如下：
　　1.录制一段自己的语音信号，并对录制的信号进行采样;画出采样后的语音信号的频谱图;
　　2. 插入噪音;画出含有噪音的频谱图;
　　3.给定滤波器的性能指标，采用窗函数法，划出滤波器的频域响应
　　4.用该滤波器对含噪音的信号进行滤波，画出滤波后信号的频谱;
　　5.回放语音信号。
　　程序如下：
　　fs=8000;
　　x1=wavread('E:\我的论文\论文\luyin.wav');
　　y1=fft(x1);
　　figure(1);
　　plot(abs(y1));
　　title('原始信号频谱');
　　xlabel('频率');
　　ylabel('频谱幅度');
　　t=0:1/8000:(length(x1)-1)/8000;
　　d=sin(2*pi*2000*t);
　　x2=x1+d';
　　sound(x2,fs);
　　y2=fft(x2);
　　figure(2);
　　title(‘含噪音的信号频谱’);
　　xlabel(‘频率’);
　　ylabel(‘频谱幅度’);
　　wp=0.25*pi;
　　ws=0.3*pi;
　　wdelta=ws-wp;
　　N=ceil(6.6*pi/wdelta); %取整
　　wn=(0.2+0.3)*pi/2;
　　b=fir1(N,wn/pi,hamming(N+1)); %选择窗函数，并归一化截止频率
　　f2=filter(bz,az,x2)
　　sound(f2,8000); %播放滤波后的语音信号
　　y3=fft(f2,);
　　figure(3)
　　y2=fft(x2);
　　subplot(3,1,1);
　　plot(abs(y2));
　　title('滤波前的频谱')
　　xlabel('频率');
　　ylabel('频谱幅度');
　　subplot(3,1,2)
　　f2=plot(abs(y3));
　　title('滤波后的频谱')
　　xlabel('频率');
　　ylabel('频谱幅度');
　　上面程序为什么在运行到plot(abs(y2));这步时,应该是原始信号和噪音信号叠加的混合信号啊?但是在上面只有噪音信号啊???????????

yangzj

对比下大小,是不是信号相对你的噪声太小了