matlab加余弦噪声[x,fs,nbits]=wavread('chi1.wav');fs=16000;M =length (x) ; t=0:1/16000:(M-1)/16000;x1=x+cos(3*pi*t);sound(x1,fs)X=fft(x1,8192);N=8192;n=0:N-1;q=n*2*pi/N;figure(2)subplot(2,1,1);plot(x1);title('add noise 信号波形');subplot(2,1,

matlab加余弦噪声[x,fs,nbits]=wavread('chi1.wav');fs=16000;M =length (x) ; t=0:1/16000:(M-1)/16000;x1=x+cos(3*pi*t);sound(x1,fs)X=fft(x1,8192);N=8192;n=0:N-1;q=n*2*pi/N;figure(2)subplot(2,1,1);plot(x1);title('add noise 信号波形');subplot(2,1,
matlab加余弦噪声
[x,fs,nbits]=wavread('chi1.wav');
fs=16000;
M =length (x) ;
t=0:1/16000:(M-1)/16000;
x1=x+cos(3*pi*t);
sound(x1,fs)
X=fft(x1,8192);
N=8192;n=0:N-1;
q=n*2*pi/N;
figure(2)
subplot(2,1,1);plot(x1);title('add noise 信号波形');
subplot(2,1,2);plot(abs(X)); title('add noise信号频谱');axis([0 9000 0 8])
到底哪错了,怎么改啊

matlab加余弦噪声[x,fs,nbits]=wavread('chi1.wav');fs=16000;M =length (x) ; t=0:1/16000:(M-1)/16000;x1=x+cos(3*pi*t);sound(x1,fs)X=fft(x1,8192);N=8192;n=0:N-1;q=n*2*pi/N;figure(2)subplot(2,1,1);plot(x1);title('add noise 信号波形');subplot(2,1,
(1)利用Matlab对MIT-BIH数据库提供的数字心电信号进行读取,并还原实际波形.
(2)对数字信号数据做一次线性插值,使其成为均匀数字信号,以便后面的信号分析.
(3)根据心电信号的频域特征（自己查阅相关资料）,设计相应的低通和高通滤波器.
(4)编程绘制实现信号处理前后的频谱,做频谱分析,得出相关结论.