chy19840208 发表于 2007-4-24 21:47

这个程序怎么改啊？

本设计步骤如下：
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是就显示＂+＂错误之类的
谁能给改一下啊


eight 发表于 2007-4-24 21:48

原帖由 chy19840208 于 2007-4-24 21:47 发表 http://forum.vibunion.com/forum/images/common/back.gif
本设计步骤如下：
1.录制一段自己的语音信号，并对录制的信号进行采样；画出采样后的语音信号的频谱图；
   2. 插入噪音；画出含有噪音的频谱图；
   3.给定滤波器的性能指标，采用窗函数法，划出滤波器的频 ...


请参考置顶贴：聚宝盆给出错误的完整信息

xjzuo 发表于 2007-4-25 09:17

注意向量长度是否匹配.
从你的程序看,x1,y1需要作一个转置．

[ 本帖最后由 xjzuo 于 2007-4-25 10:01 编辑 ]

chy19840208 发表于 2007-4-26 10:10

问题

本设计步骤如下：
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));这步时,应该是原始信号和噪音信号叠加的混合信号啊?但是在上面只有噪音信号啊???????????

eight 发表于 2007-4-26 10:22

原帖由 chy19840208 于 2007-4-26 10:10 发表 http://forum.vibunion.com/forum/images/common/back.gif
本设计步骤如下：
1.录制一段自己的语音信号，并对录制的信号进行采样；画出采样后的语音信号的频谱图；
   2. 插入噪音；画出含有噪音的频谱图；
   3.给定滤波器的性能指标，采用窗函数法，划出滤波器的频 ...


我也没有看出来，自己调试一下吧

寻找鱼的水 发表于 2007-11-22 22:19

回复 #4 chy19840208 的帖子

楼主的帖子似乎错误还比较多啊
title(‘含噪音的信号频谱’);
xlabel(‘频率’);
ylabel(‘频谱幅度’);
这里单引号应该是半角的；

y3=fft(f2,);似乎应该改成y3=fft(f2);去掉逗号

f2=filter(bz,az,x2) 这里不知道老是提示说bz ， az没有定义，这里看不懂，至于后面的就不好说了
也许我用的信号过于简单……
页: [1]
查看完整版本: 这个程序怎么改啊？