这个程序怎么改啊？
本帖最后由 wdhd 于 2016-9-6 13:15 编辑本设计步骤如下：
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是就说+错误之类
谁能该啊
很可能x1 和d一个是列向量，一个是纵向量，如果是这样，加上一句 d= d'，试试！
问题
本帖最后由 wdhd 于 2016-9-6 13:15 编辑这个程序怎么改啊?
本设计步骤如下：
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));这步时,应该是原始信号和噪音信号叠加的混合信号啊?但是在上面只有噪音信号啊???????????
对比下大小,是不是信号相对你的噪声太小了
页:
[1]