之前只有贴上部分的程式码
没注意到部分的宣告在其他地方
对不起OTZ
很抱歉造成大家的困扰
感谢labbat大大!!!
以下是全部的程式码
clear all;
close all;
y= wavread('H:\MyVoice1.wav');
y1=25*y(:,1);
t=10/423936*[0:1:423935].';
sig=[t,y1];
x=0;
Lt=length(t);
tx=10/423936*[0:1:423935].';
Ltx=length(tx);
s=sin(4000*pi*tx);
Mf=fftshift(fft(y1,Lt));
Sf=fftshift(fft(s,Ltx));
figure(1)
subplot(211)
plot(t,y1)
xlabel('time,sec')
title('My Voice');
subplot(212)
plot(tx,s);
xlabel('time,sec')
title('Sinwave');
yt=s+y1;
load H:\BWLPF.mat G SOS
[b,a] = sos2tf(SOS,G)
z1=filter(b,a,yt);
load H:\BWHPF.mat q w
[d,c] = sos2tf(w,q)
z2=filter(d,c,yt);
figure(2);
subplot(311);
plot(t,yt)
xlabel('time,sec');
title('Recieved signal')
subplot(312);
plot(tx,z1);
xlabel('time,sec');
title('Recieved signal after lowpass filter')
subplot(313);
plot(tx,z2);
xlabel('time,sec');
title('Recieved signal after highpass filter')
figure(3);
subplot(211);
plot(abs(Mf));
title('Message spectrum')
subplot(212);
plot(abs(Sf));
title('sinwave spectrum')
如题我想要将音档用低通滤波器和高通滤波器绘出图
但是结果却是两个都是白的
请问能帮我看一下程式码有错误吗?
感激不尽!