Question

在线等~!!!MATLAB中的问题!!!跪求!!!

有错误提示
Error in ==> qq1 at 11
[B,A] = buffer(n,fc/(fs/2));
??? Buffer size N must be an integer-valued scalar double.

下面是我的程序
close all;clear;
Ts=1e-3;
r=0.5;
t=(-10e-3):1e-4:10e-3;
[num,den]=rcosine(1e3,1e4,'fir/sqrt',r,10);
data=sign(rand(1,5000)-0.5);
datain=[data;zeros(9,length(data))];
datain=reshape(datain,1,10*length(data));
wavout=filter(num,den,datain);
n=10;fs=10000;fc=500;
[B,A] = buffer(n,fc/(fs/2));
noise=0.00*randn(size(wavout));
wavout=wavout+noise;
wavout1=filter(B,A,wavout);
wavout2=0.2*[zero(1,24),wavout(1:length(wavout)-24)]+1*[zero(1,34),wavout(1:length(wavout)-34)]+0.4*[zero(1,44),wavout(1:length(wavout)-44)];
wavout=[wavout1(1:20000),wavout2(20001:length(wavout))];
wavout=filter(num,den,wavout);
recpj=sign(rec);
X=rec;
L=1;
B=[zeros(1,L),1,zeros(1,L)];
adjstep=0.01;
M=length(B);
N=length(X);
bb(N-M+1,1:length(B))=0;
for k=1:N-M+1;
y(k)=sum(X(k:k+M-1)*filplr(B));
XK=X(k:k+M-1);
ek=sign(y(k))-y(k);
e(k)=ek;
dlt=fliplr(XK)*ek*adjstep;
B=B+dlt;
bb(k,:)=B;
end
BB1=[bb(1900,:);zeros(9,2*L+1)];
BB1=reshape(BB1,1,10*(2*l+1));
BBout1=filter(BB1,1,wavout);
BB2=[B;ZEROS(9,2*L+1)];
BB2=reshape(BB2,1,10*(2*L+1));
BBout2=filter(BB2,1,wavout);
figure(1);plot(e);title('误差曲线、epsilon_k');
figure(2);plot(bb);title('自适应均衡器抽头系数c_j');
eyediagram(BBout1(15000:20000),20);
eyediagram(BBout2(46001:50000),20);
eyediagram(wavout(10001:19000),20);
eyediagram(wavout(46001:50000),20);

Answer 1

1.函数调用出错.应该是butter,不是buffer.
2.wavout2=0.2*[zero(1,24),wavout(1:length(wavout)-24)]+1*[zero(1,34),wavout(1:length(wavout)-34)]+0.4*[zero(1,44),wavout(1:length(wavout)-44)];中所有的zero应改成zeros.
3.recpj=sign(rec);中rec变量之前没有用过.是不是漏掉了?
4.y(k)=sum(X(k:k+M-1)*filplr(B));中filplr函数是你自己定义的?
5.BB1=reshape(BB1,1,10*(2*l+1));中l为大写L.
6.BB2=[B;ZEROS(9,2*L+1)];中ZEROS改成小写zeros.
如果你要用百度hi,请在你留言之后注意别人是否回复.