海潮改正matlab,求改正LS和LMMSE算法的matlab程序,。
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
function output=ls_estimation(input,pilot_inter,pilot_sequence,pilot_num);
[N,NL]=size(input);
output=zeros(N,NL-pilot_num);
i=1;
count=0;
while i<=NL
Hi=input(:,i)./pilot_sequence;
count=count+1;
if count*pilot_inter<=(NL-pilot_num)
for j=((count-1)*pilot_inter+1):count*pilot_inter
output(:,j)=input(:,(i+j-(count-1)*pilot_inter))./Hi;
end
else
for j=((count-1)*pilot_inter+1):(NL-pilot_num)
output(:,j)=input(:,(i+j-(count-1)*pilot_inter))./Hi;
end
end
i=i+pilot_inter+1;
end
这是LS的
function output=lmmse_estimation(input,pilot_inter,pilot_sequence,pilot_num,trms,t_max,snr);
%trms为多经信道的平均延时,t_max为最大延时,此处所有的时间都是已经对采样间隔做了归一化后的结果
beta=17/9;
[N,NL]=size(input);
Rhh=zeros(N,N);
for k=1:N
for l=1:N
Rhh(k,l)=(1-exp((-1)*t_max*((1/trms)+j*2*pi*(k-l)/N)))./(trms*(1-exp((-1)*t_max/trms))*((1/trms)+j*2*pi*(k-l)/N));
end
end
output=zeros(N,NL-pilot_num);
i=1;
count=0;
while i<=NL
Hi=input(:,i)./pilot_sequence;
Hlmmse=Rhh*inv(Rhh+(beta/snr)*eye(N))*Hi;
count=count+1;
if count*pilot_inter<=(NL-pilot_num)
for j=((count-1)*pilot_inter+1):count*pilot_inter
output(:,j)=input(:,(i+j-(count-1)*pilot_inter))./Hlmmse;
end
else
for j=((count-1)*pilot_inter+1):(NL-pilot_num)
output(:,j)=input(:,(i+j-(count-1)*pilot_inter))./Hlmmse;
end
end
i=i+pilot_inter+1;
end
这是LMMSE的 老是程序跑不起来。。求大神gaiz
海潮改正matlab,求改正LS和LMMSE算法的matlab程序,。相关推荐
- matlab求偏微分方法解析解,偏微分数值解(2,MATLAB求解方法)学案.ppt
这部分主要讨论如何用MATLAB实现对偏微分方程的数值仿真求解.MATLAB的偏微分方程工具箱(PDE Toolbox)的出现,为偏微分方程的求解以及定性研究提供了捷径.主要步骤为: 2.1 用偏微分 ...
- matlab求分段函数的值.,如何用MATLAB求分段函数的最小值和最大值?
7.1.1 分段线性插值 所谓分段线性插值就是通过插值点用折线段连接起来逼近原曲线,这也是计算机绘制图形的基本原理.实现分段线性插值不需编制函数程序,MATLAB自身提供了内部函数interp1其主要 ...
- matlab整定串级pid,PID算法在Matlab串级控制中的应用
PID算法在Matlab串级控制中的应用 自114 1112002039 陈艳 前言:这个专题是由王娟老师给我们授课,我感觉收获挺大的,尤其是matlab仿真软件的使用,为我以后的实验课打下良好的基础 ...
- matlab 求向量的交集_从零开始的matlab学习笔记——(16)函数绘图
matlab应用--求极限,求导,求积分,解方程,函数绘图....更多内容尽在个人专栏:matlab学习 这一节开始我们正式进入函数绘图,内容比较繁杂,我尽量把它整合的简洁一些. 单一图像: 首先我们 ...
- matlab求二阶微分方程的通解,二阶常微分方程matlab的数值解和解析解分析总报告.ppt...
方程求解程序清单 a=-1,b=2,c=-1; w=1; m=2; n=1; h = 0.02; t=0:h:30; s1=dsolve('a*D2y+b*Dy+c*y=sin(w*t)','y(0) ...
- matlab求方程在X附近的根,matlab 实验03 求代数方程的近似根(解)
实验三 求代数方程的近似根(解) 求代数方程的根是最常见的数学问题之一(这里称为代数方程,主要是想和后面的微分方程区别开.为简明起见,在本实验的以下叙述中,把代数方程简称为方程),当是一次多项式时,称 ...
- 用Matlab求拐点和凹凸区间需要使用MATLAB
1.启动MATLAB,输入如下代码,求解函数曲线的一阶导数.二阶导数.拐点和拐点处的切线方程. close all; clear all; clc;format compact % MATLAB之考研 ...
- matlab求系统稳定时k的范围,Matlab大作业
一.通过举例说明运用MATLAB 判别控制系统稳定的所有方法 稳定是控制系统是否能进行工作的首要条件.一般来说,稳定性成为区分系统是否有用的标志.从实际应用的角度来看,可以认为只有稳定的系统才有用.而 ...
- matlab求不等式的方法,求解变分不等式的matlab程序我需要
2008-02-21 不等式的解法关于一元三次不等式的解法, 全日制普通高级中学教科书(试验修订本.必修)第一章第四节是"含绝对值的不等式解法",第五节是"一元二次不等式 ...
最新文章
- redis+主从复制+集群配置
- mysql 多个unique key_[MySQL]MySQL 中通过使用UNIQUE KEY 来控制字段值不重复的问题.
- 三维重建:闭环检测-相机闭环
- 基于 WPF 模块化架构下的本地化设计实践
- [转] android学习和广告平台赚钱
- Clickhouse 数据库连接工具——DBeaver
- 学习3 二维游戏动画合成(侠客行)
- 从0到10W粉,这一年我在csdn的成长之路
- 第十二节:通信之WLAN(Miracast)
- 贪心法 第1关:找零钱
- 光一个html文件能加图片吗,网站页面设计中光的特效设计
- 从零开始学C++之STL(七):剩下5种算法代码分析与使用示例(remove 、rotate 、sort、lower_bound、accumulate)
- 计算机显示器文字不清楚,笔记本电脑字体显示不清晰解决方法分辨率和显卡驱动...
- json文件怎么写注释
- AI智能化测试——各大厂质量保障实践分享汇总(下)
- mysql查询某学期开设的课程_查询KCXX表中2、3、4学期开设的课程情况
- adonis视图——nodejs应用第一季04
- Apple watch 回收如何选择
- 如何给下属安排工作?
- 联通java终面_中国联通面试经验