a.m

递归最小二乘滤波器算法

fid1=fopen('G:\研究生课程\随机信号处理\RLS\1.dat','r');
signal=fread(fid1,800,'short');
fclose(fid1);

fid2=fopen('G:\研究生课程\随机信号处理\RLS\3.dat','r');
noisy=fread(fid2,800,'short');
fclose(fid2);

fid3=fopen('G:\研究生课程\随机信号处理\RLS\out.dat','r');
out=fread(fid3,800,'short');
fclose(fid3);

fid4=fopen('G:\研究生课程\随机信号处理\RLS\mse.dat','r');
mse=fread(fid4,800,'double');
fclose(fid4);

fid5=fopen('G:\研究生课程\随机信号处理\RLS\2.dat','r');
noise=fread(fid5,800,'short');
fclose(fid5);

figure(1);
subplot(3,1,1); plot(signal); title('Clean Signal Waveform'); xlabel('Time'); ylabel('Amplitude');
subplot(3,1,2); plot(noise); title('Noise Waveform'); xlabel('Time'); ylabel('Amplitude');
subplot(3,1,3); plot(noisy); title('Noisy Signal Waveform'); xlabel('Time'); ylabel('Amplitude');

figure(3);
subplot(2,1,1); plot(signal,'b'); hold; plot(out, 'r'); title('Signal Waveform: Blue for clean signal, Red for filtered signal'); xlabel('Time'); ylabel('Amplitude');
subplot(2,1,2); plot(mse); title('Squared error');