資源簡介
用matlab做的基于小波降噪和自適應(yīng)濾波的程序,已經(jīng)調(diào)試可以使用。
代碼片段和文件信息
clear?all
clc
%%
load(‘dls202.txt‘);?
x1=dls202(1:102401);?t=[1:10240];fs=20000;
figure(1);subplot(331);plot(tx1);title(‘原始信號(大理石)‘);xlabel(‘樣本序號/?n‘);ylabel(‘幅度/?V‘);
xd1=wden(x1‘minimaxi‘‘s‘‘one‘5‘db5‘);
subplot(332);plot(xd1);title(‘經(jīng)小波降噪后輸出信號(大理石)‘);xlabel(‘樣本序號/?n‘);ylabel(‘幅度/?V‘);
%%%%%%%%%%%自適應(yīng)濾波算法實現(xiàn)%%%%
n=10240;??%輸入信號抽樣點數(shù)N
a=1;
Signal?=?xd1;?
%%%%%%%%%%%%%%%%%%自適應(yīng)濾波LMS算法仿真%%%%%%%%%%%%%%%%
g=100;??%統(tǒng)計仿真次數(shù)為g
N=n;??%輸入信號抽樣點數(shù)N
k=128;???%時域抽頭LMS算法濾波器階數(shù)
pp=zeros(gN-k);??%將每次獨立循環(huán)的誤差結(jié)果存于矩陣pp中,以便后面對其平均
u=0.001;
for?q=1:g
????%設(shè)置初值
????yn_1=zeros(1N);???%output?signal
????yn_1(1:k)=Signal(1:k);??%將輸入信號SignalAddNoise的前k個值作為輸出yn_1的前k個值
????w=zeros(1k);???%設(shè)置抽頭加權(quán)初值
????e=zeros(1N);???%誤差信號
????
????%用LMS算法迭代濾波
????for?i=(k+1):N
????????XN=Signal((i-k+
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????3879??2012-11-08?20:19??lbjz.m
-----------?---------??----------?-----??----
?????????????????3879????????????????????1
評論
共有 條評論