matlab减小数据方差,最小方差控制
各位大神好,根据matlab文件编了个相应的simulink仿真,但是输出一直发散,该怎么查原因啊,另外各位大神有自适应控制其他的simulink仿真吗?可否发我,急需!这是相应的matlab代码
%最小方差控制(MVC)
clear all; close all;
a=[1 -1.7 0.7]; b=[1 0.5]; c=[1 0.2]; d=4; %对象参数
na=length(a)-1; nb=length(b)-1; nc=length(c)-1; %na、nb、nc为多项式A、B、C阶次
nf=nb+d-1; %nf为多项式F的阶次
L=400; %控制步数
uk=zeros(d+nb,1); %输入初值:uk(i)表示u(k-i);
yk=zeros(na,1); %输出初值
yrk=zeros(nc,1); %期望输出初值
xik=zeros(nc,1); %白噪声初值
yr=10*[ones(L/4,1);-ones(L/4,1);ones(L/4,1);-ones(L/4+d,1)]; %期望输出
xi=sqrt(0.1)*randn(L,1); %白噪声序列
[e,f,g]=sindiophantine(a,b,c,d); %求解单步Diophantine方程
for k=1:L
time(k)=k;
y(k)=-a(2:na+1)*yk+b*uk(d:d+nb)+c*[xi(k);xik];%采集输出数据
u(k)=(-f(2:nf+1)*uk(1:nf)+c*[yr(k+d:-1:k+d-min(d,nc));yrk(1:nc-d)]-g*[y(k);yk(1:na-1)])/f(1);%求控制量
%更新数据
for i=d+nb:-1:2
uk(i)=uk(i-1);
end
uk(1)=u(k);
for i=na:-1:2
yk(i)=yk(i-1);
end
yk(1)=y(k);
for i=nc:-1:2
yrk(i)=yrk(i-1);
xik(i)=xik(i-1);
end
if nc>0
yrk(1)=yr(k);
xik(1)=xi(k);
end
end
MVC.png
(19.12 KB, 下载次数: 5)
2015-12-25 09:29 上传
最小方差
matlab减小数据方差,最小方差控制相关推荐
- 波束形成:最小方差无畸变响应波束形成器(MVDR)
MVDR波束形成器 1. 问题描述 假设远场空间有一个感兴趣的信号 d ( t ) d\left(t\right) d(t),称为期望信号,其波达方向为 θ d \theta_d θd,和 J ...
- 如何用MATLAB优化投资组合的风险(方差)或得到投资组合的有效前沿
如何用MATLAB优化投资组合的风险(方差)或得到投资组合的有效前沿 在金融领域,如果降低投资的风险是一个非常重要的问题,没有风险控制的投资无异于在资本的汪洋大海里裸奔. 我们一般用资产收益率的方差来 ...
- 【51nod - 1098】 最小方差(基础数学,公式化简,前缀和,积的前缀和)
题干: 若x1,x2,x3--xn的平均数为k. 则方差s^2 = 1/n * [(x1-k)^2+(x2-k)^2+--.+(xn-k)^2] . 方差即偏离平方的均值,称为标准差或均方差,方差描述 ...
- python量化策略——大类资产配置模型(最小方差模型)
最小方差模型 寻求风险最小的大类资产组合. max=XTΣX\max=X^{T}\Sigma Xmax=XTΣX s.t.ΣX=1,Xi≥0,i=1,2,3...s.t. \quad \Sigma ...
- 无偏估计和最小方差无偏估计简介
无偏估计和最小方差无偏估计 无偏估计:Unbiased Estimation 最小方差无偏估计: Minimum Variance Unbiased Estimation (MVU) 前言 在正式开始 ...
- 【转】MATLAB的数组期望和方差
设M为一个矩阵: 1.求和 sum(M):以矩阵M的每一列为对象,对每一列的数据分别求和. sum(M,2):以矩阵的每一行为对象,对每一行的数据分别求和. sum(M(:)):将矩阵中的所有元素相加 ...
- 分布估计算法java_数理统计8:点估计的有效性、一致最小方差无偏估计(UMVUE)、零无偏估计法...
在之前的学习中,主要基于充分统计量给出点估计,并且注重于点估计的无偏性与相合性.然而,仅有这两个性质是不足的,无偏性只能保证统计量的均值与待估参数一致,却无法控制统计量可能偏离待估参数的程度:相合性只 ...
- 2019-10-12 线性最小方差估计和正交定理
很多东西再看第二遍的时候,都有一种恍然大悟的感觉 参数估计方法 现代控制理论中往往要通过观测值对系统内部结构参数进行估计,类似于我们常用的线性最小二乘法,举例:二维坐标(身高.年龄)给一堆的点,找俩参 ...
- 【图论】Tree之最小方差树(bzoj 3754)
正题 bzoj 3754 题目大意 求最小方差生成树 解题思路 可以枚举数值之和,然后跑最小生成树 每次加25,这样可以提高速度,且如果平均数接近,那么所选的边和答案会相同,再暴力跑一遍即可 code ...
最新文章
- 游戏控制杆OUYA游戏开发快速入门教程
- 几款开源的数据挖掘工具
- a = 3 中3在栈还是堆_56.堆、栈、方法区(3.0)
- Python手动实现kmeans聚类和调用sklearn实现
- 【Python】手把手教你用Python做一个图像融合demo,小白可上手!
- IIS启用页面Gzip压缩实现说明
- mysql双主 绿色_mysql (双主,互主)
- 双11专刊|云原生数据仓库AnalyticDB支撑双11,大幅提升分析实时性和用户体验
- laravel构造函数和中间件执行顺序问题
- clickhouse 数据存储
- 袁承兴:【译】Async/Await(一)——多任务
- 百旺智能编码_极速开票,智能编码,一键匹配,颠覆了传统手动输入的开票模式,再也不用担心选错税收分类编码啦...
- arctan查表法_DSP正弦波查表法原理
- c#样条曲线命令_C#绘制由点定义的贝塞尔样条曲线
- python安装numpy模块教程_python安装numpy科学计算模块
- LoRa无线模块的优势
- 前端开发之百度地图使用的API和实例
- 通过天眼查查询相关企业信息
- HTML5游戏引擎(二)02-egret引擎之hello world——快速上手-清理项目 程序入口 绘制单色背景 调整屏幕的适配模式 添加文字 响应用户操作-让文字变色
- 程序员专属壁纸十七张
热门文章
- HDLBites学习笔记之Error (10200): Verilog HDL Conditional Statement error at top_module
- 手机系统更新(提示已是最新版本),怎么升级更新
- 乐普生物通过聆讯:预计年亏10亿 海通证券与阳光人寿是股东
- yolov7 打开深度摄像头 realsences
- 查查自己五行属什么命_ywyuan_新浪博客
- 华为手机怎么用云歌_华为手机功能之语音助手小艺,带你了解小艺的使用方式...
- 【大数据开发】scala——tuple、list(含高阶方法)、wordcount案例、set、并行处理数据和sorted、sortBy、sortWith的区别
- 三屏指挥调调度终端/三屏计算机/一机三屏指挥调度/三屏融合指挥调度
- 【postgresql】数据聚合之PERCENTILE_CONT
- win10可用空间变成未分配_有关如何在win10系统中对未分配的磁盘空间进行分区的详细教程...