声卡采样滤波后怎么输出来就变化呢

theason

请问高手，怎么我采样之后经声卡输出在示波器里看到的波形和采样滤波后的信号不一样呢？请高手帮忙

不知道怎么把采样滤波后的信号原样输出来，谢谢。

程序输出部分有问题吗？







ai = analoginput('winsound');
ch=addchannel(ai,1);
set(ai,'SampleRate',44100)
set(ai,'SamplesPerTrigger',44100*5);
start(ai)
[data,time] = getdata(ai);
delete(ai)



figure;
subplot(211),plot(time,data); %时域图
title('Time-Domain');
xlabel('time /s');ylabel('magnitude /V');


%cheby2型滤波
[b,a] = cheby2(9,20,300/22050);%设计9阶，带阻波纹20dB，不
                  %确定是从20开始，300/22050为归一化频率
[fil,zf]=filter(b,a,data); %zf为滤波器延迟
subplot(212),plot(time,fil);title('fitered signal');
xlabel('time /s');ylabel('magnitude /V');

                             
%进行FFT变换并做频谱图
y=fft(data);%进行fft变换
mag=abs(y);%求幅值
mag1=20*log10(mag);
f=(0:length(y)-1)'*44100/length(y);%进行对应的频率转换
figure;
subplot(211);
plot(f(1:64),mag1(1:64));%做频谱图
%axis([0,150,0,500]);
xlabel('Frequency /Hz');
ylabel('Magnitude /dB');
title('Samples FFT N=128');
grid;

%绘制功率谱
sq=mag1.^2;
subplot(212);
plot(f(1:64),sq(1:64));
xlabel('Frequency(Hz)');
ylabel('Power Spec');
title('Spec');
grid;



%处理信号输出
ao = analogoutput('winsound');
addchannel(ao,1);
set(ao,'SampleRate',44100);
set(ao,'TriggerType','Manual')
putdata(ao,fil);
start(ao);
trigger(ao);
waittilstop(ao,6)
stop(ao);
delete(ao);

[ 本帖最后由 sigma665 于 2008-1-16 10:46 编辑 ]

sigma665

版主不是万能的
耐心等待

花如月

原帖由 theason 于 2008-1-16 10:40 发表

                               
登录/注册后可看大图

请问高手，怎么我采样之后经声卡输出在示波器里看到的波形和采样滤波后的信号不一样呢？请高手帮忙

不知道怎么把采样滤波后的信号原样输出来，谢谢。

程序输出部分有问题吗？







ai = analoginput ...

如果采样滤波后的信号是通过matlab分析后在图像窗口输出的，那么先检查你采的原始信号和示波器的输入信号是否有区别。如果区别不大，说明原始信号没问题。那就检查是否是滤波器设计不合理导致的波形失真，一步一步来信号源到处理的每个步骤来排除可能的错误。

theason

谢谢两位回答。
我检查了传感器信号和滤波后的信号相差不大，但是在声卡扬声器口输出时就变了
通过改动不同的ao的SampleRate变化也很大（示波器接在声卡输出端看到），通常是ao的SampleRate越大输出也越大，失真更厉害；如果不设置SampleRate值也有输出，却不清楚这个时候的SampleRate值是多大了
而现在我的目的就是想把那滤波后的信号原样输出来，私底下认为是信号输出程序有问题，不知道哪位有这个方面的经验，请赐教，拜谢

w89986581

排出程序问题的话，只可能两个原因：1采用率不够；2滤波器阈值设置有问题。

theason

采样频率的话，因为采用声卡，只有几个档位：8000，11025，22050，44010
按照采样定理，44100足够了，通常实验室的干扰都在100Hz内，即使超出也不会超出多少。

至于所说的滤波器，我采用的cheby2型带通滤波，没有涉及到那个值的设置吧，直接运用的matlab函数了

不过还是谢谢你的意见
我估计是输出程序的问题
还有一个请教是哪位有没有用matlab实时控制声卡输出的经验呢？
现在做了一个实时采集的，但是不会实时输出，请高手指教

[ 本帖最后由 eight 于 2008-1-23 18:38 编辑 ]

songzy41

楼主是用44100Hz通过声卡采集数据，就是说信号最高频率可达22050Hz。而楼主又用了一个低通滤波器，低通滤波器的截止频率是300Hz，所以滤波后的波形当然和原始采样的波形有很大的差别，高频都被滤去了。

theason

现在问题不是这个，而是想通过MATLAB编程实时输出实时采集到的信号，不知道哪位有办法？

theason

实时采集能够完成，但实时输出没有成功，希望各位能人指导，谢谢

songzy41

在实时采集后
ai = analoginput('winsound');
ch=addchannel(ai,1);
set(ai,'SampleRate',44100)
set(ai,'SamplesPerTrigger',44100*5);
start(ai)
[data,time] = getdata(ai);
delete(ai)
加
wavplay(data,44100);
便能输出声音了.

theason

谢谢，这样能实时输出吗？我实时采集参考别人的程序
我试试

theason

好像这个方案不行。实时采集程序都是采用回调函数处理的

[ 本帖最后由 eight 于 2008-1-29 15:36 编辑 ]

theason

希望高手指教哦

theason

大家都没有做个这个方面的吗？:'( :'(