Question

Error in ==> Untitled2 at 4y=w/max(abs(w)); 错误？

T0=8000;
[w,Fs,bits]=wavread('wan.WAV');
N=length(w);
y=w/max(abs(w));
%sound(w,Fs,bits);
t=0:T0-1;
subplot(6,1,1);plot(t/Fs,real(y(1:T0)),'k');
xlabel('Time/s');ylabel('Amp');
axis([0 0.2 -0.5 1.3]);
z=fft(y);
H=z.*conj(z)/N;
f=0:N/2-1:subplot(6,1,2);plot(f*Fs/N,H(1:N/2),'k');
xlabel('Fre/s');ylabel('Amp');
axis([0 N/200 0 25]);
%%%%%%%%%%%%%
[w,Fs]=wavread('wan.WAV');
N=length(w);
SNR=10;
y=awgn(y,SNR);
sound(y,Fs,bits);
z=fft(y);
H=z.*conj(z)/N;
t=0:T0-1;
subplot(6,1,3);plot(t/Fs,real(y(1:T0)),'k');
xlabel('Time/s');ylabel('Amp');
axis([0 0.2 -0.5 1.3]);
f=0:N/2-1;subplot(6,1,4);plot(f*Fs/N,H(1:N/2),'k');
xlabel('Fre/s');ylabel('Amp');
axis([0 N/200 0 25]);
%%%%%%%%%%%%%
fp=1000;fs=1200;As=50;
Wp=2*pi*fp/Fs;Ws=2*pi*fs/Fs;
Nj=ceil(6.6*pi/(Ws-Wp));
Nj=mod(Nj+1,2)+Nj;
w=hamming(Nj);
Wc=(Wp+Ws)/2/pi;
b=firl(Nj-1,Wc,w);
[x,Fs]=wavread('wan.WAV');
N=length(x);
for i=1:N
y(i)=y(i)+(sin(2*pi*1500*i/Fs)+cos(2*pi*1600*i/Fs)+sin(2*pi*2000*i/Fs))/10;
end
YO=filter(b,l,y);
sound(YO,Fs,bits);
z=fft(YO);
H=z*conj(z)/N;
t=0:T0-1;
subplot(6,1,5);plot(t/Fs,real(YO(1:T0)),'k');
xlabel('Time/s');ylabel('Amp');
f=0:Nj/2-1;subplot(6,1,6);plot(f*Fs/Nj,H(1:Nj/2),'k');
xlabel('Fre/Hz');ylabel('Amp');
axis([0 Nj/2 0 25]);
哪里错了?

Answer 1

T0=8000;
[w,Fs,bits]=wavread('wan.WAV');
N=length(w);
y=w/max(abs(w));
%sound(w,Fs,bits);
t=0:T0-1;
subplot(6,1,1);plot(t/Fs,real(y(1:T0)),'k');
xlabel('Time/s');ylabel('Amp');
axis([0 0.2 -0.5 1.3]);
z=fft(y);
H=z.*conj(z)/N;
f=0:N/2-1:subplot(6,1,2);plot(f*Fs/N,H(1:N/2),'k');
xlabel('Fre/s');ylabel('Amp');
axis([0 N/200 0 25]);
%%%%%%%%%%%%%
[w,Fs]=wavread('wan.WAV');
N=length(w);
SNR=10;
y=awgn(y,SNR);
sound(y,Fs,bits);
z=fft(y);
H=z.*conj(z)/N;
t=0:T0-1;
subplot(6,1,3);plot(t/Fs,real(y(1:T0)),'k');
xlabel('Time/s');ylabel('Amp');
axis([0 0.2 -0.5 1.3]);
f=0:N/2-1;subplot(6,1,4);plot(f*Fs/N,H(1:N/2),'k');
xlabel('Fre/s');ylabel('Amp');
axis([0 N/200 0 25]);
%%%%%%%%%%%%%
fp=1000;fs=1200;As=50;
Wp=2*pi*fp/Fs;Ws=2*pi*fs/Fs;
Nj=ceil(6.6*pi/(Ws-Wp));
Nj=mod(Nj+1,2)+Nj;
w=hamming(Nj);
Wc=(Wp+Ws)/2/pi;
b=firl(Nj-1,Wc,w);
[x,Fs]=wavread('wan.WAV');
N=length(x);
for i=1:N
y(i)=y(i)+(sin(2*pi*1500*i/Fs)+cos(2*pi*1600*i/Fs)+sin(2*pi*2000*i/Fs))/10;
end
YO=filter(b,l,y);
sound(YO,Fs,bits);
z=fft(YO);
H=z*conj(z)/N;
t=0:T0-1;
subplot(6,1,5);plot(t/Fs,real(YO(1:T0)),'k');
xlabel('Time/s');ylabel('Amp');
f=0:Nj/2-1;subplot(6,1,6);plot(f*Fs/Nj,H(1:Nj/2),'k');
xlabel('Fre/Hz');ylabel('Amp');
axis([0 Nj/2 0 25]);

追问

亲 你改的哪里啊？ 还是复制一下的？