我如何设计使用LMS自适应滤波器算法? ?
20视图(30天)
显示旧的评论
答案(2)
sandeep URITI
2014年2月16日
clc清晰;关闭所有;
Fc = 4000000;%载波频率= 1 mhz f = Fc;Fs = 16000000;%采样频率fm = 700000;Fd =调频;传感器= 4;%的传感器数量angle_tx =π/ 4;%发射机= 45度角angle_jam =π/ 4;%干扰机= 45度角sens_wts = (0.2 - 0.4 0.6 - 0.8);%分配传感器加权c = 3 e08; lembda=c/f; % Transmitted signal wavelength
桑普= 2 * (6 * Fs / fm);%的样本数据max = (1 / Fs) * (samps-1);t = 0:1 / Fs: (max);
%调制数据modsignal =罪(2 *π*调频* t);%基带信号modsignal (modsignal > 0) = 1;modsignal (modsignal < 0) = 1;车=罪(2 *π* Fc * t);%载体t_sig = modsignal。*车;ch_samps = %调制数据长度(t_sig);%调制信号数据样本g = 0: ch_samps-1;图(2)%;%的阴谋(g, t_sig);%情节传播信号%标题(调制信号传输); grid on; % axis([0 8e6 0 1]);
d = lembda / 2;%传感器分离并= 15 e-6;%步长
%传输信号的FFT k = ch_samps;fft_samps = 2 ^ nextpow2 (k);t_fft = fft (t_sig fft_samps) / k;fprime = f / 2 * linspace (0, 1, fft_samps / 2);图(3);次要情节(1,1)%情节传输信号的频率成分的情节(fprime 2 * abs (t_fft (1: fft_samps / 2)));网格;标题(“原始传输信号频谱”);包含(Hz频率);ylabel(级); axis([0 8e6 0 1]) %Jammer Signal t=0:1/Fs:max; j_sig=1*sin(2*pi*Fc*t);
%干扰信号的FFT fft_samps = 2 ^ nextpow2 (k);j_fft = fft (j_sig fft_samps) / k;fprime = f / 2 * linspace (0, 1, fft_samps / 2);图(4)%;%的阴谋(fprime 2 * abs (j_fft (1: fft_samps / 2)));%画出干扰信号的频率成分%网格;%的标题(“干扰信号的频谱”);%包含(频率(赫兹));% ylabel(级);%轴([0 8 e6 0 1])
%阵列传播向量t2 = 1:传感器;v (t2) = exp(我* (t2-1) * 2 *π* d * sin (angle_tx) * 1 / lembda); %传播向量
结束
t3 = 1:传感器;eeta (t3) = exp(我* (t3-1) * 2 *π* d * sin (angle_jam) * 1 / lembda); %传播向量对干扰信号
%的干扰机接收传感器阵列j_rcvd = j_sig ' * eeta;%干扰信号接收@传感器
% j_sig + t_sig =数据接收后@ anteena x = t_sig ' * v + j_rcvd;
%收到信号的FFT @传感器fft_samps = 2 ^ nextpow2 (k);fft_samps x_fft = fft (x) / k;fprime = f / 2 * linspace (0, 1, fft_samps / 2);图(3);次要情节(3,1,2)情节(fprime 2 * abs (x_fft (1: fft_samps / 2)));%的阴谋Rxd信号的频率成分网格;标题(“信号在接收天线(干扰信号+期望信号)');包含(频率(赫兹));ylabel(级);轴([0 8 e6 0 1])
% LMS算法n1 = 1: ch_samps间= sens_wts * x ';E = t_sig-x_est;sens_wts = sens_wts +(并* E * x);结束
%的FFT估计fft_samps = 2 ^ nextpow2 (k);fft_samps z_fft = fft(间)/ k;fprime = f / 2 * linspace (0, 1, fft_samps / 2);图(3);次要情节(3,1,3);情节(fprime 2 * abs (z_fft (1: fft_samps / 2)));%情节估计信号的频率成分网格;标题(估计信号频谱的);包含(频率(赫兹));ylabel(级); axis([0 8e6 0 1]) %************************** %Scatter plot of estimated signal
% y_est = ddemodce (Fd间,f * 2、相移键控,2);% % y_est = circshift (y_ester, 3) %图(1);%次要情节(2,1,2);%干(y_est,“填满”),网格;%轴([2 101 0 1]);
