fs=44100;%语音信号采样频率为22050
x1=wavread('ss.wav');%读取语音信号的数据，赋给变量x1
sound(x1,44100);%播放语音信号
figure(1)
plot(x1);%做原始语音信号的时域图形
title('原始语音信号');
xlabel('时间/T');
ylabel('幅值/N');
>> fs=44100;
x1=wavread('ss.wav');
t=0:1/44100:(length(x1)-1)/44100;%将所加噪声信号的点数调整到与原始信号相同
Au=0.5;
d=Au*cos(2*pi*8000*t)';%噪声为8KHz的余弦信号
x2=x1+d;%信号的叠加
wp=0.25*pi;ws=0.3*pi;%通带和阻带的截止频率分别为0.25pi和0.3pi
Rp=1;Rs=15;%通带和阻带区的波纹系数分别为1和15
Fs=82050;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(2)
plot(W*Fs/(2*pi),abs(H))
xlabel('频率/Hz')
ylabel('频率响应幅度/N')
title('Butterworth')
??? Error using ==> plus
Matrix dimensions must agree.
谁能告诉我怎么改啊，我看不出来了！！！！
你x1长度为多少？？d的长度为多少。让两者长度一样就行了。。。。。