function?[ykalman]?=?ex1kalman（）
????load（‘Data/simdata.mat‘‘t‘‘y‘‘ynf‘‘ysin‘）;
????%?State?=?（y?ydot）.?We?only?observe?（y）.
????%?X（t+1）?=?F?X（t）?+?noise（Q）
????%?Y（t）?=?H?X（t）?+?noise（R）
????fprintf（1?‘-----------------Kalman-----------------\n‘）;
????ss?=?2;?%?state?size
????os?=?1;?%?observation?size
????F?=?[1?0;?0?1];?
????H?=?[1?0];
????Q?=?0.1*eye（ss）;
????R?=?1*eye（os）;
????initx?=?[0?0]‘;
????inity?=?1*eye（ss）;

????%[xfilt?Vfilt?VVfilt?loglik]?=?kalman_filter（y?F?H?Q?R?initx?inity）;
????ykalman?=?kalman_smoother（y?F?H?Q?R?initx?inity）;

????%?dfilt?=?x（1:）?-?xfilt（1:）;
????%?mse_filt?=?sqrt（sum（sum（dfilt.^2）））

%?????dsmooth?=?x（1:）?-?xsmooth（1:）;
%?????mse_smooth?=?sqrt（sum（sum（dsmooth.^2）））
????ykalman?=?ykalman（1:）;
????
????%?帶阻濾波，消除50Hz工頻干擾
?????%?經(jīng)頻譜分析可知，有用信號的頻率成分可能落在工頻帶，因此不能使用此方法
%?????Hd?=?bandstop;
%?????ykalman?=?filter（Hdykalman）;
????
????%計(jì)算信噪比
????eta?=?ynf?-?ykalman;?%?Residual?Signal
????snrpkalman?=?10*log（?（ynf*ynf‘）/（eta*eta‘）?）;???
%?????figure
%?????plot（tynf‘r‘tynf+ysin‘m‘tykalman‘b‘）;
%?????legend（‘ynf‘?‘ynf+ysin‘?‘ykalman‘）
%?????title（‘Kalman‘）;
????
????save（‘Data/kalman.mat‘‘ykalman‘‘snrpkalman‘）;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????69460??2013-07-11?22:15??ircbp2.1_wkp\Data\dictCA.mat

?????文件??????15987??2013-07-11?22:15??ircbp2.1_wkp\Data\realdata.mat

?????文件??????16063??2013-07-11?22:15??ircbp2.1_wkp\Data\simdata.mat

?????文件???????1337??2013-07-11?22:15??ircbp2.1_wkp\ex1kalman.m

?????文件????????748??2013-07-11?22:15??ircbp2.1_wkp\ex2sgolay.m

?????文件???????1065??2013-07-11?22:15??ircbp2.1_wkp\ex3bp.m

?????文件???????9925??1980-01-01?14:04??ircbp2.1_wkp\ex4ircbp.m

?????文件???????9955??1980-01-01?14:04??ircbp2.1_wkp\ex4ircbp.m~

?????文件???????1453??2013-07-11?22:15??ircbp2.1_wkp\getdata.m

?????文件???????2145??2013-07-11?22:15??ircbp2.1_wkp\getdictCA.m

?????文件???16393591??2013-12-07?23:32??ircbp2.1_wkp\toolbox\cvx-a64.tar.gz

-----------?---------??----------?-----??----

?????????????16521729????????????????????11