【优化算法】改进型的LMS算法-SVSLMS算法【含Matlab源码 632期】
一、获取代码方式
获取代码方式1:
完整代码已上传我的资源:【优化算法】改进型的LMS算法-SVSLMS算法【含Matlab源码 632期】
获取代码方式2:
通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码。
备注:
订阅紫极神光博客付费专栏,可免费获得1份代码(有效期为订阅日起,三天内有效);
二、简介
最小均方(LMS, Least Mean Squares)是最基本的自适应滤波算法。
LMS算法是自适应滤波器中常用的一种算法与维纳算法不同的是其系统的系数随输入序列而改变。维纳算法中截取输入序列自相关函数的一段构造系统的最佳系数。而LMS算法则是对初始化的滤波器系数依据最小均方误差准则进行不断修正来实现的。因此理论上讲LMS算法的性能在同等条件下要优于维纳。但是LMS是在初始值下逐步调整的,因此在系统稳定前,会有一段调整时间,调整时间受步长因子的控制,一定范围内,步长因子越大,调整时间越小,步长因子的最大取值为R的迹。LMS采用平方误差最小的原则代替均方误差最小的原则,信号基本关系如下:
三、部分源代码
%%------------------- 定义变量:----------------------%%
% N - # of elements in array %
% d - 阵元间距 %
% ang - theta in deg %
% thetaS - 期望用户波达角 %
% thetaI - 干扰波达角 %
% T - 期望信号的周期 %
% t - 期望信号的时间轴 %
% S - 期望信号 %
% I - 干扰信号 %
% vS,vI - 期望用户和干扰用户的转移矢量 %
% X - 总阵列因子 %
% Rxx - 总接收信号相关矩阵 %
% mu - 收敛参数 %
% w - 用LMS算法确定的均匀线阵权值 %
% x - 总接收信号 %
% y - 阵列输出 %
% e - 阵列输出和期望信号的误差 %
% theta - range of AOA's (rad) %
% AF - 加权阵列输出 %
%%-----------------------------------------------------------%%%%----- 赋值-----%%N = 8; d =0.5; thetaS = 30*pi/180; thetaI = -60*pi/180;%期望信号方向30. 干扰-60.%%%----- 期望和干扰信号-----%%T = 1E-3; t = (1:100)*T/100; it = 1:100;%1E-3表示1ms%S = cos(2*pi*t/T); %,此时,S已经是一个矩阵!1行100列I = randn(1,100);%生成1行100列的正态分布随机矩阵%%%----- 为每个用户的线性阵列信号创建的阵列因子 -----%%vS = []; vI = [];
i = 1:N;
vS = exp(1j*(i-1)*2*pi*d*sin(thetaS)).';%’表示厄米特转置,产生N行1列的阵列向量,然后后面在此基础上乘以信号就得到天线每个阵元接收到的信号
vI = exp(1j*(i-1)*2*pi*d*sin(thetaI)).';%%----- 用LMS解决权值 -----%%w = zeros(N,1);%初始天线阵权值全为0%
X = vS + vI;
Rxx = X*X';%X'为X的厄米特转置%
mu = 1/(4*abs(trace(Rxx)));%trace,迹%wi = zeros(N,max(it));%初始权值都为0%
for n = 1:length(S)%n是迭代次数?%x = S(n)*vS + I(n)*vI;y(n) = w'*x;e = conj(S(n)) - y(n); esave(n) = abs(e)^2;%conj,求复数的共轭%nu=mu*(1-exp(-(abs(e))^2));w = w + nu*conj(e)*x;wi(:,n) = w;
end %没有判断收敛,要判断是从误差收敛曲线来看
w = w/w(1); % 第一权值的规范解%%%----- 显示权值 -----%%disp(' ')
disp('%------------------------------------------------------------------------%')
disp(' ')
disp([' N = ',num2str(N),' 的权值是:'])
disp(' ')
for m = 1:length(w)disp([' w',num2str(m),' = ',num2str(w(m),3)])
end
disp(' ')%%----- 输出结果 -----%%% 1.) Plot 权重与迭代次数wi = wi.';
figure(1), plot(it,abs(wi(:,1)),'kx',it,abs(wi(:,2))...,'ko',it,abs(wi(:,3)),'ks',it,abs(wi(:,4)),'k+',...it,abs(wi(:,5)),'kd','markersize',2)
xlabel('Iteration no.'), ylabel('|weights|')
title('\bf 阵列权值')% 2.) Plot 信号采集和跟踪figure(2)
plot(it,S,'k',it,real(y),'k--')
xlabel('迭代次数'), ylabel('Signals')
title('\bf 期望信号采集与跟踪')
legend('期望信号','阵列输出')% 3.) Plot 最小均方误差figure(3), plot(it,esave,'k')
xlabel('Iteration no.'), ylabel('|e|^2')
title('\bf 均方误差与迭代次数.')% 4.) Plot 阵列因子theta = -pi/2:.01:pi/2;
AF = 0;
for i = 1:NAF = AF + w(i)'.*exp(j*(i-1)*2*pi*d*sin(theta));
end
四、运行结果
四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1] 包子阳,余继周,杨杉.智能优化算法及其MATLAB实例(第2版)[M].电子工业出版社,2016.
[2]张岩,吴水根.MATLAB优化算法源代码[M].清华大学出版社,2017.
【优化算法】改进型的LMS算法-SVSLMS算法【含Matlab源码 632期】相关推荐
- 【微电网优化】基于matlab粒子群算法求解综合能源系统优化问题【含Matlab源码 1969期】
一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[微电网优化]基于matlab粒子群算法求解综合能源系统优化问题[含Matlab源码 1969期] 点击上面蓝色字体,直接付费下载,即可. ...
- 【优化布局】基于matlab免疫算法求解充电站最优布局【含Matlab源码 2539期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[生产调度]基于matlab免疫算法求解生产调度零等待问题[含Matlab源码 1178期] 点击上面蓝色字体,直接付费下载,即可. 获取 ...
- 【数字信号去噪】基于matlab粒子群算法优化VMD分解分量选择数字信号降噪【含Matlab源码 1979期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[数字信号去噪]基于matlab粒子群算法优化VMD分解分量选择数字信号降噪[含Matlab源码 1979期] 点击上面蓝色字体,直接付费 ...
- 【优化布局】免疫算法求解充电站最优布局【含Matlab源码 2539期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[生产调度]基于matlab免疫算法求解生产调度零等待问题[含Matlab源码 1178期] 点击上面蓝色字体,直接付费下载,即可. 获取 ...
- 【雷达通信】基于matlab NCP算法SAR回波生成和成像【含Matlab源码 1185期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[雷达通信]基于matlab NCP算法SAR回波生成和成像[含Matlab源码 1185期] 点击上面蓝色字体,直接付费下载,即可. 获 ...
- 【配送路径规划】蚁群算法求解配送路径最短问题【含Matlab源码 2222期】
⛄一.VRP简介 1 VRP基本原理 车辆路径规划问题(Vehicle Routing Problem,VRP)是运筹学里重要的研究问题之一.VRP关注有一个供货商与K个销售点的路径规划的情况,可以简 ...
- 【配送路径规划】基于matlab蚁群算法求解配送路径最短问题【含Matlab源码 2222期】
⛄一.VRP简介 1 VRP基本原理 车辆路径规划问题(Vehicle Routing Problem,VRP)是运筹学里重要的研究问题之一.VRP关注有一个供货商与K个销售点的路径规划的情况,可以简 ...
- 【ABC三维路径规划】基于matlab人工蜂群算法多无人机三维路径规划【含Matlab源码 170期】
一.无人机简介 0 引言 随着现代技术的发展,飞行器种类不断变多,应用也日趋专一化.完善化,如专门用作植保的大疆PS-X625无人机,用作街景拍摄与监控巡察的宝鸡行翼航空科技的X8无人机,以及用作水下 ...
- 【BA三维路径规划】基于matlab改进的蝙蝠算法农用无人机三维路径规划【含Matlab源码 1514期】
一.无人机简介 0 引言 随着现代技术的发展,飞行器种类不断变多,应用也日趋专一化.完善化,如专门用作植保的大疆PS-X625无人机,用作街景拍摄与监控巡察的宝鸡行翼航空科技的X8无人机,以及用作水下 ...
- 【PSO三维路径规划】基于matlab粒子群算法无人机山地三维路径规划【含Matlab源码 1405期】
⛄一.无人机简介 0 引言 随着现代技术的发展,飞行器种类不断变多,应用也日趋专一化.完善化,如专门用作植保的大疆PS-X625无人机,用作街景拍摄与监控巡察的宝鸡行翼航空科技的X8无人机,以及用作水 ...
最新文章
- 干货 | 20个教程,掌握时间序列的特征分析(附代码)
- 给AI一张高清照片,分分钟还你细节满满的3D人体模型,GitHub标星3.6k | 在线可玩...
- MySQL5.7新特性——在线收缩undo表空间 (转载)
- php 二维数组 根据某个字段排序
- python之绝对导入和相对导入
- 什么从什么写短句_2020抖音文案短句:爱情、励志、伤感合集,值得收藏!
- Nagios 安装方法
- java 多态 显式隐式,Java 构造器中的显式参数和this隐式参数
- 51CTO博客——架起我与读者沟通、见面的桥梁[博友话题]
- 用互联网模式打造智慧城市的“城市云”
- click事件的执行顺序
- Centos--swoole平滑重启服务
- ExpandableListView 实现三级菜单中grou_item与child_item点击无响应
- 产品开发之小步快跑,快速迭代
- java类加载机制?双亲委派模型有可能被破坏吗
- Unity3D基础3:贴图与材质球
- 类似新浪微博中取消关注的弹出确认框
- 历尽磨难,探索未来生活的2014年
- 多个视频合成一个视频(带合成工具)
- python读取Excel中关联表格的数据(只要是同Excel中
热门文章
- ZOJ 1060 Count the Color
- Git学习笔记(2)-创建仓库
- Servlet API 中文版
- 禁止选择,右键菜单,拷贝,拖拽
- python opencv 创建滑动条调整值
- 20200203每日一句
- Atitit 软件开发体系法规大全v2.docx Atitit 软件开发体系大全 目录 1. 基本法(类似宪法)	1 2. 行政法	1 3. 流程法民商法	2 3.1. Ui提升法	2 3.2. 功
- Atitit everthing placeholder index list 目录 1.1. sumdoc 2019 zipver t1---t91	1 1.2. 00 sumdoc ever
- Atitit 云计算体系树.docx Atitit 云计算之道 attilax著 艾龙 著 serverless bomb 1. 什么才是云计算的根本特征..	2 2. 云计算体系	2 3. “云
- Atitit 获取数据库表主键功能的实现 数据库模块 艾提拉attilax总结 package com.attilax.sql; import java.sql.Connection; im