關於數位濾波器 |
|
缺席
|
|
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 發送簡訊給我 |
發表時間：2005-06-23 19:37:25
IP:203.71.xxx.xxx 未訂閱
系統時間：2024-05-04 10:01:02
聯絡我們 | Delphi K.Top討論版
| 本站聲明 |
|
1. 本論壇為無營利行為之開放平台，所有文章都是由網友自行張貼，如牽涉到法律糾紛一切與本站無關。 2. 假如網友發表之內容涉及侵權，而損及您的利益，請立即通知版主刪除。 3. 請勿批評中華民國元首及政府或批評各政黨，是藍是綠本站無權干涉，但這裡不是政治性論壇！ |