Question

MATLAB编程问题，拜托了。

具体分为以下步骤：
（1）语音信号的采集：利用Windows下的录音机，录制一段话音。然后在Matlab软件平台下，利用函数wavread对语音信号进行采样，播放语音信号，并绘制原始语音信号；
（2）对原始信号加入噪声：对原始语音信号加入
s=sin(2*pi*f*Ts*n)
的噪声，采样后可知Fs = 16000，选择f = 2500，播放加入噪声信号的语音信号，并绘制噪声信号和含噪语音信号；
（3）频谱分析：分别对原始语音信号，噪声信号和含噪声的语音信号进行频谱分析，并绘出各频谱图；
（4）设计滤波器：计算滤波器的性能指标，设计滤波器，绘制滤波器的特性曲线；
（5）滤波器滤波：用自己设计的滤波器对采集的信号进行滤波，得出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化，并回放语音信号，感觉滤波前后的声音有变化。

可以的话，讲详细点吧，真是看书看的头都大了。
分我有的是，诚心求教而已，需要分能给多少就给多少。

Answer 1

帮助有限，有个别问题我也不是很明白
(1)这个可以解决
你可能要修改FS F，
fs=22050;
%语音信号采样频率为22050
x1=wavread('h:\课程设计2\shuzi.wav'); %读取语音信号的数据，赋给变量x1
sound(x1,22050); %播放语音信号
y1=fft(x1,1024); %对信号做1024点FFT变换
f=fs*(0:511)/1024;
figure(1)
plot(x1) %做原始语音信号的时域图形
title('原始语音信号');
xlabel('time n');
ylabel('fuzhi n');
figure(2)
freqz(x1) %绘制原始语音信号的频率响应图
title('频率响应图')
figure(3)
subplot(2,1,1);
plot(abs(y1(1:512))) %做原始语音信号的FFT频谱图
title('原始语音信号FFT频谱')
subplot(2,1,2);
plot(f,abs(y1(1:512)));
title('原始语音信号频谱')
xlabel('Hz');
ylabel('fuzhi');
（2）对原始语音信号加入
s=sin(2*pi*f*Ts*n)
你可以设s=sin(2*pi*f*Ts*n)修改下面的d值

fs=22050;
x1=wavread('h:\课程设计2\shuzi.wav');
f=fs*(0:511)/1024;
t=0:1/22050:(size(x1)-1)/22050; %将所加噪声信号的点数调整到与原始信号相同
Au=0.03;
d=[Au*cos(2*pi*5000*t)]'; %噪声为5kHz的余弦信号
x2=x1+d;
sound(x2,22050); %播放加噪声后的语音信号
y2=fft(x2,1024);
figure(1)
plot(t,x2)
title('加噪后的信号');
xlabel('time n');
ylabel('fuzhi n');
figure(2)
subplot(2,1,1);
plot(f,abs(y1(1:512)));
title('原始语音信号频谱');
xlabel('Hz');
ylabel('fuzhi');
subplot(2,1,2);
plot(f,abs(y2(1:512)));
title('加噪后的信号频谱');
xlabel('Hz');
ylabel('fuzhi');
(3)上面已经有了
（4）这个你可以借鉴下面程序
程序3:双线性变换法设计Butterworth滤波器

fs=22050;
x1=wavread('h:\课程设计2\shuzi.wav');
t=0:1/22050:(size(x1)-1)/22050;
Au=0.03;
d=[Au*cos(2*pi*5000*t)]';
x2=x1+d;
wp=0.25*pi;
ws=0.3*pi;
Rp=1;
Rs=15;
Fs=22050;
Ts=1/Fs;
wp1=2/Ts*tan(wp/2); %将模拟指标转换成数字指标
ws1=2/Ts*tan(ws/2);
[N,Wn]=buttord(wp1,ws1,Rp,Rs,'s'); %选择滤波器的最小阶数
[Z,P,K]=buttap(N); %创建butterworth模拟滤波器
[Bap,Aap]=zp2tf(Z,P,K);
[b,a]=lp2lp(Bap,Aap,Wn);
[bz,az]=bilinear(b,a,Fs); %用双线性变换法实现模拟滤波器到数字滤波器的转换
[H,W]=freqz(bz,az); %绘制频率响应曲线
figure(1)
plot(W*Fs/(2*pi),abs(H))
grid
xlabel('频率／Hz')
ylabel('频率响应幅度')
title('Butterworth')
f1=filter(bz,az,x2);
figure(2)
subplot(2,1,1)
plot(t,x2) %画出滤波前的时域图
title('滤波前的时域波形');
subplot(2,1,2)
plot(t,f1); %画出滤波后的时域图
title('滤波后的时域波形');
sound(f1,22050); %播放滤波后的信号
F0=fft(f1,1024);
f=fs*(0:511)/1024;
figure(3)
y2=fft(x2,1024);
subplot(2,1,1);
plot(f,abs(y2(1:512))); %画出滤波前的频谱图
title('滤波前的频谱')
xlabel('Hz');
ylabel('fuzhi');
subplot(2,1,2)
F1=plot(f,abs(F0(1:512))); %画出滤波后的频谱图
title('滤波后的频谱')
xlabel('Hz');
ylabel('fuzhi');

程序4:窗函数法设计滤波器：

fs=22050;
x1=wavread('h:\课程设计2\shuzi.wav');
t=0:1/22050:(size(x1)-1)/22050;
Au=0.03;
d=[Au*cos(2*pi*5000*t)]';
x2=x1+d;
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)); %选择窗函数，并归一化截止频率
figure(1)
freqz(b,1,512)
f2=filter(bz,az,x2)
figure(2)
subplot(2,1,1)
plot(t,x2)
title('滤波前的时域波形');
subplot(2,1,2)
plot(t,f2);
title('滤波后的时域波形');
sound(f2,22050); %播放滤波后的语音信号
F0=fft(f2,1024);
f=fs*(0:511)/1024;
figure(3)
y2=fft(x2,1024);
subplot(2,1,1);
plot(f,abs(y2(1:512)));
title('滤波前的频谱')
xlabel('Hz');
ylabel('fuzhi');
subplot(2,1,2)
F2=plot(f,abs(F0(1:512)));
title('滤波后的频谱')
xlabel('Hz');
ylabel('fuzhi');
这些程序中应该有一个小错误，我整理的时候发现不能直接使用公式x2=x1+d
（5）综合上述，你可以借鉴下面的程序
fs=22050;
%语音信号采样频率为22050
x1=wavread('Windows Critical Stop.wav'); %读取语音信号的数据，赋给变量x1
sound(x1,22050); %播放语音信号
y1=fft(x1,1024); %对信号做1024点FFT变换
f=fs*(0:511)/1024;

figure(1)
plot(x1) %做原始语音信号的时域图形
title('原始语音信号');
xlabel('time n');
ylabel('fuzhi n');

x2=randn(1,length(x1)); %产生一与x长度一致的随机信号
sound(x2,22050);
figure(2)
plot(x2) %做原始语音信号的时域图形
title('高斯随机信号');
xlabel('time n');
ylabel('fuzhi n');

figure(3)
x3=x1+x2;
y3=fft(x3,1024);
f=fs*(0:511)/1024;
plot(t,x3)
title('加噪后的信号');
xlabel('time n');
ylabel('fuzhi n');

figure(4)
subplot(4,1,1);
plot(f,abs(y1(1:512)));
title('原始语音信号频谱');
xlabel('Hz');
ylabel('fuzhi');
subplot(4,1,2);
plot(f,abs(y3(1:512)));
title('加噪后的信号频谱');
xlabel('Hz');
ylabel('fuzhi');
这里还有上述出现的问题，公式x3=x2+x1不可用，有人说是类型不一样，不能直接相加，但我感觉是函数公式没定义，你找个人看看吧，我也只是整到这里。

Answer 2

100分解决这个问题……太少了吧？这个账号都有点少呢

Answer 3

%(1)
%例如录下了一段声音：voice.wav放在程序同一个文件夹下
[y,f,b]=wavread('voice.wav');
sound(f,y)
l=length(y);
t=(0:l-1)./f;
plot(t,y(:,1),'r',t,y(:,1),'g')
%(2)
F=2500;
s=sin(2*pi*F*Ts*n)%Ts=?????????????

Answer 4

用simulink可以