matlab中卡尔曼滤波函数中变量sys,Qn,Rn,Nn是什么

卡尔曼滤波器的MATLAB程序，实现对一个输入信号的滤波，哪位大神给提供一下，急用！！！~

clear
N=200;
w(1)=0;
x(1)=5;
a=1;
c=1;
Q1 = randn(1,N)*1;%过程噪声
Q2 = randn(1,N);%测量噪声
for k=2:N;x(k)=a*x(k-1)+Q1(k-1);end%状态矩阵
for k=1:N;Y(k)=c*x(k)+Q2(k);end
p(1)=10;
s(1)=1;
for t=2:N;
Rww=cov(Q1(1:t));
Rvv=cov(Q2(1:t));
p1(t)=a.^2*p(t-1)+Rww;
b(t)=c*p1(t)/(c.^2*p1(t)+Rvv);%kalman增益
s(t)=a*s(t-1)+b(t)*(Y(t)-a*c*s(t-1)); p(t)=p1(t)-c*b(t)*p1(t); end
t=1:N;
plot(t,s,'r',t,Y,'g',t,x,'b');%红色卡尔曼，绿色观测值，蓝色状态值 legend('kalman estimate','ovservations','truth');
下面的链接地址就有，希望对您有帮助。
http://wenku.baidu.com/link?url=VExBsM4W1bg9adNs7I6H3EqAST9Jk7GYLpS0XRR0qOyqLucNxa6Jr4240kNsB6CSwVOfLQxPTIJ5I4Q-JM3Uib99uW0gpF8cGvemX7VXPxq

作为一个计量学中不太常用的方法，专门介绍的书应该不多吧
在主要的计量学中应该都有涉及，梁逸增白灰黑复杂多组分分析体系及其化学计量学算法中有算法介绍，许禄的化学计量学有详细介绍还有例子，建议参考。

[kest,L,P] = kalman(sys,Qn,Rn,Nn)
卡尔曼滤波器信号模型
x(k) = A * x(k-1) + w(k)
y(k) = C * x(k) + v(k)

下边的w和v就应该是上边这两个w和v了
E{ww'} = QN,这个是系统噪声协方差阵；
E{vv'} = RN,这个是观测噪声协方差阵；
E{wv'} = NN,这个看字面应该是系统噪声与观测噪声的互协方差阵；
白噪声均值为0，所以上面的几个值也可以认为是自相关函数和互相关函数；
sys给定系统模型；

---

楚雄彝族自治州17860218452： MATLAB中kalman函数里的参数怎么确定 - ？
捷媚十五： 卡尔曼滤波器信号模型 x(k) = A * x(k-1) + w(k) y(k) = C * x(k) + v(k) 下边的w和v就应该是上边这两个w和v了 E{ww'} = QN,这个是系统噪声协方差阵;E{vv'} = RN,这个是观测噪声协方差阵;E{wv'} = NN,这个看字面应该是系统噪声与观测噪声的互协方差阵;这个值的话我当初是自己设的值,w和v的都是高斯噪声,则R和Q应该都是只 有对角线上有值的非奇异矩阵,我假设噪声为单位高斯白噪声,则R和Q都为 单位阵,阶数根据你的系统确定,而NN阵应该为0,因为一般都设两个噪声不相关 这只是我自己的一点见解,如有高人看到欢迎拍砖

楚雄彝族自治州17860218452： 求:卡尔曼滤波在变形监测数据处理中应用时,用matlab实现的程序,需要具体实例程序!其中数据为沉降量. - ？
捷媚十五： [r,c] = find(R == max(R(:))); 检索R中最大元素所在的位置(行标r和列标c) thetap = theta(c(1)); theta()是自定义函数

楚雄彝族自治州17860218452： matlab中的局部变量和全局变量是什么意思 - ？
捷媚十五： global 是定义全局变量的 局部变量是单独某个函数下的变量 在其他函数中不能使用该数值...

楚雄彝族自治州17860218452： 请问,卡尔曼滤波方程中Z(k)=H(k)X(k)+V(k)中的观测矩阵H(k)是怎么得到的呢? - ？
捷媚十五： 为了推广卡尔曼滤波对多个变量的计算,矩阵是一个同时刻的一组测量.观测矩阵H(k)就是你的测量矩阵.

楚雄彝族自治州17860218452： MATLAB中实现扩展卡尔曼和无味卡尔曼滤波,想对一条简单的曲线和直线进行滤波, - ？
捷媚十五： % INPUTS:% y(:,t) - the observation at time t% A - the system matrix% C - the observation matrix % Q - the system covariance % R - the observation covariance% init_x - the initial state (column) vector % init_V - the initial state covariance

楚雄彝族自治州17860218452： 卡尔曼滤波中的过程噪声与观测噪声选取问题 - ？
捷媚十五： Q表征了系统模型的统计特性,增加Q中元素的值,等价于增加系统噪声或增加系统参数的不确定性,从而使得增益矩阵Gf增大,加大了系统校正权值,提高了系统动态性能和稳态值. 矩阵 R 表征了测量噪声, 增加 R 中元素的值, 意味加大测量噪声的影响, 同时使 Gf 减小, 减弱系统校正权值, 降低了系统瞬态响应和稳态值. 所以矩阵 Q 和 R 的选取具有矛盾性, 二者进行选取时, 应视具体情况权衡而定.

楚雄彝族自治州17860218452： matlab中怎么在函数中输出一个变量 - ？
捷媚十五： 1、很简单啊. 2、比如我们输出一个函数f=A*sin(x)-B*exp(C*x)这种表达式,A,B,C是你输入的任何参数. (1)你写上如下代码的M文件 function f=dispf(A,B,C) syms x; f=A*sin(x)-B*exp(C*x); end (2)现在你直接在命令窗口输入命令dispf就可以

楚雄彝族自治州17860218452： 在matlab中如何查看并获得Embedded Matlab Function 中的变量值？
捷媚十五： 可能你的程序中使用了动态数组,例如: y=1; y(2)=2; 这样的语句在matlab中运行当然没问题,数组会动态扩展,但在embedded matlab function模块中不允许这样用.如果还查不出,可以把m文件代码贴出来帮你分析.

楚雄彝族自治州17860218452： 话说matlab 中spear函数可以求出两个向量或矩阵的spearman 相关系数,但是数理统计工具箱中并没有这个函数 - ？
捷媚十五： 的确没有spear函数,我们是用corr来计算的: [RHO,PVAL] = corr(X,Y,'name',value) 其中name可以是type,rows,tail,而value分别如下: type: 'Pearson' (the default) computes Pearson's linear correlation coefficient'Kendall' computes Kendall's tau...

楚雄彝族自治州17860218452： matlab数字滤波器系数如何量化 - ？
捷媚十五： 一、量化方法:先求出浮点系数,然后用Matlab里的量化函数进行量化,也就是点左下方第三个图标就是量化键就可.二、matlab的使用技巧:1、. Ctrl+C 中断正在执行的操作如果程序不小心进入死循环,或者计算时间太长,可以在命令窗...