作者CatEar (CatEar)
看板MATLAB
标题[问题] 将音档用滤波器绘图(更新)
时间Mon Apr 7 18:30:16 2014
之前只有贴上部分的程式码
没注意到部分的宣告在其他地方
对不起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')
如题我想要将音档用低通滤波器和高通滤波器绘出图
但是结果却是两个都是白的
请问能帮我看一下程式码有错误吗?
感激不尽!
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.135.32.184
※ 文章网址: http://webptt.com/cn.aspx?n=bbs/MATLAB/M.1396866618.A.7D6.html
1F：推 labbat:z1 z2没有定义 04/07 21:17
※ 编辑: CatEar (140.135.32.184), 04/07/2014 23:28:50
※ 编辑: CatEar (140.135.32.184), 04/07/2014 23:31:47