關於數位濾波器 |
缺席
|
sharkneil
一般會員 發表:27 回覆:33 積分:12 註冊:2005-01-21 發送簡訊給我 |
|
㊣
版主 發表:261 回覆:2302 積分:1667 註冊:2005-01-04 發送簡訊給我 |
用示波器抓下來的波形應該會有兩組,一組time,一組data.
直接設計一2階濾波器,轉為轉移函數型態.再以示波器的data為輸入,得到的輸出就是慮波過的結果了.
也可以用simulink去做.輸入的部分改為變數(示波器data)就可以了.
記得simulink block 輸入型式改為matrix.
------
------------------------------------------------------------------------- 走是為了到另一境界,停是為了欣賞人生;未走過千山萬水,怎知生命的虛實與輕重!? |
sharkneil
一般會員 發表:27 回覆:33 積分:12 註冊:2005-01-21 發送簡訊給我 |
以下程式是利用MATLAB寫出來的~~可以參考看看: clear;
fid=0;%initialize file identifier
while fid<1
filename='v60.txt';
[fid,message]=fopen(filename,'r');%open file from disk
if fid==-1
disp(message)
end
end
data=fscanf(fid,'%g');%read data from file
%get the time point reference
data=data/16;
for u=1:length(data);
x(u)=u/192;
end; %subplot(2,1,1);
%plot(x,data)
%axis([0 1 0 1200]);
%grid status=fclose(fid) fs=192;
t=0:(1/fs):length(data)*(1/fs);
mL0=[1 1 0 0];
fHz0=[0 1 5 96];
f0=fHz0/(fs/2);
[bl,al]=yulewalk(6,f0,mL0);
fHz1=linspace(0,fs/2,50);
om1=2*pi*fHz1;
z=exp(sqrt(-1)*om1/fs);
mL=abs(polyval(bl,z)./polyval(al,z));
%plot(fHz0,mL0,fHz1,mL)
yl=filter(bl,al,data)
subplot(2,1,1);
plot(x,yl*23)%,x,data)
ylabel('battery votage(V)')
xlabel('time(sec)')
axis([0 60 5 8]);
%plot(data)
|
diroex
一般會員 發表:38 回覆:51 積分:22 註冊:2005-04-04 發送簡訊給我 |
|
sharkneil
一般會員 發表:27 回覆:33 積分:12 註冊:2005-01-21 發送簡訊給我 |
本站聲明 |
1. 本論壇為無營利行為之開放平台,所有文章都是由網友自行張貼,如牽涉到法律糾紛一切與本站無關。 2. 假如網友發表之內容涉及侵權,而損及您的利益,請立即通知版主刪除。 3. 請勿批評中華民國元首及政府或批評各政黨,是藍是綠本站無權干涉,但這裡不是政治性論壇! |